Molecular recognition of DNA by ligands

Roughness and complexity of the free energy profile

Wenwei Zheng, Attlio Vittorio Vargiu, Mary A. Rohrdanz, Paolo Carloni, Cecilia Clementi

Research output: Contribution to journalArticle

12 Citations (Scopus)

Abstract

Understanding the molecular mechanism by which probes and chemotherapeutic agents bind to nucleic acids is a fundamental issue in modern drug design. From a computational perspective, valuable insights are gained by the estimation of free energy landscapes as a function of some collective variables (CVs), which are associated with the molecular recognition event. Unfortunately the choice of CVs is highly non-trivial because of DNA's high flexibility and the presence of multiple association-dissociation events at different locations and/or sliding within the grooves. Here we have applied a modified version of Locally-Scaled Diffusion Map (LSDMap), a nonlinear dimensionality reduction technique for decoupling multiple-timescale dynamics in macromolecular systems, to a metadynamics-based free energy landscape calculated using a set of intuitive CVs. We investigated the binding of the organic drug anthramycin to a DNA 14-mer duplex. By performing an extensive set of metadynamics simulations, we observed sliding of anthramycin along the full-length DNA minor groove, as well as several detachments from multiple sites, including the one identified by X-ray crystallography. As in the case of equilibrium processes, the LSDMap analysis is able to extract the most relevant collective motions, which are associated with the slow processes within the system, i.e., ligand diffusion along the minor groove and dissociation from it. Thus, LSDMap in combination with metadynamics (and possibly every equivalent method) emerges as a powerful method to describe the energetics of ligand binding to DNA without resorting to intuitive ad hoc reaction coordinates.

Original languageEnglish (US)
Article number145102
JournalJournal of Chemical Physics
Volume139
Issue number14
DOIs
StatePublished - Oct 14 2013
Externally publishedYes

Fingerprint

Molecular recognition
Free energy
Anthramycin
roughness
deoxyribonucleic acid
Surface roughness
free energy
grooves
Ligands
ligands
drugs
DNA
profiles
sliding
dissociation
X ray crystallography
nucleic acids
detachment
Pharmaceutical Preparations
decoupling

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry

Cite this

Molecular recognition of DNA by ligands : Roughness and complexity of the free energy profile. / Zheng, Wenwei; Vargiu, Attlio Vittorio; Rohrdanz, Mary A.; Carloni, Paolo; Clementi, Cecilia.

In: Journal of Chemical Physics, Vol. 139, No. 14, 145102, 14.10.2013.

Research output: Contribution to journalArticle

Zheng, Wenwei ; Vargiu, Attlio Vittorio ; Rohrdanz, Mary A. ; Carloni, Paolo ; Clementi, Cecilia. / Molecular recognition of DNA by ligands : Roughness and complexity of the free energy profile. In: Journal of Chemical Physics. 2013 ; Vol. 139, No. 14.
@article{6fc89aeb84514b2aa537be5e544255c8,
title = "Molecular recognition of DNA by ligands: Roughness and complexity of the free energy profile",
abstract = "Understanding the molecular mechanism by which probes and chemotherapeutic agents bind to nucleic acids is a fundamental issue in modern drug design. From a computational perspective, valuable insights are gained by the estimation of free energy landscapes as a function of some collective variables (CVs), which are associated with the molecular recognition event. Unfortunately the choice of CVs is highly non-trivial because of DNA's high flexibility and the presence of multiple association-dissociation events at different locations and/or sliding within the grooves. Here we have applied a modified version of Locally-Scaled Diffusion Map (LSDMap), a nonlinear dimensionality reduction technique for decoupling multiple-timescale dynamics in macromolecular systems, to a metadynamics-based free energy landscape calculated using a set of intuitive CVs. We investigated the binding of the organic drug anthramycin to a DNA 14-mer duplex. By performing an extensive set of metadynamics simulations, we observed sliding of anthramycin along the full-length DNA minor groove, as well as several detachments from multiple sites, including the one identified by X-ray crystallography. As in the case of equilibrium processes, the LSDMap analysis is able to extract the most relevant collective motions, which are associated with the slow processes within the system, i.e., ligand diffusion along the minor groove and dissociation from it. Thus, LSDMap in combination with metadynamics (and possibly every equivalent method) emerges as a powerful method to describe the energetics of ligand binding to DNA without resorting to intuitive ad hoc reaction coordinates.",
author = "Wenwei Zheng and Vargiu, {Attlio Vittorio} and Rohrdanz, {Mary A.} and Paolo Carloni and Cecilia Clementi",
year = "2013",
month = "10",
day = "14",
doi = "10.1063/1.4824106",
language = "English (US)",
volume = "139",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "14",

}

