Thermodynamics of protein denatured states

In Situ Observation of Chemically Induced Protein Denaturation at Solvated Interfaces

Proteins unfold in chaotropic salt solutions, a process that is difficult to observe at the single protein level. The work presented here demonstrates that a liquid-based atomic force microscope and graphene liquid-cell-based scanning transmission electron microscope make it possible to observe chemically induced protein unfolding. To illustrate this capability, ferritin proteins were deposited on a graphene surface, and the concentration-dependent urea- or guanidinium-induced changes of morphology were monitored for holo-ferritin with its ferrihydrite core as well as apo-ferritin without this core. Depending on the chaotropic agent the liquid-based imaging setup captured an unexpected transformation of natively folded holo-ferritin proteins into rings after urea treatment but not after guanidinium treatment. Urea treatment of apo-ferritin did not result in nanorings, confirming that nanorings are a specific signature of denaturation of holo-ferritins after exposture to sufficiently high urea concentrations. Mapping the in situ images with molecular dynamics simulations of ferritin subunits in urea solutions suggests that electrostatic destabilization triggers denaturation of ferritin as urea makes direct contact with the protein and also disrupts the water H-bonding network in the ferritin solvation shell. Our findings deepen the understanding of protein denaturation studied using label-free techniques operating at the solid-liquid interface.

Influence of the Interdomain Interface on Structural and Redox Properties of Multiheme Proteins

Multiheme proteins are important in energy conversion and biogeochemical cycles of nitrogen and sulfur. A diheme cytochrome c4 (c4) was used as a model to elucidate roles of the interdomain interface on properties of iron centers in its hemes A and B. Isolated monoheme domains c4-A and c4-B, together with the full-length diheme c4 and its Met-to-His ligand variants, were characterized by a variety of spectroscopic and stability measurements. In both isolated domains, the heme iron is Met/His-ligated at pH 5.0, as in the full-length c4, but becomes His/His-ligated in c4-B at higher pH. Intradomain contacts in c4-A are minimally affected by the separation of c4-A and c4-B domains, and isolated c4-A is folded. In contrast, the isolated c4-B is partially unfolded, and the interface with c4-A guides folding of this domain. The c4-A and c4-B domains have the propensity to interact even without the polypeptide linker. Thermodynamic cycles have revealed properties of monomeric folded isolated domains, suggesting that ferrous (FeII), but not ferric (FeIII) c4-A and c4-B, is stabilized by the interface. This study illustrates the effects of the interface on tuning structural and redox properties of multiheme proteins and enriches our understanding of redox-dependent complexation.

Denatured State Conformational Biases in Three-Helix Bundles Containing Divergent Sequences Localize near Turns and Helix Capping Residues

Rhodopseudomonas palustris cytochrome c', a four-helix bundle, and the second ubiquitin-associated domain, UBA(2), a three-helix bundle from the human homologue of yeast Rad23, HHR23A, deviate from random coil behavior under denaturing conditions in a fold-specific manner. The random coil deviations in each of these folds occur near interhelical turns and loops in their tertiary structures. Here, we examine an additional three-helix bundle with an identical fold to UBA(2), but a highly divergent sequence, the first ubiquitin-associated domain, UBA(1), of HHR23A. We use histidine-heme loop formation methods, employing eight single histidine variants, to probe for denatured state conformational bias of a UBA(1) domain fused to the N-terminus of iso-1-cytochrome c (iso-1-Cytc). Guanidine hydrochloride (GuHCl) denaturation shows that the iso-1-Cytc domain unfolds first, followed by the UBA(1) domain. Denatured state (4 and 6 M GuHCl) histidine-heme loop formation studies show that as the size of the histidine-heme loop increases, loop stability decreases, as expected for the Jacobson-Stockmayer relationship. However, loops formed with His35, His31, and His15, of UBA(1), are 0.6-1.1 kcal/mol more stable than expected from the Jacobson-Stockmayer relationship, confirming the importance of deviations of the denatured state from random coil behavior near interhelical turns of helical domains for facilitating folding to the correct topology. For UBA(1) and UBA(2), hydrophobic clusters on either side of the turns partially explain deviations from random coil behavior; however, helix capping also appears to be important.

A pH-dependent protein stability switch coupled to the perturbed pKa of a single ionizable residue

The contribution of electrostatic interactions in protein stability has not been fully understood. Burial of an ionizable amino acid inside the hydrophobic protein core can affect its ionization equilibrium and shift its pKa differentially in the native (N) and unfolded (U) states of a protein and this coupling between the folding/unfolding cycle and the ionization equilibria of the ionizable residue can substantially influence the protein stability. Here, we studied the coupling of the folding/unfolding cycle with the ionization of a buried ionizable residue in a multi-domain protein, Human Serum Albumin (HSA) using fluorescence spectroscopy. A pH-dependent change in the stability of HSA was observed in the near native pH range (pH 6.0-9.0). The protonation-deprotonation equilibrium of a single thiol residue that is buried in the protein structure was identified to give rise to the pH-dependent protein stability. We quantified the pKa of the thiol residue in the N and the U states. The mean pKa of the thiol in the N state was upshifted by 0.5 units to 8.7 due to the burial of the thiol in the protein structure. Surprisingly, the mean pKa of the thiol in the U state was observed to be downshifted by 1.3 units to 6.9. These results indicate that some charged residues are spatially proximal to the thiol group in the U state. Our results suggest that, in addition to the N state, electrostatic interactions in the U state are important determinants of protein stability.

Histological and Biomechanical Characteristics of Human Decellularized Allograft Heart Valves After Eighteen Months of Storage in Saline Solution

