Thermal Unfolding Curves Research Articles

Determining the Conformational Stability of a Protein from Urea and Thermal Unfolding Curves

This article contains protocols for determining the conformational stability of a globular protein from either urea or thermal unfolding curves. Circular dichroism is the optical spectroscopic technique most commonly used to monitor protein unfolding. These protocols describe how to analyze data from an unfolding curve to obtain the thermodynamic parameters necessary to calculate conformational stability, and how to determine differences in stability between protein variants. Curr. Protoc. Protein Sci. 71:28.4.1-28.4.14. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Determining protein conformational stability from urea-induced unfolding curves Support Protocol 1: Preparing a urea stock solution Support Protocol 2: Analyzing urea unfolding curves Basic Protocol 2: Determining the conformational stability of a protein from thermal unfolding curves Support Protocol 3: Analyzing thermal unfolding curves Support Protocol 4: Determining differences in conformational stability for protein variants.

Residue-Dependent Transition Temperatures and Denaturant Midpoints in the Folding of a Multidomain Protein.

As a consequence of the finite size of globular proteins, it is expected that there should be dispersions in the global melting temperature (Tm) and the denaturation midpoint (Cm). Thermodynamic considerations dictate that the dispersions, ΔTm in Tm, and ΔCm in Cm, should decrease with N, the number of residues in the protein. We performed coarse-grained simulations of the self-organized polymer (SOP) model of the multidomain protein adenylate kinase (ADK) with N = 214 in order to calculate thermal and denaturation unfolding titration curves. The results show that ΔTm/Tm and ΔCm/Cm are nonzero and follow the previously established ( Phys. Rev. Lett. 2004, 93, 268107) thermodynamic 1/N scaling for proteins accurately. For ADK, the dispersions are small (≈0.004), which implies that the melting temperature is more or less unique, which is unlike in BBL (N = 40) where ΔTm/Tm ≈ 0.03.

Sequence-specific destabilization of azurin by tetramethylguanidinium-dipeptide ionic liquids

The thermal unfolding of the copper redox protein azurin was studied in the presence of four different dipeptide-based ionic liquids (ILs) utilizing tetramethylguanidinium as the cation. The four dipeptides have different sequences including the amino acids Ser and Asp: TMG-AspAsp, TMG-SerSer, TMG-SerAsp, and TMG-AspSer. Thermal unfolding curves generated from temperature-dependent fluorescence spectroscopy experiments showed that TMG-AspAsp and TMG-SerSer have minor destabilizing effects on the protein while TMG-AspSer and TMG-SerAsp strongly destabilize azurin. Red-shifted fluorescence signatures in the 25 °C correlate with the observed protein destabilization in the solutions with TMG-AspSer and TMG-SerAsp. These signals could correspond to interactions between the Asp residue in the dipeptide and the azurin Trp residue in the unfolded state. These results, supported by appropriate control experiments, suggest that dipeptide sequence-specific interactions lead to selective protein destabilization and motivate further studies of TMG-dipeptide ILs.

Protein crosslinking improves the thermal resistance of plastocyanin immobilized on a modified gold electrode

Increasing the thermal stability of immobilized proteins is a motivating goal for improving the performance of electrochemical biodevices. In this work, we propose the immobilization of crosslinked plastocyanin from the thermophilic cyanobacterium Phormidium laminosum by simultaneous incubation of a mixture of plastocyanin and the coupling reagents. The thermal stability of the so built covalently immobilized protein films has been assessed by cyclic voltammetry in the 0–90 °C temperature range and has been compared to that of physisorbed films. It is shown that the protein loss along a thermal cycle is significantly reduced in the case of the crosslinked films, whose redox properties remain unaltered along a cyclic heating-cooling thermal scan, and can withstand the contact with 70 °C solutions for four hours. Comparison of thermal unfolding curves obtained by circular dichroism spectroscopy of both free and crosslinked protein confirms the improved thermic resistance of the crosslinked plastocyanin. Notably, the electron transfer thermodynamics of physisorbed and crosslinked plastocyanin films are quite similar, suggesting that the formation of intra- and inter-protein amide bonds do not affect the integrity and functionality of the copper redox centers. UV–Vis absorption and circular dichroism measurements corroborate that protein crosslinking does not alter the coordination geometry of the metal center.

