Solvation Free Energies of Alanine Peptides: The Effect of Flexibility

The electrostatic (Delta G(el)), van der Waals cavity-formation (Delta G(vdw)), and total (Delta G) solvation free energies for 10 alanine peptides ranging in length (n) from 1 to 10 monomers were calculated. The free energies were computed both with fixed, extended conformations of the peptides and again for some of the peptides without constraints. The solvation free energies, Delta G(el), and components Delta G(vdw), and Delta G, were found to be linear in n, with the slopes of the best-fit lines being gamma(el), gamma(vdw), and gamma, respectively. Both gamma(el) and gamma were negative for fixed and flexible peptides, and gamma(vdw) was negative for fixed peptides. That gamma(vdw) was negative was surprising, as experimental data on alkanes, theoretical models, and MD computations on small molecules and model systems generally suggest that gamma(vdw) be positive. A negative gamma(vdw) seemingly contradicts the notion that Delta G(vdw) drives the initial collapse of the protein when it folds by favoring conformations with small surface areas. When we computed Delta G(vdw) for the flexible peptides, thereby allowing the peptides to assume natural ensembles of more compact conformations, gamma(vdw) was positive. Because most proteins do not assume extended conformations, a Delta G(vdw) that increases with increasing surface area may be typical for globular proteins. An alternative hypothesis is that the collapse is driven by intramolecular interactions. We find few intramolecular H-bonds but show that the intramolecular van der Waals interaction energy is more favorable for the flexible than for the extended peptides, seemingly favoring this hypothesis. The large fluctuations in the vdw energy may make attributing the collapse of the peptide to this intramolecular energy difficult.