A refined hydrogen bond potential for flexible protein models

Abstract

One of the major disadvantages of coarse-grained hydrogen bond potentials, for their use in protein folding simulations, is the appearance of abnormal structures when these potentials are used in flexible chain models, and no other geometrical restrictions or energetic contributions are defined into the system.We have efficiently overcome this problem, for chains of adequate size in a relevant temperature range, with a refined coarse-grained hydrogen bond potential. With it, we have been able to obtain nativelike alpha-helices and beta-sheets in peptidic systems, and successfully reproduced the competition between the populations of these secondary structure elements by the effect of temperature and concentration changes. In this manuscript we detail the design of the interaction potential and thoroughly examine its applicability in energetic and structural terms, considering factors such as chain length, concentration, and temperature.