Background: The best storage preservation method for maintaining the quality and safety of human decellularized allograft heart valves is yet to be established. Objective: The aim of the present study was to evaluate the stability in terms of extracellular matrix (ECM) integrity of human heart valve allografts decellularized using sodium dodecyl sulfate-ethylenediaminetetraacetic acid (SDS-EDTA) and stored for 6, 12, and 18 months. Methods: A total of 70 decellularized aortic and pulmonary valves were analyzed across different storage times (0, 6, 12, and 18 months) for solution pH measurements, histological findings, cytotoxicity assay results, biomechanical test results, and microbiological suitability test results. Continuous data were analyzed using one-way analysis of variance comparing the follow-up times. Results: The pH of the stock solution did not change during the different time points, and no microbial growth occurred up to 18 months. Histological analysis showed that the decellularized allografts did not present deleterious outcomes or signs of structural degeneration in the ECM up to 12 months. The biomechanical properties showed changes over time in different aspects. Allografts stored for 18 months presented lower tensile strength and elasticity than those stored for 12 months (p < 0.05). The microbiological suitability test suggested no residual antimicrobial effects. Conclusion: Changes in the structure and functionality of SDS-EDTA decellularized heart valve allografts occur after 12 months of storage.

Two-step self-assembly of a spider silk molecular clamp

Web spiders synthesize silk fibers of unique strength and extensibility through the controlled self-assembly of protein building blocks, so-called spidroins. The spidroin C-terminal domain is highly conserved and connects two polypeptide chains through formation of an all-helical, intertwined dimer. Here we use contact-induced fluorescence self-quenching and resonance energy transfer in combination with far-UV circular dichroism spectroscopy as three orthogonal structural probes to dissect the mechanism of folding and dimerization of a spidroin C-terminal domain from the major ampullate gland of the nursery web spider Euprosthenops australis. We show that helices forming the dimer core assemble very rapidly and fold on association. Subsequently, peripheral helices fold and dock slowly onto the preformed core. Lability of outer helices facilitates formation of a highly expanded, partially folded dimer. The high end-to-end distance of chain termini in the partially folded dimer suggests an extensibility module that contributes to elasticity of spider silk.

Probing Denatured State Conformational Bias in a Three-Helix Bundle, UBA(2), Using a Cytochrome c Fusion Protein

Previous work with the four-helix-bundle protein cytochrome c' from Rhodopseudomonas palustris using histidine-heme loop formation methods revealed fold-specific deviations from random coil behavior in its denatured state ensemble. To examine the generality of this finding, we extend this work to a three-helix-bundle polypeptide, the second ubiquitin-associated domain, UBA(2), of the human DNA excision repair protein. We use yeast iso-1-cytochrome c as a scaffold, fusing the UBA(2) domain at the N-terminus of iso-1-cytochrome c. We have engineered histidine into highly solvent accessible positions of UBA(2), creating six single histidine variants. Guanidine hydrochloride denaturation studies show that the UBA(2)-cytochrome c fusion protein unfolds in a three-state process with iso-1-cytochrome c unfolding first. Furthermore, engineered histidine residues in UBA(2) strongly destabilize the iso-1-cytochrome c domain. Equilibrium and kinetic histidine-heme loop formation measurements in the denatured state at 4 and 6 M guanidine hydrochloride show that loop stability decreases as the size of the histidine-heme loop increases, in accord with the Jacobson-Stockmayer equation. However, we observe that the His27-heme loop is both more stable than expected from the Jacobson-Stockmayer relationship and breaks more slowly than expected. These results show that the sequence near His27, which is in the reverse turn between helices 2 and 3 of UBA(2), is prone to persistent interactions in the denatured state. Therefore, consistent with our results for cytochrome c', this reverse turn sequence may help to establish the topology of this fold by biasing the conformational distribution of the denatured state.

Helical Propensity Affects the Conformational Properties of the Denatured State of Cytochrome c'

Changing the helical propensity of a polypeptide sequence might be expected to affect the conformational properties of the denatured state of a protein. To test this hypothesis, alanines at positions 83 and 87 near the center of helix 3 of cytochrome c' from Rhodopseudomonas palustris were mutated to serine to decrease the stability of this helix. A set of 13 single histidine variants in the A83S/A87S background were prepared to permit assessment of the conformational properties of the denatured state using histidine-loop formation in 3 M guanidine hydrochloride. The data are compared with previous histidine-heme loop formation data for wild-type cytochrome c'. As expected, destabilization of helix 3 decreases the global stabilities of the histidine variants in the A83S/A87S background relative to the wild-type background. Loop stability versus loop size data yields a scaling exponent of 2.1 ± 0.2, similar to the value of 2.3 ± 0.2 obtained for wild-type cytochrome c'. However, the stabilities of all histidine-heme loops, which contain the helix 3 sequence segment, are increased in the A83S/A87S background compared to the wild-type background. Rate constants for histidine-heme loop breakage are similar for the wild-type and A83S/A87S variants. However, for histidine-heme loops that contain the helix 3 sequence segment, the rate constants for loop formation increase in the A83S/A87S background compared to the wild-type background. Thus, residual helical structure appears to stiffen the polypeptide chain slowing loop formation in the denatured state. The implications of these results for protein folding mechanisms are discussed.

Identification and characterization of a cyclosporin binding cyclophilin from Staphylococcus aureus Newman

