Molecular spectroscopic and molecular simulation studies on the interaction of oral contraceptive drug Ormeloxifene with CT-DNA

The interaction between oral contraceptive drug Ormeloxifene (ORM) and calf thymus DNA (CT-DNA) was studied using UV-Vis, fluorescence, circular dichroism (CD) and ¹H NMR spectral techniques under physiological buffer (pH 7.4). Competitive binding assays with ethidium bromide (EB) and Hoechst 33258, viscosity measurements, KI quenching studies, molecular docking and metadynamics simulation studies were also substantiated the spectroscopic results. ORM is found to binds in the minor groove of CT-DNA as evidenced by: (1) non-displacement of EB from EB/CT-DNA complex; (2) appreciable displacement of Hoechst 33258 from its CT-DNA complex; (3) slight alteration in the CD signal; (4) small shifts (Δδ < 0.033 ppm) without broadening in ¹H NMR signals and (5) the nearly equal extent of quenching of fluorescence of ORM by KI in the absence and presence of CT-DNA. Negative values of both enthalpy and entropy changes pointed out that the interaction between ORM and CT-DNA is governed mainly by H-bonding and van der Waals forces. Negative free energy change suggested a spontaneous interaction between ORM and CT-DNA. The free energy landscape of the binding process was computed using metadynamics simulation. The simulation study results disclosed that ORM binds to the minor groove of DNA through H-bonding and π-π stacking interactions. The results of molecular docking and simulation studies corroborate the available experimental data.

Water-Controlled Keto-Enol Tautomerization of a Prebiotic Nucleobase

Barbituric acid is believed to be a proto-RNA nucleobase that was used for biological information transfer on prebiotic earth before DNA and RNA in their present forms evolved. Nucleobases have various tautomeric forms and the relative stability of these forms is critical to their biological function. It has been shown that barbituric acid has a tri-keto form in the gas phase and an enol form in the solid state. However, its dominant tautomeric form in aqueous medium that is most relevant for biology has been investigated only to a limited extent and the findings are inconclusive. We have used multiple approaches, namely, molecular dynamics, quantum chemistry, NMR, and IR spectroscopy to determine the most stable tautomer of barbituric acid in solution. We find a delicate balance in the stability of the two tautomers, tri-keto and enol, which is tipped toward the enol as the extent of solvation by water increases.

Influence of conventional hydrogen bonds in the intercalation of phenanthroline derivatives with DNA: The important role of the sugar and phosphate backbone

The influence of hydrogen bonds in model intercalated systems between guanine-cytosine and adenine-thymine DNA base pairs (bps) was analyzed with the popular intercalator 1,10-phenanthroline (phen) and derivatives obtained by substitution with OH and NH2 groups in positions 4 and 7. Semiempirical and Density Functional Theory (DFT) methods were used both including dispersion effects: PM6-DH2, M06-2X and B3LYP-D3 along with the recently developed near linear-scaling coupled cluster method DLPNO-CCSD(T) for benchmark calculations. Our results given by QTAIM and non-covalent interaction analysis confirmed the existence of hydrogen bonds created by OH and NH2 . The trends in the energy decomposition analysis for the interaction energy, ΔEint , showed that the ΔEelstat contributions are equal or even a little bit higher than the values for ΔEdisp . Such important ΔEelstat attractive contribution comes mainly from the conventional hydrogen bonds formed by OH and NH2 functional groups with DNA not only with bps but specially with the sugar and phosphate backbone. This behavior is very different from that of phen and other classical intercalators that cannot form conventional hydrogen bonds, where the ΔEdisp is the most important attractive contribution to the ΔEint . The inclusion of explicit water molecules in molecular dynamics simulations showed, as a general trend, that the hydrogen bonds with the bps disappear during the simulations but those with the sugar and phosphate backbone remain in time, which highlights the important role of the sugar and phosphate backbone in the stabilization of these systems.

