Free-energy relationship

A free energy relationship, also known as a Gibbs energy relationship, relates the Gibbs free energy change ΔG of one series of chemical reactions to the ΔG of the same or another, related series of reactions [1]. Here, ΔG represents the difference in Gibbs energy between two states of a reaction.

Such reaction series are typically generated by systematic structural modifications, such as changing substituents on reactants or substrates, or by introducing amino acid mutations into enzyme proteins. Additionally, reaction conditions such as the type of solvent may also be changed. In many cases, the relationship between two Gibbs free energy changes can be approximated by a linear correlation. This linear relationship is known as a linear free energy relationship [2]. It is often abbreviated as LFER.

In a broad sense, free-energy relationships include correlations between two ΔG° (the standard free energy difference between two stable states) values and also correlations between two ΔG (the activation free energy) values. However, the term is most commonly used to describe correlations between ΔG and ΔG°. When expressed in terms of experimentally accessible quantities, this corresponds to a correlation between the logarithm of a reaction rate constant (log k) and the logarithm of a equilibrium constant (log K). This form is known as the rate–equilibrium free energy relationship (REFER).

Using REFER, a linear free energy relationship can be written as

or relative to a reference r,

which may also be expressed as

where Δlog k represents the difference in log k value relative to the reference log kr. In this form, the relationship is a proportional one that passes through the origin, emphasizing the linear dependence between the two free energy changes.

The equilibrium constant K is determined solely by the free energy difference between the initial and final stable states and remains unchanged regardless of the reaction path. In contrast, the rate constant k depends on the reaction path and the transition state. Consequently, a simple correlation between k and K is not generally expected. Thus, a linear correlation in the REFER plot highlights the significance of linear free energy relationships.

Linear free energy relationships are often used to estimate rate or equilibrium constants since they are experimentally difficult to determine [3].

In the classical sense, free energy relationships establish the extent at which bond formation and breakage happen in the transition state of a reaction. Establishing free energy relationships allows us to understand the reaction mechanism. For this purpose, it is effective when combined with kinetic isotope experiments.

Although the International Union of Pure and Applied Chemistry (IUPAC) recommends the term linear Gibbs energy relationship in place of linear free energy relationship, the latter remains widely used in the literature [1].

Chemistry

Well-known examples of linear free energy relationships in chemistry include the Brønsted relation [4] and Hammett relation [5] used in chemical kinetics. The Brønsted catalysis equation describes the relationship between the ionization constant of a series of catalysts and the reaction rate constant for a reaction on which the catalyst operates. The Hammett equation predicts the equilibrium constant or reaction rate of a reaction from a substituent constant and a reaction type constant.

Beyond these classical examples, the Edwards equation relates nucleophilicity to both polarizability and basicity [6][7], while the Marcus theory explains the rates of electron transfer reactions – the rate at which an electron can move or jump from one chemical species (called the electron donor) to another (called the electron acceptor) [8][9]. The Marcus theory provides an example of a quadratic free energy relationship (QFER), in which the reaction barrier depends quadratically on the reaction free energy.

Physical Chemistry (partition between two phases)

A typical LFER relation for predicting the equilibrium concentration of a compound or solute in the vapor phase to a condensed (or solvent) phase can be defined as follows (following M.H. Abraham and co-workers):[10][11]

where SP is some free-energy related property, such as an adsorption or absorption constant, log K, anesthetic potency, etc. The lowercase letters (e, s, a, b, l) are system constants describing the contribution of the aerosol phase to the sorption process.[12] The capital letters (E, S, A, B, L) are solute descriptors representing the complementary properties of the compounds. Specifically,

  • L is the gas–liquid partition constant on n-hexadecane at 298 K;
  • E = the excess molar refraction (E = 0 for n-alkanes).
  • S = the ability of a solute to stabilize a neighbouring dipole by virtue of its capacity for orientation and induction interactions;
  • A = the solute's effective hydrogen bond acidity; and
  • B = the solute's effective hydrogen-bond basicity.

