Aarhus University Seal

Accurate Methyl Group Dynamics in Protein Simulations with AMBER Force Fields

Research output: Contribution to journal/Conference contribution in journal/Contribution to newspaperJournal articleResearchpeer-review

Standard

Accurate Methyl Group Dynamics in Protein Simulations with AMBER Force Fields. / Hoffmann, Falk; Mulder, Frans A. A.; Schaefer, Lars V.

In: Journal of Physical Chemistry B, Vol. 122, No. 19, 17.05.2018, p. 5038-5048.

Research output: Contribution to journal/Conference contribution in journal/Contribution to newspaperJournal articleResearchpeer-review

Harvard

Hoffmann, F, Mulder, FAA & Schaefer, LV 2018, 'Accurate Methyl Group Dynamics in Protein Simulations with AMBER Force Fields', Journal of Physical Chemistry B, vol. 122, no. 19, pp. 5038-5048. https://doi.org/10.1021/acs.jpcb.8b02769

APA

Hoffmann, F., Mulder, F. A. A., & Schaefer, L. V. (2018). Accurate Methyl Group Dynamics in Protein Simulations with AMBER Force Fields. Journal of Physical Chemistry B, 122(19), 5038-5048. https://doi.org/10.1021/acs.jpcb.8b02769

CBE

MLA

Hoffmann, Falk, Frans A. A. Mulder and Lars V. Schaefer. "Accurate Methyl Group Dynamics in Protein Simulations with AMBER Force Fields". Journal of Physical Chemistry B. 2018, 122(19). 5038-5048. https://doi.org/10.1021/acs.jpcb.8b02769

Vancouver

Hoffmann F, Mulder FAA, Schaefer LV. Accurate Methyl Group Dynamics in Protein Simulations with AMBER Force Fields. Journal of Physical Chemistry B. 2018 May 17;122(19):5038-5048. doi: 10.1021/acs.jpcb.8b02769

Author

Hoffmann, Falk ; Mulder, Frans A. A. ; Schaefer, Lars V. / Accurate Methyl Group Dynamics in Protein Simulations with AMBER Force Fields. In: Journal of Physical Chemistry B. 2018 ; Vol. 122, No. 19. pp. 5038-5048.

Bibtex

@article{7c2b7227cd5342b9855bfa09c5f2b6b4,
title = "Accurate Methyl Group Dynamics in Protein Simulations with AMBER Force Fields",
abstract = "An approach is presented to directly simulate the dynamics of methyl groups in protein side-chains, as accessible via NMR spin relaxation measurements, by all-atom MD simulations. The method, which does not rely on NMR information or any system-specific adjustable parameters, is based on calculating the time-correlation functions (TCFs) of the C-H bonds in methyl groups and explicitly takes the truncation of the TCFs due to overall tumbling of the molecule into account. Using ubiquitin as a model protein, we show (i) that an accurate description of the methyl dynamics requires reparametrization of the potential energy barriers of methyl group rotation in the AMBER ff99SB*-ILDN force field (and related parameter sets), which was done with CCSD(T) coupled cluster calculations of frequency isolated dipeptides as reference, and (ii) that the TIP4P/2005 solvation model yields overall tumbling correlation times that are in close agreement with experimental data. The methyl axis squared order parameters S-axis(2) and associated correlation times tau(f), obtained within the Lipari-Szabo formalism, are in good agreement with the values derived from NMR deuterium relaxation experiments. Importantly, the relaxation rates and spectral densities derived from MD and NMR agree as well, enabling a direct comparison without assumptions inherent to simplified motional models.",
keywords = "SIDE-CHAIN DYNAMICS, NMR RELAXATION DATA, MAGNETIC-RESONANCE RELAXATION, ORDER-PARAMETER ANALYSIS, DEUTERIUM SPIN PROBES, MODEL-FREE APPROACH, MOLECULAR-DYNAMICS, ROTATIONAL DIFFUSION, CONFORMATIONAL ENTROPY, COMPUTER-SIMULATION",
author = "Falk Hoffmann and Mulder, {Frans A. A.} and Schaefer, {Lars V.}",
year = "2018",
month = may,
day = "17",
doi = "10.1021/acs.jpcb.8b02769",
language = "English",
volume = "122",
pages = "5038--5048",
journal = "Journal of Physical Chemistry Part B: Condensed Matter, Materials, Surfaces, Interfaces & Biophysical",
issn = "1520-6106",
publisher = "American Chemical Society",
number = "19",

}

