Molecular dynamics and quantum chemical approaches in the study of the hydration of protonated cyclohexyldiamines

Explicit hydration of the neutral and charged cyclohexylamine and of the cyclohexyldiamine isomers in their mono- or diprotonated forms is investigated through classical molecular dynamics (MD) simulations in aqueous solutions combined with DFT calculations in amine–water complexes. The MD studies performed in the monoamines reveal that the structure of the hydration shell around the neutral amino group (NH₂) is quite distinct from that around the charged one (NH₃⁺). On average, the number of water molecules surrounding the two groups is calculated to be ～2 and 3–4, respectively. The variation of the hydration structure prompted by the groups’ proximity is discussed based on the data found for the mono- and diprotonated diamines. To have a more detailed picture of the water molecules’ arrangement around the amino groups and of the amine–water hydrogen bonds, geometry optimisations in hydrates with up to six water molecules are carried out at the B3LYP/aug-cc-pVDZ level. Complexation energies are also computed. The main findings emerging from these calculations are found to be very helpful to rationalise the mutual influence of the amino groups and therefore to better elucidate the MD findings. The complementary nature of the two research methods is emphasised as an excellent tool in order to closely examine the hydration of polyamines, as exemplified for the cyclohexyldiamines.