Development of CHARMM additive potential energy parameters for -methyl amino acids

Potential energy parameters for α-methyl amino acids were generated with ab initio calculations on α-methyl-N-acetylalanyl-N′-methylamide (the α-methyl “alanine dipeptide”) which served as an input to a grid-based correction to the backbone torsional potential (known as CMAP) consistent with the CHARMM36m additive protein force field. The new parameters were validated by comparison with experimentally determined helicities of the 22 residue C-terminal peptide (H10) from apolipoprotein A1 and five α-methylated variants in water and 0.3:0.7 trifluoroethanol (TFE)/water. Conventional molecular dynamics simulation totaling 30 μs for each peptide is in overall good agreement with the experiment, including the increased helicity in 30% TFE. An additional 500 ns of simulation using two-dimensional dihedral biasing (bpCMAP) replica exchange reduced left-handed conformations, increased right-handed helices, and thereby mostly decreased agreement with the experiment. Analysis of side chain–side chain salt bridges suggests that the overestimation of the helical content may be, in part, due to such interactions. The increased helicity of the peptides in 30% TFE arises from decreased hydrogen bonding of the backbone atoms to water and a concomitant increase in intramolecular backbone hydrogen bonds. The present inclusion of parameters for amino acids methylated at the α-carbon can be extended to all amino acids that contain two β-carbons such as isovaline and stapled peptides.