Standard Gibbs Energy of Metabolic Reactions: I. Hexokinase

Sep 21, 2016 - Florian Meurer, Maria Bobrownik, Gabriele Sadowski, and Christoph Held ... The standard Gibbs energy of reaction enables calculation of...
0 downloads 0 Views 742KB Size
Subscriber access provided by Northern Illinois University

Article

Standard Gibbs energy of metabolic reactions: I. Hexokinase reaction Florian Meurer, Maria Bobrownik, Gabriele Sadowski, and Christoph Held Biochemistry, Just Accepted Manuscript • DOI: 10.1021/acs.biochem.6b00471 • Publication Date (Web): 21 Sep 2016 Downloaded from http://pubs.acs.org on September 22, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Biochemistry is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

1 2

Standard Gibbs energy of metabolic reactions: I. Hexokinase

3

reaction

4

Florian Meurer, Maria Bobrownik, Gabriele Sadowski, Christoph Held*

5 6 7

Laboratory of Thermodynamics, Department of Biochemical and Chemical Engineering,

8

Technische Universität Dortmund, Emil-Figge-Str. 70, 44227 Dortmund, Germany

9 10 11 12 13 14 15 16 17 18 19 20

* corresponding author: [email protected]

21

ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

22

Symbols

23

Greek letters Symbol r

Page 2 of 37

Property

Unit

relative dielectric permittivity

-

association-energy parameter     

 ∗,

activity coefficient of component i on molality base

-

rational activity coefficient of component i on molality base

-

  

activity coefficient of component i at infinite dilution on molality base

-

association-volume parameter

-

stoichiometric coefficient of component i

-

segment diameter of component i

Å

osmotic coefficient of aqueous solution

-

 , 



24 25

K



Latin letters Symbol

Property

Unit

res

residual Helmholtz energy

J

hc

hard-chain contribution to Helmholtz energy

J

disp

dispersive contribution to Helmholtz energy

J

association contribution to Helmholtz energy

J

ionic contribution to Helmholtz energy

J

f-function of component i

-

chemical standard Gibbs energy of chemical reaction

J·mol

biochemical standard Gibbs energy of reaction

J·mol

Boltzmann constant

J·K

chemical thermodynamic equilibrium constant

-

biochemical thermodynamic equilibrium constant

-

A

A A

assoc

A

A

ion



∆  

∆ ′ kB

 

-1

-1

-1

ACS Paragon Plus Environment

Page 3 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

∗

chemical activity-coefficient ratio

-

biochemical activity-coefficient ratio

-



chemical molality-based equilibrium constant

-



biochemical molality-based equilibrium constant

-

molality of component i

mol·kg

molality of component i at thermodynamic equilibrium

mol·kg

segment number

-

molar mass of component i

g·mol

number of association sites of component i

-

 ∗  





mi

seg



assoc

Ni

!

osmolality, ! =   

-1

-1

-1

-1

osmol kg -1

-1

-1

-1

R

ideal gas constant, 8.31446 J·mol ·K

J·mol ·K

T

temperature

K

dispersion-energy parameter of component i

K

ui/kB

26 27

ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 37

28

Abstract

29

The standard Gibbs energy of reaction enables calculation of the driving force of a

30

(bio)chemical reaction. Gibbs energies of reaction are required in thermodynamic

31

approaches to determine fluxes as well as single reaction conversions of metabolic

32

bioreactions. The hexokinase reaction (phosphorylation of glucose) is the entrance step

33 34

of glycolysis, and thus its standard Gibbs energy of reaction (∆ ) is of great impact.

∆  is accessible from equilibrium measurements, and the very small concentrations of

35

the reacting agents cause usually high error bars in data reduction steps. Even worse,

36

works from literature do not account for the non-ideal behavior of the reacting agents

37 38

(activity coefficients were assumed to be unity), thus published ∆  values are not standard data. Consistent treatment of activity coefficients of reacting agents is crucial

39

for

40

measurements.

41

In this work equilibrium molalities of hexokinase reaction were measured with an

42

enzyme kit. These results were combined with reacting agents’ activity coefficients

43

obtained with the thermodynamic model ePC-SAFT. Pure-component parameters for

44

adenosine triphosphate (ATP) and adenosine diphosphate (ADP) were fitted to

45

the

accurate

determination

of

standard

Gibbs

energy

from

equilibrium