Structural analyses combined with small-angle X-ray scattering reveals that the retention of heme is critical for maintaining the structure of horseradish peroxidase under denaturing conditions.

We analyzed the structure of horseradish peroxidase (HRP) under denaturing conditions of 9M urea or 6M guanidine hydrochloride (GdnHCl). Far-UV circular dichroism (CD) spectra indicated the existence of native-like secondary structure of holo-HRP in 9M urea. In addition, slight changes in near-UV and Soret region CD spectra of holo-HRP in 9M urea suggest that the tertiary structure of holo-HRP and the binding of heme remain partially intact in this condition. A transition in the thermal unfolding transition curve of holo-HRP in 9M urea indicated the existence of a considerable amount of secondary structure. However, no secondary structure, tertiary structure, or interaction between heme and HRP were observed in holo-HRP in 6M GdnHCl. Small-angle X-ray scattering indicated that although distal and proximal domains of holo-HRP in 9M urea might be partially unfolded, the central region that contains the heme might maintain its tertiary structure. Our results suggest that retention of the heme is essential for maintenance of the structure of HRP under highly denaturing conditions.

Microscopic nucleation and propagation rates of an alanine-based α-helix.

An infrared temperature-jump (T-jump) study by Huang et al. (Proc. Natl. Acad. Sci. U. S. A., 2002, 99, 2788-2793) showed that the conformational relaxation kinetics of an alanine-based α-helical peptide depend not only on the final temperature (Tf) but also on the initial temperature (Ti) when Tf is fixed. Their finding indicates that the folding free energy landscape of this peptide is non-two-state like, allowing for the population of conformational ensembles with different helical lengths and relaxation times in the temperature range of the experiment. Because α-helix folding involves two fundamental events, nucleation and propagation, the results of Huang et al. thus present a unique opportunity to determine their rate constants - a long-sought goal in the study of the helix-coil transition dynamics. Herein, we capitalize on this notion and develop a coarse-grained kinetic model to globally fit the thermal unfolding curve and T-jump kinetic traces of this peptide. Using this strategy, we are able to explicitly determine the microscopic rate constants of the kinetic steps encountered in the nucleation and propagation processes. Our results reveal that the time taken to form one α-helical turn (i.e., an α-helical segment with one helical hydrogen bond) is about 315 ns, whereas the time taken to elongate this nucleus by one residue (or backbone unit) is 5.9 ns, depending on the position of the residue.

Adaptive thermostability of light-harvesting complexes in marine picocyanobacteria.

Marine Synechococcus play a key role in global oceanic primary productivity. Their wide latitudinal distribution has been attributed to the occurrence of lineages adapted to distinct thermal niches, but the physiological and molecular bases of this ecotypic differentiation remain largely unknown. By comparing six strains isolated from different latitudes, we showed that the thermostability of their light-harvesting complexes, called phycobilisomes (PBS), varied according to the average sea surface temperature at strain isolation site. Comparative analyses of thermal unfolding curves of the three phycobiliproteins (PBP) constituting PBS rods suggested that the differences in thermostability observed on whole PBSs relied on the distinct molecular flexibility and stability of their individual components. Phycocyanin was the least thermostable of all rod PBP, constituting a fragility point of the PBS under heat stress. Amino-acid composition analyses and structural homology modeling notably revealed the occurrence of two amino-acid substitutions, which might have a role in the observed differential thermotolerance of this phycobiliprotein among temperature ecotypes. We hypothesize that marine Synechococcus ancestors occurred first in warm niches and that during the colonization of cold, high latitude thermal niches, their descendants have increased the molecular flexibility of PBP to maintain optimal light absorption capacities, this phenomenon likely resulting in a decreased stability of these proteins. This apparent thermoadaptability of marine Synechococcus has most probably contributed to the remarkable ubiquity of these picocyanobacteria in the ocean.

ProteoPlex: stability optimization of macromolecular complexes by sparse-matrix screening of chemical space.

