Dynamics and deformability of α-, 310- and πhelices

Abstract

Protein structures are often represented as seen in crystals as (i) rigid macromolecules (ii) with helices, sheets and coils. However, both definitions are partial because (i) proteins are highly dynamic macromolecules and (ii) the description of protein structures could be more precise. With regard to these two points, we analyzed and quantified the stability of helices by considering α-helices as well as 310- and πhelices. Molecular dynamic (MD) simulations were performed on a large set of 169 representative protein domains. The local protein conformations were followed during each simulation and analyzed. The classical flexibility index (B-factor) was confronted with the MD root mean square flexibility (RMSF) index. Helical regions were classified according to their level of helicity from high to none. For the first time, a precise quantification showed the percentage of rigid and flexible helices that underlie unexpected behaviors. Only 76.4% of the residues associated with α-helices retain the conformation, while this tendency drops to 40.5% for 310-helices and is never observed for πhelices. α-helix residues that do not remain as an α-helix have a higher tendency to assume β-turn conformations than 310- or πhelices. The 310-helices that switch to the α-helix conformation have a higher B-factor and RMSF values than the average 310-helix but are associated with a lower accessibility. Rare π-helices assume a β-turn, bend and coil conformations, but not α- or 310-helices. The view on πhelices drastically changes with the new DSSP (Dictionary of Secondary Structure of Proteins) assignment approach, leading to behavior similar to 310-helices, thus underlining the importance of secondary structure assignment methods.