Group additivity calculation of the standard molal thermodynamic properties of aqueous amino acids, polypeptides and unfolded proteins as a function of temperature, pressure and ionization state

Abstract

Abstract. Thermodynamic calculation of the chemical speciation of proteins and the limits of protein metastability affords a quantitative understanding of the biogeochemical constraints on the distribution of proteins within and among different organisms and chemical environments. These calculations depend on accurate determination of the ionization states and standard molal Gibbs free energies of proteins as a function of temperature and pressure, which are not generally available. Hence, to aid predictions of the standard molal thermodynamic properties of ionized proteins as a function of temperature and pressure, calculated values are given below of the standard molal thermodynamic properties at 25°C and 1 bar and the revised Helgeson-Kirkham-Flowers equations of state parameters of the structural groups comprising amino acids, polypeptides and unfolded proteins. Group additivity and correlation algorithms were used to calculate contributions by ionized and neutral sidechain and backbone groups to the standard molal Gibbs free energy (Δ G°), enthalpy (Δ H°), entropy (S°), isobaric heat capacity (C°P), volume (V°) and isothermal compressibility (κ°T) of multiple reference model compounds. Experimental values of C°P, V° and κ°T at high temperature were taken from the recent literature, which ensures an internally consistent revision of the thermodynamic properties and equations of state parameters of the sidechain and backbone groups of proteins, as well as organic groups. As a result, Δ G°, Δ H°, S° C°P, V° and κ°T of unfolded proteins in any ionization state can be calculated up to T~-300°C and P~-5000 bars. In addition, the ionization states of unfolded proteins as a function of not only pH, but also temperature and pressure can be calculated by taking account of the degree of ionization of the sidechain and backbone groups present in the sequence. Calculations of this kind represent a first step in the prediction of chemical affinities of many biogeochemical reactions, as well as of the relative stabilities of proteins as a function of temperature, pressure, composition and intra- and extracellular chemical potentials of O2 and H2, NH3, H2PO4 and CO2.