Hans V. Westerhoff
Sessions
The history of Biology abounds in discoveries where the integral of precise theory and precise experimentation, the trade mark of the ISGSB, has made the difference. These discoveries include no less than the structure of DNA and the mechanisms underlying biological free-energy transduction and the multiplicity of oncogenes. This opening lecture of the 2024 International Study Group on Systems Biology (ISGSB) will sketch the essence and history of ISGSB/BTK. The essence will be identified as ISGSB’s intensive and informal discussions of controversial issues in biology. The latter will be described not in terms of numbers, but in terms of leaps in (my) understanding, produced by ISGSB’s core methodologies.
Metabolic networks need to meet more requirements than single enzymes, in order to be functional. Aware of the heterogeneity of (tumor) cell populations, we went after this principle and engaged in network based drug design. We thought that by identifying the metabolic potential of individual cells, we could identify which targets could be used to most effectively incapacitate most individuals of a tumor cell population.
We projected mRNA sequence counts obtained for >3000 cells out of a growing tumor-cell population onto the genome-wide metabolic map after converting the numbers to Vmax’s. We used Flux Balance Analysis to predict the pathways the individual cells could be using and thereby their vulnerabilities to potential metabolic drugs. That is at least what we thought we would do.
Much to our surprise however, none of the cells was predicted to be able to grow.
We then considered whether this could be due to the cells being social metabolically, i.e. massively exchanging metabolites, with some cells taking care of the upper part of glycolysis, others the TCA cycle, yet others the lower part of glycolysis. This was a nice and social idea, but apparently not realistic: subdividing the cells into subpopulations and offering metabolites synthesized by one subpopulation as substrates to the others, did not lead to growth of either subpopulation.
In this presentation we shall discuss what explanation of the growth in the absence of mRNA for the metabolic pathways, we did come up with.
And, we discuss how the actual resulting model did identify cholesterol and asparagine synthesis pathways as relevant, though complex, drug targets.
In physical chemistry, Gibbs free energies and chemical potentials (not energies) are the descriptors for the energetics of compounds and processes. They integrate energy, entropy and volume-work. Changes in Gibbs energy and between chemical potentials equal the useful work at constant temperature and pressure. However, this systematics requires one to be explicit about the protonation, Mg-complexation and hydration states of the molecules, and to correct of activity coefficients differing from 1. Moreover the standard chemical potentials are defined for biologically irrelevant situations such as pH=pMg=0, concentrations of 1 Molar, hydrogen gas at 1 atmosphere and crystalline phosphorous. In addition, one must monitor the number of protons and Mg ions liberated or consumed by the reaction. Rather than a single ATP synthesis reaction, this methodology requires describing some 20 reaction variants, explicating the various protonation, Mg-complexation and hydration states of the three molecules involved. Indeed, multiple literature studies have presented highly complex ways of calculating the standard Gibbs energy of this highly important reaction, which then typically differed by 17 kJ/mol from its actual work potential.
We here present new ‘metabolic’ energies and potentials. These should replace the Gibbs free energy and the chemical potential for life processes and metabolites under in-vivo like conditions. We also present a new ‘Thermotable’, which contains the standard metabolic potentials for up to a thousand metabolites of interest to systems biology. The Thermotable enables the direct computation of standard (i.e., concentration-independent) reaction energies by simple subtraction. These are immediately relevant for the in vivo reference conditions of pH=7, pMg=3, ionic strength 0.15 M, and T=310 K, and concentrations of 1 mM. There is no need for correcting concentrations to activities. The new standard reaction energies are very good approximations towards the actual reaction energies, for when concentrations are unknown (yet on the order of 1 mM rather than 1 M).
The metabolic potentials given in the Thermotable are much better descriptors of the energetic potential of biochemicals than the chemical potentials were, because they are taken relative to a growth medium of biological significance, i.e., 1 mM of (total) bicarbonate, ammonium, phosphate, and sulfate, as well as liquid water, H+ at pH7, and Mg2+ at pMg3. This will be illustrated by plotting the metabolic energy landscape for travelling down major metabolic pathways. Whilst the maps of chemical potentials were rugged and irregular, the magnitudes of the potentials are realistic (in terms of numbers of ATP energies) and the new maps are smoothly down-hill except for steps up of about 50 kJ/mol where ATP energy is invested (and down where ATP is made). With the new metabolic energies, the calculation of free energy differences and equilibrium constants becomes facile, omitting the many usual points of confusion in calculating from standard chemical potentials.
The new thermodynamics will be useful for turning metabolic network maps into metabolic energy landscapes and for warning against perpetua mobilia proposed by Flux Balance Analysis.