Cyclophilins, a class of peptidyl-prolyl cis-trans isomerase (PPIase) enzymes, are inhibited by cyclosporin A (CsA), an immunosuppressive drug. Staphylococcus aureus Newman, a pathogenic bacterium, carries a gene for encoding a putative cyclophilin (SaCyp). SaCyp shows significant homology with other cyclophilins at the sequence level. A three-dimensional model structure of SaCyp harbors a binding site for CsA. To verify whether SaCyp possesses both the PPIase activity and the CsA binding ability, we have purified and investigated a recombinant SaCyp (rCyp) using various in vitro tools. Our RNase T1 refolding assay indicates that rCyp has a substantial extent of PPIase activity. rCyp that exists as a monomer in the aqueous solution is truly a cyclophilin as its catalytic activity specifically shows sensitivity to CsA. rCyp appears to bind CsA with a reasonably high affinity. Additional investigations reveal that binding of CsA to rCyp alters its structure and shape to some extent. Both rCyp and rCyp-CsA are unfolded via the formation of at least one intermediate in the presence of guanidine hydrochloride. Unfolding study also indicates that there is substantial extent of thermodynamic stabilization of rCyp in the presence of CsA as well. The data suggest that rCyp may be exploited to screen the new antimicrobial agents in the future.

Molecular determinant of the effects of hydrostatic pressure on protein folding stability

Hydrostatic pressure is an important environmental variable that plays an essential role in biological adaptation for many extremophilic organisms (for example, piezophiles). Increase in hydrostatic pressure, much like increase in temperature, perturbs the thermodynamic equilibrium between native and unfolded states of proteins. Experimentally, it has been observed that increase in hydrostatic pressure can both increase and decrease protein stability. These observations suggest that volume changes upon protein unfolding can be both positive and negative. The molecular details of this difference in sign of volume changes have been puzzling the field for the past 50 years. Here we present a comprehensive thermodynamic model that provides in-depth analysis of the contribution of various molecular determinants to the volume changes upon protein unfolding. Comparison with experimental data shows that the model allows quantitative predictions of volume changes upon protein unfolding, thus paving the way to proteome-wide computational comparison of proteins from different extremophilic organisms.

Proteins can be both stabilized and destabilized by pressure. Here the authors analyse the factors contributing to both negative and positive protein volume change upon denaturation, and shed light on the molecular determinants allowing proteins to be stable at high pressures.

Effects of Visible-Light Irradiation of Protoporphyrin IX on the Self-Assembly of Tubulin Heterodimers

The formation and the effects of laser irradiation of the complex formed by protoporphyrin IX (PPIX) and tubulin was investigated. We have used tubulin as a model protein to investigate whether docked photoactive ligands can affect the structure and function of polypeptides upon exposure to visible light. We observed that laser irradiation in the Soret band prompts bleaching of the PPIX, which is accompanied by a sharp decrease in the intensity and average fluorescence lifetime of the protein (dominated by the four tryptophan residues of the tubulin monomer). The kinetics indicate non-trivial effects and suggest that the photosensitization of the PPIX bound to tubulin prompts structural alterations of the protein. These modifications were also observed through changes in the energy transfer between Trp residues and PPIX. The results suggest that laser irradiation produces localized partial unfolding of tubulin and that the changes prompt modification of the formation of microtubules in vitro. Measurements of singlet oxygen formation were inconclusive in determining whether the changes are prompted by reactive oxygen species or other excited state mechanisms.

Inhibitor-induced conformational stabilization and structural alteration of a mip-like peptidyl prolyl cis-trans isomerase and its C-terminal domain

FKBP22, an Escherichia coli-encoded PPIase (peptidyl-prolyl cis-trans isomerase) enzyme, shares substantial identity with the Mip-like pathogenic factors, caries two domains, exists as a dimer in solution and binds some immunosuppressive drugs (such as FK506 and rapamycin) using its C-terminal domain (CTD). To understand the effects of these drugs on the structure and stability of the Mip-like proteins, rFKBP22 (a chimeric FKBP22) and CTD⁺ (a CTD variant) have been studied in the presence and absence of rapamycin using different probes. We demonstrated that rapamycin binding causes minor structural alterations of rFKBP22 and CTD⁺. Both the proteins (equilibrated with rapamycin) were unfolded via the formation of intermediates in the presence of urea. Further study revealed that thermal unfolding of both rFKBP22 and rapamycin-saturated rFKBP22 occurred by a three-state mechanism with the synthesis of intermediates. Intermediate from the rapamycin-equilibrated rFKBP22 was formed at a comparatively higher temperature. All intermediates carried substantial extents of secondary and tertiary structures. Intermediate resulted from the thermal unfolding of rFKBP22 existed as the dimers in solution, carried an increased extent of hydrophobic surface and possessed relatively higher rapamycin binding activity. Despite the formation of intermediates, both the thermal and urea-induced unfolding reactions were reversible in nature. Unfolding studies also indicated the considerable stabilization of both proteins by rapamycin binding. The data suggest that rFKBP22 or CTD⁺ could be exploited to screen the rapamycin-like inhibitors in the future.

Denatured state ensembles with the same radii of gyration can form significantly different long-range contacts

