Improved Modeling of Peptidic Foldamers Using a Quantum Chemical Parametrization Based on Torsional Minimum Energy Path Matching

Abstract

The increasing interest in novel foldamer constructs demands an accurate computational treatment on an extensive timescale. However, it is still a challenge to derive a force field (FF) that can reproduce the experimentally known fold while also allowing the spontaneous exploration of other structures. Here, aiming at a realistic reproduction of backbone torsional barriers, the relevant proper dihedrals of acyclic β2-, β3- and β2,3-amino acids were added to the CHARMM FF and optimized using a novel, self-consistent iterative procedure based on quantum chemical relaxed scans. The new FF was validated by molecular dynamics simulations on three acyclic peptides. While they resided most of the time in their preferred fold (>80 % in helices and >50 % in hairpin), they also visited other conformations. Owing to the CHARMM36m-consistent parametrization, the proposed extension is suitable for exploring new foldamer structures and assemblies, and their interactions with diverse biomolecules.