Molecular Dynamics Simulations of Protein RNA Complexes by Using an Advanced Electrostatic Model

Area of Science:

• Computational Biology
• Structural Biology
• Biophysics

Background:

• Protein-RNA interactions are fundamental to RNA biological functions.
• Simulating protein-RNA complexes with molecular dynamics (MD) is challenging due to strong electrostatic forces and complex physical interactions.
• Previous MD simulations often used fixed-charge force fields, with limitations in accurately reproducing some protein-RNA structures.

Purpose of the Study:

• To evaluate the efficacy of the AMOEBA polarizable force field for MD simulations of protein-RNA complexes.
• To investigate if refined force field parameters can improve the stability and accuracy of these simulations.
• To compare simulation results with experimental data, including crystal structures.

Main Methods:

• Performed MD simulations on two representative protein-RNA complexes.
• Utilized the AMOEBA polarizable force field.
• Refined van der Waals parameters to match quantum-mechanical data for base-base and base-amino acid interactions.

Main Results:

• Refined parameters enhanced the stability of the hydrogen-bond network at the protein-RNA interface.
• One complex, previously unstable in simulations, remained stable using the AMOEBA force field.
• Observed reversible hydrogen bond dynamics consistent with crystal structure data, suggesting potential solution vs. crystal structure differences.

Conclusions:

• Polarizable force fields, like AMOEBA, show significant promise for simulating protein-RNA complexes.
• Force field refinement is crucial for improving the accuracy of molecular dynamics simulations in this area.
• Further validation and refinement of polarizable force fields are necessary for broader application.