Heat capacity of hydrogen-bonded networks: an alternative view of protein folding thermodynamics

Large changes in heat capacity (deltaCp) have long been regarded as the characteristic thermodynamic signature of hydrophobic interactions. However, similar effects arise quite generally in order-disorder transitions in homogeneous systems, particularly those comprising hydrogen-bonded networks, and this may have significance for our understanding of protein folding and other biomolecular processes. The positive deltaCp associated with unfolding of globular proteins in water, thought to be due to hydrophobic interactions, is also typical of the values found for the melting of crystalline solids, where the effect is greatest for the melting of polar compounds, including pure water. This suggests an alternative model of protein folding based on the thermodynamics of phase transitions in hydrogen-bonded networks. Folded proteins may be viewed as islands of cooperatively-ordered hydrogen-bonded structure, floating in an aqueous network of less-well-ordered H-bonds in which the degree of hydrogen bonding decreases with increasing temperature. The enthalpy of melting of the protein consequently increases with temperature. A simple algebraic model, based on the overall number of protein and solvent hydrogen bonds in folded and unfolded states, shows how deltaCp from this source could match the hydrophobic contribution. This confirms the growing view that the thermodynamics of protein folding, and other interactions in aqueous systems, are best described in terms of a mixture of polar and non-polar effects in which no one contribution is necessarily dominant.     

Heat capacity of protein folding.

We construct a Hamiltonian for a single domain protein where the contact enthalpy and the chain entropy decrease linearly with the number of native contacts. The hydration effect upon protein unfolding is included by modeling water as ideal dipoles that are ordered around the unfolded surfaces, where the influence of these surfaces, covered with an "ice-like" shell of water, is represented by an effective field that directs the water dipoles. An intermolecular pair interaction between water molecules is also introduced. The heat capacity of the model exhibits, the common feature of small globular proteins, two peaks corresponding to cold and warm unfolding, respectively. By introducing ad hoc vibrational modes, we obtain quantitatively good accordance with experiments on myoglobin.     

Advances in free-energy-based simulations of protein folding and ligand binding

1. Download : Download high-res image (122KB)
2. Download : Download full-size image

     

Confinement effects on the thermodynamics of protein folding: Monte Carlo simulations

The effects of chaperonin-like cage-induced confinement on protein stability have been studied for molecules of varying sizes and topologies. Minimalist models based on Gō-like interactions are employed for the proteins, and density-of-states-based Monte Carlo simulations are performed to accurately characterize the thermodynamic transitions. This method permits efficient sampling of conformational space and yields precise estimates of free energy and entropic changes associated with protein folding. We find that confinement-driven stabilization is not only dependent on protein size and cage radius, but also on the specific topology. The choice of the confining potential is also shown to have an effect on the observed stabilization and the scaling behavior of the stabilization with respect to the cage size.     

The role of water in ligand binding

1. Download : Download high-res image (132KB)
2. Download : Download full-size image

Exploration of the complex modulatory role of water in ligand–target binding is a current challenge of drug design. This review reports on recent advances of prediction of water structure and function in the context of ligand engineering. The surveyed theoretical approaches showed remarkable progress in the past years. Beyond complementing experiments, they also supplied unmeasurable data. For example, thermodynamic calculations improved ligand binding by the selection of certain water molecules for structural replacement. Molecular dynamics and explicit solvent models remained indispensable to achieve precise results. Topographical analyses of hydration networks proved useful for the prediction of the stabilizing role of interconnected water clusters mediating target–ligand interactions.     

Disease-causing SAP mutants are defective in ligand binding and protein folding

The X-linked lymphoproliferative (XLP) syndrome is caused by mutations or deletions in the SH2D1A gene that encodes an SH2 domain protein named SH2D1A or SAP. The identification of a number of missense mutations within the protein's SH2 domain, each of which can directly cause disease, provides a unique opportunity to investigate the function of an interaction protein module, SH2, in the pathogenesis of XLP. We show here that SAP mutants found in XLP patients are defective in binding its physiological ligands signaling lymphocyte activating molecule (SLAM), a co-receptor in T cell activation, and Fyn, a Src family protein tyrosine kinase. Consequently, these mutants are deficient in signaling through the SLAM receptor. This is reflected by compromised abilities for the mutants to recruit Fyn to SLAM and to activate Fyn, by reduced phosphorylation of the receptor, and by deficiencies for the mutants in blocking binding of SHP-2 to SLAM. Furthermore, all mutants examined are defective in protein folding as manifested by their significantly reduced melting temperatures upon thermal denaturation, compared to that of SAP. Taken together, these results suggest that defects in ligand binding, receptor signaling, and protein folding collectively contribute to the loss of function for disease-causing SAP mutants.     

