Predicting NMR relaxation of proteins from molecular dynamics simulations with accurate methyl rotation barriers

2020 ◽  
Vol 152 (8) ◽  
pp. 084102 ◽  
Author(s):  
Falk Hoffmann ◽  
Frans A. A. Mulder ◽  
Lars V. Schäfer
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Nmr Relaxation ◽  
Dynamics Simulations ◽  
Methyl Rotation

Structure and Dynamics of Spherical and Rodlike Alkyl Ethoxylate Surfactant Micelles Investigated Using NMR Relaxation and Atomistic Molecular Dynamics Simulations

2019 ◽  
Author(s):  
Allison Edwards ◽  
Abdolreza Javidialesaadi ◽  
Katie Weigandt ◽  
George Stan ◽  
Charles Eads
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Nmr Relaxation ◽  
Cross Relaxation ◽  
Relaxation Rates ◽  
Surfactant Micelles ◽  
Cylindrical Micelles ◽  
Spherical Micelles ◽  
Dynamics Simulations ◽  
Relative Motions

We study molecular arrangements and dynamics in alkyl ethoxylate nonionic surfactant micelles by combining high field (600 and 700 MHz) NMR relaxation measurements with large-scale atomistic molecular dynamics simulations. For spherical micelles, but not for cylindrical micelles, cross relaxation rates are positive only for surfactant alkyl tail atoms connected to the hydrophilic head group. All cross relaxation rates are negative for cylindrical micelles. This effect is reproducible either by changing composition (ratios of the nonionic surfactants) or changing temperature of a single surfactant in order to change the micelle shape. We validate the micelle shape by SANS and use the results as a guide for our simulations. We calculate parameters that determine relaxation rates directly from simulated trajectories, without introducing specific functional forms. Results indicate that relative motions of nearby atoms are liquid-like, in agreement with 13C T1 measurements, though constrained by micelle morphology. Relative motions of distant atoms have slower components because the relative changes in distances and angles are smaller when the moving atoms are further apart. The slow, long-range motions appear to be responsible for the predominantly negative cross relaxation rates observed in NOESY spectra. The densities of atoms from positions 1 and 2 in the boundary region are lower in spherical micelles compared to cylindrical micelles. Correspondingly, motions in this region are less constrained by micelle morphology in the spherical compared to the cylindrical cases. The two effects of morphology lead to the unusual occurrence of positive cross relaxation involving positions 1 and 2 for spheres.

Molecular dynamics simulations of NMR relaxation and diffusion of hydrocarbons

2018 ◽  
Author(s):  
Dilip Asthagiri ◽  
Philip M. Singer ◽  
Arjun Valiya Parambathu ◽  
Zeliang Chen ◽  
George J. Hirasaki ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Nmr Relaxation ◽  
Dynamics Simulations ◽  
And Diffusion

Structure and Dynamics of Spherical and Rodlike Alkyl Ethoxylate Surfactant Micelles Investigated Using NMR Relaxation and Atomistic Molecular Dynamics Simulations

2019 ◽  
Author(s):  
Allison Edwards ◽  
Abdolreza Javidialesaadi ◽  
Katie Weigandt ◽  
George Stan ◽  
Charles Eads
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Nmr Relaxation ◽  
Cross Relaxation ◽  
Relaxation Rates ◽  
Surfactant Micelles ◽  
Cylindrical Micelles ◽  
Spherical Micelles ◽  
Dynamics Simulations ◽  
Relative Motions

We study molecular arrangements and dynamics in alkyl ethoxylate nonionic surfactant micelles by combining high field (600 and 700 MHz) NMR relaxation measurements with large-scale atomistic molecular dynamics simulations. For spherical micelles, but not for cylindrical micelles, cross relaxation rates are positive only for surfactant alkyl tail atoms connected to the hydrophilic head group. All cross relaxation rates are negative for cylindrical micelles. This effect is reproducible either by changing composition (ratios of the nonionic surfactants) or changing temperature of a single surfactant in order to change the micelle shape. We validate the micelle shape by SANS and use the results as a guide for our simulations. We calculate parameters that determine relaxation rates directly from simulated trajectories, without introducing specific functional forms. Results indicate that relative motions of nearby atoms are liquid-like, in agreement with 13C T1 measurements, though constrained by micelle morphology. Relative motions of distant atoms have slower components because the relative changes in distances and angles are smaller when the moving atoms are further apart. The slow, long-range motions appear to be responsible for the predominantly negative cross relaxation rates observed in NOESY spectra. The densities of atoms from positions 1 and 2 in the boundary region are lower in spherical micelles compared to cylindrical micelles. Correspondingly, motions in this region are less constrained by micelle morphology in the spherical compared to the cylindrical cases. The two effects of morphology lead to the unusual occurrence of positive cross relaxation involving positions 1 and 2 for spheres.

