Unfolding Transition Research Articles

The deconvolution of differential scanning calorimetry unfolding transitions

This paper is a review of a process for deconvolution of unfolding thermal transitions measured by differential scanning calorimetry. The mathematical background is presented along with illustrations of how the unfolding data is processed to resolve the number of sequential transitions needed to describe an unfolding mechanism and to determine thermodynamic properties of the intermediate states. Examples of data obtained for a simple two-state unfolding of a G-quadruplex DNA structure derived from the basic human telomere sequence, (TTAGGG)4TT are used to present some of the basic issues in treating the DSC data. A more complex unfolding mechanism is also presented that requires deconvolution of a multistate transition, the unfolding of a related human telomere structure, (TTAGGG)12 TT. The intent of the discussion is to show the steps in deconvolution, and to present the data at each step to help clarify how the information is derived from the various mathematical manipulations.

Dissecting cooperative communications in a protein with a high-throughput single-molecule scalpel.

Miscued communication often leads to misfolding and aggregation of the proteins involved in many diseases. Owing to the ensemble average property of conventional techniques, detailed communication diagrams are difficult to obtain. Mechanical unfolding affords an unprecedented perspective on cooperative transitions by observing a protein along a trajectory defined by two mutated cysteine residues. Nevertheless, this approach requires tedious sample preparation at the risk of altering native protein conformations. To address these issues, we applied click chemistry to tether a protein to the two dsDNA handles through primary amines in lysine residues as well as at the N terminus. As a proof of concept, we used laser tweezers to mechanically unfold and refold calmodulin along 36 trajectories, maximally allowed by this strategy in a single batch of protein preparation. Without a priori knowledge of the particular residues to which the double-stranded DNA handles attach, we used hierarchical cluster analysis to identify 20 major trajectories, according to the size and the pattern of unfolding transitions. We dissected the cooperativity into all-or-none and partially cooperative events, which represent strong and weak high-order interactions in proteins, respectively. Although the overall cooperativity is higher within the N or C lobe than that between the lobes, the all-or-none cooperativity is anisotropic among different the unfolding trajectories and becomes relatively more predominant when the size of the protein segments increases. The average cooperativity for all-or-none transitions falls within the expected range observed by ensemble techniques, which supports the hypothesis that unfolding of a free protein can be reconstituted from individual trajectories.

Interactions of urea with native and unfolded proteins: a volumetric study.

We describe a statistical thermodynamic approach to analyzing urea-dependent volumetric properties of proteins. We use this approach to analyze our urea-dependent data on the partial molar volume and adiabatic compressibility of lysozyme, apocytochrome c, ribonuclease A, and α-chymotrypsinogen A. The analysis produces the thermodynamic properties of elementary urea-protein association reactions while also yielding estimates of the effective solvent-accessible surface areas of the native and unfolded protein states. Lysozyme and apocytochrome c do not undergo urea-induced transitions. The former remains folded, while the latter is unfolded between 0 and 8 M urea. In contrast, ribonuclease A and α-chymotrypsinogen A exhibit urea-induced unfolding transitions. Thus, our data permit us to characterize urea-protein interactions in both the native and unfolded states. We interpreted the urea-dependent volumetric properties of the proteins in terms of the equilibrium constant, k, and changes in volume, ΔV0, and compressibility, ΔKT0, for a reaction in which urea binds to a protein with a concomitant release of two waters of hydration to the bulk. Comparison of the values of k, ΔV0, and ΔKT0 with the similar data obtained on small molecules mimicking protein groups reveals lack of cooperative effects involved in urea-protein interactions. In general, the volumetric approach, while providing a unique characterization of cosolvent-protein interactions, offers a practical way for evaluating the effective solvent accessible surface area of biologically significant fully or partially unfolded polypeptides.

PH-Induced conformational isomerization of leghemoglobin from Arachis hypogea.