The complementary system constants are identified as

  • l = the contribution from cavity formation and dispersion interactions;
  • e = the contribution from interactions with solute n-electrons and pi electrons;
  • s = the contribution from dipole-type interactions;
  • a = the contribution from hydrogen-bond basicity (because a basic sorbent will interact with an acidic solute); and
  • b = the contribution from hydrogen-bond acidity to the transfer of the solute from air to the aerosol phase.

Similarly, the correlation of solvent–solvent partition coefficients as log SP, is given by

where V is McGowan's characteristic molecular volume in cubic centimeters per mole divided by 100.

Biochemistry

In biochemistry, linear free energy relationships are widely used to analyze enzyme reaction mechanisms using series of structurally related substrates [13] or systematic amino-acid substitutions within enzyme proteins [14].

Biophysics

LFER analysis is also applied to the folding and unfolding reactions of protein molecules. In the presence of high concentrations of denaturants such as urea and guanidine hydrochloride, protein molecules are denatured and do not form stable structures. When the denaturant is removed, the protein molecules return to their native structures in a short timescale of microseconds to milliseconds. This refolding process is called protein refolding.

The φ-value analysis, based on single–amino acid mutations, is a well-established method for probing the transition state of protein refolding reactions [15]. The φ value represents the extent to which the native-state interactions around the side chain of a mutated residue are formed in the transition state ensemble [16]. In an idealized case, φ values are either 0 or 1, corresponding to completely unformed or fully formed native interactions, respectively; however, experimentally determined values typically fall between 0 and 1.

Residue-based LFER is a special form of linear free-energy relationship in which multiple amino acid residues within a single polypeptide chain are monitored as local probes under single experimental conditions [17]. This approach is applicable to a wide range of protein conformational changes, from small-scale fluctuations to large structural changes, such as coupled binding and folding of intrinsically disordered proteins (IDPs) to target proteins [18].

Residue-based LFER can be derived mathematically from the consistency principle of protein folding [19], originally proposed by Nobuhiro Gō in 1983 [20]. Both residue-based LFER and its generalized quadratic free energy relationship (QFER) provide residue-specific information about the transition states of protein structural changes. The information obtained is equivalent to that obtained from the φ-value analysis [21]. Combining local probe measurements with amino-acid substitutions further extends the applicability of residue-based LFER and QFER approaches [21].