Defining the structural, dynamic, and energetic properties of the unfolded state of proteins is critical for an in-depth understanding of protein folding, protein thermodynamics, and protein aggregation. Here we analyze long-range contacts and compactness in two apparently fully unfolded ensembles of the same protein: the acid unfolded state of the C-terminal domain of ribosomal protein L9 in the absence of high concentrations of urea as well as the urea unfolded state at low pH. Small angle X-ray scattering reveals that the two states are expanded with values of Rg differing by <7%. Paramagnetic relaxation enhancement (PRE) nuclear magnetic resonance studies, however, reveal that the acid unfolded state samples conformations that facilitate contacts between residues that are distant in sequence while the urea unfolded state ensemble does not. The experimental PRE profiles for the acid unfolded state differ significantly from these predicted using an excluded volume limit ensemble, but these long-range contacts are largely eliminated by the addition of 8 M urea. The work shows that expanded unfolded states can sample very different distributions of long-range contacts yet still have similar radii of gyration. The implications for protein folding and for the characterization of unfolded states are discussed.

Secondary structure, a missing component of sequence-based minimotif definitions

Minimotifs are short contiguous segments of proteins that have a known biological function. The hundreds of thousands of minimotifs discovered thus far are an important part of the theoretical understanding of the specificity of protein-protein interactions, posttranslational modifications, and signal transduction that occur in cells. However, a longstanding problem is that the different abstractions of the sequence definitions do not accurately capture the specificity, despite decades of effort by many labs. We present evidence that structure is an essential component of minimotif specificity, yet is not used in minimotif definitions. Our analysis of several known minimotifs as case studies, analysis of occurrences of minimotifs in structured and disordered regions of proteins, and review of the literature support a new model for minimotif definitions that includes sequence, structure, and function.

Tryptophan stabilizes His-heme loops in the denatured state only when it is near a loop end

We use a host-guest approach to evaluate the effect of Trp guest residues relative to Ala on the kinetics and thermodynamics of formation of His-heme loops in the denatured state of iso-1-cytochrome c at 1.5, 3.0, and 6.0 M guanidine hydrochloride (GdnHCl). Trp guest residues are inserted into an alanine-rich segment placed after a unique His near the N-terminus of iso-1-cytochrome c. Trp guest residues are either 4 or 10 residues from the His end of the 28-residue loop. We find the guest Trp stabilizes the His-heme loop at all GdnHCl concentrations when it is the 4th, but not the 10th, residue from the His end of the loop. Thus, residues near loop ends are most important in developing topological constraints in the denatured state that affect protein folding. In 1.5 M GdnHCl, the loop stabilization is ~0.7 kcal/mol, providing a thermodynamic rationale for the observation that Trp often mediates residual structure in the denatured state. Measurement of loop breakage rate constants, k(b,His), indicates that loop stabilization by the Trp guest residues occurs completely after the transition state for loop formation in 6.0 M GdnHCl. Under poorer solvent conditions, approximately half of the stabilization of the loop develops in the transition state, consistent with contacts in the denatured state being energetically downhill and providing evidence for funneling even near the rim of the folding funnel.

Highly perturbed pKa values in the unfolded state of hen egg white lysozyme

The majority of pK(a) values in protein unfolded states are close to the amino acid model pK(a) values, thus reflecting the weak intramolecular interactions present in the unfolded ensemble of most proteins. We have carried out thermal denaturation measurements on the WT and eight mutants of HEWL from pH 1.5 to pH 11.0 to examine the unfolded state pK(a) values and the pH dependence of protein stability for this enzyme. The availability of accurate pK(a) values for the folded state of HEWL and separate measurements of mutant-induced effects on the folded state pK(a) values, allows us to estimate the pK(a) values of seven acidic residues in the unfolded state of HEWL. Asp-48 and Asp-66 display pK(a) values of 2.9 and 3.1 in our analysis, thus representing the most depressed unfolded state pK(a) values observed to date. We observe a strong correlation between the folded state pK(a) values and the unfolded state pK(a) values of HEWL, thus suggesting that the unfolded state of HEWL possesses a large degree of native state characteristics.

Residual structure in unfolded proteins

The denatured state ensemble (DSE) of unfolded proteins, once considered to be well-modeled by an energetically featureless random coil, is now well-known to contain flickering elements of residual structure. The position and nature of DSE residual structure may provide clues toward deciphering the protein folding code. This review focuses on recent advances in our understanding of the nature of DSE collapse under folding conditions, the quantification of the stability of residual structure in the DSE, the determination of the location and types of residues involved in thermodynamically significant residual structure and advances in detection of long-range interactions in the DSE.

Analysis of electrostatic interactions in the denatured state ensemble of the N-terminal domain of L9 under native conditions

The pH dependence of protein stability is defined by the difference in the number of protons bound to the folded state and to the denatured state ensemble (DSE) as a function of pH. In many cases, the protonation behavior can be described as the sum of a set of independently titrating residues; in this case, the pH dependence of stability reflects differences in folded and DSE pK(a)'s. pH dependent stability studies have shown that there are energetically important interactions involving charged residues in the DSE of the N-terminal domain of L9 (NTL9), which affect significantly the stability of the protein. The DSE of wild type NTL9 cannot be directly characterized under native conditions because of its high stability. A destabilized double mutant of NTL9, V3AI4A, significantly populates the folded state and the DSE in the absence of denaturant. The two states are in slow exchange on the nuclear magnetic resonance time scale, and diffusion measurements indicate that the DSE is compact. The DSE pK(a)'s of all of the acidic residues were directly determined. The DSE pK(a) of Asp8 and Asp23 are depressed relative to model compounds values. Use of the mutant DSE pK(a)'s together with known native state pK(a)'s leads to a significantly improved agreement between the measured pH dependent stability and that predicted by the Tanford-Wyman linkage relationship. An analysis of the literature suggests that DSE interactions involving charged residues are relatively common and should be considered in discussions of protein stability.

How do thermophilic proteins and proteomes withstand high temperature?

