Simulation protocols for coherent femtosecond vibrational spectra of peptides.
Two algorithms for simulating the response of peptides to sequences of IR pulses are developed and applied to N-methyl acetamide (NMA) and a 17 residue alpha-helical peptide (YKKKH17) in D(2)O. A fluctuating vibrational-exciton Hamiltonian for the amide I mode is constructed from molecular dynamics trajectories. Coupling with the environment is described using a density functional theory electrostatic map. The cumulant expansion of Gaussian fluctuation incorporates motional narrowing due to fast frequency fluctuations and is adequate for NMA and for isotopically labeled bands in large peptides. Real-space truncation of the scattering matrix of the nonlinear exciton equations significantly reduces the computational cost, making it particularly attractive for slow fluctuations in large globular proteins.
Zhuang
Wei
W
Department of Chemistry, University of California at Irvine, California 92697-2025, USA.
Abramavicius
Darius
D
Hayashi
Tomoyuki
T
Mukamel
Shaul
S
D002627:Chemistry, Physical / Q000379:methods*; D003198:Computer Simulation; D017666:Deuterium Oxide / Q000737:chemistry; D008956:Models, Chemical; D008958:Models, Molecular; D015233:Models, Statistical; D008968:Molecular Conformation; D016011:Normal Distribution; D010455:Peptides / Q000737:chemistry*; D011487:Protein Conformation; D017510:Protein Folding; D017433:Protein Structure, Secondary
D010455:Peptides; D017666:Deuterium Oxide