Analysis of the loop-helix interaction in bundle motif protein structures.

Molecular dynamics simulations and energy analysis have been carried out to study the structural mobility and stability of the four alpha-helix bundle motifs. The simulation results as well as the X-ray data show that the atomic RMS fluctuation is larger at the loop region for four representative proteins investigated: methemerythrin, cytochrome b-562, cytochrome c', and bovine somatotropin. The loop-loop, helix-helix, and loop-helix interactions are computed for the unfolded and folded proteins. In the folded and solvated protein structures the loop-helix interaction is stronger than the helix-helix interaction, especially in the electrostatic component. But the stabilization energies of both the loop-helix and the helix-helix interactions relative to those of an unfolded structure are of the same order of magnitude. The stabilization due to protein-solvent interaction is greater in the helix region than in the loop region. The percentage of hydrophilic solvent accessible area for the four proteins studied was calculated with the method of Eisenberg and McLachlan. The percentage of the hydrophilic area is greater in the loops than in the helices. A Poisson-Boltzmann calculation shows that the potential from the loops acting on a helix is generally more negative than that from other helices.