RNA/peptide binding driven by electrostatics -- Insight from bi-directional pulling simulations
arXiv:1307.5565 · doi:10.1021/ct3009914
Abstract
RNA/protein interactions play crucial roles in controlling gene expression. They are becoming important targets for pharmaceutical applications. Due to RNA flexibility and to the strength of electrostatic interactions, standard docking methods are insufficient. We here present a computational method which allows studying the binding of RNA molecules and charged peptides with atomistic, explicit-solvent molecular dynamics. In our method, a suitable estimate of the electrostatic interaction is used as an order parameter (collective variable) which is then accelerated using bi-directional pulling simulations. Since the electrostatic interaction is only used to enhance the sampling, the approximations used to compute it do not affect the final accuracy. The method is employed to characterize the binding of TAR RNA from HIV-1 and a small cyclic peptide. Our simulation protocol allows blindly predicting the binding pocket and pose as well as the binding affinity. The method is general and could be applied to study other electrostatics-driven binding events.
Reprinted (adapted) with permission from J. Chem. Theory Comput., 2013, 9 (3) 1720 (2013). Copyright (2013) American Chemical Society
References in corpus (5)
- Canonical sampling through velocity-rescaling
- Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
- PLUMED: a portable plugin for free-energy calculations with molecular dynamics
- Optimized free energies from bidirectional single-molecule force spectroscopy
- RNA unwinding from reweighted pulling simulations
Cited by in corpus (8)
- PLUMED 2: New feathers for an old bird
- Enhanced Sampling in Molecular Dynamics Using Metadynamics, Replica-Exchange, and Temperature-Acceleration
- Analyzing and biasing simulations with PLUMED
- Ligand-induced stabilization of the aptamer terminal helix in the add adenine riboswitch
- Accurate multiple time step in biased molecular simulations
- ATP dependent NS3 helicase interaction with RNA: insights from molecular simulations
- Kissing loop interaction in adenine riboswitch: insights from umbrella sampling simulations
- maze: Heterogeneous Ligand Unbinding along Transient Protein Tunnels