Unit 3: Thermodynamics of biomolecules - Subjective Questions

Enthalpy (): The heat content of a system at constant pressure. In biomolecules it reflects the energy of bond formation/breaking, including hydrogen bonds, van der Waals interactions, and electrostatic interactions.

Entropy (): A measure of disorder or the number of accessible microstates. In biological systems it includes conformational entropy of the chain and the entropy of surrounding water molecules.

Gibbs Free Energy (): The thermodynamic potential that determines spontaneity at constant temperature and pressure.

The fundamental relationship is:

Significance:

• If : the process is spontaneous (favorable).
• If : the process is non-spontaneous.
• If : the system is at equilibrium.

In biomolecular processes, the balance between enthalpic contributions (favorable interactions) and entropic contributions (chain disorder vs. hydrophobic effect) governs whether structures form spontaneously.

The hydrophobic effect refers to the tendency of nonpolar molecules or side chains to aggregate in aqueous solution, minimizing their contact with water.

Molecular Basis:

• Water molecules around a nonpolar solute form ordered, cage-like structures (clathrates), reducing water entropy.
• When nonpolar groups cluster together in the protein core, ordered water is released into the bulk solvent.

Thermodynamic Contribution:

• The release of ordered water increases the entropy of the solvent ().
• This makes strongly negative, contributing favorably to .

Role in Folding:

• The hydrophobic effect drives the burial of nonpolar residues (Leu, Ile, Val, Phe) in the protein interior.
• It is considered the dominant force stabilizing the folded (native) state.
• It explains why proteins fold into compact globular structures with hydrophobic cores and hydrophilic surfaces.

Interestingly, the hydrophobic effect is largely entropy-driven at physiological temperatures, distinguishing it from ordinary enthalpic bonding.

Two-State Model:
For many small, single-domain proteins, folding can be approximated as a reversible transition between only two significantly populated states:

where = unfolded state and = native (folded) state. Intermediate states are negligibly populated.

Equilibrium Constant:

Derivation of Free Energy Relationship:
The standard free energy of folding is related to the equilibrium constant by:

where:

• = universal gas constant ()
• = absolute temperature (K)

Fraction Folded:
If is the fraction folded and :

Key Features of the Model:

• Only and are populated; the transition is cooperative.
• Folding curves are sigmoidal.
• The stability of a protein is typically small, to (marginal stability).

This marginal stability allows proteins to be flexible enough for function yet stable enough to maintain structure.

Thermodynamics of Folding:

• Concerns the equilibrium and stability of the folded state.
• Determines whether folding is favorable, defined by .
• Described by Anfinsen's hypothesis: the native structure corresponds to the global free energy minimum.
• Independent of the pathway taken.

Kinetics of Folding:

• Concerns the rate and pathway by which folding occurs.
• Determines how fast and by what route a protein reaches its native state.
• Described by transition states, folding intermediates, and energy barriers.

Comparison Table:

Why Both Matter:

• Thermodynamics explains the stability of the final structure.
• Kinetics explains Levinthal's paradox — that proteins fold rapidly (ms–s) rather than searching all conformations randomly (which would take longer than the age of the universe).
• Together they explain the folding funnel concept, where the energy landscape is biased toward the native state, providing both a thermodynamic sink and kinetically accessible routes.

Anfinsen's Thermodynamic Hypothesis:
The hypothesis states that the native three-dimensional structure of a protein is determined solely by its amino acid sequence, and that this structure corresponds to the thermodynamically most stable conformation (global free energy minimum) under physiological conditions.

Key Implications:

• All information needed for folding is encoded in the primary sequence.
• Folding is a spontaneous, reversible process.
• The native state is at the global minimum of free energy, not merely a kinetic trap.

Experimental Evidence (Ribonuclease A experiments):

• Anfinsen denatured ribonuclease A using urea (denaturant) and -mercaptoethanol (reduces disulfide bonds), completely unfolding and inactivating the enzyme.
• Upon removal of denaturant and reducing agent, the protein spontaneously refolded.
• The refolded protein regained full enzymatic activity and the correct disulfide bonds reformed.

Significance of the Result:

• Refolding to the correct native state without any external template proved that the sequence alone dictates structure.
• If disulfides reoxidized under denaturing conditions, incorrect ("scrambled") bonds formed, giving inactive protein — showing native structure is the thermodynamic minimum only under proper conditions.

Limitations:

• Some proteins require chaperones or specific conditions to fold correctly in vivo, and some are kinetically trapped, so the hypothesis is a valuable but not universal principle.

Levinthal's Paradox:
Proposed by Cyrus Levinthal, this paradox highlights a contradiction in protein folding:

• Consider a protein of residues, each with several possible conformations.
• If a protein samples conformations randomly, the number of possible states is astronomically large (e.g., for 100 residues).
• A random search would require a time longer than the age of the universe.
• Yet real proteins fold within milliseconds to seconds.

Conclusion: Folding cannot be a random search; it must be directed by some mechanism.

Resolution — The Folding Funnel:
The modern energy landscape theory resolves this paradox:

• The free energy landscape is shaped like a funnel, not a flat golf course.
• The width of the funnel represents conformational entropy (number of states).
• The depth represents decreasing free energy toward the native state.

Key Features:

• As folding proceeds, both energy and entropy decrease, guiding the chain downhill.
• Multiple pathways lead to the native state at the bottom of the funnel.
• Partially folded intermediates and local minima may exist but are biased toward the native state.
• Folding is a biased, parallel search, not a random one.

Ruggedness of the Landscape:

• A smooth funnel gives fast folding.
• A rugged funnel with deep local minima can cause kinetic traps and misfolding.

Thus, the funnel model shows folding is directed by a downhill free energy gradient, enabling rapid folding without exhaustive searching.

Protein structure is stabilized by several weak, non-covalent interactions that collectively determine the native fold.

1. Hydrogen Bonds:

• Formed between donor (N–H, O–H) and acceptor (C=O) groups.
• Crucial for secondary structure (-helices, -sheets).
• Energy: – each. Net stabilization is small because unfolded state also H-bonds with water.

2. Van der Waals (London Dispersion) Forces:

• Weak, short-range attractions between transiently polarized atoms.
• Individually weak but numerous in the tightly packed core.
• Important for close packing and shape complementarity.

3. Electrostatic (Ionic) Interactions / Salt Bridges:

• Between oppositely charged side chains (e.g., Lys and Glu).
• Energy depends on distance and dielectric environment.
• Stronger in the protein interior (low dielectric).

4. Hydrophobic Interactions:

• Aggregation of nonpolar groups driven by the entropy of released water.
• Considered the dominant stabilizing force.

5. Disulfide Bonds (covalent, but stabilizing):

• Covalent S–S bonds between cysteine residues.
• Reduce conformational entropy of the unfolded state, stabilizing the fold.

Relative Contributions:

• Hydrophobic effect provides the largest net stabilization.
• Hydrogen bonds contribute to specificity more than net stability.
• Van der Waals forces stabilize tight packing.
• The net stability is a small difference between large opposing forces, giving proteins marginal stability ( to ).

Cold Denaturation:
Most proteins denature (unfold) not only at high temperatures but also at low temperatures, a counterintuitive phenomenon known as cold denaturation.

Thermodynamic Basis:
The stability of a protein is given by the free energy of folding:

Both and are temperature-dependent because of the large heat capacity change () on unfolding:

Stability Curve:

• Plotting against temperature gives an inverted parabola (the protein stability curve).
• The curve crosses at two temperatures:
  • (high) = heat denaturation temperature.
  • (low) = cold denaturation temperature.

Role of the Hydrophobic Effect:

• The large positive arises mainly from exposure of hydrophobic groups to water upon unfolding.
• At low temperatures, the entropic penalty of the hydrophobic effect decreases, so the hydrophobic effect weakens.
• This makes the folded state less stable, leading to unfolding at low temperature.

Significance:

• Cold denaturation is usually observed below 0 °C (often requiring pressure or additives to observe experimentally).
• It confirms the temperature-dependent, entropy-driven nature of hydrophobic stabilization.

Melting Temperature ():
The melting temperature is the temperature at which half of the protein molecules are unfolded and half are folded. At :

At this point:

A higher indicates greater thermal stability.

Differential Scanning Calorimetry (DSC):
DSC directly measures the heat capacity () of a protein solution as a function of temperature.