Molecular machines or macromolecular complexes are supramolecular assemblies of biomolecules with a variety of functions. Structure determination of these complexes in a purified state is often tedious owing to their compositional complexity and the associated relative structural instability. To improve the stability of macromolecular complexes in vitro, we present a generic method that optimizes the stability, homogeneity and solubility of macromolecular complexes by sparse-matrix screening of their thermal unfolding behavior in the presence of various buffers and small molecules. The method includes the automated analysis of thermal unfolding curves based on a biophysical unfolding model for complexes. We found that under stabilizing conditions, even large multicomponent complexes reveal an almost ideal two-state unfolding behavior. We envisage an improved biochemical understanding of purified macromolecules as well as a substantial boost in successful macromolecular complex structure determination by both X-ray crystallography and cryo-electron microscopy.

A comprehensive alanine-scanning mutagenesis study reveals roles for salt bridges in the structure and activity of Pseudomonas aeruginosa elastase.

The relationship between salt bridges and stability/enzymatic activity is unclear. We studied this relationship by systematic alanine-scanning mutation analysis using the typical M4 family metalloprotease Pseudomonas aeruginosa elastase (PAE, also known as pseudolysin) as a model. Structural analysis revealed seven salt bridges in the PAE structure. We constructed ten mutants for six salt bridges. Among these mutants, six (Asp189Ala, Arg179Ala, Asp201Ala, Arg205Ala, Arg245Ala and Glu249Ala) were active and four (Asp168Ala, Arg198Ala, Arg253Ala, and Arg279Ala) were inactive. Five mutants were purified, and their catalytic efficiencies (k cat/K m), half-lives (t 1/2) and thermal unfolding curves were compared with those of PAE. Mutants Asp189Ala and Arg179Ala both showed decreased thermal stabilities and increased activities, suggesting that the salt bridge Asp189-Arg179 stabilizes the protein at the expense of catalytic efficiency. In contrast, mutants Asp201Ala and Arg205Ala both showed slightly increased thermal stability and slightly decreased activity, suggesting that the salt bridge Asp201-Arg205 destabilizes the protein. Mutant Glu249Ala is related to a C-terminal salt bridge network and showed both decreased thermal stability and decreased activity. Furthermore, Glu249Ala showed a thermal unfolding curve with three discernable states [the native state (N), the partially unfolded state (I) and the unfolded state (U)]. In comparison, there were only two discernable states (N and U) in the thermal unfolding curve of PAE. These results suggest that Glu249 is important for catalytic efficiency, stability and unfolding cooperativity. This study represents a systematic mutational analyses of salt bridges in the model metalloprotease PAE and provides important insights into the structure-function relationship of enzymes.

Determination of oligomeric states of peptide complexes using thermal unfolding curves.

In their native form peptides are often found as oligomeric complexes, meaning they consist of more than one peptide chain. Coiled coils and helical bundles are common examples of such complexes. Their oligomeric state needs to be known precisely as this tremendously influences their biochemical and biophysical properties. The extensive analysis of circular dichroism spectroscopic data is commonly used to investigate the thermodynamics of binding and folding of these complexes. Here we present FitDis! an easy-to-use programme, which fits the most common two-state unfolding transition to the measured thermal unfolding curves of any oligomer of any stoichiometry. We demonstrate, with simulated and real examples, that the comparison of different stoichiometric models fitted to the same dataset reveals the oligomeric states of these complexes along with detailed thermodynamic information. This method will significantly ease the analysis of and increase the amount of information gained from, the thermal unfolding curves of peptide complexes.

Factor defining the effects of glycine betaine on the thermodynamic stability and internal dynamics of horse cytochrome C.