experimental osmotic coefficients (water+Na2ATP, water+NaADP). ∆  of the

46

hexokinase reaction at 298.15 K and pH 7 was found to be -17.83 ± 0.52 kJ·mol-1. This

47

value was compared with experimental literature data; very good agreement between

48 49

the different ∆  values was obtained by accounting for pH, pMg, and the activity coefficients of the reacting agents.

50

ACS Paragon Plus Environment

Page 5 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

51

Introduction

52

The activation of glucose by phosphorylation to glucose-6-phosphate (G6P) is the entry

53

step of glycolysis in biological cells. This reaction is catalyzed by the enzyme

54

hexokinase, and it contains adenosine triphosphate (ATP) as reactant that is reduced to

55

adenosine diphosphate (ADP) according to: Glucose + ATP ↔ G6P + ADP

(1)

Glucose + ATP3- ↔ G6P- + ADP2-

(2)

56

Eq. (1) shows the biochemical reaction while Eq. (2) displays the chemical expression

57

that contains the highest ionized species of each reacting agent present at biological

58

standard conditions (298.15 K and pH 7). Due to its importance for metabolism the

59

hexokinase reaction was subject to different investigations in the last decades. While

60

Fromm et al. in 1964 (1) and later Blázquez et al. (2) worked on enzyme kinetics and

61

inhibition of the hexokinase reaction using different hexokinases, Bianconi et al. (3)

62

studied the kinetics as well as reaction enthalpy of reaction (1), and they used two

63

isozymes of hexokinase for that.

64

The thermodynamic equilibrium of this hexokinase reaction has been subject to some

65

investigations (4, 5). The reaction equilibrium was found to be far on the product side,

66

and the reaction therefore has usually been presented as an irreversible reaction.

67

Nevertheless, hexokinase reaction was proven to be reversible. Studies on (standard)

68

enthalpy of reaction by calorimeter experiments were published by different authors (3,

69

5, 6). They applied titration calorimetry to determine enthalpies of reaction at different

70

conditions considering the influence of pH, pMg (6) and the influence of some inert

71

sugars (5). The conditions under investigation in those works are non-ideal from the

72

thermodynamic point-of-view with respect to deviations from the biological standard ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

73

state, which is a one molal hypothetically ideal solution. Deviations from this standard

74

state are accounted for by the activity coefficients of the reacting agents. However, all

75

previous works on hexokinase have assumed that the activity coefficients of the reacting

76

agents are unity. As demonstrated by Hoffmann et al. (7, 8) the equilibrium of biological

77

reactions depend on the reacting agents’ concentrations, even if the molalities of

78

reacting agents are very low (in the mmolal range). All data available on literature on the

79

hexokinase reaction were measured at conditions that yielded equilibrium molalities of

80

the reacting agents up to 47.21 mmol·kg-1. One exception are the conditions applied by

81

Robbins and Boyer (4) who determined equilibrium-composition data at very low

82

equilibrium molalities (max. 7.6 mmol·kg-1) of the reacting agents. Only at these highly

83

diluted conditions activity coefficients might be close to unity, whereas activity

84

coefficients are required for investigations at higher molalities.

85

Activity coefficients might be accessible from measurements (e.g. osmometry).

86

However, measuring activity coefficients of the reacting agents of the hexokinase

87

reaction depending on a broad variety of conditions (reacting agents’ molalities, ionic

88

strength, pH, temperature, presence of additives and cofactors) is not reasonable.

89

Application of a thermodynamic model is more convenient, especially since

90

thermodynamic models for biomolecules improved significantly. To give some examples,

91

gE-models like the Pitzer equation (9) as well as equations of state (EOS) like Peng-

92

Robinson EOS and (electrolyte) Perturbed-Chain Statistical Associating Fluid Theory

93

(ePC-SAFT) (10, 11) were applied to describe thermodynamic properties of biological

94

solutions. Using ePC-SAFT, excellent results were obtained on modeling activity

95

coefficients in solutions containing urea, amino acids, or adenine (12-14). Further, ePC-

96

SAFT was successfully applied to model thermodynamic equilibrium of bioreactions (7, ACS Paragon Plus Environment

Page 6 of 37

Page 7 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biochemistry

97

8, 15, 16). ePC-SAFT requires pure-component parameters and binary interaction

98

parameters. In (7, 8, 15, 16), the pure-component ePC-SAFT parameters for the

99