References

  1. ^ a b IUPAC, Compendium of Chemical Terminology, 5th ed. (the "Gold Book") (2025). Online version: (2006–) "linear free-energy relation". doi:10.1351/goldbook.L03551
  2. ^ Jaffeé H. "A reëxamination of the Hammett equation". Chem Rev 53, 191-261 (1953). doi:10.1021/cr60165a003
  3. ^ Lassila JK, Zalatan JG, Herschlag D (2011-06-15). "Biological phosphoryl-transfer reactions: understanding mechanism and catalysis". Annual Review of Biochemistry. 80 (1): 669–702. doi:10.1146/annurev-biochem-060409-092741. PMC 3418923. PMID 21513457.
  4. ^ J. N. Brönsted and K. Pedersen. "Die katalytische Zersetzung des Nitramids und ihre physikalisch-chemische Bedeutung". Z. Phys Chem 108, 185 (1924). doi:10.1515/zpch-1924-10814
  5. ^ Hammett, L.P. "Some relations between reaction rates and equilibrium constants". Chem Rev 17, 125-136 (1935). doi:10.1021/cr60056a010
  6. ^ Edwards, J.O. (1954). "Correlation of Relative Rates and Equilibria with a Double Basicity Scale". Journal of the American Chemical Society. 76 (6): 1540–1547. doi:10.1021/ja01635a021.
  7. ^ Edwards, J.O. (1956). "Polarizability, Basicity and Nucleophilic Character". Journal of the American Chemical Society. 78 (9): 1819–1820. doi:10.1021/ja01590a012.
  8. ^ Marcus, R.A. "On the Theory of Oxidation-Reduction Reactions Involving Electron Transfer I" J Chem Phys 24, 966 (1956). doi:10.1063/1.1742723
  9. ^ Marcus.R.A. "Electrostatic Free Energy and Other Properties of States Having Nonequilibrium Polarization I. J Chem Phys 24, 979. doi:10.1063/1.1742724
  10. ^ Abraham MH, Ibrahim A, Zissimos AM, Zhao YH, Comer J, Reynolds DP (October 2002). "Application of hydrogen bonding calculations in property based drug design". Drug Discovery Today. 7 (20): 1056–63. doi:10.1016/s1359-6446(02)02478-9. PMID 12546895.
  11. ^ Poole CF, Atapattu SN, Poole SK, Bell AK (October 2009). "Determination of solute descriptors by chromatographic methods". Analytica Chimica Acta. 652 (1–2): 32–53. doi:10.1016/j.aca.2009.04.038. PMID 19786169.
  12. ^ Bradley JC, Abraham MH, Acree WE, Lang AS (2015). "Predicting Abraham model solvent coefficients". Chemistry Central Journal. 9: 12. doi:10.1186/s13065-015-0085-4. PMC 4369285. PMID 25798192.
  13. ^ Alnajjar, K. S., Wang, K., Alvarado-Cruz, I., Chavira, C., Negahbani, A., Nakhjiri, M., et al. "Modifying the Basicity of the dNTP Leaving Group Modulates Precatalytic Conformational Changes of DNA Polymerase β". Biochemistry 63, 1412–1422 (2024) doi:10.1021/acs.biochem.4c00065
  14. ^ Michałowski, M. A., Terejko, K., Gos, M., Iżykowska, I., Czyżewska, M. M., Kłopotowski, K., et al. "Φ value analysis underscores strong functional and structural compactness of the GABA A receptor". Proc Natl Acad Sci, U.S.A. 122, (2025) doi:10.1073/pnas.2512278122
  15. ^ Fersht, A. R. "From covalent transition states in chemistry to noncovalent in biology: from β- to Φ- value analysis of protein folding". Q Rev Biophys 57, (2024). doi:10.1017/S0033583523000045
  16. ^ Fersht, A. R., Matouschek, A., Serrano, L. "The folding of an enzyme. I. Theory of protein engineering analysis of stability and pathway of protein folding". J Mol Biol 224, 771–782 (1992). doi:10.1016/0022-2836(92)90561-W
  17. ^ Fujinami, D., Hayashi, S., Kohda, D. "Residue-Specific Kinetic Insights into the Transition State in Slow Polypeptide Topological Isomerization by NMR Exchange Spectroscopy". J Phys Chem Lett 12, 10551–10557 (2021). doi:10.1021/acs.jpclett.1c02387
  18. ^ Fujinami, D., Hayashi, S., Kohda, D. "Retrospective study for the universal applicability of the residue-based linear free energy relationship in the two-state exchange of protein molecules". Sci Rep 12, 16843. (2022). doi:10.1038/s41598-022-21226-z
  19. ^ Kohda, D., Hayashi, S., Fujinami, D. "Residue-based correlation between equilibrium and rate constants is an experimental formulation of the consistency principle for smooth structural changes of proteins". Biophys Physicobiol 20, e200046 (2023). doi:10.2142/biophysico.bppb-v20.0046
  20. ^ Gō, N. "Theoretical Studies of Protein Folding". Annu Rev Biophys Bioeng 12, 183–210 (1983). doi:10.1146/annurev.bb.12.060183.001151
  21. ^ a b Kohda, D., Hayashi, S., Fujinami, D. "NMR evidence for residue-based LFER relationship in two-state folding-unfolding equilibrium of the spectrin SH3 domain in acidic solutions". Biophys Physicobiol 23, e230011 (2026). doi:10.2142/biophysico.bppb-v23.0011

See also