Dielectric constant of aqueous solutions of proteins and organic polymers from molecular dynamics simulations

Abstract

The dielectric constant of water/oligomer mixtures, spanning the range from pure water to pure oligomeric melts, is investigated using molecular dynamics (MD) simulations. As prototypical water-soluble organic substances, we consider neutral poly-glycine, poly-ethylene glycol, and charged monomeric propionic acid. As the water content is reduced, the dielectric constant decreases but does not follow an ideal mixing behavior. The deviations from ideal mixing originate primarily in the non-linear relation between the oligomer mass fraction and collective polarization effects. We find that the dielectric constant is dominated by water polarization, even if the oligomer mass fraction exceeds 50%. By a double extrapolation of the MD simulation results to the limit of vanishing water fraction and to the limit of infinite oligomeric chain length, we estimate the orientational contribution to the dielectric constant of the pure polymeric melts. By this procedure, we obtain ϵ = 17 ± 2 for polyglycine and ϵ = 1 ± 0.3 for polyethylene glycol. The large difference is rationalized by polarization correlations of glycine units. Interestingly, we find constant temperature simulations to outperform replica exchange simulations in terms of equilibration speed.