The pH dependence of proteins is related to the thermodynamic stability and electrostatic interactions in the native state of a protein. Here we report the pH-induced conformational transition of the heme protein leghemoglobin (Lb) isolated from root nodules of the leguminous plant Arachis hypogea. Unlike the other heme proteins myoglobin, hemoglobin, and cytochrome c, the structural characteristics and interactions of Lb is almost unknown, though its functional importance is already established since it binds oxygen to maintain the environment for N2 fixation. We investigated pH-induced unfolding of this protein and identified a number of conformational isomers using multiple fluorescence observables as a function of pH titration. We have characterized the acid- and base-induced conformational transitions among the structural states over the pH range 2-11. Depending on the solution conditions, Lb can exist in one of three phases: pH 2, 3, 4; pH 5, 6, 7; pH 8, 9, 10. The secondary structure as revealed by CD spectroscopy indicated the maximum percentage of α-helix to be present at pH 7, where the structure of Lb is also most rigid according to fluorescence anisotropy experiments. The fluorescence lifetime of tryptophan was observed to be maximum at pH 10 and minimum at pH 6, suggesting unfolding transitions of Lb. Thus, alteration of the microenvironment of the globin moiety during pH transition ultimately leads to the conformational change of this monomeric protein Lb.

A role for the Parkinson's disease protein DJ-1 as a chaperone and antioxidant in the anhydrobiotic nematode Panagrolaimus superbus.

Mutations in the human DJ-1/PARK7 gene are associated with familial Parkinson's disease. DJ-1 belongs to a large, functionally diverse family with homologues in all biological kingdoms. Several activities have been demonstrated for DJ-1: an antioxidant protein, a redox-regulated molecular chaperone and a modulator of multiple cellular signalling pathways. The majority of functional studies have focussed on human DJ-1 (hDJ-1), but studies on DJ-1 homologues in Drosophila melanogaster, Caenorhabditis elegans, Dugesia japonica and Escherichia coli also provide evidence of a role for DJ-1 as an antioxidant. Here, we show that dehydration is a potent inducer of a dj-1 gene in the anhydrobiotic nematode Panagrolaimus superbus. Our secondary structure and homology modelling analyses shows that recombinant DJ-1 protein from P. superbus (PsuDJ-1.1) is a well-folded protein, which is similar in structure to the hDJ-1. PsuDJ-1.1 is a heat stable protein; with T1/2 unfolding transition values of 76 and 70°C obtained from both circular dichroism (CD) and Fourier transform infrared spectroscopy (FTIR) measurements respectively. We found that PsuDJ-1.1 is an efficient antioxidant that also functions as a 'holdase' molecular chaperone that can maintain its chaperone function in a reducing environment. In addition to its chaperone activity, PsuDJ-1.1 may also be an important non-enzymatic antioxidant, capable of providing protection to P. superbus from oxidative damage when the nematodes are in a desiccated, anhydrobiotic state.

Unfolding properties of recombinant human serum albumin products are due to bioprocessing steps.

We have used differential scanning calorimetry (DSC) to determine the unfolding properties of commercial products of human serum albumin (HSA) prepared from pooled human blood, transgenic yeast, and transgenic rice. The initial melting temperatures (Tm1 ) for the unfolding transitions of the HSA products varied from 62°C to 75°C. We characterized the samples for purity, fatty acid content, and molecular weight. The effects of adding fatty acids, heat pasteurization, and a low pH defatting technique on the transition temperatures were measured. Defatted HSA has a structure with the lowest stability (Tm of ∼62°C). When fatty acids are bound to HSA, the structure is stabilized (Tm of ∼64-72°C), and prolonged heating (pasteurization at 60°C) results in a heat-stabilized structural form containing fatty acids (Tm of ∼75-80°C). This process was shown to be reversible by a low pH defatting step. This study shows that the fatty acid composition and bioprocessing history of the HSA commercial products results in the large differences in the thermal stability.

Residues Phe103 and Phe149 are critical for the co-chaperone activity of Bacillus licheniformis GrpE