Advances in heterocycles as DNA intercalating cancer drugs

The insertion of a molecule between the bases of DNA is known as intercalation. A molecule is able to interact with DNA in different ways. DNA intercalators are generally aromatic, planar, and polycyclic. In chemotherapeutic treatment, to suppress DNA replication in cancer cells, intercalators are used. In this article, we discuss the anticancer activity of 10 intensively studied DNA intercalators as drugs. The list includes proflavine, ethidium bromide, doxorubicin, dactinomycin, bleomycin, epirubicin, mitoxantrone, ellipticine, elinafide, and echinomycin. Considerable structural diversities are seen in these molecules. Besides, some examples of the metallo-intercalators are presented at the end of the chapter. These molecules have other crucial properties that are also useful in the treatment of cancers. The successes and limitations of these molecules are also presented.

From groove binding to intercalation: unravelling the weak interactions and other factors modulating the modes of interaction between methylated phenanthroline-based drugs and duplex DNA

Several antitumor drugs base their cytotoxicity on their capacity to intercalate between base pairs of DNA. Nevertheless, it has been established that the mechanism of intercalation of drugs in DNA starts with the prior groove binding mode of interaction of the drug with DNA. Sometimes, for some kind of flat small molecules, groove binding does not produce any cytotoxic effect and the fast transition of such flat small molecules to the cytotoxic intercalation mode is desirable. This is the case of methylated phenanthroline (phen) derivatives, where, changes in the substitution in the position and number of methyl groups determine their capability as cytotoxic compounds and, therefore, it is a way for the modulation of cytotoxic effects. In this work, we studied this modulation by means of the interaction of the [Pt(en)(phen)]²⁺ complex and several derivatives by methylation of phen in different number and position and the d(GTCGAC)₂ DNA hexamer via groove binding using PM6-DH2 and DFT-D methods. The analysis of the geometries, electronic structure and energetics of the studied systems was compared to experimental works to gain insight into the relation structure-interaction for the studied systems with cytotoxicity. The trends are explained by means of the Non-Covalent Interaction (NCI) index, the Energy Decomposition Analysis (EDA) and solvation contributions. Our results are in agreement with the experiments, in which the methylation of position 4 of phen seems to favour the interaction via groove binding thus making the transition to the intercalation cytotoxic mode difficult. Looking at the NCI results, these interactions come not only from the CH/π and CH/n interactions of the methyl group in position 4 but also from the ethylenediamine (en) ligand, whose orientation in the Pt complex was found in such a way that it produces a high number of weak interactions with DNA, especially with the sugar and phosphate backbone.

New Insights on the Interaction of Phenanthroline Based Ligands and Metal Complexes and Polyoxometalates with Duplex DNA and G-Quadruplexes

This work provides new insights from our team regarding advances in targeting canonical and non-canonical nucleic acid structures. This modality of medical treatment is used as a form of molecular medicine specifically against the growth of cancer cells. Nevertheless, because of increasing concerns about bacterial antibiotic resistance, this medical strategy is also being explored in this field. Up to three strategies for the use of DNA as target have been studied in our research lines during the last few years: (1) the intercalation of phenanthroline derivatives with duplex DNA; (2) the interaction of metal complexes containing phenanthroline with G-quadruplexes; and (3) the activity of Mo polyoxometalates and other Mo-oxo species as artificial phosphoesterases to catalyze the hydrolysis of phosphoester bonds in DNA. We demonstrate some promising computational results concerning the favorable interaction of these small molecules with DNA that could correspond to cytotoxic effects against tumoral cells and microorganisms. Therefore, our results open the door for the pharmaceutical and medical applications of the compounds we propose.

Role of intercalation in the electrical properties of nucleic acids for use in molecular electronics