We attempt to understand the origin of enhanced stability in thermophilic proteins by analyzing thermodynamic data for 116 proteins, the largest data set achieved to date. We compute changes in entropy and enthalpy at the convergence temperature where different driving forces are maximally decoupled, in contrast to the majority of previous studies that were performed at the melting temperature. We find, on average, that the gain in enthalpy upon folding is lower in thermophiles than in mesophiles, whereas the loss in entropy upon folding is higher in mesophiles than in thermophiles. This implies that entropic stabilization may be responsible for the high melting temperature, and hints at residual structure or compactness of the denatured state in thermophiles. We find a similar trend by analyzing a homologous set of proteins classified based only on the optimum growth temperature of the organisms from which they were extracted. We find that the folding free energy at the temperature of maximal stability is significantly more favorable in thermophiles than in mesophiles, whereas the maximal stability temperature itself is similar between these two classes. Furthermore, we extend the thermodynamic analysis to model the entire proteome. The results explain the high optimal growth temperature in thermophilic organisms and are in excellent quantitative agreement with full thermal growth rate data obtained in a dozen thermophilic and mesophilic organisms.

Manifestations of native topology in the denatured state ensemble of Rhodopseudomonas palustris cytochrome c'

To provide insight into the role of local sequence in the nonrandom coil behavior of the denatured state, we have extended our measurements of histidine-heme loop formation equilibria for cytochrome c' to 6 M guanidine hydrochloride. We observe that there is some reduction in the scatter about the best fit line of loop stability versus loop size data in 6 M versus 3 M guanidine hydrochloride, but the scatter is not eliminated. The scaling exponent, ν(3), of 2.5 ± 0.2 is also similar to that found previously in 3 M guanidine hydrochloride (2.6 ± 0.3). Rates of histidine-heme loop breakage in the denatured state of cytochrome c' show that some histidine-heme loops are significantly more persistent than others at both 3 and 6 M guanidine hydrochloride. Rates of histidine-heme loop formation more closely approximate random coil behavior. This observation indicates that heterogeneity in the denatured state ensemble results mainly from contact persistence. When mapped onto the structure of cytochrome c', the histidine-heme loops with slow breakage rates coincide with chain reversals between helices 1 and 2 and between helices 2 and 3. Molecular dynamics simulations of the unfolding of cytochrome c' at 498 K show that these reverse turns persist in the unfolded state. Thus, these portions of the primary structure of cytochrome c' set up the topology of cytochrome c' in the denatured state, predisposing the protein to fold efficiently to its native structure.

Identification of a novel functional domain of ricin responsible for its potent toxicity

Ribosome-inactivating proteins (RIPs) are toxic N-glycosidases that depurinate the universally conserved α-sarcin loop of large rRNAs. They have received attention in biological and biomedical research because of their unique biological activities toward animals and human cells as cell-killing agents. A better understanding of the depurination mechanism of RIPs could allow us to develop potent neutralizing antibodies and to design efficient immunotoxins for clinical use. Among these RIPs, ricin exhibited remarkable efficacy in depurination activity and highly conserved tertiary structure with other RIPs. It can be considered as a prototype to investigate the depurination mechanism of RIPs. In the present study, we successfully identified a novel functional domain responsible for controlling the depurination activity of ricin, which is located far from the enzymatic active site reported previously. Our study indicated that ricin A-chain mAbs binding to this domain (an α-helix comprising the residues 99-106) exhibited an unusual potent neutralizing ability against ricin in vivo. To further investigate the potential role of the α-helix in regulating the catalytic activity of ricin, ricin A-chain variants with different flexibility of the α-helix were rationally designed. Our data clearly demonstrated that the flexibility of the α-helix is responsible for controlling the depurination activity of ricin and determining the extent of protein synthesis inhibition, suggesting that the conserved α-helix might be considered as a potential target for the prevention and treatment of RIP poisoning.

Propensities of aromatic amino acids versus leucine and proline to induce residual structure in the denatured-state ensemble of iso-1-cytochrome c

Histidine-heme loop formation in the denatured state of a protein is a sensitive means for probing residual structure under unfolding conditions. In this study, we use a host-guest approach to investigate the relative tendencies of different amino acids to promote residual structure under denaturing conditions. The host for this work is a 6-amino-acid insert of five alanines, followed by a lysine engineered immediately following a unique histidine near the N-terminus of yeast iso-1-cytochrome c. We substitute the fourth alanine in this sequence HAAAXAK (with X=Trp, Phe, Tyr, and Leu). The effects of proline are tested with substitutions at positions 1 and 5 in the insert (HPAAAAK and HAAAAPK, respectively). Thermodynamic studies on His-heme loop formation in 3 M guanidine hydrochloride reveal significant stabilization of residual structure by aromatic amino acids, particularly Trp and Phe, and minimal stabilization of residual structure by Leu. Prolines slightly disfavor His-heme loop formation, presumably due to enhanced chain stiffness. Kinetic studies reveal that much of the change in His-heme loop stability for the aromatic amino acids is caused by a slowdown in the rate of His-heme loop breakage, indicating that residual structure is preferentially stabilized in the closed-loop form of the denatured state.

Denatured states of low-complexity polypeptide sequences differ dramatically from those of foldable sequences