A tryptophan-free Bacillus licheniformis nucleotide exchange factor (BlGrpE) and its Trp mutants (F70W, F103W, F149W, F70/103W, F70/149W, F103/149W and F70/103/149W) were over-expressed and purified to near homogeneity. Simultaneous addition of B. licheniformis DnaJ, NR-peptide and individual variants synergistically stimulated the ATPase activity of a recombinant DnaK (BlDnaK) from the same bacterium by 3.1–14.7-fold, which are significantly lower than the synergistic stimulation (18.9-fold) of BlGrpE. Protein–protein interaction analysis revealed that Trp mutants relevant to amino acid positions 103 and 149 lost the ability to bind BlDnaK. Circular dichroism measurements indicate that F70W displayed a comparable level of secondary structure to that of BlGrpE, and the wild-type protein and the Trp mutants as well all experienced a reversible behavior of thermal denaturation. Guanidine hydrochloride (GdnHCl)-induced unfolding transition for BlGrpE was calculated to be 1.25M corresponding to ΔGN–U of 4.29kcal/mol, whereas the unfolding transitions of mutant proteins were in the range of 0.77–1.31M equivalent to ΔGN–U of 2.41–4.14kcal/mol. Taken together, the introduction of tryptophan residue, especially at positions 103 and 149, into the primary structure of BlGrpE has been proven to be detrimental to structural integrity and proper function of the protein.

Circuit topology of self-interacting chains: implications for folding and unfolding dynamics.

Understanding the relationship between molecular structure and folding is a central problem in disciplines ranging from biology to polymer physics and DNA origami. Topology can be a powerful tool to address this question. For a folded linear chain, the arrangement of intra-chain contacts is a topological property because rearranging the contacts requires discontinuous deformations. Conversely, the topology is preserved when continuously stretching the chain while maintaining the contact arrangement. Here we investigate how the folding and unfolding of linear chains with binary contacts is guided by the topology of contact arrangements. We formalize the topology by describing the relations between any two contacts in the structure, which for a linear chain can either be in parallel, in series, or crossing each other. We show that even when other determinants of folding rate such as contact order and size are kept constant, this 'circuit' topology determines folding kinetics. In particular, we find that the folding rate increases with the fractions of parallel and crossed relations. Moreover, we show how circuit topology constrains the conformational phase space explored during folding and unfolding: the number of forbidden unfolding transitions is found to increase with the fraction of parallel relations and to decrease with the fraction of series relations. Finally, we find that circuit topology influences whether distinct intermediate states are present, with crossed contacts being the key factor. The approach presented here can be more generally applied to questions on molecular dynamics, evolutionary biology, molecular engineering, and single-molecule biophysics.

Stability and Activity of Porcine Lipase Against Temperature and Chemical Denaturants

Lipases are the class of hydrolases with wide industrial applications. The present study analyses the stability of porcine pancreatic lipase (PPL) against urea, guanidine hydrochloride (Gdn), sodium dodecyl sulphate (SDS), and temperature using different spectroscopic techniques. Interestingly, this two-domain protein shows a two-state unfolding transition against urea and Gdn. The free energy of unfolding of PPL calculated from global analysis of the unfolding transitions obtained from different spectroscopic techniques is ~2.2 kcal/mol. In the presence of SDS, PPL shows a cooperative loss of secondary and tertiary structures above 0.2 mM of SDS. At above 2 mM of SDS, PPL forms irreversible, non-native, thermally stable structure. PPL loses its activity even at lower concentrations of urea (3 M), Gdn (0.5 M), and SDS (0.8 mM). Thermal denaturation of PPL shows an irreversible unfolding, and the protein lost its activity even by increasing the temperature to 45 °C. Though PPL in higher concentrations of SDS (>5 mM) shows stable conformation against temperature, its activity is completely lost. The results suggest that the structure and activity of PPL are more sensitive against chemical denaturants and temperature, and forms irreversible, non-native (in SDS) or completely unfolded (in urea, Gdn, and at higher temperature) conformations in different denaturing conditions.

Thermodynamic description of Hofmeister effects on the LCST of thermosensitive polymers.

Cosolvent effects on protein or polymer collapse transitions are typically discussed in terms of a two-state free energy change that is strictly linear in cosolute concentration. Here we investigate in detail the nonlinear thermodynamic changes of the collapse transition occurring at the lower critical solution temperature (LCST) of the role-model polymer poly(N-isopropylacrylamide) [PNIPAM] induced by Hofmeister salts. First, we establish an equation, based on the second-order expansion of the two-state free energy in concentration and temperature space, which excellently fits the experimental LCST curves and enables us to directly extract the corresponding thermodynamic parameters. Linear free energy changes, grounded on generic excluded-volume mechanisms, are indeed found for strongly hydrated kosmotropes. In contrast, for weakly hydrated chaotropes, we find significant nonlinear changes related to higher order thermodynamic derivatives of the preferential interaction parameter between salts and polymer. The observed non-monotonic behavior of the LCST can then be understood from a not yet recognized sign change of the preferential interaction parameter with salt concentration. Finally, we find that solute partitioning models can possibly predict the linear free energy changes for the kosmotropes, but fail for chaotropes. Our findings cast strong doubt on their general applicability to protein unfolding transitions induced by chaotropes.

