Revisiting the association of cationic groove-binding drugs to DNA using a Poisson-Boltzmann approach
Related Experiment Video
Updated: Jun 10, 2026

Incorporating Target Protein Structure Flexibility and Dynamics in Computational Drug Discovery Using Ensemble-Based Docking Analysis
Published on: June 20, 2025
Marcia O Fenley1, Robert C Harris, B Jayaram
1Department of Physics and Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida, USA. mfenley@sb.fsu.edu
Abstract:
Proper modeling of nonspecific salt-mediated electrostatic interactions is essential to understanding the binding of charged ligands to nucleic acids. Because the linear Poisson-Boltzmann equation (PBE) and the more approximate generalized Born approach are applied routinely to nucleic acids and their interactions with charged ligands, the reliability of these methods is examined vis-à-vis an efficient nonlinear PBE method. For moderate salt concentrations, the negative derivative, SK(pred), of the electrostatic binding free energy, DeltaG(el), with respect to the logarithm of the 1:1 salt concentration, [M(+)], for 33 cationic minor groove drugs binding to AT-rich DNA sequences is shown to be consistently negative and virtually constant over the salt range considered (0.1-0.4 M NaCl). The magnitude of SK(pred) is approximately equal to the charge on the drug, as predicted by counterion condensation theory (CCT) and observed in thermodynamic binding studies. The linear PBE is shown to overestimate the magnitude of SK(pred), whereas the nonlinear PBE closely matches the experimental results. The PBE predictions of SK(pred) were not correlated with DeltaG(el) in the presence of a dielectric discontinuity, as would be expected from the CCT. Because this correlation does not hold, parameterizing the PBE predictions of DeltaG(el) against the reported experimental data is not possible. Moreover, the common practice of extracting the electrostatic and nonelectrostatic contributions to the binding of charged ligands to biopolyelectrolytes based on the simple relation between experimental SK values and the electrostatic binding free energy that is based on CCT is called into question by the results presented here. Although the rigid-docking nonlinear PB calculations provide reliable predictions of SK(pred), at least for the charged ligand-nucleic acid complexes studied here, accurate estimates of DeltaG(el) will require further development in theoretical and experimental approaches.
More Related Videos
09:17Structure-Based Simulation and Sampling of Transcription Factor Protein Movements along DNA from Atomic-Scale Stepping to Coarse-Grained Diffusion
Published on: March 1, 2022
10:05Genome-wide Mapping of Drug-DNA Interactions in Cells with COSMIC (Crosslinking of Small Molecules to Isolate Chromatin)
Published on: January 20, 2016
Related Concept Videos
Protein-Drug Binding: Mechanism and Kinetics
Various forces drive these interactions, including hydrogen bonds, hydrophobic interactions, ionic bonds, electrostatic interactions, and van der Waals forces. These bonds enable drugs to bind to specific sites on proteins,...
Pore Transport and Ion-Pair Transport
Pore transport, also known as convective transport, is a process where small molecules like urea, water, and sugars rapidly cross cell membranes as though there were channels or pores in the membrane. Although direct microscopic evidence is limited but the concept of pores or channels is widely accepted based on physiological evidence. Despite the lack of direct microscopic...
Drug-Receptor Bonds
In...
Maxam-Gilbert Sequencing
Challenges of the Maxam-Gilbert Method
The...