Intercalating ds-DNA/RNA with small molecules can play an essential role in controlling the electron transmission probability for molecular electronics applications such as biosensors, single-molecule transistors, and data storage. However, its applications are limited due to a lack of understanding of the nature of intercalation and electron transport mechanisms. We addressed this long-standing problem by studying the effect of intercalation on both the molecular structure and charge transport along the nucleic acids using molecular dynamics simulations and first-principles calculations coupled with the Green's function method, respectively. The study on anthraquinone and anthraquinone-neomycin conjugate intercalation into short nucleic acids reveals some universal features: (1) the intercalation affects the transmission by two mechanisms: (a) inducing energy levels within the bandgap and (b) shifting the location of the Fermi energy with respect to the molecular orbitals of the nucleic acid, (2) the effect of intercalation was found to be dependent on the redox state of the intercalator: while oxidized anthraquinone decreases, reduced anthraquinone increases the conductance, and (3) the sequence of the intercalated nucleic acid further affects the transmission: lowering the AT-region length was found to enhance the electronic coupling of the intercalator with GC bases, hence yielding an increase of more than four times in conductance. We anticipate our study to inspire designing intercalator-nucleic acid complexes for potential use in molecular electronics via creating a multi-level gating effect.

Atomistic De-novo Inhibitor Generation-Guided Drug Repurposing for SARS-CoV-2 Spike Protein with Free-Energy Validation by Well-Tempered Metadynamics

Computational drug design is increasingly becoming important with new and unforeseen diseases like COVID-19. In this study, we present a new computational de novo drug design and repurposing method and applied it to find plausible drug candidates for the receptor binding domain (RBD) of SARS-CoV-2 (COVID-19). Our study comprises three steps: atom-by-atom generation of new molecules around a receptor, structural similarity mapping to existing approved and investigational drugs, and validation of their binding strengths to the viral spike proteins based on rigorous all-atom, explicit-water well-tempered metadynamics free energy calculations. By choosing the receptor binding domain of the viral spike protein, we showed that some of our new molecules and some of the repurposable drugs have stronger binding to RBD than hACE2. To validate our approach, we also calculated the free energy of hACE2 and RBD, and found it to be in an excellent agreement with experiments. These pool of drugs will allow strategic repurposing against COVID-19 for a particular prevailing conditions.

Concentration-dependent mode of binding of drug oxatomide with DNA: multi-spectroscopic, voltammetric and metadynamics simulation analysis

The interaction between antihistaminic drug oxatomide (OXT) and calf-thymus DNA (CT-DNA) has been investigated in a physiological buffer (pH 7.4) using UV-Vis, fluorescence, ¹H NMR and circular dichroism spectral techniques coupled with viscosity measurements, KI quenching, voltammetry and in silico molecular modeling studies. OXT binds with CT-DNA in a concentration-dependent manner. At a lower [Drug]/[CT-DNA] molar ratio (0.6-0.1), OXT intercalates into the base pairs of CT-DNA, while at a higher [Drug]/[CT-DNA] molar ratio (13-6), the drug binds in the minor grooves of CT-DNA. The binding constants for the interaction are found to be in the order of 10³-10⁵ M^-1, and the groove binding mode of interaction exhibits a slightly higher binding constant than that of intercalative mode. Thermodynamic analysis of binding constants at three different temperatures suggests that both these modes of binding are mainly driven by hydrophobic interactions (ΔH^o > 0 and ΔS^o > 0). Voltammetric investigations indicate that the electro-reduction of OXT is an adsorption controlled process and shifts in reduction peak potentials reiterate the concentration-dependent mode of binding of the drug with CT-DNA. The free energy landscape obtained at the all-atom level, using metadynamics simulation studies, revealed two major binding forces: partial intercalation and minor groove binding, which corroborate well with the experimental results.Communicated by Ramaswamy H. Sarma.

A computational study of intercalation of streptozotocin (STZ) into DNA base pairs

