Strong electrostatic loop-helix interactions in bundle motif protein structures

Based on CHARMM potential (Brooks et al., 1983) an energetic analysis has been carried out for four typical 4-alpha-helix bundle proteins, i.e., methemerythrin, cytochrome b-562, cytochrome c', and bovine somatotropin. The bovine somatotropin possesses long loops, but all the other three proteins have short loops. It was found that in all these four 4-alpha-helix bundle motif structures the interaction between loops and helices was much stronger than the interaction among the four helices themselves. Particularly for the electrostatic interaction energy, the loop-helix interaction is overwhelmingly stronger than the interhelix interaction although the latter involves the favorable helix dipole interaction due to the antiparallel arrangement of neighboring alpha-helices. The present study indicates that such a conclusion holds true regardless of what loops, long or short, are in the 4-alpha-helix bundle protein, and also regardless of which empirical potential, ECEPP or CHARMM, is used for calculations although in CHARMM the electrostatic energy is much more heavily emphasized than in ECEPP. Therefore, no appropriate conclusion can be drawn in arguing whether the dipole interaction among the four alpha-helices play a stabilizing role or destabilizing role for a 4-alpha-helix bundle protein without taking into consideration the effect of interaction between helices and loops. The calculated results reported here provide, from a different point of view, insights that might be useful for revealing the essence of the driving forces during the folding of proteins.