Full Text PDF:

PDF icon paoli09.pdf

SI PDF:

PDF icon paoli09_s.pdf

Authors:
B. Paoli; M. Seeber; E.H.G. Backus; J.A. Ihalainen; P. Hamm; A. Caflisch

Journal: J. Phys. Chem. B
Year: 2009
Volume: 113
Issue: 13
Pages: 4435-4442
DOI: 10.1021/jp810431s
Type of Publication: Journal Article

Keywords:
Computer Simulation; Kinetics; Models, Molecular; Mutation; Peptides; Protein Folding; Protein Structure, Secondary; Protein Structure, Tertiary; Salts; Spectrophotometry, Infrared; Time Factors

Abstract:

Multiple 4-µs molecular dynamics (MD) simulations are used to study the folding process of the cross-linked α-helical peptide Ac-EACAR5EAAAR10EAACR15Q-NH2 (EAAAR peptide). The folding kinetics are single exponential at 330 K, while they are complex at 281 K with a clear deviation from single-exponential behavior, in agreement with time-resolved infrared (IR) spectroscopy measurements. Network analysis of the conformation space sampled by the MD simulations reveals four main folding channels which start from conformations with partially formed helical structure and non-native salt-bridges in a kinetically partitioned unfolded state. The independent folding pathways explain the comparable quality of models based on stretched exponential and multiexponential fitting of the kinetic traces at low temperature. The rearrangement of bulky side chains, and in particular their reorientation with respect to the cross-linker, makes the EAAAR peptide a slower folder at 281 K than a similar peptide devoid of the three glutamate side chains. On the basis of this simulation result, extracted from a total MD sampling of 1.0 ms, a mutant with additional bulky side chains (three methionines replacing alanines at positions 2, 7, and 12) is suggested to fold slower than the EAAAR peptide. This prediction is confirmed by time-resolved IR spectroscopy.