We report over 30 μs of unrestrained molecular dynamics simulations of six protein-RNA complexes in explicit solvent. We utilize the AMBER ff99bsc0χ(OL3) RNA force field combined with the ff99SB protein force field and its more recent ff12SB version with reparametrized side-chain dihedrals. The simulations show variable behavior, ranging from systems that are essentially stable to systems with progressive deviations from the experimental structure, which we could not stabilize anywhere close to the starting experimental structure. For some systems, microsecond-scale simulations are necessary to achieve stabilization after initial sizable structural perturbations. The results show that simulations of protein-RNA complexes are challenging and every system should be treated individually. The simulations are affected by numerous factors, including properties of the starting structures (the initially high force field potential energy, resolution limits, conformational averaging, crystal packing, etc.), force field imbalances, and real flexibility of the studied systems. These factors, and thus the simulation behavior, differ from system to system. The structural stability of simulated systems does not correlate with the size of buried interaction surface or experimentally determined binding affinities but reflects the type of protein-RNA recognition. Protein-RNA interfaces involving shape-specific recognition of RNA are more stable than those relying on sequence-specific RNA recognition. The differences between the protein force fields are considerably smaller than the uncertainties caused by sampling and starting structures. The ff12SB improves description of the tyrosine side-chain group, which eliminates some problems associated with tyrosine dynamics.

Institute of Biophysics Academy of Sciences of the Czech Republic Královopolská 135 612 65 Brno Czech Republic

Regional Centre of Advanced Technologies and Materials Department of Physical Chemistry Faculty of Science Palacký University Tř 17 Listopadu 12 771 46 Olomouc Czech Republic

Citace poskytuje Crossref.org

How Binding Site Flexibility Promotes RNA Scanning by TbRGG2 RRM: A Molecular Dynamics Simulation Study

Spontaneous binding of single-stranded RNAs to RRM proteins visualized by unbiased atomistic simulations with a rescaled RNA force field

MD simulations reveal the basis for dynamic assembly of Hfq-RNA complexes

RuvC uses dynamic probing of the Holliday junction to achieve sequence specificity and efficient resolution

Improving the Performance of the Amber RNA Force Field by Tuning the Hydrogen-Bonding Interactions

Molecular basis for the increased affinity of an RNA recognition motif with re-engineered specificity: A molecular dynamics and enhanced sampling simulations study

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

Aromatic side-chain conformational switch on the surface of the RNA Recognition Motif enables RNA discrimination

Structural study of the Fox-1 RRM protein hydration reveals a role for key water molecules in RRM-RNA recognition

Coordination between the polymerase and RNase H activity of HIV-1 reverse transcriptase

Synergy between NMR measurements and MD simulations of protein/RNA complexes: application to the RRMs, the most common RNA recognition motifs