Structural stability of E. coli trigger factor studied by synchrotron small-angle X-ray scattering

Solution small-angle X-ray scattering (SAXS) is an effective technique for quantitatively measuring the compactness and shape of proteins. We use SAXS to study the structural characteristics and unfolding transitions induced by urea for full length Escherichia coli trigger factor (TF) and a series of truncation mutants, obtaining and comparing the radiuses of gyration (Rg), the distance-distribution function (P(r) function) and integrated intensity of TF variants in native and unfolding states. The C-terminal 72-residue truncated mutant TF360 exhibited dramatic structural differences and reduced stability compared with the whole TF molecule, while the N-domain truncated mutant MC maintained its compact structure with reduced stability. These results indicate that the C-terminal region of TF plays an important role in the structural and conformational stabilities of the TF molecule, while the N-domain is relatively independent.

Collisional and Coulombic Unfolding of Gas‐Phase Proteins: High Correlation to Their Domain Structures in Solution

AbstractThe three‐dimensional structures adopted by proteins are predicated by their many biological functions. Mass spectrometry has played a rapidly expanding role in protein structure discovery, enabling the generation of models for both proteins and their higher‐order assemblies. While important coursed‐grained insights have been generated, relatively few examples exist where mass spectrometry has been successfully applied to the characterization of protein tertiary structure. Here, we demonstrate that gas‐phase unfolding can be used to determine the number of autonomously folded domains within monomeric proteins. Our ion mobility‐mass spectrometry data highlight a strong, positive correlation between the number of protein unfolding transitions observed in the gas phase and the number of known domains within a group of sixteen proteins ranging from 8–78 kDa. This correlation and its potential uses for structural biology is discussed.

Collisional and Coulombic unfolding of gas-phase proteins: high correlation to their domain structures in solution.

The three-dimensional structures adopted by proteins are predicated by their many biological functions. Mass spectrometry has played a rapidly expanding role in protein structure discovery, enabling the generation of models for both proteins and their higher-order assemblies. While important coursed-grained insights have been generated, relatively few examples exist where mass spectrometry has been successfully applied to the characterization of protein tertiary structure. Here, we demonstrate that gas-phase unfolding can be used to determine the number of autonomously folded domains within monomeric proteins. Our ion mobility-mass spectrometry data highlight a strong, positive correlation between the number of protein unfolding transitions observed in the gas phase and the number of known domains within a group of sixteen proteins ranging from 8-78 kDa. This correlation and its potential uses for structural biology is discussed.

Structure and functional investigation of ligand binding by a family 35 carbohydrate binding module (CtCBM35) of β-mannanase of family 26 glycoside hydrolase from Clostridium thermocellum.

Functional attributes of recombinant CtCBM35 (family 35 carbohydrate binding module) of β-mannanase of family 26 Glycoside Hydrolase from Clostridium thermocellum were deduced by biochemical and in silico approaches. Ligand-binding analysis of expressed CtCBM35 analyzed by affinity-gel electrophoresis and fluorescence spectroscopy exhibited association constants Ka ~ 1.2·10(5)and 3.0·10(5)M(-1) with locust bean galactomannan and mannotriose, respectively. However, CtCBM35 showed low ligand-binding affinity with insoluble ivory nut mannan with Ka of 5.0·10(-5)M(-1). Unfolding transition analysis by fluorescence spectroscopy explained the conformational changes of CtCBM35 in the presence of guanidine hydrochloride (5 M) and urea (6.25M). This explained that CtCBM35 has good conformational stability and requires higher free energy of denaturation to invoke unfolding. The three-dimensional (3-D) model of CtCBM35 from C. thermocellum generated by Modeller9v8 displayed predominance of β-sheets arranged as β-jelly-roll fold. The secondary structure of CtCBM35 by PredictProtein showed the presence of two α-helices (3%), 12 β-sheets (45%), and 15 random coils (52%). Secondary structural element analysis of cloned, expressed, and purified recombinant CtCBM35 by circular dichroism also corroborated the in silico predicted secondary structure. Multiple sequence alignment of CtCBM35 showed conserved residues (Tyr123, Gly124, and Phe125), which are commonly observed in mannan specific CBMs. Docking analysis of CtCBM35 with manno-oligosaccharide displayed the involvement of Tyr26, Gln29, Asn43, Trp66, Tyr68, Leu69, Arg76, and Leu127 residues, making polar contact with the ligand molecules. Ligand docking analysis of CtCBM35 exhibiting higher binding affinity with mannotriose and galactomannan (Man-Gal-Man moiety) substantiated the affinity binding and fluorescence results, displaying similar values of Ka.