The binding of phosphatidylcholine to the phosphatidylcholine transfer protein: affinity and role in folding

Bovine liver phosphatidylcholine transfer protein (PC-TP) has been expressed in Escherichia coli and purified to homogeneity from the cytosol fraction at a yield of 0.45 mg PC-TP per 10 mg total cytosolic protein. In addition, active PC-TP was obtained from inclusion bodies. An essential factor in the activation of PC-TP was phosphatidylcholine (PC) present in the folding buffer. PC-TP from the cytosol contains phosphatidylethanolamine (PE) and phosphatidylglycerol (PG) with a preference for the di-monounsaturated species over the saturated species as determined by fast atom bombardment mass spectrometry (FAB-MS). By incubation with microsomal membranes the endogenous PE and PG were replaced by PC. Relative to the microsomal PC species composition, PC-TP bound preferentially C16:0/C20:4-PC and C16:0/C18:2-PC (twofold enriched) whereas the major microsomal species C18:0/C18:1-PC and C18:0/C18:2-PC were distinctly less bound. PC-TP is structurally homologous to the lipid-binding domain of the steroidogenic acute regulatory protein (Nat. Struct. Biol. 7 (2000) 408). Replacement of Lys(55) present in one of the beta-strands forming the lipid-binding site, with an isoleucine residue yielded an inactive protein. This suggests that Lys(55) be involved in the binding of the PC molecule.     

Crystal structure of a berenil-dodecanucleotide complex: the role of water in sequence-specific ligand binding.

The three-dimensional structure of a complex between the dodecanucleotide d(CGCGAATTCGCG) and the anti-trypanocidal drug berenil, has been determined to a resolution of 2.5 A. The structure has been solved by molecular replacement and refined to an R factor of 0.177. A total of 49 water molecules have been located. The drug is bound at the 5'-AAT-3' region of the oligonucleotide. At one end of the drug the amidinium group is in hydrogen-bonded contact with N3 of the adenine base complementary to the thymine of the AAT. The other amidinium group does not make direct interactions with the DNA. Instead, a water molecule mediates between them. This is in hydrogen-bonded contact with an amidinium nitrogen atom, N3 of the 5' end adenine base and the ring oxygen atom of an adjacent deoxyribose. Molecular mechanics calculations have been performed on this complex, with the drug at various positions along the sequence. These show that the observed position is only 0.8 kcal/mol higher in energy than the best position. It is suggested that there is a broad energy well in the AATT region for this drug, and that water molecules as well as the neighbouring sequence, will determine precise positioning. More general aspects of minor groove binding are discussed.     

Heat capacity effects of water molecules and ions at a protein-DNA interface

The interaction of the TATA-box binding protein from the thermophilic and halophilic archaea Pyrococcus woesei (PwTBP) with an oligonucleotide containing a specific binding site is stable over a very broad range of temperatures and ionic strengths, and is consequently an outstanding system for characterising general features of protein-DNA thermodynamics. In common with other specific protein-DNA recognition events, the PwTBP-TATA box interaction is accompanied by a large negative change in heat capacity (ΔC_p) arising from the total change in solvation that occurs upon binding, which in this case involves a net uptake of cations. Contrary to previous hypotheses, we find no overall effect of ionic strength on this heat capacity change. We investigate the local contributions of site-specific ion and water binding to the overall change in heat capacity by means of a series of site-directed mutations of PwTBP. We find that although changes in the local ion binding capacity affect the enthalpic and entropic contributions to the free energy of the interaction, they do not affect the change in heat capacity. In contrast, we find remarkably large heat capacity effects arising from two particular symmetry-related mutations. The great magnitude of these effects is not explicable in terms of current semi-empirical models of heat capacity change. Previously reported X-ray crystal structures show that these mutated residues are at the centre of an evolutionarily conserved network of water-mediated hydrogen bonds between the protein and the DNA backbone. Consequently, we conclude that, in addition to water molecules buried in the protein-DNA interface that have been previously shown to influence heat capacity, bridging water molecules in a highly polar surface environment can also contribute substantially to negative heat capacity change on formation of a protein-DNA complex.     

Genomic determinants of protein folding thermodynamics in prokaryotic organisms