RIS

TY - JOUR

T1 - Accurate Methyl Group Dynamics in Protein Simulations with AMBER Force Fields

AU - Hoffmann, Falk

AU - Mulder, Frans A. A.

AU - Schaefer, Lars V.

PY - 2018/5/17

Y1 - 2018/5/17

N2 - An approach is presented to directly simulate the dynamics of methyl groups in protein side-chains, as accessible via NMR spin relaxation measurements, by all-atom MD simulations. The method, which does not rely on NMR information or any system-specific adjustable parameters, is based on calculating the time-correlation functions (TCFs) of the C-H bonds in methyl groups and explicitly takes the truncation of the TCFs due to overall tumbling of the molecule into account. Using ubiquitin as a model protein, we show (i) that an accurate description of the methyl dynamics requires reparametrization of the potential energy barriers of methyl group rotation in the AMBER ff99SB*-ILDN force field (and related parameter sets), which was done with CCSD(T) coupled cluster calculations of frequency isolated dipeptides as reference, and (ii) that the TIP4P/2005 solvation model yields overall tumbling correlation times that are in close agreement with experimental data. The methyl axis squared order parameters S-axis(2) and associated correlation times tau(f), obtained within the Lipari-Szabo formalism, are in good agreement with the values derived from NMR deuterium relaxation experiments. Importantly, the relaxation rates and spectral densities derived from MD and NMR agree as well, enabling a direct comparison without assumptions inherent to simplified motional models.

AB - An approach is presented to directly simulate the dynamics of methyl groups in protein side-chains, as accessible via NMR spin relaxation measurements, by all-atom MD simulations. The method, which does not rely on NMR information or any system-specific adjustable parameters, is based on calculating the time-correlation functions (TCFs) of the C-H bonds in methyl groups and explicitly takes the truncation of the TCFs due to overall tumbling of the molecule into account. Using ubiquitin as a model protein, we show (i) that an accurate description of the methyl dynamics requires reparametrization of the potential energy barriers of methyl group rotation in the AMBER ff99SB*-ILDN force field (and related parameter sets), which was done with CCSD(T) coupled cluster calculations of frequency isolated dipeptides as reference, and (ii) that the TIP4P/2005 solvation model yields overall tumbling correlation times that are in close agreement with experimental data. The methyl axis squared order parameters S-axis(2) and associated correlation times tau(f), obtained within the Lipari-Szabo formalism, are in good agreement with the values derived from NMR deuterium relaxation experiments. Importantly, the relaxation rates and spectral densities derived from MD and NMR agree as well, enabling a direct comparison without assumptions inherent to simplified motional models.

KW - SIDE-CHAIN DYNAMICS

KW - NMR RELAXATION DATA

KW - MAGNETIC-RESONANCE RELAXATION

KW - ORDER-PARAMETER ANALYSIS

KW - DEUTERIUM SPIN PROBES

KW - MODEL-FREE APPROACH

KW - MOLECULAR-DYNAMICS

KW - ROTATIONAL DIFFUSION

KW - CONFORMATIONAL ENTROPY

KW - COMPUTER-SIMULATION

U2 - 10.1021/acs.jpcb.8b02769

DO - 10.1021/acs.jpcb.8b02769

M3 - Journal article

C2 - 29695158

VL - 122

SP - 5038

EP - 5048

JO - Journal of Physical Chemistry Part B: Condensed Matter, Materials, Surfaces, Interfaces & Biophysical

JF - Journal of Physical Chemistry Part B: Condensed Matter, Materials, Surfaces, Interfaces & Biophysical

SN - 1520-6106

IS - 19

ER -