The nature of phonons and solitary waves in alpha-helical proteins

A parametric study of the Davydov model of energy transduction in alpha-helical proteins is described. Previous investigations have shown that the Davydov model predicts that nonlinear interactions between phonons and amide-I excitations can stabilize the latter and produce a long-lived combined excitation (the so-called Davydov soliton), which propagates along the helix. The dynamics of this solitary wave are approximately those of solitons described using the nonlinear Schrödinger equation. The present study extends these previous investigations by analyzing the effect of helix length and nonlinear coupling efficiency on the phonon spectrum in short and medium length alpha-helical segments. The phonon energy accompanying amide-I excitation shows periodic variation in time with fluctuations that follow three different time scales. The phonon spectrum is highly dependent upon chain length but a majority of the energy remains localized in normal mode vibrations even in the long chain alpha-helices. Variation of the phonon-exciton coupling coefficient changes the amplitudes but not the frequencies of the phonon spectrum. The computed spectra contain frequencies ranging from 200 GHz to 6 THz, and as the chain length is increased, the long period oscillations increase in amplitude. The most important prediction of this study, however, is that the dynamics predicted by the numerical calculations have more in common with dynamics described by using the Frohlich polaron model than by using the Davydov soliton. Accordingly, the relevance of the Davydov soliton model was applied to energy transduction in alpha-helical proteins is questionable. We conclude that the Raman lines that have been assigned to solitons in E. coli are either associated with low frequency normal modes or are instrumental- or fluorescence-induced artifacts.