Hydrophobicity of an isobutane dimer in water, methanol and acetonitrile as solvents — A classical molecular dynamics study

Abstract

The potential of mean forces (PMFs) has been determined for an isobutane dimer in various solvent environments such as water, methanol and acetonitrile at a temperature of 298 K and pressure of 1 bar using GROMACS software. All the molecular dynamics (MD) simulations are carried out using a TIP3P water model under a CHARMM36 forcefield. Following Umbrella Sampling technique, PMFs are calculated and analyzed using Weighted Histogram Analysis Method (WHAM) and coordination number of first solvation shell is extracted for all solvents using radial distribution function. The shape of PMFs contains contact minima, solvent-separated minima and desolvation maxima. The values of contact minima are not affected much by solvent environment and found to be at 0.5377, 0.5480 and 0.5495 nm for water, methanol and acetonitrile respectively. The corresponding energy depths are found −0.9134, −0.7080 and −0.5295 kcalmol[Formula: see text]. The variation observed at solvent-separated minima is noticeable and found at 0.9012, 0.9721 and 0.9151 nm for water, methanol and acetonitrile, respectively. The coordination number of the first solvation shell by taking an isobutane molecule as a reference from their center of mass is found to be 28.1, 16.9 and 14.8 for water, methanol and acetonitrile, respectively. There is a soft hydrophobic interaction between isobutane dimer and solvents like methanol and acetonitrile relative to water, might be due to the presence of competitive methyl group of methanol and acetonitrile in the solvent medium.