reacting agents were fitted to reaction-independent data e.g., thermodynamic water

100

activity data (7). This yielded good prediction results of activity coefficients in biological

101

multi-component systems. In contrast, the prediction of activity coefficients in multi-

102

component mixtures failed using classical models like (Pitzer)-Debye-Hückel (17). The

103

reason for this is the complex interaction behavior among the components present in the

104

mixtures. This yields a strongly non-ideal solution with activity coefficients distinctly

105

different from unity, which can only be predicted by state-of-the-art models (e.g. ePC-

106

SAFT) that account for short-range forces and long-range interactions.

107

The availability of equilibrium molalities and activity coefficients of the reacting agents

108 109 110 111

allows determining the standard Gibbs energy of the hexokinase reaction ∆ . Since

the entrance in the glycolytic cycle is an important step in catabolic activity of the cell, a thermodynamically-consistent ∆  value is of great impact for other scientists working on strand and production optimization in terms of metabolic engineering and

112

thermodynamic metabolic flux analysis (18, 19). All previous equilibrium studies on

113

hexokinase have assumed that the activity coefficients of the reacting agents are unity

114

despite that fact that the reactions were carried out at non-standard conditions. Thus,

115 116 117 118

∆ 

values calculated from these results are just observed values but not consistent

∆  values at standard state. Thus, in this work, the reaction equilibrium of the

hexokinase reaction was measured, activity coefficients of the reacting agents were

modeled, and ∆  of the hexokinase reaction at 298.15 K and pH 7 is derived as the

119

main result. These values will further be compared to experimental reaction data from

120

literature in the very last section of the results chapter. This will show the importance of ACS Paragon Plus Environment

Biochemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 37

121

consistent treatment of activity coefficients of reacting agents for the accurate

122

determination of Gibbs energy at standard state.

123

124 125 126 127 128

129 130 131 132 133 134 135 136

137 138

Thermodynamic formalism for hexokinase reaction

The Gibbs energy of reaction ∆ g is a fundamental property for the characterization of

reactions. It allows predicting whether or not a reaction will occur at defined conditions. At thermodynamic equilibrium ∆ g will be zero. Therefore at equilibrium Eq. (3) will be

fulfilled.

∆  = −RT ln* +

(3)

Here, R is the ideal gas constant with value of 8.31446 J·mol-1·K-1 and T is the absolute temperature and  is the thermodynamic equilibrium constant. Many compounds in biological reactions are electrolytes. They dissociate to a certain extent depending on solution conditions and may even form complexes with metal cations e.g., with Mg2+.

The chemical definition of  is formulated in terms of species and would therefore

require single species activities. Thus, the equilibrium constant  belongs to Eq. (2). In contrast, the biochemical equilibrium constant  refers to Eq. (1). In this work,  was

investigated, yielding the biochemical standard Gibbs energy of reaction ∆ ′ ∆ ′ = −RT ln * ′+

(4)

The thermodynamic equilibrium constant  ′ for the hexokinase reaction (Eq. (1)) is thus defined via the sum of all species activities Σaeq at reaction equilibrium according to:   ∑ -./0 ∙ ∑ -20  ′ =   ∑ -345678 ∙ ∑ -90

ACS Paragon Plus Environment

(5)

Page 9 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

139 140

Biochemistry

As the aeq values are defined as product of molality and activity coefficient on molality   basis  ∗, ,  ′ can be expressed as product of 

and ∗ via

 ′ = ′ ∙

141

 ∗

  ∑ ./0 ∙ ∑ 20





./0 ∙ 20 =   ∙ ∗



∙ 90 ∑ 345678 ∙ ∑ 90 345678

(6)

In Eq. (6)  ∗, values have to be calculated at equilibrium molalities   . The

142

assumption behind Eq. (6) is that the activity coefficients of all species of one

143

component are equal. The reference state in biological systems is the one molal

144 145 146

hypothetically ideal solution, in which the activity coefficients  ∗, are assumed to be

one. At conditions other than the biological standard state,  ∗, differ from unity. The rational activity coefficients can be obtained according to:  ∗,

147 148 149

150 151 152 153 154 155



= :, 

(7)

where  and  , are the molality-based activity coefficients of component i at any

given molality and at infinite dilution in the reaction medium. With this definition the formulation of the thermodynamically-consistent calculation of ∆ ′ is expressed as:  +  ∆ ′ = −RT ln*

− RT ln;∗