Thermodynamic Basis of Molecular Diffusion through Cyanobacterial Septal Junctions.

LETTER

A recent paper published by Nieves-Morión and colleagues proposes that the intercellular transfer of small molecules in heterocyst-forming cyanobacteria (calcein, 5-carboxyfluorescein [5-CF], and esculin) is mediated by simple diffusion based on the Q10 values of the kinetics (1). In this letter, we show our analysis of their results to provide a thermodynamic basis of the kinetics.

We obtained rate constant values of the intercellular transfer for those molecules from Fig. 1 of the original paper (1) using the ImageJ program (https://imagej.nih.gov/ij/). We calculated molar free energy of activation (ΔG°) for each molecule at each temperature (T) using the Eyring equation (equation 1) from the corresponding rate constant (k): where R, h, and k are the ideal gas constant (8.3145 J/K mol), Planck’s constant (6.6261 × 10−34 J s), and the Boltzmann constant (1.3807 × 10−23 J/K), respectively (2). ΔG° for each molecule shows a strong linear relationship with temperature (Fig. 1A). From this relationship, we obtained molar enthalpy and entropy of activation (ΔH° and ΔS°) for the process using the following equation:

FIG 1 

(A) Temperature-dependent standard molar free energy of activation for the intercellular transfer of calcein, 5-carboxyfluorescein (5-CF), and esculin. R2 of the fitting is shown in the figure (three values in parentheses). (B) Relationship between enthalpy and entropy of activation (R2 = 0.9907). SigmaPlot (version 11; Systat Software, Inc., San Jose, CA) was used for graph preparation and statistical analysis.

Figure 1B shows there is a correlation between ΔH° and ΔS° for these three molecules. This correlation is called enthalpy-entropy compensation, suggesting that a common mechanism is responsible for the intercellular transfer of these three molecules (3). Furthermore, ΔS° was found to be the major component in ΔG° for all three molecules, as it explains at least 95% of the value of ΔG° at 37°C. This indicates that the intercellular transfer process is highly unfavorable in terms of entropy. This feature correlates with the suggestion made by the original paper in that the simple diffusion of molecules through a restricted space reduces the entropy of the molecules because entropy is proportional to the volume available (2).

Membrane transport kinetics is typically analyzed in terms of activation energy (E). Here we obtained the E values for all three molecules from ΔH° using the following equation:

The E values at 37°C are 6.1, 3.3, and 1.4 kJ/mol for calcein, 5-CF, and esculin, respectively. These values are much smaller than that of typical facilitated diffusion. For example, facilitated diffusion for glucose uptake has an E of 64 kJ/mol (4). The values of E are comparable to the hydrogen bonding energy (4 to 40 kJ/mol) (2), suggesting that dissociation of water molecules from the molecules is involved in the diffusion. Our thermodynamic analysis confirms that the intercellular transfer of those molecules is mediated by simple diffusion.

After data analyses were completed, it was found that there was a relationship between both molecular weight and charge of a molecule with the molecule’s entropy and enthalpy values (Table 1). We were unable to conclude which of the structural parameters is responsible for the correlation of the thermodynamic parameters because of the inherent correlation between those two structural parameters themselves. Calcein is a relatively large molecule possessing a charge of −5, while esculin is a small molecule with no charge and 5-CF is of moderate size and possess a charge of −1. Therefore, we suggest a future experiment in which intercellular transfer is measured with small molecules possessing a high negative charge and/or large molecules with a low negative charge to identify which factor is responsible for the correlation to the thermodynamic parameters. This will greatly enhance understanding of the process.

TABLE 1 

Pearson correlation coefficient values between thermodynamic and structural parameters

Parameter	Pearson correlation coefficient value between parameter:
	ΔHa°	ΔSa°	Mol wt	Negative charge
ΔHa°	1	0.995	0.957	0.976
ΔSa°		1	0.981	0.992
Mol wt			1	0.997
Negative charge				1