Principle:

• The sample cell (with protein) and reference cell (buffer only) are heated at a constant rate.
• The instrument measures the excess heat needed to keep both cells at the same temperature.

Information Obtained:

• Peak position gives .
• Area under the peak gives the calorimetric enthalpy () of unfolding:
  
• Van't Hoff enthalpy () is obtained from the peak shape.
• The ratio indicates a two-state transition.
• The shift in baseline gives the heat capacity change ().

Advantages:

• Model-independent, direct measurement of .
• Provides , , , and in a single experiment.
• Useful for studying protein stability, ligand binding, and folding cooperativity.

The molten globule is a partially folded intermediate state observed during protein folding and under mildly denaturing conditions.

Comparison Table:

Key Features of the Molten Globule:

• Retains native-like secondary structure (-helices, -sheets).
• Lacks the specific tight tertiary packing of the native state.
• Binds hydrophobic dyes like ANS due to exposed hydrophobic surfaces.

Significance:

• Considered an important kinetic intermediate on the folding pathway.
• Provides insight into the sequence of events during folding — secondary structure often forms before final tertiary packing.
• Relevant to understanding misfolding and aggregation.

Conformational Entropy:
Conformational entropy refers to the entropy associated with the large number of possible conformations (backbone , angles and side-chain rotamers) available to a polypeptide chain.

Opposition to Folding:

• In the unfolded state, the chain can adopt an enormous number of conformations → high entropy.
• In the folded (native) state, the chain is restricted to essentially one conformation → low entropy.
• Therefore folding involves a large decrease in conformational entropy ().

This contributes unfavorably to because:

A negative makes positive, opposing folding.

How the Penalty is Overcome:
Folding still occurs because favorable terms outweigh the entropic penalty:

• Hydrophobic effect: Release of ordered water increases solvent entropy, offsetting the loss of chain entropy.
• Favorable enthalpy: Formation of hydrogen bonds, van der Waals contacts, and salt bridges provides negative .
• Disulfide bonds: Reduce the entropy of the unfolded state, lowering the entropic cost of folding.

Net Result:
The balance produces a small negative (marginal stability), reflecting the fine tuning between opposing forces.

Transition State Theory in Folding:
Protein folding kinetics can be described using an energy barrier separating the unfolded () and folded () states, with a high-energy transition state (TS or ) at the top of the barrier.

The folding rate depends on the height of the free energy barrier ():

• A lower barrier gives faster folding.
• The transition state represents the critical, partially structured folding nucleus.

-Value Analysis:
Developed by Alan Fersht, this experimental method probes the structure of the transition state at the residue level using site-directed mutagenesis.

Definition:

where values are the changes in free energy caused by a mutation.

Interpretation:

• : The residue is fully folded (native-like) in the transition state — its interactions are already formed at the TS.
• : The residue is fully unfolded in the transition state — it forms interactions only after the TS.
• : The residue is partially structured in the transition state.

Significance:

• Maps out which parts of the protein form structure early (folding nucleus) vs. late.
• Provides a residue-by-residue picture of the folding pathway.

Functional Design of Proteins:
Proteins are evolutionarily optimized not only for stability but for function, requiring a balance between structural integrity and dynamic flexibility.

Key Principles of Functional Design:

1. Marginal Stability:

• Proteins have only small net stability, to (equivalent to a few hydrogen bonds).
• This is by design, not accident.

2. Why Marginal Stability is Functionally Important:

• Flexibility: Enzymes must undergo conformational changes (induced fit, allostery) to bind substrates and catalyze reactions. Overly rigid proteins cannot do this.
• Regulation & Turnover: Marginally stable proteins can be more easily degraded and recycled, allowing regulation of protein levels.
• Dynamic motions: Breathing motions and local unfolding enable ligand binding, signal transduction, and catalysis.

3. Structure–Function Relationships:

• Active sites are often located in clefts formed by specific folding, using precisely positioned catalytic residues.
• Allosteric sites allow regulation distant from the active site.
• Binding specificity arises from shape and chemical complementarity.

4. Trade-offs:

• Increasing stability (e.g., extra disulfides) may reduce flexibility and impair function.
• Evolution tunes stability to the level needed for function under physiological conditions.