Ligand-binding pocket bridges DNA-binding and dimerization domains of the urate-responsive MarR homologue MftR from Burkholderia thailandensis.

Members of the multiple antibiotic resistance regulator (MarR) family often regulate gene activity by responding to a specific ligand. In the absence of ligand, most MarR proteins function as repressors, while ligand binding causes attenuated DNA binding and therefore increased gene expression. Previously, we have shown that urate is a ligand for MftR (major facilitator transport regulator), which is encoded by the soil bacterium Burkholderia thailandensis. We show here that both mftR and the divergently oriented gene mftP encoding a major facilitator transport protein are upregulated in the presence of urate. MftR binds two cognate sites in the mftR-mftP intergenic region with equivalent affinity and sensitivity to urate. Mutagenesis of four conserved residues previously reported to be involved in urate binding to Deinococcus radiodurans HucR and Rhizobium radiobacter PecS significantly reduced protein stability and DNA binding affinity but not ligand binding. These data suggest that residues equivalent to those implicated in ligand binding to HucR and PecS serve structural roles and that MftR relies on distinct residues for ligand binding. MftR exhibits a two-step melting transition suggesting independent unfolding of the dimerization and DNA-binding regions; urate binding or mutations in the predicted ligand-binding sites result in one-step unfolding transitions. We suggest that MftR binds the ligand in a cleft between the DNA-binding lobes and the dimer interface but that the mechanism of ligand-mediated attenuation of DNA binding differs from that proposed for other urate-responsive MarR homologues. Since DNA binding by MftR is attenuated at 37 °C, our data also suggest that MftR responds to both ligand and a thermal upshift by attenuated DNA binding and upregulation of the genes under its control.

Distinguishing unfolding and functional conformational transitions of calmodulin using ultraviolet resonance Raman spectroscopy.

Calmodulin (CaM) is a ubiquitous moderator protein for calcium signaling in all eukaryotic cells. This small calcium-binding protein exhibits a broad range of structural transitions, including domain opening and folding-unfolding, that allow it to recognize a wide variety of binding partners in vivo. While the static structures of CaM associated with its various binding activities are fairly well-known, it has been challenging to examine the dynamics of transition between these structures in real-time, due to a lack of suitable spectroscopic probes of CaM structure. In this article, we examine the potential of ultraviolet resonance Raman (UVRR) spectroscopy for clarifying the nature of structural transitions in CaM. We find that the UVRR spectral change (with 229 nm excitation) due to thermal unfolding of CaM is qualitatively different from that associated with opening of the C-terminal domain in response to Ca(2+) binding. This spectral difference is entirely due to differences in tertiary contacts at the interdomain tyrosine residue Tyr138, toward which other spectroscopic methods are not sensitive. We conclude that UVRR is ideally suited to identifying the different types of structural transitions in CaM and other proteins with conformation-sensitive tyrosine residues, opening a path to time-resolved studies of CaM dynamics using Raman spectroscopy.

Effects of crowding on the stability of a surface-tethered biopolymer: an experimental study of folding in a highly crowded regime.

