Skip to content

Electrochemical and theoretical studies of the interaction between anticancer drug ponatinib and dsDNA

MetadataDetails
Publication Date2024-01-27
JournalScientific Reports
AuthorsSylwia Smarzewska, Anna Ignaczak, Kamila Koszelska
InstitutionsUniversity of ƁódĆș
Citations9
AnalysisFull AI Review Included
  • Core Objective: Investigated the interaction mechanism between the anticancer drug ponatinib (PNT), a third-generation tyrosine kinase inhibitor, and double-stranded DNA (dsDNA).
  • Methodology: Combined electrochemical techniques (Square Wave Voltammetry, SWV) using a Boron-Doped Diamond Electrode (BDDE) with hierarchical quantum-chemical calculations (PM7 and DFT).
  • Interaction Type: Voltammetric and theoretical results strongly suggest that the primary binding mechanism is groove binding, with intercalation largely excluded due to minimal peak potential shifts observed.
  • Binding Site Preference (Theoretical): DFT calculations identified the dsDNA major groove (MaG) as the energetically preferred site for PNT complexation, regardless of the PNT conformer (stretched or lowest energy).
  • Nucleobase Specificity: PNT interacts preferentially with electrochemically detectable nucleobases, specifically guanine (dGua) and adenine (dAdo) residues, confirmed by both pH-dependent voltammetry and hydrogen bonding analysis.
  • Binding Strength: Computational analysis showed that hydrogen bonds formed between PNT and guanine residues are geometrically stronger than those formed with adenine.
  • Electrode Kinetics: PNT oxidation on the BDDE was determined to be an irreversible process, controlled by adsorption in acetate buffer (pH 4.7) and a mixed diffusion-adsorption process in PBS (pH 7.4).
ParameterValueUnitContext
Working ElectrodeBoron-Doped Diamond (BDDE)3 mm diameterConventional three-electrode system
Physiological pH7.4pHPhosphate-Buffered Saline (PBS)
Acidic pH4.7pHAcetate Buffer
SWV Amplitude30mVVoltammetric interaction studies
SWV Frequency25HzVoltammetric interaction studies
SWV Step Potential4mVVoltammetric interaction studies
CV Scan Rate Range50 to 500mV s-1Kinetics study
PNT Oxidation Peak 1 (Ep-pH slope)13mV pH-1Complex, multi-step reaction
PNT Oxidation Peak 2 (Ep-pH slope)51mV pH-1Suggests equal proton/electron participation
Log Ip - Log v Slope (PBS)0.71N/AMixed diffusion-adsorption control
Log Ip - Log v Slope (Acetate)0.97N/AAdsorption control (close to 1.0)
Most Stable Ecompl (119D:PNT_ST)-65.6kcal mol-1DFT, Major Groove (MaG) binding
Most Stable Ecompl (1BNA:PNT_LE)-56.6kcal mol-1DFT, Major Groove (MaG) binding
Guanine O
H-N H-Bond Distance1.96 to 2.37AngstromStronger hydrogen bonding
Adenine H
O H-Bond Distance2.20AngstromWeaker hydrogen bonding

The study employed a rigorous, three-stage hierarchical approach for computational modeling, coupled with advanced electrochemical analysis using a BDDE.

  1. Electrode Preparation: BDDE surface was polished using alumina slurry, followed by thorough rinsing with distilled/deionized water to ensure a clean, active surface.
  2. PNT Characterization: Cyclic Voltammetry (CV) and Square Wave Voltammetry (SWV) were used to determine the irreversible oxidation behavior of PNT across a wide pH range (1.7-9.0).
  3. Interaction Study (Incubation): PNT (5.0 ”mol L-1) and dsDNA (80 mg L-1) were mixed in acetate buffer (pH 4.7) or PBS (pH 7.4) and incubated at room temperature for varying periods (10, 30, 60, 90 min).
  4. Concentration Variation Study: SWV was performed immediately after adding increasing amounts of dsDNA (10-80 mg L-1) to a fixed concentration of PNT, analyzing changes in peak current and potential shift.
  1. Step 1: Molecular Mechanics (HyperChem):
    • Initial models of dsDNA:PNT complexes were generated using two B-DNA dodecamers (1BNA and 119D).
    • PNT conformers (PNT_LE and PNT_ST) were systematically rotated (30° increments around X, Y, Z axes) across four binding sites: External Binding (ExB), Major Groove (MaG), Minor Groove (MiG), and Intercalation (InC).
  2. Step 2: Semiempirical Optimization (MOPAC):
    • All generated structures were optimized using the PM7 method, incorporating solvent effects via the Conductor-like Screening Model (COSMO) in water.
    • Complexation enthalpies (Hcompl) were calculated to identify the most stable structure for each site.
  3. Step 3: DFT Calculation (Gaussian 16):
    • The most stable structures from Step 2 underwent partial re-optimization (PNT relaxed, dsDNA frozen) using the high-accuracy M062X-GD3 functional with the 6-31G(d,p) basis set.
    • Solvent effects were modeled using the Polarizable Continuum Model (PCM) in water.
    • Final complexation energies (Ecompl) were computed to confirm binding stability and identify preferred sites.
  • Electrochemical Drug Sensing: Utilizing the high sensitivity and stability of the BDDE for the trace analysis and quantification of tyrosine kinase inhibitors (TKIs) like PNT in complex biological matrices (plasma, urine), offering a “green chemistry” alternative to traditional chromatographic methods.
  • Genotoxicity Screening: Applying electrochemical biosensors (BDDE/DNA) to rapidly screen new drug candidates for potential DNA damage or interaction mechanisms (groove binding vs. intercalation), crucial for early-stage pharmaceutical safety assessment.
  • Personalized Medicine & TKI Monitoring: Development of robust, portable electrochemical devices for clinical monitoring of TKI drug concentrations in cancer patients (e.g., CML, ALL), ensuring optimal therapeutic windows and minimizing toxic side effects.
  • Advanced Biosensor Design: Leveraging the detailed structural and energetic data from DFT calculations to rationally design and optimize DNA-based biosensors that selectively recognize and bind specific drug molecules or metabolites based on groove geometry and nucleobase preference.
  • Material Science for Electroanalysis: Confirms the utility of the BDDE as a stable, wide-potential-window platform for studying complex, irreversible redox processes involving large organic molecules and biomolecules, suitable for harsh or physiological environments.