Conclusion:
Proteins are designed at the edge of stability — stable enough to maintain a defined structure yet flexible enough to perform dynamic biological functions. This delicate balance is central to their functional design.

Heat Capacity Change ():
Protein unfolding is accompanied by a large positive change in heat capacity ().

Molecular Origin:

• Upon unfolding, buried hydrophobic groups become exposed to water.
• Water forms ordered structures around these nonpolar groups, and this structured water has a high heat capacity.
• Hence, exposure of hydrophobic surface area is the primary contributor to positive .

Temperature Dependence of Thermodynamic Functions:
Because , both and vary with temperature:

Combining these gives the Gibbs–Helmholtz equation for stability:

Significance:

• The nonzero produces the characteristic curved (parabolic) stability curve of vs. .
• It is responsible for the existence of both heat and cold denaturation.
• The magnitude of is proportional to the hydrophobic surface area buried in the native state, linking thermodynamics to structure.

Molecular Chaperones:
Molecular chaperones are specialized proteins that assist in the correct folding and assembly of other proteins, preventing misfolding and aggregation, especially in the crowded cellular environment.

Major Classes:

• Hsp70 (DnaK): Binds exposed hydrophobic segments of nascent or unfolded chains, preventing premature aggregation. ATP-dependent.
• Chaperonins (Hsp60 / GroEL–GroES): Form barrel-shaped complexes providing an isolated chamber where a single polypeptide can fold protected from aggregation.
• Hsp90: Assists folding of specific signaling proteins.
• Small heat shock proteins: Hold unfolded proteins until refolding is possible.

Mechanism (e.g., GroEL/GroES):

1. Unfolded protein binds inside the GroEL cavity via hydrophobic interactions.
2. ATP and the GroES "cap" bind, enclosing the protein in a hydrophilic chamber.
3. The protein folds in isolation over a folding cycle (~10–15 s).
4. ATP hydrolysis releases the folded (or partially folded) protein.

Why Anfinsen's Principle is Not Violated:

• Chaperones do not provide structural information for the final fold.
• They do not dictate the native conformation — the sequence still determines the structure.
• They simply prevent unproductive interactions (aggregation, misfolding) and provide a favorable environment.
• They lower the kinetic barriers and increase the yield of correctly folded protein.

Conclusion:
Chaperones are kinetic assistants, not templates. The native state still corresponds to the thermodynamic minimum encoded in the sequence, fully consistent with Anfinsen's hypothesis.

Starting Point:
The standard Gibbs free energy change is related to the equilibrium constant by:

We also know:

Combining the two:

Dividing throughout by :

Differentiating with respect to (assuming and are approximately constant over the range):

This is the van't Hoff equation (differential form).

Integrated Form:

Application to Protein Unfolding:

• Measure the equilibrium constant (from fraction folded/unfolded) at several temperatures using spectroscopy or calorimetry.
• Plot vs. (a van't Hoff plot).
• The slope , giving the van't Hoff enthalpy .
• The intercept gives .

Comparison with Calorimetry:

• If , the transition is two-state.
• A ratio suggests intermediates; suggests intermolecular cooperativity/aggregation.

Protein Misfolding:
Misfolding occurs when a protein fails to reach its correct native conformation and instead adopts an alternative, non-functional structure, often becoming trapped in a local energy minimum on the folding landscape.

Causes:

• Mutations altering the sequence.
• Errors in translation or post-translational modification.
• Environmental stress (temperature, pH).
• Failure of chaperone systems.

Consequences:

• Loss of biological function.
• Exposure of hydrophobic surfaces normally buried in the core.
• Tendency to aggregate with other misfolded molecules.

Aggregation and Amyloid Formation:

• Misfolded proteins can associate into ordered, insoluble aggregates called amyloid fibrils.
• Amyloids are characterized by a cross- sheet structure, where -strands run perpendicular to the fibril axis.
• These structures are highly stable and resistant to degradation.

Thermodynamic Perspective:

• Amyloid states can be more thermodynamically stable than the native state in some cases (deep off-pathway minima).
• The native functional state may be only kinetically favored, not the global minimum.

Associated Diseases (Amyloidoses / Conformational Diseases):

• Alzheimer's disease — amyloid- plaques and tau tangles.
• Parkinson's disease — -synuclein aggregates.
• Prion diseases (Creutzfeldt–Jakob) — misfolded PrP propagates its conformation.
• Type II diabetes, Huntington's disease, and others.

Conclusion:
Protein misfolding illustrates the delicate balance of the folding landscape; disruption leads to aggregation and a range of serious degenerative diseases.

Several models have been proposed to describe the sequence of events during protein folding. They differ in the order in which secondary and tertiary structures form.

1. Framework Model (Diffusion-Collision Model):

• Secondary structure forms first, independently and locally (-helices, -hairpins).
• These preformed elements then diffuse and collide, coalescing into the tertiary structure.
• Folding is hierarchical: local → global.

2. Hydrophobic Collapse Model:

• The protein first undergoes a rapid, non-specific collapse driven by the hydrophobic effect, expelling water and forming a compact globule.
• Tertiary contacts (collapse) precede the full development of secondary structure.
• Secondary and tertiary structures then rearrange within the collapsed state (molten globule).

3. Nucleation-Condensation Model:

• A combined/intermediate mechanism.
• A small set of key residues form a folding nucleus where secondary and tertiary interactions develop simultaneously and cooperatively.
• The rest of the structure rapidly condenses around this nucleus.
• Secondary structure is stabilized by forming tertiary contacts at the same time.

Comparison Table:

Current View:
The nucleation-condensation model is most widely supported for many small proteins, with framework and collapse mechanisms representing limiting cases of a general unified mechanism.

Chemical Denaturants:

• Urea and guanidinium chloride (GdmCl) are common denaturants that unfold proteins by weakening the hydrophobic effect and preferentially solvating the unfolded state (interacting favorably with the exposed peptide backbone).

Denaturation Curve:

• The protein is exposed to increasing concentrations of denaturant, and the fraction unfolded is monitored using a spectroscopic probe (e.g., fluorescence, circular dichroism, absorbance).
• The result is a sigmoidal unfolding transition with pre-transition (folded), transition, and post-transition (unfolded) regions.

Determining Fraction Folded:
From the observed signal ():

and the equilibrium constant:

Linear Extrapolation Method (LEM):

• The free energy of unfolding varies linearly with denaturant concentration :
  

where:

• = stability in the absence of denaturant (the desired quantity).
• = the m-value, reflecting the change in solvent-accessible surface area upon unfolding.

Procedure:

1. Calculate at each denaturant concentration within the transition region.
2. Plot vs. .
3. Extrapolate the straight line to to obtain .

Midpoint ():
At the midpoint concentration, :

This method provides a robust measure of intrinsic protein stability.

Thermodynamics of Nucleic Acid Structures:
The stability of the DNA double helix and RNA secondary structures is governed by a balance of enthalpic and entropic contributions.

1. Hydrogen Bonding (Base Pairing):

• A–T (A–U) pairs form 2 hydrogen bonds; G–C pairs form 3 hydrogen bonds.
• Provides specificity of base pairing.
• Contributes to stability, but like proteins, the net contribution is modest because bases also H-bond with water in the single-stranded state.

2. Base Stacking Interactions:

• The dominant stabilizing force in nucleic acids.
• Arises from van der Waals and hydrophobic/dispersion interactions between adjacent stacked bases.
• Stacking is favorable and depends on sequence — G–C rich regions stack more strongly.

3. Electrostatic Effects:

• The negatively charged phosphate backbone causes repulsion, destabilizing the helix.
• Counterions (Na, Mg) shield these charges, stabilizing the duplex. Hence stability increases with ionic strength.

Melting Temperature () of DNA:

• is the temperature at which half of the double-stranded DNA has separated into single strands (denatured).
• Monitored by the hyperchromic effect — increased UV absorbance at 260 nm upon melting.

Factors Affecting :

• GC content: Higher GC → higher (3 H-bonds + stronger stacking).
  
• Ionic strength: Higher salt → higher .
• Length: Longer duplexes → higher .
• pH and denaturants also affect stability.

Thermodynamic Relation:

Conclusion:
Nucleic acid stability is primarily a balance between favorable base stacking and hydrogen bonding versus unfavorable electrostatic repulsion, with serving as a key experimental measure of structural stability.