Helix motion in protein C12A-p8MTCP1: Comparison of molecular dynamics simulations and multifield NMR relaxation data

2002 ◽  
Vol 23 (16) ◽  
pp. 1577-1586 ◽  
Author(s):  
Philippe Barthe ◽  
Christian Roumestand ◽  
Hélène Déméné ◽  
Laurent Chiche
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Nmr Relaxation ◽  
Relaxation Data ◽  
Dynamics Simulations

scholarly journals Predicting NMR Relaxation of Proteins from Molecular Dynamics Simulations with Accurate Methyl Rotation Barriers

2020 ◽  
Author(s):  
Falk Hoffmann ◽  
Frans Mulder ◽  
Lars V. Schäfer
Keyword(s):  
Molecular Dynamics ◽  
Force Field ◽  
Force Fields ◽  
Nmr Relaxation ◽  
Md Simulations ◽  
Quantitative Agreement ◽  
Relaxation Rates ◽  
Protein Force Fields ◽  
Dynamics Simulations ◽  
Methyl Rotation

The internal dynamics of proteins occurring on time scales from picoseconds to nanoseconds can be sensitively probed by nuclear magnetic resonance (NMR) spin relaxation experiments, as well as by molecular dynamics (MD) simulations. This complementarity offers unique opportunities, provided that the two methods are compared at a suitable level. Recently, several groups have used MD simulations to compute the spectral density of backbone and side-chain molecular motions, and to predict NMR relaxation rates from these. Unfortunately, in the case of methyl groups in protein side-chains, inaccurate energy barriers to methyl rotation were responsible for a systematic discrepancy in the computed relaxation rates, as demonstrated for the AMBER ff99SB*-ILDN force field (and related parameter sets), impairing quantitative agreement between simulations and experiments. However, correspondence could be regained by emending the MD force field with accurate coupled cluster quantum chemical calculations. Spurred by this positive result, we tested whether this approach could be generally applicable, in spite of the fact that different MD force fields employ different water models. Improved methyl group rotation barriers for the CHARMM36 and AMBER ff15ipq protein force fields were derived, such that the NMR relaxation data obtained from the MD simulations now also display very good agreement with experiment. Results herein showcase the performance of present-day MD force fields, and manifest their refined ability to accurately describe internal protein dynamics.

scholarly journals Predicting NMR Relaxation of Proteins from Molecular Dynamics Simulations with Accurate Methyl Rotation Barriers

2019 ◽  
Author(s):  
Falk Hoffmann ◽  
Frans Mulder ◽  
Lars V. Schäfer
Keyword(s):  
Molecular Dynamics ◽  
Force Field ◽  
Force Fields ◽  
Nmr Relaxation ◽  
Md Simulations ◽  
Quantitative Agreement ◽  
Relaxation Rates ◽  
Protein Force Fields ◽  
Dynamics Simulations ◽  
Methyl Rotation

The internal dynamics of proteins occurring on time scales from picoseconds to nanoseconds can be sensitively probed by nuclear magnetic resonance (NMR) spin relaxation experiments, as well as by molecular dynamics (MD) simulations. This complementarity offers unique opportunities, provided that the two methods are compared at a suitable level. Recently, several groups have used MD simulations to compute the spectral density of backbone and side-chain molecular motions, and to predict NMR relaxation rates from these. Unfortunately, in the case of methyl groups in protein side-chains, inaccurate energy barriers to methyl rotation were responsible for a systematic discrepancy in the computed relaxation rates, as demonstrated for the AMBER ff99SB*-ILDN force field (and related parameter sets), impairing quantitative agreement between simulations and experiments. However, correspondence could be regained by emending the MD force field with accurate coupled cluster quantum chemical calculations. Spurred by this positive result, we tested whether this approach could be generally applicable, in spite of the fact that different MD force fields employ different water models. Improved methyl group rotation barriers for the CHARMM36 and AMBER ff15ipq protein force fields were derived, such that the NMR relaxation data obtained from the MD simulations now also display very good agreement with experiment. Results herein showcase the performance of present-day MD force fields, and manifest their refined ability to accurately describe internal protein dynamics.

Sign in / Sign up

Export Citation Format

Share Document