A compatible osmolyte such as glycine betaine (GB) and low concentrations of a denaturant constrain the internal dynamics of natively folded carbonmonoxycytochrome c (NCO) at pH 7.0. GB and subdenaturing concentrations of guanidine hydrochloride (GdnHCl) or urea have a cumulative effect on the constrained dynamics of NCO. At higher denaturant concentrations, large-scale unfolding fluctuations dominate the dynamics and inclusion of GB opposes the structural fluctuations that cause unfolding of the protein. These deductions are made from kinetic and thermodynamic parameters measured for the CO dissociation reaction of NCO at varying concentrations of denaturant in the absence and presence of 1.0 M GB. Intermolecular docking between horse ferrocytochrome c and a denaturant or GB reveals that the denaturant-mediated constrained dynamics of the protein is due to polyfunctional interactions between the denaturant and different groups of protein while the GB-mediated restricted dynamics of the protein arises from both the direct interactions of GB with different side chains of Lys or Arg residues of the protein and indirect interactions of GB with the protein surface. Thermodynamic analysis of the thermal and GdnHCl-induced unfolding curves of ferrocytochrome c measured in the absence and presence of 1.0 M GB at pH 7.0 indicates that GB increases the thermodynamic stability of ferrocytochrome c at neutral pH. Analysis of thermal and urea-induced unfolding curves of ferricytochrome c measured at different GdnHCl concentrations in the absence and presence of 0.5-1.0 M GB at pH 7.0 and 3.8 suggests that GB counteracts the destabilizing effect of the denaturant at pH 7.0 but exhibits an additive effect on the destabilizing effect of the denaturant at pH 3.8.

Role of the Cytoplasmic N-terminal Cap and Per-Arnt-Sim (PAS) Domain in Trafficking and Stabilization of Kv11.1 Channels

The N-terminal cytoplasmic region of the Kv11.1a potassium channel contains a Per-Arnt-Sim (PAS) domain that is essential for the unique slow deactivation gating kinetics of the channel. The PAS domain has also been implicated in the assembly and stabilization of the assembled tetrameric channel, with many clinical mutants in the PAS domain resulting in reduced stability of the domain and reduced trafficking. Here, we use quantitative Western blotting to show that the PAS domain is not required for normal channel trafficking nor for subunit-subunit interactions, and it is not necessary for stabilizing assembled channels. However, when the PAS domain is present, the N-Cap amphipathic helix must also be present for channels to traffic to the cell membrane. Serine scan mutagenesis of the N-Cap amphipathic helix identified Leu-15, Ile-18, and Ile-19 as residues critical for the stabilization of full-length proteins when the PAS domain is present. Furthermore, mutant cycle analysis experiments support recent crystallography studies, indicating that the hydrophobic face of the N-Cap amphipathic helix interacts with a surface-exposed hydrophobic patch on the core of the PAS domain to stabilize the structure of this critical gating domain. Our data demonstrate that the N-Cap amphipathic helix is critical for channel stability and trafficking.

Ca 2 + and Mg 2 + binding induce conformational stability of Calfumirin-1 from Dictyostelium discoideum

The apo-Calfumirin-1 (CAF-1) binds to Ca2+ with high affinity and also to Mg2+ with high positive cooperativity. The thermal unfolding curves of wtCAF-1 monitored at neutral pH by CD spectroscopy are reversible and show different thermal stabilities in the absence or presence of Ca2+ and Mg2+ ions. Metal-free wtCAF-1 shows greater thermal stability than EF-IV mutant protein. We observed that GdnHCl-induced unfolding of apo-wtCAF-1 monitored by CD and fluorescence spectroscopies increases co-operative folding with approximately same Cm values. Binding of Ca2+ and Mg2+ ions to CAF-1 dramatically altered the fluorescence and CD spectra, indicating metal ion-induced conformational changes both in the wild-type and mutant proteins. The hydrophobic probe, ANS is used to observe alteration in surface hydrophobicity of the protein in different ligation states. In apo-wtCAF-1, the exposed hydrophobic surfaces are able to bind ANS which is in contrast to the unfolded or the metal ions ligated conformations. Isothermal titration calorimetry (ITC) results show two possible independent binding sites of comparable affinity for the metal ions. However, their binding to the EF-IV E helix-loop-F helix mutant apo-protein happens with different affinities. The present study demonstrates that Ca2+ or Mg2+ binding plays a possible role in the conformational stability of the protein. This study aims to understand the mechanisms of Ca2+ and Mg2+ binding in CAF-1 using biophysical methods and correlate its improvement in D. discoideum cell differentiation. Results clearly indicate that wtCAF-1 is more stable than the mutant and retains its activity in Ca2+ bound form.