The high packing densities and fixed geometries with which biomolecules can be attached to macroscopic surfaces suggest that crowding effects may be particularly significant under these often densely packed conditions. Exploring this question experimentally, we report here the effects of crowding on the stability of a simple, surface-attached DNA stem-loop. We find that crowding by densely packed, folded biomolecules destabilizes our test-bed biomolecule by ∼2 kJ/mol relative to the dilute (noninteracting) regime, an effect that presumably occurs due to steric and electrostatic repulsion arising from compact neighbors. Crowding by a dense brush of unfolded biomolecules, in contrast, enhances its stability by ∼6 kJ/mol, presumably due to excluded volume and electrostatic effects that reduce the entropy of the unfolded state. Finally, crowding by like copies of the same biomolecule produces a significantly broader unfolding transition, likely because, under these circumstances, the stabilizing effects of crowding by unfolded molecules increase (and the destabilizing effects of neighboring folded molecules decrease) as more and more neighbors unfold. The crowding of surface-attached biomolecules may thus be a richer, more complex phenomenon than that seen in homogeneous solution.

Free energy landscape of lysozyme: multiple near-native conformational states and rollover in the urea dependence of folding energy.

Deviation from linearity of the equilibrium folding free energy (ΔG) of proteins along the reaction coordinate is scarcely known. Optical spectroscopic observables and NMR-measured average molecular dimensional property of lysozyme with urea at pH 5 reveal that ΔG rolls over from linearity under mild to strongly native-like conditions. The urea dependence of ΔG is graphed in the 0-7 M range of the denaturant by employing a series of guanidine hydrochloride (GdnHCl)-induced equilibrium unfolding transitions, each in the presence of a fixed level of urea. The observed linear dependence of ΔG on urea under denaturing conditions begins to deviate as moderately native-like conditions are approached and eventually rolls over under strongly native-like conditions. This is atypical of the upward curvature in the ΔG vs denaturant plot predicted by the denaturant binding model. On increasing the denaturant concentration from 0 to 5 M, the hydrodynamic radius of lysozyme shrinks by ∼2 Å. We suggest subdenaturing levels of urea affect the population distribution among multiple near-native isoenergetic conformational states so as to promote them sequentially with increments of the denaturant. We use a multiple-state sequential model to show that the keel over of ΔG occurs due to these near-native alternative states in the native ensemble used for defining the unfolding equilibrium constant (KU), which we assume to vary linearly with urea. The results and the model appear to indicate a rugged flat bottom in the free energy landscape wherein population distribution of native-like states is modulated by urea-affected interstate motions.

Quantitative thermodynamic model for globular protein folding

We present a statistical mechanics formalism for theoretical description of the process of protein folding ↔ unfolding transition in water environment. The formalism is based on the construction of the partition function of a protein obeying two-stage-like folding kinetics. Using the statistical mechanics model of solvation of hydrophobic hydrocarbons we obtain the partition function of infinitely diluted solution of proteins in water environment. The calculated dependencies of the protein heat capacities upon temperature are compared with the corresponding results of experimental measurements for staphylococcal nuclease and metmyoglobin.

Evaluation of protein formulation and its viscosity with DSC, DLS, and microviscometer

The viscosity of highly concentrated protein solutions was evaluated using lysozyme as model protein. Viscosity profiles of lysozyme were examined with the effect of buffer and pH-value at various concentrations. The viscosity of lysozyme dissolved in water increased continuously with the concentration as the slope of shear stress against shear rate increased with the concentration. In addition, the viscosity of lysozyme was higher in histidine buffer than in acetate buffer at selected pH ranges. The effect of various excipient concentrations was also investigated in means of unfolding transition temperature (Tm), viscosity, hydrodynamic size and zeta potential by using differential scanning calorimetry (DSC), microviscometer and dynamic light scattering (DLS). The selected excipients except surfactants increased the viscosity of protein solution with their concentration. Carbohydrates increased the viscosity relatively higher than amino acids and also they increased the conformational stability (Tm) by enhancing the protein molecule more in compact form. Also amino acids increased the viscosity but decreased the conformational stability since they seemed to be only dispersed in the solution avoiding protein–protein interactions, resulting in a decrease of zeta potential. Consequently, the applied methods—DSC, DLS and microviscometer demonstrated the potential to develop a highly concentrated protein formulation to decrease the high viscosity effect with acceptable conformational stability.