Langevin dynamics of peptides: The frictional dependence of isomerization rates of N‐acetylalanyl‐N′‐methylamide

Abstract

The rate constant for the transition between the equatorial and axial conformations of N‐acetylalanyl‐N′‐methylamide has been determined from Langevin dynamics (LD) simulations with no explicit solvent. The isomerization rate is maximum at collision frequency γ = 2 ps⁻¹, shows diffusive character for γ ≥ 10 ps⁻¹, but does not approach zero even at γ = 0.01 ps⁻¹. This behavior differs from that found for a one‐dimensional bistable potential and indicates that both collisional energy transfer with solvent and vibrational energy transfer between internal modes are important in the dynamics of barrier crossing for this system. It is suggested that conformational searches of peptides be carried out using LD with a collision frequency that maximizes the isomerization rate (i.e., γ ≈ 2 ps⁻¹). This method is expected to be more efficient than either molecular dynamics in vacuo (which corresponds to LD with γ = 0) or molecular dynamics in solvent (where dynamics is largely diffusive).