How the primary sequence of a protein encodes conformational preferences that operate early in folding to promote efficient formation of the correct native topology is still poorly understood. To address this issue, we have prepared a set of yeast iso-1-cytochrome c variants that contain polyalanine inserts ranging from 6 to 30 residues in length near the N terminus of the protein. We study the thermodynamics and kinetics of His-heme loop formation in the denatured state at 3 and 6 M guanidine-HCl concentration. We find that polyalanine closely approximates a random coil with excluded volume giving scaling exponents, nu(3), for equilibrium loop formation of 2.26 +/- 0.13 and 1.97 +/- 0.04 in 3 and 6 M guanidine-HCl, respectively. The rate of loop breakage initially decreases and then becomes independent of loop size as would be expected for a random coil. Comparison with previously reported data for denatured state His-heme loop formation for iso-1-cytochrome c and Rhodopseudomonas palustris cytochrome c', shows that foldable sequences deviate significantly from random coil behavior and that the deviation is fold-dependent.

Thermodynamics of loop formation in the denatured state of rhodopseudomonas palustris cytochrome c': scaling exponents and the reconciliation problem

The observation that denatured proteins yield scaling exponents, nu, consistent with random-coil behavior and yet can also have pockets of residual or nonrandom structure has been termed the "reconciliation problem". To provide greater insight into the denatured state of a foldable sequence, we have measured histidine-heme loop formation equilibria in the denatured state of a class II c-type cytochrome, cytochrome c' from Rhodopseudomonas palustris. We have prepared a series of variants that provide His-heme loop stabilities, pK(loop)(His), for loop sizes ranging from 10 to 111 residues at intervals of 7 to 11 residues along the sequence of the protein. We observe a scaling exponent for loop formation, nu(3), of 2.5+/-0.3. Theoretical values for nu(3) range from 1.8 to 2.4; thus, the observed nu(3) is consistent with random-coil behavior. However, in contrast to data for loop formation as a function of loop size obtained with peptides of homogeneous sequence, we observe considerable scatter about the linear dependence of loop stability on loop size. Thus, foldable sequences behave very differently from homogeneous peptide sequences. The observed scatter suggests that there is considerable variation in the conformational properties along the backbone of a foldable sequence, consistent with alternating compact and extended regions. With regard to the reconciliation problem, it is evident that a scaling exponent consistent with a random coil is necessary but not sufficient to demonstrate random-coil behavior.

The unfolded state of the C-terminal domain of the ribosomal protein L9 contains both native and non-native structure

Interest in the structural and dynamic properties of unfolded proteins has increased in recent years owing to continued interest in protein folding and misfolding. Knowledge of the unfolded state under native conditions is particularly important for obtaining a complete picture of the protein folding process. The C-terminal domain of protein L9 is a globular alpha, beta protein with an unusual mixed parallel and antiparallel beta-strand structure. The folding kinetics and equilibrium unfolding of CTL9 strongly depend on pH, and follow a simple two state model. Both the native and the unfolded state can be significantly populated at pH 3.8 in the absence of denaturant, allowing the native state and the unfolded state to be characterized under identical conditions. Backbone (15)N, (13)C, (1)H and side-chain (13)C(beta), (1)H(beta) chemical shifts, amide proton NOEs, and (15)N R(2) relaxation rates were obtained for the two conformational states at pH 3.8. All the data indicate that the pH 3.8 native state is well folded and is similar to the native state at neutral pH. There is significant residual structure in the pH 3.8 unfolded state. The regions corresponding to the two native state alpha-helices show strong preference to populate helical phi and psi angles. The segment that connects alpha-helix 2 and beta-strand 2 has a significant tendency to form non-native alpha-helical structure. Comparison with the pH 2.0 unfolded state and the urea unfolded state indicates that the tendency to adopt both native and non-native helical structure is stronger at pH 3.8, demonstrating that the unfolded state of CTL9 under native-like conditions is more structured. The implications for the folding of CTL9 are discussed.

Importance of contact persistence in denatured state loop formation: kinetic insights into sequence effects on nucleation early in folding

Protein folding is dependent on the formation and persistence of simple loops early in folding. Ease of loop formation and persistence is believed to be dependent on the steric interactions of the residues involved in loop formation. We have previously investigated this factor in the denatured state of iso-1-cytochrome c using a five-amino-acid insert in front of a unique histidine in the N-terminal region of the protein. Previously, we reported that the apparent pK(a) values of loop formation for the most flexible (all Gly) and least flexible (all Ala) insert were, within error, the same. We evaluate whether this observation is due to differences in the persistence of loop contacts or due to effects of local sequence sterics and main-chain hydration on the persistence length of the chain. We also test whether sequence order affects loop formation. Here, we report kinetic results coupled to further mutagenesis of the insert to discern between these possibilities. We find that the amino acid-glycine versus alanine-next to the loop forming histidine has a dominant effect on loop kinetics and equilibria. A glycine in this position speeds loop breakage relative to alanine, resulting in less stable loops. At high percentage of Gly in the insert, rates of loop formation and breakage exactly compensate, leading to a leveling out in loop stability. Loop formation rates also increase with glycine content, inconsistent with poly-Gly segments being more extended than previously suspected due to main-chain hydration or local sterics. Unlike loop breakage rates, loop formation rates are insensitive to local sequence. Together, these observations suggest that contact persistence plays a more important role in defining the "folding code" than rates of loop formation.

NMR analysis of native-state protein conformational flexibility by hydrogen exchange

The rate of hydrogen exchange for the most protected amides of a protein is widely used to provide an estimate of global conformational stability by analyzing the exchange kinetics in the unfolded state in terms of model peptide exchange rates. The exchange behavior of the other amides of the protein which do not exchange via a global unfolding mechanism can provide insight into the smaller-scale conformational transitions that facilitate access to solvent as required for the exchange reaction. However, since the residual tertiary structure in the exchange-competent conformation can modulate the chemistry of the exchange reaction, equilibrium values estimated from normalization with model peptide rates are open to question. To overcome this limitation, the most robust approaches utilize differential analyses as a function of experimental variables such as denaturant concentration, temperature, pH, and mutational variation. Practical aspects of these various differential analysis techniques are considered with illustrations drawn from the literature.

