Orientation restraints in molecular dynamics simulations using time and ensemble averaging

Area of Science:

• Biophysics
• Structural Biology
• Computational Chemistry

Background:

• Nuclear Magnetic Resonance (NMR) spectroscopy provides valuable data on molecular structure and dynamics.
• Partially ordered molecules present challenges for traditional molecular dynamics (MD) simulations.
• Residual dipolar couplings (RDCs) and chemical-shift anisotropies (CSAs) are key NMR-derived restraints.

Purpose of the Study:

• To develop and validate a methodology for simulating protein dynamics with NMR-derived restraints.
• To investigate the impact of averaging restraints on molecular fluctuations and order parameters.
• To provide a robust computational tool for studying partially ordered biological molecules.

Main Methods:

• Definition of a restraint potential compatible with molecular dynamics and energy minimization.
• Implementation of time and ensemble averaging for restraint application.
• Extensive MD simulations of histidine containing phosphocarrier protein using RDC-derived restraints.

Main Results:

• Simulations demonstrate that time- or ensemble-averaged restraints preserve the natural fluctuations of restrained vectors.
• Unaveraged restraints significantly reduce the fluctuations of backbone N-H vectors.
• Averaging is essential for accurate determination of the molecular order-parameter tensor.

Conclusions:

• The presented methodology enables accurate simulation of protein dynamics with NMR restraints.
• Averaging restraints is critical for maintaining realistic molecular motion and order parameters.
• This approach enhances the utility of NMR data in computational structural biology.