Turn-Directed α-β Conformational Transition of α-syn12 Peptide at Different pH Revealed by Unbiased Molecular Dynamics Simulations

The transition from α-helical to β-hairpin conformations of α-syn12 peptide is characterized here using long timescale, unbiased molecular dynamics (MD) simulations in explicit solvent models at physiological and acidic pH values. Four independent normal MD trajectories, each 2500 ns, are performed at 300 K using the GROMOS 43A1 force field and SPC water model. The most clustered structures at both pH values are β-hairpin but with different turns and hydrogen bonds. Turn(9-6) and four hydrogen bonds (HB(9-6), HB(6-9), HB(11-4) and HB(4-11)) are formed at physiological pH; turn(8-5) and five hydrogen bonds (HB(8-5), HB(5-8), HB(10-3), HB(3-10) and HB(12-1)) are formed at acidic pH. A common folding mechanism is observed: the formation of the turn is always before the formation of the hydrogen bonds, which means the turn is always found to be the major determinant in initiating the transition process. Furthermore, two transition paths are observed at physiological pH. One of the transition paths tends to form the most-clustered turn and improper hydrogen bonds at the beginning, and then form the most-clustered hydrogen bonds. Another transition path tends to form the most-clustered turn, and turn(5-2) firstly, followed by the formation of part hydrogen bonds, then turn(5-2) is extended and more hydrogen bonds are formed. The transition path at acidic pH is as the same as the first path described at physiological pH.