A theoretical analysis of chlorophyll a (Chla) hydration processes in aqueous organic solvents has been carried out by means of quantum chemistry calculations. A detailed knowledge of the thermodynamics of these processes is fundamental in order to better understand the organization of chlorophyll molecules in vivo, specifically the structure of chlorophyll pairs in photosystems I and II. In the present work, we assumed a Chla model in which the phytyl chain is replaced by a methyl group. Calculations were performed at the B3LYP/6-31G(d) level corrected for basis set superposition errors and dispersion interaction energy. This computational scheme was previously shown to provide data close to MP2/6-311++(2d,2p) results. Solvents effects were taken into account using either continuum (for nonpolar solvents) or discrete-continuum (for polar coordinating solvents) methods. In the latter case, we first examined the structure of Chla in rigorously dry solutions. Two types of solvents were characterized according to Mg-atom coordination: In type I solvents (acetone, acetonitrile, DMSO), Mg exhibits five-coordination, whereas in type II solvents (THF, pyridine), Mg exhibits six-coordination. Hydration processes are quite dependent on solvent nature. In nonpolar or low-polarity solvents such as cyclohexane or chloroform, hydration is always exothermic and exergonic, despite a large entropy term that strongly opposes hydration. In polar solvents of type II, hydration is quite unfavorable, and essentially no hydrates are expected in these media, except perhaps at very large water concentrations (although, in such a case, the medium cannot be simply described as an organic solvent). In polar solvents of type I, the situation is intermediate, and dihydration is favorable in some cases (acetone, acetonitrile) and unfavorable in others (DMSO). It is interesting to note that first hydration processes in coordinating solvents (of either type I or type II), where a water molecule must displace a solvent molecule coordinated to Mg, exhibit values of DeltaH > 0 and DeltaS > 0, in sharp contrast to first hydration processes in nonpolar media. The present results represent the first theoretical attempt to rationalize the large amount of experimental data on hydration and aggregation of Chla in aqueous organic media that have been accumulated over the past four decades. The data stress, in particular, the key role of Chla dihydrates, a point that has been the object of intense debate in the literature. Clearly, dihydrates are found to be more stable than monohydrates owing to a particular structure in which cooperative interactions occur between the water molecules and Chla. The calculations also explain the irregular behavior observed for Chla in aqueous THF or pyridine: In these media, Chla remains basically unhydrated because the Chla-solvent adducts are stabilized by strong dispersion interactions.