Thermodynamic analysis of osmolyte effect on thermal stability of ribonuclease A in terms of water activity

Thermal unfolding of ribonuclease A (RNase) was analyzed in various osmolyte solutions of glycine, proline, sarcosine, N,N-dimethylglycine, betaine, myo-inositol, taurine, and trimethylamine-N-oxide (TMAO). All the osmolytes tested stabilized the protein. The thermal unfolding curve was described well by the van't Hoff equation and the melting temperature and the enthalpy of protein unfolding were obtained. The Wyman–Tanford equation, which describes the unfolded-to-folded protein ratio as a function of water activity, was successfully applied to obtain a linear plot. In consideration of this experimentally obtained linearity, the Wyman–Tanford plot could be integrated to calculate the stabilization free energy of the protein (∆∆G) in the solution. The ∆∆G was proved to be described by the property of the microstructure around the protein surface, which is composed of the protein hydration, the cosolute-binding, and the preferential exclusion, and the property of the bulk solution; water activity. The m-values of osmolytes for protein unfolding were obtained from ∆∆G calculated. Among the osmolytes tested, myo-inositol showed the highest m-value.

Determining the Conformational Stability of a Protein from Urea and Thermal Unfolding Curves

This unit contains basic protocols for determining the conformational stability of a globular protein from either urea or thermal unfolding curves. Circular dichroism is the optical spectroscopic technique most commonly used to monitor protein unfolding. The protocols describe how to analyze data from an unfolding curve to obtain the thermodynamic parameters necessary to calculate conformational stability, and how to determine differences in stability between protein variants.

A Survey of Antiprion Compounds Reveals the Prevalence of Non-PrP Molecular Targets

Prion diseases are fatal neurodegenerative diseases caused by the accumulation of the misfolded isoform (PrP(Sc)) of the prion protein (PrP(C)). Cell-based screens have identified several compounds that induce a reduction in PrP(Sc) levels in infected cultured cells. However, the molecular targets of most antiprion compounds remain unknown. We undertook a large-scale, unbiased, cell-based screen for antiprion compounds and then investigated whether a representative subset of the active molecules had measurable affinity for PrP, increased the susceptibility of PrP(Sc) to proteolysis, or altered the cellular localization or expression level of PrP(C). None of the antiprion compounds showed in vitro affinity for PrP or had the ability to disaggregate PrP(Sc) in infected brain homogenates. These observations suggest that most antiprion compounds identified in cell-based screens deploy their activity via non-PrP targets in the cell. Our findings indicate that in comparison to PrP conformers themselves, proteins that play auxiliary roles in prion propagation may be more effective targets for future drug discovery efforts.

Structural and Functional Relationships between the Lectin and Arm Domains of Calreticulin

Calreticulin and calnexin are key components in maintaining the quality control of glycoprotein folding within the endoplasmic reticulum. Although their lectin function of binding monoglucosylated sugar moieties of glycoproteins is well documented, their chaperone activity in suppressing protein aggregation is less well understood. Here, we use a series of deletion mutants of calreticulin to demonstrate that its aggregation suppression function resides primarily within its lectin domain. Using hydrophobic peptides as substrate mimetics, we show that aggregation suppression is mediated through a single polypeptide binding site that exhibits a K(d) for peptides of 0.5-1 μM. This site is distinct from the oligosaccharide binding site and differs from previously identified sites of binding to thrombospondin and GABARAP (4-aminobutyrate type A receptor-associated protein). Although the arm domain of calreticulin was incapable of suppressing aggregation or binding hydrophobic peptides on its own, it did contribute to aggregation suppression in the context of the whole molecule. The high resolution x-ray crystal structure of calreticulin with a partially truncated arm domain reveals a marked difference in the relative orientations of the arm and lectin domains when compared with calnexin. Furthermore, a hydrophobic patch was detected on the arm domain that mediates crystal packing and may contribute to calreticulin chaperone function.

When a module is not a domain: the case of the REJ module and the redefinition of the architecture of polycystin-1.