Experimental characterization of the denatured state ensemble of proteins

The traditional view of the denatured state ensemble of proteins is that it behaves as a classic random coil. This model has important implications for the analysis of protein stability, protein folding, and cooperativity; namely that the effects of mutations on the free energy of the denatured state ensemble can be ignored. This assumption, which is still routinely made, at least at the implicit level, greatly simplifies the analysis of such experiments. However it has long been recognized that the denatured state ensemble (DSE) of real proteins is often quite different from a random coil and can exhibit significant structural preferences. In some cases parts of the chain can even adopt relatively well-defined conformations, particularly under native conditions. Well-studied examples of DSE interactions include elements of hydrogen-bonded secondary structure, particularly helices or turns, as well hydrophobic clusters, hydrophobic aromatic clusters, and more recently interactions involving charged residues. Deviations from random-coil behavior are of practical importance if they influence protein folding, stability, or function, or if they compromise our analysis and interpretation of experiments. The existence of residual structure in the DSE naturally leads to the question of its role in protein folding and stability, and raises the possibility that some mutations could exert a significant part of their effect by altering the DSE. Much of our understanding of the interactions governing protein stability and the folding process have been generated by mutational studies; thus, a detailed understanding of the denatured state ensemble is critical.

Liaison amid disorder: non-native interactions may underpin long-range coupling in proteins

A lattice-model study of double-mutant cycles published in BMC Structural Biology underscores how interactions in non-native conformations can lead to thermodynamic coupling between distant residues in globular proteins, adding to recent advances in delineating the often crucial roles played by disordered conformational ensembles in protein behavior.

Insight into "insoluble proteins" with pure water

Many proteins are not refoldable and also insoluble. Previously no general method was available to solubilize them and consequently their structural properties remained unknown. Surprisingly, we recently discovered that all insoluble proteins in our laboratory, which are highly diverse, can be solubilized in pure water. Structural characterization by CD and NMR led to their classification into three groups, all of which appear trapped in the highly disordered or partially-folded states with a substantial exposure of hydrophobic side chains. In this review, I discuss our results in a wide context and subsequently propose a model to rationalize the discovery. The potential applications are also explored in studying protein folding, design and membrane proteins.

Identification and characterization of folding inhibitors of hen egg lysozyme: an example of a new paradigm of drug design

Studies of protein folding indicate the presence of native contacts in the denatured state, giving rise to folding elements which contribute to the accomplishment of the native state. The possibility of finding molecules which can interact with specific folding elements of a target protein preventing it from reaching its native state, and hence from becoming biologically active, is particularly attractive. The notion that folding elements not only provide molecular recognition directing the folding process, but also have conserved sequence, implies that targeting such elements will make protein folding inhibitors less susceptible to mutations which, in many cases, abrogate drug effects. The folding-inhibition strategy can lead to a truly novel and rational approach to drug design, aside from providing new insight into folding. This is illustrated in the case of hen egg lysozyme.

Competition between reversible aggregation and loop formation in denatured iso-1-cytochrome c

The competition between intramolecular histidine-heme loop formation and ligand-mediated oligomer formation in the denatured state is investigated for two yeast iso-1-cytochrome c variants, AcH26I52 and AcA25H26I52. Besides the native His 18 heme ligand, both variants contain a single His at position 26. The AcA25H26I52 variant has Pro 25 mutated to Ala. The concentration dependence of the apparent pK(a) for His 26-heme binding in 3 M guanidine hydrochloride indicates that the P25A mutation disfavors oligomerization mediated by intermolecular heme ligation by 10-fold. Single- and double-pH-jump stopped-flow experiments with the AcH26I52 variant show that fast phases for His-heme bond formation and breakage are due to intramolecular loop formation and slow phases for His-heme bond formation and breakage are due to intermolecular aggregation. The presence of two closely spaced slow phases in the kinetics of loop formation for both variants suggests that intermolecular His 26-heme ligation results in both dimers and higher-order aggregates. The P25A mutation slows formation and speeds breakdown of an initial dimer, demonstrating a strong effect of local sequence on aggregation. Analysis of the kinetic data yields equilibrium constants for intramolecular loop formation and intermolecular dimerization at pH 7.1 and indicates that the rate constant for intermolecular aggregation is very fast at this pH (10(7)-10(8) M(-1) s(-1)). In light of the very fast rates of aggregation in the denatured state, comparison of models involving reversible or irreversible oligomerization steps suggests that equilibrium control of the partitioning between folding and aggregation is advantageous for productive protein folding in vivo.

Stability of the heme Fe-N-terminal amino group coordination bond in denatured cytochrome c

