You can search for a compound or enzyme by entering their common name into the search box. For example, try "Succinyl-CoA" or "pyruvate formate-lyase." As you type, you should see a list of potential matches which should give you a sense of what eQuilibrator thinks you are looking for. To the right of each suggestion you will see they type of the match: "Compound" or "Enzyme". Complete your search by selecting a compound or enzyme from the list or typing the name in full. When you've finished typing your query, click 'Search' or press 'Enter' to see search results.
The search results page shows the top matches for your query. Results are shown in two sections: one for compounds and another for enzymes. If you're interested in a compound, click on it to see more details and to find ΔfG. If you're interested in an enzyme, click on the enzyme name to see detailed information about it. If you want to explore the thermodynamics of an enzyme-catalyzed reaction, choose a reaction listed under the enzyme name to get a ΔrG estimation.
You can search for a reaction by entering an (almost) free-text string of the reaction. For example, you can type "glucose = 2 ethanol + 2 co2".
In general, write your reaction using " = " or " <=> " to separate the substrates from the products. The spaces around " = " are required. Always use spaces before and after compound names, "+" and stoichiometric coefficients.
Also note that it is possible to find reactions through the enzymes that catalyze them. For example, instead of typing in the full reaction "Fructose 1,6-bisphosphate = Glycerone phosphate + Glyceraldehyde 3-phosphate" you can instead search for "aldolase" or "fructose bisphosphate aldolase."
"Standard conditions" implies that all substrates and products of a reaction are at 1 M concentration. If their concentrations are different, it will affect the ΔrG'.
For gases the standard condition is defined as 1 atmosphere (bar) partial pressure. However, if one knows the soluble concentration of the gas of interest it should be specified by choosing "custom" concentrations. Alternatively, if you know the partial pressure of a reactant in the gas phase of the reaction chamber, and assume there is an equilibrium with the solution phase, then you can select a "custom" concentration for the gas (in units of mbar). You can also indicate that you want to use the standard gas phase for the ΔG'° by adding (g) to the end of the compound name. For example, try to search for: CO2(g) + PEP + H2O = Oxaloacetate + Pi. This will work only for compounds for which the formation energy in gas phase is in our database, namely O2, N2, H2, CO2, and CO.
G is the Gibbs free energy of a physical system. It quantifies the amount of energy available for work in a system that is at constant temperature and pressure. Biological systems do not often experience radical changes in temperature or pressure, so G is appropriate for their analysis. On eQuilibrator we assume that the temperature is 25 °C (298.15 Kelvin) and the pressue is 1 bar. As a consequence of the second law of thermodynamics, spontaneous changes in systems at constant temperature and pressure must lead to a reduction in G. Chemical reactions change biochemical systems, so the same must be true of them: in order to occur, a reaction must lead to a reduction in the Gibbs free energy - ΔrG' < 0.
ΔrG° is the change in Gibbs free energy due to a chemical reaction in standard conditions and without accounting for pH, ionic strength or any other cellular factors (i.e. they are taken as 0). ΔrG'° is the change in Gibbs free energy due to a chemical reaction due a reaction at a particular pH and ionic strength. These are convenient reference points, but neither ΔrG° nor ΔrG'° account for the effect of reactant concentrations on the free energy.
In cells there are certain prevailing metabolite concentrations. Concentrations affect the entropic component of the free energy: if substrate concentrations are higher than product concentrations then forward flow of the reaction will balance the concentrations and increase the entropy of the system. ΔrG' accounts for the effect of concentrations. It is very important to use ΔrG' values derived from physiologically plausible concentrations in thermodynamic analysis because reactant concentrations have a substantial effect on ΔrG' and because the laws of thermodynamics constrain only ΔrG', not ΔrG° or ΔrG'°.
The formation energy of a compound is the change in free energy due to the compound's formation reaction - the reaction forming it from elements in their standard form. For example, the formation reaction of carbon dioxide is O2 + C ⇌ CO2. ΔfG° is the formation energy of a compound in standard conditions and without accounting for pH, ionic strength or any other cellular factors (i.e. they are taken as 0).
In contrast, ΔfG'° is the formation energy of a compound at a particular set of cellular conditions. When we show a ΔfG'° on eQuilibrator, as we do on compound and reaction pages, it has been transformed to the pH and ionic strength values. The default pH is 7 and the default ionic strength is 0.1 M.
As explained above, the reaction Gibbs energy depends on the concentrations of the reactants (and similarly for the formation energy ΔfG' and the reduction potential E'). It is useful to have a standard for these concentrations in order to compare reactions without specifying a concentration explicitly. It is standard to use 1 M concentrations to compare reactions in solution. However, in the context of metabolic reactions inside living cells, 1 mM is a much more appropriate standard concentration because metabolite concentrations typically range from 1 nM to 10 mM and never exceed 100 mM. As such, 1 M standard concentrations are entirely non-physiological.
The notation for the reaction Gibbs energy in the standard 1M concentrations is ΔrG'°, i.e. the degree sign (°) represents 1M. For our physiological standard we use the superscript 'm' to mark the 1mM concentration that is used for all reactants: ΔrG'm, ΔfG'm, and E'm.
There are several reasons that a ΔG value might not be available.
In the case of a reaction, ΔrG is not available when ΔfG is missing for one of the reactants. ΔfG is very difficult to measure experimentally. We use the group contribution method to estimate ΔfG for compounds that have not been characterized experimentally [2,10]. However, we need certain data about a compound in order to approximate ΔfG.
Our list of compounds is taken from the Kyoto Encyclopedia of Genomes and Genes (KEGG) . Occasionally, KEGG compound entries don't contain an InChI identifier. Without an InChI identifier we cannot compute the structure of the compound and without structural information we can't estimate ΔfG.
In some cases, a compound may contain a chemical group that was not observed in any of the compounds that have experimental ΔfG measurements. We cannot produce estimates for these compounds as we cannot estimate the contribution of unknown groups. In other cases, our algorithm for analyzing compound structure fails to decompose the compound into groups. If we are unable to decompose the compound then we cannot use the group contribution method.
Moreover, we can't estimate ΔrG for any reaction that contains a compound that we cannot estimate ΔfG for.
In order to fully understand how to calculate estimation uncertainties, you'll probably need to read our paper on the Component Contribution method . The short answer would be that we ran a cross-validation benchmark using a set of reactions for which the ΔrG' has been measured. Any reaction that you type in, is decomposed into compounds and these compounds are decomposed into groups. By comparing this decomposition vector with the ones from our database, we can estimate the ΔrG'. Along the way, we can also evaluate how good our estimation is, by checking how good we were for similar reactions in our benchmark.
eQuilibrator uses the "Alberty method" for biochemical thermodynamics. In the Alberty method, H+ is defined to have 0 free energy [6,8]. Instead of correcting for H+ concentration, a pH correction accounts for the abundance of H+. You can use the pH "slider" to see the effect of H+ concentration on your reaction.
Biochemical systems are generally assumed to be aqueous environments [6,8]. Therefore, the concentration of water is fixed.
The temperature is fixed at 25 °C (298.15 K) for all ΔG values given. The group contribution method enables us to approximate ΔfG of compounds at a particular temperature (the temperature at which they were measured) . As the change in free energy is defined as ΔG = ΔH - TΔS and we don't know the value of ΔS in most cases, we cannot predict how changes in temperature will affect ΔfG.
CO2 in solution gives rise to several chemical species. It can be quite confusing to think about the equilibrium between these species - doing so requires care. CO2(aq) is dissolved CO2. CO2(aq) undergoes a spontaneous hydration reaction to form carbonic acid:
Several hydrated species form in water through addition or release of protons: carbonic acid (H2CO3), bicarbonate (HCO3-) and carbonate (CO32-). In thermodynamics of biochemical reactions, different ionic states (known as pseudo-isomers) are lumped together. If you search for any of these hydrated species, eQuilibrator will use their lumped form - HCO3-(aq). Sometimes it is easier to measure or analyze the sum of CO2(aq) and its three hydrated forms (H2CO3, HCO3- and CO32-). This sum of species is termed CO2(total). Note that in the chemical formula of CO2(total) there are actually 3 oxygen atoms because it also includes the hydrating water molecule.
When only the total concentration is known, it is assumed that there is equilibrium among the four species and one uses CO2(total) in place of CO2(aq). If, however, you know or can measure the concentration of CO2(aq) alone, then it is reasonable to use CO2(aq) as a reactant. The concentration of CO2(aq) is usually straightforward to calculate based on Henry’s law dictating, for example, that under standard atmospheric conditions of about 400ppm CO2(g) the concentration of CO2(aq) is about 10 uM. The concentrations of bicarbonate and CO2(total), however, depend strongly on pH. More information is supplied in the figure below and in this link.
The uncatalyzed hydration reaction (CO2(aq) + H2O ⇌ H2CO3) takes minutes to equilibrate. In many organisms, however, this reaction is catalyzed by the enzyme carbonic anhydrase, which speeds up the reaction by several orders of magnitude . In cells, therefore, CO2(aq) is generally considered to be in equilibrium with its hydrated forms (carbonic acid, bicarbonate and carbonate) save in some special cases such as in cyanobacterial carbon concentrating mechanisms where carbonic anhydrase is absent from some parts of the cell.
We note that anaplerotic reactions use bicarbonate as their substrate whereas decarboxylation reactions release CO2 but in finding the ΔG' they can be written using either CO2(aq), HCO3-(aq) or CO2(total) as long as the concentrations used are accurate. This is true because of the equilibrium among these species.
We know that this whole issue is quite confusing. You are welcome to suggest how to explain this better.
A half-reaction can be the oxidation or reduction component of a redox reaction, without the other component. When you search for such a reaction, eQuilibrator recognizes that the number of electrons is not balanced and automatically switches to 'half-reaction' mode. Without knowing the other half, the change in Gibbs energy is not well defined. The parameter that is used to describe the potential difference (in Volts) between the products and substrates of a half-reaction is called redox potential and is marked by E'°. The redox potential is equal to the voltage at equilibrium under standard conditions of an electrochemical cell in which the cathode reaction is the half-reaction considered, and the anode is a standard hydrogen electrode where hydrogen is oxidized: ½ H2 ⇌ H+ + e-.
Assuming you do want to get the Gibbs energy of a reaction, you have two options. You can either balance the electrons in the half-reaction by completing the other half (we make it easy for you by providing a link for balancing with NAD+/NADH). Or, you could use the bottom panel in the results table to adjust the potential of the electrons in the reaction (i.e. change the value of e- potential [mV]). This is useful in cases where the half-reaction is coupled to an electron transfer for which the redox potential is known, but the compounds themselves are not, or they are too complex to use in eQuilibrator (e.g. protein complexes or membranal electron transport chains).
The reduction or oxidation of the pair Fe(III)/Fe(II) is ubiquitous in biology, for example in the iron-sulfur clusters of ferredoxins. However, the chemical environment of the iron atom can have a large effect on the reduction potential of the Fe(III)/Fe(II) pair with the measured reduction potentials of natural ferredoxins varying by more than 350 mV . That is, variation in the measured reduction potential of ferredoxins equals to reduction potential of NAD/NADH!
Similarly, in dissimilatory iron reduction the specific chemical form of Fe(III) can drastically affect the reduction potential. For example, a half reaction with a well-characterized crystalline form Goethite has a redox potential of about -300 mV while y-FeOOH, (Lepidocrocite), which can be treated as having the same empirical formula, has a redox potential of about -100 mV at pH 7 . As a result we strongly suggest that you enter the second half-reaction that interests you (e.g. reduction of pyruvate to acetyl-CoA) and use the bottom panel to adjust the potential of the electrons in the reaction to match the iron donor-acceptor pair that interests you.
The ΔG'° of the ATP hydrolysis reaction is affected by many factors, notably also by the concentration of free Mg2+ ions. The value cited in  and used in the original version of eQuilibrator (-36.4 kJ/mol) assumes no magnesium ([Mg2+] = 0). In the current version of eQuilibrator2.0 we use the component contribution method  that uses measurements collected in the NIST thermodynamic database for enzyme-catalyzed reactions  that were performed under varying levels of Mg2+. This is also the more relevant situation in vivo. As noted in many studies, when taking into account [Mg2+], the value changes and is observed to be in the range -26 to -32 kJ/mol depending on the reference. A clear discussion can be found at .
We define the driving force of a reaction or pathway as -ΔG' - i.e. a favorable reaction has a positive driving force. The total driving force for a pathway is the driving force associated with the pathway net reaction given particular metabolite and cofactor concentrations. The minimum total driving force, which we calculate in analyzing a pathway, is the smallest driving force associated with that pathway given the limits assumed for metabolite an cofactor concentrations. Similarly, the maximum total driving force is the largest driving force associated with pathway given those same limits.
The MDF of a pathway is a metric of how thermodynamically favorable a pathway can be in physiological conditions. The value of the MDF is smallest -ΔG' obtained by any pathway reaction when metabolite concentrations are chosen to make all pathway reactions as favorable as possible (-ΔG' as positive as possible).
You can read more about the MDF in PLoS Computational Biology.
If you have questions about eQuilibrator, please consult the Google Group to see if your question has been answered before. Please also post your questions to the group so that all eQuilibrator users can benefit from your experience. If you have questions about the data and methods on which eQuilibrator is based, please consult these references. If you'd like to contact us directly, reach out to the Milo Lab, which maintains eQuilibrator.