Stabilization of intermediate density states in globular proteins by homogeneous intramolecular attractive interactions

On-lattice simulations of two-dimensional self-avoiding chains subject to homogeneous intramolecular attractive interactions were performed as a model for studying various density regimes in globular proteins. For short chains of less than 15 units, all conformations were generated and classified by density. The range of intramolecular interactions was found to increase uniformly with density, and the average number of topological contacts is directly proportional to density. The uniform interaction energy increases the probability of high density states but does not necessarily lead to dominance of the highest density state. Typically, several large peaks appear in the probability distribution of packing densities, their location and amplitude being determined by the balance between entropic effects enhancing more expanded conformations and attractive interactions favoring compact forms. Also, the homogeneous interaction energy affects the distribution of most probable interacting points in favor of the longer range interactions over the short range ones, but in addition it introduces some more detailed preferences even among short range interactions. There are some implications about the characteristics of the intermediate density states and also for the likelihood that the native state does not correspond completely to the lowest energy conformation.

Cooperative structural transitions induced by non-homogeneous intramolecular interactions in compact globular proteins

The role played by non-homogeneous interactions in stabilizing cooperative structural changes in proteins was investigated by exhaustive simulations of all compact conformations compatible with several well-defined globule-like shapes in three dimensions. Conformational free energies corresponding to the association of residues i and j were computed both for the unperturbed system, all subject to identical intramolecular interactions, and for the perturbed system in which a single pair of residues is probed by changing its interactions with an attractive or repulsive interaction. The high packing density leads to strong coupling between residues so that specific interactions between a given pair of residues are accompanied by considerable enthalpy changes. Relatively weak, about 1-2 kcal/mol, attractive interactions can exert a dramatic effect on the free energy distribution. Usually, central positions in the sequence most affect the conformational characteristics. Some of these interaction pairs appear to be capable of effecting major conformation transitions because of the high level of cooperativity in the dense state. Effects of repulsive interactions, however, do not depend so strongly on residue pair and cause more localized structural changes. This approach can suggest more, or less, sensitive loci for amino acid substitution.