In the denatured states of Hydrogenobacter thermophilus cytochrome c(552) (HT) and Pseudomonas aeruginosa cytochrome c(551) (PA), and their mutants, the N-terminal amino group of the polypeptide chain is coordinated to heme Fe in place of the axial Met, the His-N(term) form being formed. The coordination of the N-terminal amino group to heme Fe leads to loop formation by the N-terminal stretch preceding the first Cys residue bound to the heme, and the N-terminal stretches of HT and PA are different from each other in terms of both the sequence and the number of constituent amino acid residues. The His-N(term) form was shown to be rather stable, and hence it can influence the stability of the denatured state. We have investigated the heme Fe coordination structures and stabilities of the His-N(term) forms emerging upon guanidine hydrochloric acid-induced unfolding of the oxidized forms of the proteins. The Fe-N(term) coordination bond in the His-N(term) form with a 9-residue N-terminal stretch of HT proteins was found to be tilted to some extent away from the heme normal, as reflected by the great heme methyl proton shift spread. On the other hand, the small heme methyl proton shift spread of the His-N(term) form with an 11-residue stretch of PA proteins indicated that its Fe-N(term) bond is nearly parallel with the heme normal. The stability of the His-N(term) form was found to be affected by the structural properties of the N-terminal stretch, such as its length and the N-terminal residue. With a given N-terminal residue, the stability of the His-N(term) form is higher for a 9-residue N-terminal stretch than an 11-residue one. In addition, with a given length of the N-terminal stretch, the His-N(term) form with an N-terminal Glu is stabilized by a few kJ mol(-1) relative to that with an N-terminal Asn. These results provide a novel insight into the stabilizing interactions in the denatured cyts c that will facilitate elucidation of the folding/unfolding mechanisms of the proteins.

The extremely slow-exchanging core and acid-denatured state of green fluorescent protein

Green fluorescent protein (GFP) is a large protein with a complex eleven-stranded beta-barrel structure. Previous studies have shown that it has a complex energy landscape for folding on which there are several intermediate states and a denatured state with significant residual structure. Here, we use two different types of HD exchange measurement and nuclear magnetic resonance (NMR) techniques to probe the energy landscape for folding of GFP in further detail. HD exchange experiments were performed over a wide range of conditions including different concentrations of denaturant. Results show that the penetration model dominates the exchange mechanism, consistent with the known stability and slow unfolding kinetics of GFP. HD exchange experiments at high pH establish that there is an extremely slow-exchanging superstable core of amide protons in GFP that are clustered and located in beta-strands 1, 2, 4, 5, and 6. These residues form part of a mini-beta-sheet which we propose constitutes a folding nucleus. Using a pulsed-labeling strategy, the acid-denatured state has been investigated and the residual structure observed in earlier studies shown to locate to beta-strands 1 and 3. There is some evidence that this residual structure is stabilized by a localized hydrophobic collapse of the polypeptide chain.

Structure and energetics of the hydrogen-bonded backbone in protein folding

We seek to understand the link between protein thermodynamics and protein structure in molecular detail. A classical approach to this problem involves assessing changes in protein stability resulting from added cosolvents. Under any given conditions, protein molecules in aqueous buffer are in equilibrium between unfolded and folded states, U(nfolded) <==> N(ative). Addition of organic osmolytes, small uncharged compounds found throughout nature, shift this equilibrium. Urea, a denaturing osmolyte, shifts the equilibrium toward U; trimethylamine N-oxide (TMAO), a protecting osmolyte, shifts the equilibrium toward N. Using the Tanford Transfer Model, the thermodynamic response to many such osmolytes has been dissected into groupwise free energy contributions. It is found that the energetics involving backbone hydrogen bonding controls these shifts in protein stability almost entirely, with osmolyte cosolvents simply dialing between solvent-backbone versus backbone-backbone hydrogen bonds, as a function of solvent quality. This reciprocal relationship establishes the essential link between protein thermodynamics and the protein's hydrogen-bonded backbone structure.

Effects of excluded volume upon protein stability in covalently cross-linked proteins with variable linker lengths

To explore the effects of molecular crowding and excluded volume upon protein stability, we used a series of cross-linking reagents with nine different single-cysteine mutants of staphylococcal nuclease to make covalently linked dimers. These cross-linkers ranged in length from 10.5 to 21.3 A, compelling separations which would normally be found only in the most concentrated protein solutions. The stabilities of the dimeric proteins and monomeric controls were determined by guanidine hydrochloride and thermal denaturation. Dimers with short linkers tend to exhibit pronounced three-state denaturation behavior, as opposed to the two-state behavior of the monomeric controls. Increasing linker length leads to less pronounced three-state behavior. The three-state behavior is interpreted in a three-state model where cross-linked native protein dimer, N-N, interconverts in a two-state transition with a dimer where one protein subunit is denatured, N-D. The remaining native protein in turn can denature in another two-state transition to a state, D-D, in which both tethered proteins are denatured. Three-state behavior is best explained by excluded volume effects in the denatured state. For many dimers, linkers longer than 17 A removed most three-state character. This sets a limit on the flexibility and size of the denatured state. Notably, in contradiction to theoretical predictions, these cross-linked dimers were not stabilized. The failure of these predictions is possibly due to neglect of the alteration in hydrophobic exposure that accompanies any significant reduction in the conformational space of the denatured state.

Molecular processes in biological thermosensation

Since thermal gradients are almost everywhere, thermosensation could represent one of the oldest sensory transduction processes that evolved in organisms. There are many examples of temperature changes affecting the physiology of living cells. Almost all classes of biological macromolecules in a cell (nucleic acids, lipids, proteins) can present a target of the temperature-related stimuli. This review discusses some features of different classes of temperature-sensing molecules as well as molecular and biological processes that involve thermosensation. Biochemical, structural, and thermodynamic approaches are applied in the paper to organize the existing knowledge on molecular mechanisms of thermosensation. Special attention is paid to the fact that thermosensitive function cannot be assigned to any particular functional group or spatial structure but is rather of universal nature. For instance, the complex of thermodynamic, structural, and functional features of hemoglobin family proteins suggests their possible accessory role as “molecular thermometers”.