Here we investigate how thermodynamic properties of orthologous proteins are influenced by the genomic environment in which they evolve. We performed a comparative computational study of 21 protein families in 73 prokaryotic species and obtained the following main results. (i) Protein stability with respect to the unfolded state and with respect to misfolding are anticorrelated. There appears to be a trade-off between these two properties, which cannot be optimized simultaneously. (ii) Folding thermodynamic parameters are strongly correlated with two genomic features, genome size and G+C composition. In particular, the normalized energy gap, an indicator of folding efficiency in statistical mechanical models of protein folding, is smaller in proteins of organisms with a small genome size and a compositional bias towards A+T. Such genomic features are characteristic for bacteria with an intracellular lifestyle. We interpret these correlations in light of mutation pressure and natural selection. A mutational bias toward A+T at the DNA level translates into a mutational bias toward more hydrophobic (and in general more interactive) proteins, a consequence of the structure of the genetic code. Increased hydrophobicity renders proteins more stable against unfolding but less stable against misfolding. Proteins with high hydrophobicity and low stability against misfolding occur in organisms with reduced genomes, like obligate intracellular bacteria. We argue that they are fixed because these organisms experience weaker purifying selection due to their small effective population sizes. This interpretation is supported by the observation of a high expression level of chaperones in these bacteria. Our results indicate that the mutational spectrum of a genome and the strength of selection significantly influence protein folding thermodynamics.     

Kinetics, thermodynamics and evolution of non-native interactions in a protein folding nucleus

A lattice model with side chains was used to investigate protein folding with computer simulations. In this model, we rigorously demonstrate the existence of a specific folding nucleus. This nucleus contains specific interactions not present in the native state that, when weakened, slow folding but do not change protein stability. Such a decoupling of folding kinetics from thermodynamics has been observed experimentally for real proteins. From our results, we conclude that specific non-native interactions in the transition state would give rise to straight phi-values that are negative or larger than unity. Furthermore, we demonstrate that residue Ile 34 in src SH3, which has been shown to be kinetically, but not thermodynamically, important, is universally conserved in proteins with the SH3 fold. This is a clear example of evolution optimizing the folding rate of a protein independent of its stability and function.     

Alpha-synuclein multistate folding thermodynamics: implications for protein misfolding and aggregation

Alpha-synuclein aggregation has been tightly linked with the pathogenesis of Parkinson's disease and other neurodegenerative disorders. Despite the protein's putative function in presynaptic vesicle regulation, the roles of lipid binding in modulating alpha-synuclein conformations and the aggregation process remain to be fully understood. This study focuses on a detailed thermodynamic characterization of monomeric alpha-synuclein folding in the presence of SDS, a well-studied lipid mimetic. Far-UV CD spectroscopy was employed for detection of conformational transitions induced by SDS, temperature, and pH. The data we present here clearly demonstrate the multistate nature of alpha-synuclein folding, which involves two predominantly alpha-helical partially folded thermodynamic intermediates that we designate as F (most folded) and I (intermediately folded) states. Likely structures of these alpha-synuclein conformational states are also discussed. These partially folded forms can exist in the presence of either monomeric or micellar forms of SDS, which suggests that alpha-synuclein has an intrinsic propensity for adopting multiple alpha-helical structures even in the absence of micelle or membrane binding, a feature that may have implications for its biological activity and toxicity. Additionally, we discuss the relation between alpha-synuclein three-state folding and its aggregation, within the context of isothermal titration calorimetry and transmission electron microscopy measurements of SDS-initiated oligomer formation.     

The intestinal fatty acid binding protein: the role of turns in fast and slow folding processes

The intestinal fatty acid binding protein is one of a family of proteins that are composed of two beta-sheets surrounding a large interior cavity into which the ligand binds. Glycine residues occur in many of the turns between adjacent antiparallel beta-strands. In previous work, the effect of replacing these glycine residues with valine has been examined with stopped flow instrumentation using intrinsic tryptophan fluorescence spectroscopy [Kim and Frieden (1998) Protein Sci. 7, 1821-1828]. To resolve the burst phase missing in the stopped flow measurements, these valine mutants have been reexamined with sub-millisecond continuous flow instrumentation. Some of the glycine residues have also been replaced with proline, and the folding reactions of these proline mutants have been compared with those of their valine counterparts. In all cases, the stability of the protein is decreased, but some turns appear to be more critical for final structure stabilization than others. Surprisingly, the rate constants observed for all the mutants measured by sub-millisecond continuous flow methods are quite similar (1400-3000 s(-1)), and in all the mutants, there is a shift in the fluorescence emission maximum from that of the unfolded protein to lower wavelengths, suggesting some collapse of the unfolded state within 200 micros. In contrast to the rate constants observed for the initial folding events measured by the sub-millisecond continuous flow method, the rate constants for the slower phase observed in the stopped flow instrument vary widely for the different mutants. The latter step appears to be related to side chain stabilization rather than secondary structure formation. It is also shown that the ligand binds tightly only to the native protein and not to any intermediate forms.     