Deoxyribonucleic acid (DNA) drug intercalation is a well-known phenomenon for the treatment of cancer. Streptozotocin (STZ) is a drug agent containing toxic properties that make it good in the pancreatic cancer. The main objective of this study is the intercalation of the anticancer drug into the stacked base pair of DNA sequence with ATGC using a density functional theory (DFT) code named as ADF-Molecule. ADF code implements DFT using the Slater-type orbitals (STO) for computational analysis of atomic and molecular structures. All the calculations were carried out with the GGA and hybrid exchange correlation functional with TZ2P basis sets. It was captivatingly studied that during the intercalation process, the bonds between the DNA base pairs broken. Moreover, during the process of intercalation, the free radicals are considered responsible for disturbance in the base configurations. It was determined that the disturbances that occurred in the base pairs lead to discontinuity in the replication of that particular sequence in the DNA strand.

Chiral molecular imprinted sensor for highly selective determination of D-carnitine in enantiomers via dsDNA-assisted conformation immobilization

In this paper, a novel approach was established on the basis of a molecularly imprinted technique with the aid of double-stranded deoxyribonucleic acid (dsDNA) embedded in a molecularly imprinted polymer (MIP) membrane as a new functional unit with chiral recognition for highly specific chiral recognition. The chiral molecules were immobilized and anchored in the cavities of the MIP membrane on the basis of the three-dimensional structure of a molecule determined by the functional groups, spatial characterization of the cavities of MIPs, and the spatial orientation with dsDNA embedded in MIPs. D-carnitine was selected as an example of a chiral molecular template, which intercalated into dsDNA immobilized on the gold electrode surface to form dsDNA-D-carnitine complex, and then the complex was embedded in the MIP during electropolymerization. After elution, the stereo-selective cavities were obtained. Our findings have shown that AAAA-TTTT base sequence had high affinity for D-carnitine intercalation. Combined with the electrochemical detection method, MIP sensor was prepared. The selectivity of the MIP sensor to ultratrace D-carnitine was significantly improved; the sensor had remarkable stereo-selectivity and highly chiral specific recognition to D-carnitine, and L-carnitine with a concentration of 10,000 times D-carnitine did not interfere with the detection of D-carnitine in the assay of raceme. The sensor also exhibited high sensitivity to ultratrace D-carnitine determination with a linear response to the concentration of D-carnitine in the range of 3.0 × 10^-16 mol/L to 4.0 × 10^-13 mol/L, with a detection limit of 2.24 × 10^-16 mol/L. The mechanism of chiral recognition was studied, and result showed that apart from the recognition effect of imprinted cavities, dsDNA provided chiral selectivity to the spatial orientation of chiral molecules via the intercalation of chiral molecules with dsDNA and electrostatic interaction with groups of DNA base.

Dynamical Recrossing in the Intercalation Process of the Anticancer Agent Proflavine into DNA

Intercalation into DNA is the interaction mode of some anthracycline antibiotics. Recently, the molecular mechanism of this process was explored using the static free energy landscape. Here we explore the dynamical effects in the intercalation of proflavine into DNA by calculating the transmission coefficient κ-providing a measure of the departure from transition state theory for the reaction rate constant-by examination of the recrossing events at the transition state. For that purpose, we first found the accurate transition state of this complex system-as judged by a committor analysis-using a set of all-atom simulations of total length 6.3 ms. In a subsequent calculation of the transmission coefficient κ in another extensive set of simulations the small value κ = 0.1 was found, indicating a significant departure from TST. Comparison of this result with Grote-Hynes and Kramers theories shows that neither theory is able to capture this complex system's recrossing events; the source of this striking failure is discussed, as are related aspects of the mechanism. This study suggests that, for biomolecular processes similar to this, dynamical effects essential for the process are complex in nature and require novel approaches for their elucidation.

Medicinal chemistry of acridine and its analogues

