Extension of the AMBER force field to cyclic alpha,alpha dialkylated peptides
Abstract
The popular biomolecular AMBER (ff99SB) force field (FF) has been extended with new parameters for the simulations of peptides containing alpha, alpha dialkylated residues with cyclic side chains. Together with the recent set of nitroxide parameters [E. Stendardo, A. Pedone, P. Cimino, M. C. Menziani, O. Crescenzi and V. Barone, Phys. Chem. Chem. Phys., 2010, 12, 11697] this extension allows treating the TOAC residue (TOAC, 2,2,6,6-tetramethylpiperidine-1-oxyl-4-amino-4-carboxylic acid) widely used as a spin label in protein studies. All the conformational minima of the Ac-Ac6C-NMe (Ac = acetyl, Ac6C = 1-aminocyclohexaneacetic acid, NMe = methylamino) and Ac-TOAC-NMe dipeptides have been examined in terms of geometry and relative energy stability by Quantum Mechanical (QM) computations employing an hybrid density functional (PBE0) for an extended training set of conformers with various folds. A very good agreement between QM and MM (molecular mechanics) data has been obtained in most o...f the investigated properties, including solvent effects. Finally, the new set of parameters has been validated by comparing the conformational and dynamical behavior of TOAC-labeled polypeptides investigated by means of classical molecular dynamics (MD) simulations with QM data and experimental evidence. The new FF accurately describes the tuning of conformational and dynamical behavior of the Ac-TOAC-NMe dipeptide and double spin-labeled heptapeptide Fmoc-(Aib-Aib-TOAC)(2)-Aib-OMe (Fmoc, fluorenyl-9-methoxycarbonyl; Aib, alpha-aminoisobutyric acid; OMe, methoxy) by solvents with different polarity. In particular, we found that the 3(10) helical structure of heptapeptide is the most stable one in vacuo, with a geometry very similar to the X-ray crystallographic structure, whereas a conformational equilibrium between the 3(10)- and alpha-helical structures is established in aqueous solution, in agreement with EPR data.
Source:
Physical Chemistry Chemical Physics, 2012, 14, 44, 15308-15320Publisher:
- Royal Soc Chemistry, Cambridge
Funding / projects:
- Italian Ministry for Education, University and Research (MIUR)
- INFN
- COST Action "COnvergent Distributed Environment for Computational Spectroscopy (CODECS)'' - CM1002
DOI: 10.1039/c2cp42713c
ISSN: 1463-9076
PubMed: 23051698
WoS: 000310153300008
Scopus: 2-s2.0-84867928399
Collections
Institution/Community
IHTMTY - JOUR AU - Grubišić, Sonja AU - Brancato, Giuseppe AU - Pedone, Alfonso AU - Barone, Vincenzo PY - 2012 UR - https://cer.ihtm.bg.ac.rs/handle/123456789/1127 AB - The popular biomolecular AMBER (ff99SB) force field (FF) has been extended with new parameters for the simulations of peptides containing alpha, alpha dialkylated residues with cyclic side chains. Together with the recent set of nitroxide parameters [E. Stendardo, A. Pedone, P. Cimino, M. C. Menziani, O. Crescenzi and V. Barone, Phys. Chem. Chem. Phys., 2010, 12, 11697] this extension allows treating the TOAC residue (TOAC, 2,2,6,6-tetramethylpiperidine-1-oxyl-4-amino-4-carboxylic acid) widely used as a spin label in protein studies. All the conformational minima of the Ac-Ac6C-NMe (Ac = acetyl, Ac6C = 1-aminocyclohexaneacetic acid, NMe = methylamino) and Ac-TOAC-NMe dipeptides have been examined in terms of geometry and relative energy stability by Quantum Mechanical (QM) computations employing an hybrid density functional (PBE0) for an extended training set of conformers with various folds. A very good agreement between QM and MM (molecular mechanics) data has been obtained in most of the investigated properties, including solvent effects. Finally, the new set of parameters has been validated by comparing the conformational and dynamical behavior of TOAC-labeled polypeptides investigated by means of classical molecular dynamics (MD) simulations with QM data and experimental evidence. The new FF accurately describes the tuning of conformational and dynamical behavior of the Ac-TOAC-NMe dipeptide and double spin-labeled heptapeptide Fmoc-(Aib-Aib-TOAC)(2)-Aib-OMe (Fmoc, fluorenyl-9-methoxycarbonyl; Aib, alpha-aminoisobutyric acid; OMe, methoxy) by solvents with different polarity. In particular, we found that the 3(10) helical structure of heptapeptide is the most stable one in vacuo, with a geometry very similar to the X-ray crystallographic structure, whereas a conformational equilibrium between the 3(10)- and alpha-helical structures is established in aqueous solution, in agreement with EPR data. PB - Royal Soc Chemistry, Cambridge T2 - Physical Chemistry Chemical Physics T1 - Extension of the AMBER force field to cyclic alpha,alpha dialkylated peptides VL - 14 IS - 44 SP - 15308 EP - 15320 DO - 10.1039/c2cp42713c ER -
@article{ author = "Grubišić, Sonja and Brancato, Giuseppe and Pedone, Alfonso and Barone, Vincenzo", year = "2012", abstract = "The popular biomolecular AMBER (ff99SB) force field (FF) has been extended with new parameters for the simulations of peptides containing alpha, alpha dialkylated residues with cyclic side chains. Together with the recent set of nitroxide parameters [E. Stendardo, A. Pedone, P. Cimino, M. C. Menziani, O. Crescenzi and V. Barone, Phys. Chem. Chem. Phys., 2010, 12, 11697] this extension allows treating the TOAC residue (TOAC, 2,2,6,6-tetramethylpiperidine-1-oxyl-4-amino-4-carboxylic acid) widely used as a spin label in protein studies. All the conformational minima of the Ac-Ac6C-NMe (Ac = acetyl, Ac6C = 1-aminocyclohexaneacetic acid, NMe = methylamino) and Ac-TOAC-NMe dipeptides have been examined in terms of geometry and relative energy stability by Quantum Mechanical (QM) computations employing an hybrid density functional (PBE0) for an extended training set of conformers with various folds. A very good agreement between QM and MM (molecular mechanics) data has been obtained in most of the investigated properties, including solvent effects. Finally, the new set of parameters has been validated by comparing the conformational and dynamical behavior of TOAC-labeled polypeptides investigated by means of classical molecular dynamics (MD) simulations with QM data and experimental evidence. The new FF accurately describes the tuning of conformational and dynamical behavior of the Ac-TOAC-NMe dipeptide and double spin-labeled heptapeptide Fmoc-(Aib-Aib-TOAC)(2)-Aib-OMe (Fmoc, fluorenyl-9-methoxycarbonyl; Aib, alpha-aminoisobutyric acid; OMe, methoxy) by solvents with different polarity. In particular, we found that the 3(10) helical structure of heptapeptide is the most stable one in vacuo, with a geometry very similar to the X-ray crystallographic structure, whereas a conformational equilibrium between the 3(10)- and alpha-helical structures is established in aqueous solution, in agreement with EPR data.", publisher = "Royal Soc Chemistry, Cambridge", journal = "Physical Chemistry Chemical Physics", title = "Extension of the AMBER force field to cyclic alpha,alpha dialkylated peptides", volume = "14", number = "44", pages = "15308-15320", doi = "10.1039/c2cp42713c" }
Grubišić, S., Brancato, G., Pedone, A.,& Barone, V.. (2012). Extension of the AMBER force field to cyclic alpha,alpha dialkylated peptides. in Physical Chemistry Chemical Physics Royal Soc Chemistry, Cambridge., 14(44), 15308-15320. https://doi.org/10.1039/c2cp42713c
Grubišić S, Brancato G, Pedone A, Barone V. Extension of the AMBER force field to cyclic alpha,alpha dialkylated peptides. in Physical Chemistry Chemical Physics. 2012;14(44):15308-15320. doi:10.1039/c2cp42713c .
Grubišić, Sonja, Brancato, Giuseppe, Pedone, Alfonso, Barone, Vincenzo, "Extension of the AMBER force field to cyclic alpha,alpha dialkylated peptides" in Physical Chemistry Chemical Physics, 14, no. 44 (2012):15308-15320, https://doi.org/10.1039/c2cp42713c . .