The extracellular region of a group of cell-surface receptors known as the polycystic kidney disease 1 family, containing, among others, polycystin-1, has been controversially described as containing four FNIII (fibronectin typeIII) domains or one REJ (receptor of egg jelly protein) module in the same portion of polypeptide. Stimulated by recent atomic force microscopy work, we re-examined the similarity of these four domains with a FNIII sequence profile showing the evolutionary relationship. Two of the predicted domains could be expressed in bacteria and refolded to give a protein suitable for biophysical study, and one of these expressed solubly. CD spectroscopy showed that both domains contain a significant amount of β-sheet, in good agreement with theoretical predictions. Confirmation of independent folding as a domain is obtained from highly co-operative thermal and urea unfolding curves. Excellent dispersion of peaks in the high-field region of one-dimensional NMR spectra confirms the presence of a hydrophobic core. Analytical ultracentrifugation and analytical gel filtration agree very well with the narrow linewidths in the NMR spectra that at least one of the domains is monomeric. On the basis of this combined theoretical and experimental analysis, we show that the extracellular portion of polycystin-1 does indeed contain β-sheet domains, probably FNIII, and that, consequently, the REJ module is not a single domain.

A Critical Assessment of Putative Gatekeeper Interactions in the Villin Headpiece Helical Subdomain

The helical subdomain of the villin headpiece (HP36) is one of the smallest naturally occurring proteins that folds cooperatively. Its small size, rapid folding, and simple three-helix topology have made it an extraordinary popular model system for computational, theoretical, and experimental studies of protein folding. Aromatic–proline interactions involving Trp64 and Pro62 have been proposed to play a critical role in specifying the subdomain fold by acting as gatekeeper residues. Note that the numbering corresponds to full-length headpiece. Mutation of Pro62 has been shown to lead to a protein that does not fold, but this may arise for two different reasons: The residue may make interactions that are critical for the specificity of the fold or the mutation may simply destabilize the domain. In the first case, the protein cannot fold, while in the second, the small fraction of molecules that do fold adopt the correct structure. The modest stability of the wild type prevents a critical analysis of these interactions because even moderately destabilizing mutations lead to a very small folded state population. Using a hyperstable variant of HP36, denoted DM HP36, as our new wild type, we characterized a set of mutants designed to assess the role of the putative gatekeeper interactions. Four single mutants, DM Pro62Ala, DM Trp64Leu, DM Trp64Lys, and DM Trp64Ala, and a double mutant, DM Pro62Ala Trp64Leu, were prepared. All mutants are less stable than DM HP36, but all are well folded as judged by CD and 1H NMR. All of the mutants display sigmoidal thermal unfolding and urea-induced unfolding curves. Double-mutant cycle analysis shows that the interactions between Pro62 and Trp64 are weak but favorable. Interactions involving Pro62 and proline–aromatic interactions are, thus, not required for specifying the subdomain fold. The implications for the design and thermodynamics of miniature proteins are discussed.

α-Hemoglobin Stabilizing Protein (AHSP), a Kinetic Scheme of the Action of a Human Mutant, AHSPV56G

A kinetic analysis has been made of the interaction of alpha-Hb chains with a mutant alpha-hemoglobin stabilizing protein, AHSP(V56G), which is the first case of an AHSP mutation associated with clinical symptoms of mild thalassemia syndrome. The chaperone AHSP is thought to protect nascent alpha chains until final binding to the partner beta-Hb. Rather than protecting alpha chains, the mutant chaperone is partially unfolded but recovers its secondary structure via interaction with alpha-Hb. For both AHSP(WT) and AHSP(V56G), the binding to alpha-Hb is quite rapid relative to the alpha-beta reaction, as expected because the chaperone binding must be quite competitive to complete the alpha chain folding process before alpha-Hb binds irreversibly to beta-Hb. The main kinetic difference is a dissociation rate of AHSP(V56G).alpha-Hb some four times faster relative to AHSP.alpha-Hb. Considering a role of protein folding, the AHSP(V56G) apparently does not bind long enough (0.5 s versus 2 s for the WT) to complete the structural modifications. The overall replacement reaction (AHSP.alpha-Hb + beta-Hb --> AHSP + alphabeta) can be quite long, especially if there is an excess of AHSP relative to beta-Hb monomers.