The role of Trp-82 in the folding of intestinal fatty acid binding protein

Multiple phases have been observed during the folding and unfolding of intestinal fatty acid binding protein (WT-IFABP) by stopped-flow fluorescence. Site-directed mutagenesis has been used to examine the role of each of the two tryptophans of this protein in these processes. The unfolding and refolding kinetics of the mutant protein containing only tryptophan 82 (W6Y-IFABP) showed that the tryptophan at this location was critical to the fluorescence signal changes observed throughout the unfolding reaction and early in the refolding reaction. However, the kinetic patterns of the mutant protein containing only tryptophan 6 (W82Y-IFABP) indicated that the tryptophan at this location participated in the fluorescence signal changes observed early in the unfolding reaction and late in the refolding reaction. Together, these data suggest that native-like structure was formed first in the vicinity of tryptophan 82, near the center of the hydrophobic core of this beta-sheet protein, prior to formation of native-like structure in the periphery of the protein.     

Gatekeepers in the ribosomal protein s6: thermodynamics, kinetics, and folding pathways revealed by a minimalist protein model

We investigate the effect of structural gatekeepers on the folding of the ribosomal protein S6. Folding thermodynamics and early refolding kinetics are studied for this system utilizing computer simulations of a minimalist protein model. When gatekeepers are eliminated, the thermodynamic signature of a folding intermediate emerges, and a marked decrease in folding efficiency is observed. We explain the prerequisites that determine the "strength" of a given gatekeeper. The investigated gatekeepers are found to have distinct functions, and to guide the folding and time-dependent packing of non-overlapping secondary structure elements in the protein. Gatekeepers avoid kinetic traps during folding by favoring the formation of "productive topologies" on the way to the native state. The trends in folding rates in the presence/absence of gatekeepers observed for our minimalist model of S6 are in very good agreement with experimental data on this protein.     

Structural energetics of protein folding and binding

Structural energetics is a method for calculating the energetics of protein folding and binding reactions as a function of temperature. This approach allows measured energetics to be interpreted with regards to the protein structure and the prediction of energetics from known structures. Recent advances include improvements in the parameterization of enthalpy, entropy and heat capacity terms and new applications, especially with regards to understanding dynamic properties of proteins and how these are affected by ligand binding.     

The role of atomic level steric effects and attractive forces in protein folding

Protein folding into tertiary structures is controlled by an interplay of attractive contact interactions and steric effects. We investigate the balance between these contributions using structure‐based models using an all‐atom representation of the structure combined with a coarse‐grained contact potential. Tertiary contact interactions between atoms are collected into a single broad attractive well between the C^β atoms between each residue pair in a native contact. Through the width of these contact potentials we control their tolerance for deviations from the ideal structure and the spatial range of attractive interactions. In the compact native state dominant packing constraints limit the effects of a coarse‐grained contact potential. During folding, however, the broad attractive potentials allow an early collapse that starts before the native local structure is completely adopted. As a consequence the folding transition is broadened and the free energy barrier is decreased. Eventually two‐state folding behavior is lost completely for systems with very broad attractive potentials. The stabilization of native‐like residue interactions in non‐perfect geometries early in the folding process frequently leads to structural traps. Global mirror images are a notable example. These traps are penalized by the details of the repulsive interactions only after further collapse. Successful folding to the native state requires simultaneous guidance from both attractive and repulsive interactions. Proteins 2012. © 2011 Wiley Periodicals, Inc.     

Design of beta-sheet systems for understanding the thermodynamics and kinetics of protein folding

Peptide beta-sheet systems have emerged as context-independent models for probing secondary structure propensities, the nature and magnitude of stabilizing weak interactions, and aspects of cooperativity both parallel and perpendicular to the strand direction. These systems have allowed fundamental advances in understanding non-covalent interactions relevant to both chemical and biological systems, and in describing the protein folding energy landscape.