TY - JOUR

T1 - Molecular recognition of DNA by ligands

T2 - Roughness and complexity of the free energy profile

AU - Zheng, Wenwei

AU - Vargiu, Attlio Vittorio

AU - Rohrdanz, Mary A.

AU - Carloni, Paolo

AU - Clementi, Cecilia

PY - 2013/10/14

Y1 - 2013/10/14

N2 - Understanding the molecular mechanism by which probes and chemotherapeutic agents bind to nucleic acids is a fundamental issue in modern drug design. From a computational perspective, valuable insights are gained by the estimation of free energy landscapes as a function of some collective variables (CVs), which are associated with the molecular recognition event. Unfortunately the choice of CVs is highly non-trivial because of DNA's high flexibility and the presence of multiple association-dissociation events at different locations and/or sliding within the grooves. Here we have applied a modified version of Locally-Scaled Diffusion Map (LSDMap), a nonlinear dimensionality reduction technique for decoupling multiple-timescale dynamics in macromolecular systems, to a metadynamics-based free energy landscape calculated using a set of intuitive CVs. We investigated the binding of the organic drug anthramycin to a DNA 14-mer duplex. By performing an extensive set of metadynamics simulations, we observed sliding of anthramycin along the full-length DNA minor groove, as well as several detachments from multiple sites, including the one identified by X-ray crystallography. As in the case of equilibrium processes, the LSDMap analysis is able to extract the most relevant collective motions, which are associated with the slow processes within the system, i.e., ligand diffusion along the minor groove and dissociation from it. Thus, LSDMap in combination with metadynamics (and possibly every equivalent method) emerges as a powerful method to describe the energetics of ligand binding to DNA without resorting to intuitive ad hoc reaction coordinates.

AB - Understanding the molecular mechanism by which probes and chemotherapeutic agents bind to nucleic acids is a fundamental issue in modern drug design. From a computational perspective, valuable insights are gained by the estimation of free energy landscapes as a function of some collective variables (CVs), which are associated with the molecular recognition event. Unfortunately the choice of CVs is highly non-trivial because of DNA's high flexibility and the presence of multiple association-dissociation events at different locations and/or sliding within the grooves. Here we have applied a modified version of Locally-Scaled Diffusion Map (LSDMap), a nonlinear dimensionality reduction technique for decoupling multiple-timescale dynamics in macromolecular systems, to a metadynamics-based free energy landscape calculated using a set of intuitive CVs. We investigated the binding of the organic drug anthramycin to a DNA 14-mer duplex. By performing an extensive set of metadynamics simulations, we observed sliding of anthramycin along the full-length DNA minor groove, as well as several detachments from multiple sites, including the one identified by X-ray crystallography. As in the case of equilibrium processes, the LSDMap analysis is able to extract the most relevant collective motions, which are associated with the slow processes within the system, i.e., ligand diffusion along the minor groove and dissociation from it. Thus, LSDMap in combination with metadynamics (and possibly every equivalent method) emerges as a powerful method to describe the energetics of ligand binding to DNA without resorting to intuitive ad hoc reaction coordinates.

UR - http://www.scopus.com/inward/record.url?scp=84886891303&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84886891303&partnerID=8YFLogxK

U2 - 10.1063/1.4824106

DO - 10.1063/1.4824106

M3 - Article

VL - 139

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 14

M1 - 145102

ER -