'Acridine' along with its functional analogue 'Acridone' is the most privileged pharmacophore in medicinal chemistry with diverse applications ranging from DNA intercalators, endonuclease mimics, ratiometric selective ion sensors, and P-glycoprotein inhibitors in countering the multi-drug resistance, enzyme inhibitors, and reversals of neurodegenerative disorders. Their interaction with DNA and ability of selectively identifying numerous biologically useful ions has cemented exploitability of the acridone nucleus in modern day therapeutics. Additionally, most derivatives and salts of acridine are planar, crystalline, and stable displaying a strong fluorescence which, when coupled with their marked bio selectivity and low cytotoxicity, enables the studying and monitoring of several biochemical, metabolic, and pharmacological processes. In this review, a detailed picture covering the important therapeutic aspects of the acridone nucleus and its functional analogues is discussed.

Anomalous water dynamics at surfaces and interfaces: synergistic effects of confinement and surface interactions

In nature, water is often found in contact with surfaces that are extended on the scale of molecule size but small on a macroscopic scale. Examples include lipid bilayers and reverse micelles as well as biomolecules like proteins, DNA and zeolites, to name a few. While the presence of surfaces and interfaces interrupts the continuous hydrogen bond network of liquid water, confinement on a mesoscopic scale introduces new features. Even when extended on a molecular scale, natural and biological surfaces often have features (like charge, hydrophobicity) that vary on the scale of the molecular diameter of water. As a result, many new and exotic features, which are not seen in the bulk, appear in the dynamics of water close to the surface. These different behaviors bear the signature of both water-surface interactions and of confinement. In other words, the altered properties are the result of the synergistic effects of surface-water interactions and confinement. Ultrafast spectroscopy, theoretical modeling and computer simulations together form powerful synergistic approaches towards an understanding of the properties of confined water in such systems as nanocavities, reverse micelles (RMs), water inside and outside biomolecules like proteins and DNA, and also between two hydrophobic walls. We shall review the experimental results and place them in the context of theory and simulations. For water confined within RMs, we discuss the possible interference effects propagating from opposite surfaces. Similar interference is found to give rise to an effective attractive force between two hydrophobic surfaces immersed and kept fixed at a separation of d, with the force showing an exponential dependence on this distance. For protein and DNA hydration, we shall examine a multitude of timescales that arise from frustration effects due to the inherent heterogeneity of these surfaces. We pay particular attention to the role of orientational correlations and modification of the same due to interaction with the surfaces.

New insights into the mechanism of action of pyrazolo[1,2-a]benzo[1,2,3,4]tetrazin-3-one derivatives endowed with anticancer potential

Due to the scarce biological profile, the pyrazolo[1,2-a]benzo[1,2,3,4]tetrazine-3-one scaffold (PBT) has been recently explored as promising core for potential anticancer candidates. Several suitably decorated derivatives (PBTs) exhibited antiproliferative activity in the low-micromolar range associated with apoptosis induction and cell cycle arrest on S phase. Herein, we selected the most active derivatives and submitted them to further biological explorations to deepen the mechanism of action. At first, a DNA targeting is approached by means of flow Linear Dichroism experiments so as to evaluate how small planar molecules might interact with DNA, including the interference with the catalytic cycle of topoisomerase II and the influence on the cleavable complex stabilization (poisoning effect). In support of the experimental data, in silico studies have been achieved to better understand the chemical space of the interactions. Interestingly some meaningful structural features, useful for further developments, were found. The 8,9-di-Cl substituted derivative revealed as the most effective in the intercalative process, as well as on the inhibition of catalytic activity of topoisomerase II. Predicted ADME studies confirm that PBTs are promising as potential drug candidates.

Effects of oxygenation on the intercalation of 1,10-phenanthroline-5,6/4,7-dione between DNA base pairs: a computational study

The effects of oxygen in positions 4,7 and 5,6 of phenanthroline have been studied computationally when this ligand intercalates between DNA base pairs. Our results indicate that solvation energy could be the driving force of the process and thus, it can be also related with the cytotoxicity of the drug.

Structure and dynamics of proflavine association around DNA

