Sequence specificity in DNA-drug intercalation: MD simulation and density functional theory approaches

Area of Science:

• Biophysics
• Computational Chemistry
• Molecular Biology

Background:

• DNA is a critical target for cancer therapy.
• Drug-DNA intercalation is a key mechanism but challenging to study due to DNA deformation.
• Daunomycin is an anticancer drug that interacts with DNA.

Purpose of the Study:

• To theoretically investigate the intercalation process of daunomycin into DNA.
• To determine the binding preference of daunomycin for specific DNA sequences.
• To understand the biophysical and energetic factors governing daunomycin-DNA interactions.

Main Methods:

• Molecular dynamics simulations of daunomycin-DNA complexes across various sequences.
• Free-energy analysis to assess binding affinity.
• Density Functional Theory (DFT) calculations with an extended counterpoise method for Basis Set Superposition Error (BSSE) correction.

Main Results:

• Classical energy analyses suggest daunomycin favors binding to TC/GA DNA sequences.
• DFT-based calculations, after BSSE correction, confirm TC/GA as the preferred binding sequence.
• The energy penalty associated with base pair un-stacking also supports TC/GA preference.

Conclusions:

• Daunomycin exhibits a sequence-specific binding preference for TC/GA sites within the DNA double helix.
• The study provides a refined computational approach for analyzing drug-DNA interactions, accounting for BSSE.
• Understanding these binding preferences can inform the design of more effective DNA-targeting cancer therapies.