Proflavine is a small molecule that intercalates into DNA and, thereby, acts as an anticancer agent. Intercalation of proflavine is shown to be a two-step process in which the first step is believed to be the formation of a pre-intercalative outside bound state. Experimental studies so far have been unable to capture the nature of the outside bound state. However, the sub-millisecond timescale observed in fluorescence kinetic experiments is often attributed to the binding of proflavine outside of DNA. Here, we have performed molecular dynamics simulations with multiple proflavine molecules to study the structure and dynamics of the formation of the outside bound state of DNA at different ion concentrations. We observed that the timescale of the outside bound state formation is, at least, five orders of magnitude faster (in nanoseconds) than the experimentally reported timescale (sub-milliseconds) attributed to binding outside DNA. Moreover, we also observed the stacked arrangement of proflavine all around DNA, which is different from the experimentally predicted stacking arrangement perpendicular to the helical axis of DNA in the close vicinity of the phosphate groups. This study, therefore, provides insight into the molecular structure and dynamics of the pre-intercalative outside bound state and will help in understanding the overall intercalation mechanism.

A theoretical study of methylation and CH/π interactions in DNA intercalation: methylated 1,10-phenanthroline in adenine–thymine base pairs

This work shows that quality is better that quantity to estabilize the intercalation of methylated phen.

Molecular Docking of Biologically Active Substances to Double Helical Nucleic Acids

Molecular docking of ligands to DNA-targets is of great importance for the design of new anticancer drugs. Unfortunately, most docking programs were developed for protein-ligand docking which raises a question about their applicability for the DNA-ligand docking. In this study, the popular docking programs AutoDock Vina, AutoDock4 and AutoDock3 were compared for a test set of 50 DNA-ligand complexes taken from the Nucleic Acid Database. It was shown that the version 3.05 of the AutoDock program was the most successful in reproducing the structures of intercalation and minor-groove complexes. The program AutoDock4 was able to re-dock to within 2 Å RMSD most of the intercalation complexes of the test set, but showed poor performance for minor groove binders. While Vina, on the contrary, failed to construct six intercalation complexes of the test set, but showed satisfactory results for DNA-ligand minor-groove complexes when small search space was used.

DNA intercalation studies and antimicrobial activity of Ag@ZrO2 core-shell nanoparticles in vitro

Ag@ZrO2 core-shell nanoparticles were prepared by one pot simultaneous reduction of AgNO3 and hydrolysis of zirconium (IV) isopropoxide. The formation of core-shell nanoparticles was confirmed by absorption, XRD, and HR-TEM techniques. The antibacterial activity of Ag@ZrO2 core-shell nanoparticles against Escherichia coli and Staphylococcus aureus and the antifungal properties against Candida albicans, Candida glabrata, Aspergillus niger and Aspergillus flavus were examined by the agar diffusion method. DNA intercalation studies were carried out in CT-DNA. As a result ZrO2 supported on the surface of AgNPs not only prevented aggregation, but also proved to have enhanced antimicrobial activity and DNA intercalation than the Ag nanoparticles.

Proflavine Hemisulfate as a Fluorescent Contrast Agent for Point-of-Care Cytology

Proflavine hemisulfate, an acridine-derived fluorescent dye, can be used as a rapid stain for cytologic examination of biological specimens. Proflavine fluorescently stains cell nuclei and cytoplasmic structures, owing to its small amphipathic structure and ability to intercalate DNA. In this manuscript, we demonstrated the use of proflavine as a rapid cytologic dye on a number of specimens, including normal exfoliated oral squamous cells, cultured human oral squamous carcinoma cells, and leukocytes derived from whole blood specimens using a custom-built, portable, LED-illuminated fluorescence microscope. No incubation time was needed after suspending cells in 0.01% (w/v) proflavine diluted in saline. Images of proflavine stained oral cells had clearly visible nuclei as well as granular cytoplasm, while stained leukocytes exhibited bright nuclei, and highlighted the multilobar nature of nuclei in neutrophils. We also demonstrated the utility of quantitative analysis of digital images of proflavine stained cells, which can be used to detect significant morphological differences between different cell types. Proflavine stained oral cells have well-defined nuclei and cell membranes which allowed for quantitative analysis of nuclear to cytoplasmic ratios, as well as image texture analysis to extract quantitative image features.

Atomistic account of structural and dynamical changes induced by small binders in the double helix of a short DNA

Nucleic acids are flexible molecules and their dynamical properties play a key role in molecular recognition events. Small binders interacting with DNA fragments induce both structural and dynamical changes in the double helix. We study the dynamics of a DNA dodecamer and of its complexes with Hoechst 33258, which is a minor groove binder, and with the ethidium cation, which is an intercalator, by molecular dynamics simulation. The thermodynamics of DNA-drug interaction is evaluated in connection with the structure and the dynamics of the resulting complexes. We identify and characterize the relevant changes in the configurational distribution of the DNA helix and relate them to the corresponding entropic contributions to the binding free energy. The binder Hoechst locks the breathing motion of the minor groove inducing a reduction of the configurational entropy of the helix, which amounts to 20 kcal mol(-1). In contrast, intercalations with the ethidium cation enhance the flexibility of the double helix. We show that the balance between the energy required to deform the helix for the intercalation and the gain in configurational entropy is the origin of cooperativity in the binding of a second ethidium and of anti-cooperativity in the binding of a third one. The results of our study provide an understanding of the relation between structure, dynamics and energetics in the interaction between DNA fragments and small binders, highlighting the role of dynamical changes and consequent variation of the configurational entropy of the DNA double helix for both types of binders.

Recent advances in small organic molecules as DNA intercalating agents: synthesis, activity, and modeling

The interaction of small molecules with DNA plays an essential role in many biological processes. As DNA is often the target for majority of anticancer and antibiotic drugs, study about the interaction of drug and DNA has a key role in pharmacology. Moreover, understanding the interactions of small molecules with DNA is of prime significance in the rational design of more powerful and selective anticancer agents. Two of the most important and promising targets in cancer chemotherapy include DNA alkylating agents and DNA intercalators. For these last the DNA recognition is a critical step in their anti-tumor action and the intercalation is not only one kind of the interactions in DNA recognition but also a pivotal step of several clinically used anti-tumor drugs such as anthracyclines, acridines and anthraquinones. To push clinical cancer therapy, the discovery of new DNA intercalators has been considered a practical approach and a number of intercalators have been recently reported. The intercalative binding properties of such molecules can also be harnessed as diagnostic probes for DNA structure in addition to DNA-directed therapeutics. Moreover, the problem of intercalation site formation in the undistorted B-DNA of different length and sequence is matter of tremendous importance in molecular modeling studies and, nowadays, three models of DNA intercalation targets have been proposed that account for the binding features of intercalators. Finally, despite DNA being an important target for several drugs, most of the docking programs are validated only for proteins and their ligands. Therefore, a default protocol to identify DNA binding modes which uses a modified canonical DNA as receptor is needed.

Drug-DNA intercalation: from discovery to the molecular mechanism

The ability of small molecules to perturb the natural structure and dynamics of nucleic acids is intriguing and has potential applications in cancer therapeutics. Intercalation is a special binding mode where the planar aromatic moiety of a small molecule is inserted between a pair of base pairs, causing structural changes in the DNA and leading to its functional arrest. Enormous progress has been made to understand the nature of the intercalation process since its idealistic conception five decades ago. However, the biological functions were detected even earlier. In this review, we focus mainly on the acridine and anthracycline types of drugs and provide a brief overview of the development in the field through various experimental methods that led to our present understanding of the subject. Subsequently, we discuss the molecular mechanism of the intercalation process, free-energy landscapes, and kinetics that was revealed recently through detailed and rigorous computational studies.