Molecular Modeling Study of the Norfloxacin-DNA Complex

Abstract

Molecular modeling and molecular dynamics were performed to investigate the interaction of norfloxacin with the DNA oligonucleotide 5′-d(ATACGTAT)₂. Eight quinolone-DNA binding structures were built by molecular modeling on the basis of experimental results. A 100ps molecular dynamics calculation was carried out on two groove binding models and six partially intercalating models. The resulting average structures were compared with each other and to free DNA structure as a reference. The favorable binding mode of norfloxacin to a DNA substrate was pursued by structural assess including steric hindrance, presence of hydrogen-bonding, non-bonding energies of the complex and presence of abnormal structural distortion. Although two of the intercalative models showed the highest binding energy and the lowest non-bonding interaction energy, they presented structural features which contrast with experimental results. On the other hand, one groove binding model demonstrated the most acceptable structure when the experimental observation was accounted. In this model, hydrogen bonding of the carbonyl and carboxyl group of the norfloxacin rings with the DNA bases was present, and norfloxacin binds to the amine group of the guanine base which protrudes toward the minor groove of B-DNA.     

Binding free energy of selected anticancer compounds to DNA--theoretical calculations

Many studies have elucidated structures and thermodynamics of complexes formed by different ligands with DNA. However, in most cases structural and free energy binding studies were not correlated with each other because of the problem of identifying which experimental free energy of binding corresponds to which experimental DNA-ligand structure. In the present work, Poisson-Boltzmann and solvent-accessible surface area methods were used to predict unknown modes of interaction between DNA and three different ligands: mitoxantrone and two pyrimidoacridine derivatives. In parallel, experimental measurements of binding free energy for the studied complexes were performed to compare experimental and calculated values. Our studies showed that the calculated values of free energy are only close to experimental data for some models of interaction between ligands and DNA. Based on this correlation, the most likely models of DNA-ligand complexes were postulated: (i) mitoxantrone and one derivative of pyrimidoacridine, both with two charged side chains, intercalate from the minor groove of DNA and bind with both chains in this groove; (ii) pyrimidoacridine, with only one side chain, very likely does not intercalate into DNA at all. Additionally, the non-electrostatic and electrostatic parts of the calculated binding free energy for the DNA-ligands studied are discussed.     

On the competition between water,sodium ions,and spermine in binding to DNA: a molecular dynamics computer simulation study

The interaction of DNA with the polyamine spermine(4+) (Spm(4+)), sodium ions, and water molecules has been studied using molecular dynamics computer simulations in a system modeling a DNA crystal. The simulation model consisted of three B-DNA decamers in a periodic hexagonal cell, containing 1200 water molecules, 8 Spm(4+), 32 Na(+), and 4 Cl(-) ions. The present paper gives a more detailed account of a recently published report of this system and compares results on this mixed Spm(4+)/Na(+)-cation system with an molecular dynamics simulation carried out for the same DNA decamer under similar conditions with only sodium counterions (Korolev et al., J. Mol. Biol. 308:907). The presence of Spm(4+) makes significant influence on the DNA hydration and on the interaction of the sodium ions with DNA. Spermine pushes water molecules out of the minor groove, whereas Na(+) attracts and organizes water around DNA. The major binding site of the Spm(4+) amino groups and the Na(+) ions is the phosphate group of DNA. The flexible polyamine spermine displays a high presence in the minor groove but does not form long-lived and structurally defined complexes. Sodium ions compete with Spm(4+) for binding to the DNA bases in the minor groove. Sodium ions also have several strong binding sites in the major groove. The ability of water molecules, Spm(4+), and Na(+) to modulate the local structure of the DNA double helix is discussed.     

On the importance of van der Waals interaction in the groove binding of DNA with ligands: restrained molecular dynamics study

Comparison of interaction energy between an oligonucleotide and a DNA-binding ligand in the minor and major groove modes was made by use of restrained molecular dynamics. Distortion in DNA was found for the major groove mode whereas less significant changes for both ligand and DNA were detected for the minor groove binding after molecular dynamics simulation. The conformation of the ligand obtained from the major groove mode resembles that computed with the ligand soaked in water. The van der Waals contact energy was found to be as significant as electrostatic energy and more important for difference in binding energy between these two binding modes. The importance of van der Waals force in groove binding was supported by computations on the complex formed by the repressor peptide fragment from the bacteriophage 434 and its operator oligonucleotide.     

Minor groove deformability of DNA: a molecular dynamics free energy simulation study

The conformational deformability of nucleic acids can influence their function and recognition by proteins. A class of DNA binding proteins including the TATA box binding protein binds to the DNA minor groove, resulting in an opening of the minor groove and DNA bending toward the major groove. Explicit solvent molecular dynamics simulations in combination with the umbrella sampling approach have been performed to investigate the molecular mechanism of DNA minor groove deformations and the indirect energetic contribution to protein binding. As a reaction coordinate, the distance between backbone segments on opposite strands was used. The resulting deformed structures showed close agreement with experimental DNA structures in complex with minor groove-binding proteins. The calculated free energy of minor groove deformation was approximately 4-6 kcal mol(-1) in the case of a central TATATA sequence. A smaller equilibrium minor groove width and more restricted minor groove mobility was found for the central AAATTT and also a significantly ( approximately 2 times) larger free energy change for opening the minor groove. The helical parameter analysis of trajectories indicates that an easier partial unstacking of a central TA versus AT basepair step is a likely reason for the larger groove flexibility of the central TATATA case.     

Methylene blue binding to DNA with alternating AT base sequence: minor groove binding is favored over intercalation

The results presented in this paper on methylene blue (MB) binding to DNA with AT alternating base sequence complement the data obtained in two former modeling studies of MB binding to GC alternating DNA. In the light of the large amount of experimental data for both systems, this theoretical study is focused on a detailed energetic analysis and comparison in order to understand their different behavior. Since experimental high-resolution structures of the complexes are not available, the analysis is based on energy minimized structural models of the complexes in different binding modes. For both sequences, four different intercalation structures and two models for MB binding in the minor and major groove have been proposed. Solvent electrostatic effects were included in the energetic analysis by using electrostatic continuum theory, and the dependence of MB binding on salt concentration was investigated by solving the non-linear Poisson-Boltzmann equation. We find that the relative stability of the different complexes is similar for the two sequences, in agreement with the interpretation of spectroscopic data. Subtle differences, however, are seen in energy decompositions and can be attributed to the change from symmetric 5'-YpR-3' intercalation to minor groove binding with increasing salt concentration, which is experimentally observed for the AT sequence at lower salt concentration than for the GC sequence. According to our results, this difference is due to the significantly lower non-electrostatic energy for the minor groove complex with AT alternating DNA, whereas the slightly lower binding energy to this sequence is caused by a higher deformation energy of DNA. The energetic data are in agreement with the conclusions derived from different spectroscopic studies and can also be structurally interpreted on the basis of the modeled complexes. The simple static modeling technique and the neglect of entropy terms and of non-electrostatic solute-solvent interactions, which are assumed to be nearly constant for the compared complexes of MB with DNA, seem to be justified by the results.     

Methylene Blue Binding to DNA with Alternating AT Base Sequence: Minor Groove Binding is Favored over Intercalation

Abstract

The results presented in this paper on methylene blue (MB) binding to DNA with AT alternating base sequence complement the data obtained in two former modeling studies of MB binding to GC alternating DNA. In the light of the large amount of experimental data for both systems, this theoretical study is focused on a detailed energetic analysis and comparison in order to understand their different behavior. Since experimental high-resolution structures of the complexes are not available, the analysis is based on energy minimized structural models of the complexes in different binding modes. For both sequences, four different intercalation structures and two models for MB binding in the minor and major groove have been proposed. Solvent electrostatic effects were included in the energetic analysis by using electrostatic continuum theory, and the dependence of MB binding on salt concentration was investigated by solving the non-linear Poisson-Boltzmann equation. We find that the relative stability of the different complexes is similar for the two sequences, in agreement with the interpretation of spectroscopic data. Subtle differences, however, are seen in energy decompositions and can be attributed to the change from symmetric 5′-YpR-3′ intercalation to minor groove binding with increasing salt concentration, which is experimentally observed for the AT sequence at lower salt concentration than for the GC sequence. According to our results, this difference is due to the significantly lower non-electrostatic energy for the minor groove complex with AT alternating DNA, whereas the slightly lower binding energy to this sequence is caused by a higher deformation energy of DNA. The energetic data are in agreement with the conclusions derived from different spectroscopic studies and can also be structurally interpreted on the basis of the modeled complexes. The simple static modeling technique and the neglect of entropy terms and of non-electrostatic solute-solvent interactions, which are assumed to be nearly constant for the compared complexes of MB with DNA, seem to be justified by the results.     

Molecular docking and dynamics simulations on the interaction of cationic porphyrin–anthraquinone hybrids with DNA G-quadruplexes

A series of cationic porphyrin–anthraquinone hybrids bearing either pyridine, imidazole, or pyrazole rings at the meso-positions have been investigated for their interaction with DNA G-quadruplexes by employing molecular docking and molecular dynamics simulations. Three types of DNA G-quadruplexes were utilized, which comprise parallel, antiparallel, and mixed hybrid topologies. The porphyrin hybrids have a preference to bind with parallel and mixed hybrid structures compared to the antiparallel structure. This preference arises from the end stacking of porphyrin moiety following G-stem and loop binding of anthraquinone tail, which is not found in the antiparallel due to the presence of diagonal and lateral loops that crowd the G-quartet. The binding to the antiparallel, instead, occurred with poorer affinity through both the loop and wide groove. All sites of porphyrin binding were confirmed by 6 ns molecular dynamics simulation, as well as by the negative value of the total binding free energies that were calculated using the MMPBSA method. Free energy analysis shows that the favorable contribution came from the electrostatic term, which supposedly originated from the interaction of either cationic pyridinium, pyrazole, or imidazole groups and the anionic phosphate backbone, and also from the van der Waals energy, which primarily contributed through end stacking interaction.     

DNA-binding study of anticancer drug cytarabine by spectroscopic and molecular docking techniques

The interaction of anticancer drug cytarabine with calf thymus DNA (CT-DNA) was investigated in vitro under simulated physiological conditions by multispectroscopic techniques and molecular modeling study. The fluorescence spectroscopy and UV absorption spectroscopy indicated drug interacted with CT-DNA in a groove-binding mode, while the binding constant of UV-vis and the number of binding sites were 4.0 ± 0.2 × 10⁴ L mol^?1 and 1.39, respectively. The fluorimetric studies showed that the reaction between the drugs with CT-DNA is exothermic. Circular dichroism spectroscopy was employed to measure the conformational change of DNA in the presence of cytarabine. Furthermore, the drug induces detectable changes in its viscosity for DNA interaction. The molecular modeling results illustrated that cytarabine strongly binds to groove of DNA by relative binding energy of docked structure ?20.61 KJ mol^?1. This combination of multiple spectroscopic techniques and molecular modeling methods can be widely used in the investigation on the interaction of small molecular pollutants and drugs with biomacromolecules for clarifying the molecular mechanism of toxicity or side effect in vivo.     

DNA binding studies of antibiotic drug cephalexin using spectroscopic and molecular docking techniques

The purpose of this study was to explore an accurate characterization of the binding interaction of antibiotic drug cephalexin with calf thymus DNA (CT-DNA) as a relevant biological target by using UV absorption, fluorescence spectroscopy and circular dichroism (CD) in vitro under simulated physiological conditions (pH = 7.4) and also through a molecular modeling study. The results showed that the drug interacts with the DNA helix via a minor groove binding mode. The thermodynamic parameters were calculated and showed that the reaction between the drug and CT-DNA was exothermic. In addition, the drug enforced traceable changes in the viscosity of DNA. The molecular modeling results indicated that cephalexin forcefully binds to the minor groove of DNA with a relative binding energy of ?21.02?kJ mol^?1. The obtained theoretical results were in good agreement with those obtained from experimental studies.     

The effect of dimerization on the interaction of ibuprofen drug with calf thymus DNA: Molecularmodeling and spectroscopic investigation

The interaction between the dimer structure of ibuprofen drug (D-IB) and calf thymus DNA under simulative physiological conditions was investigated with the use of Hoechst 33258 and methylene blue dye as spectral probes by the methods of UV-visible absorption, fluorescence spectroscopy, circular dichroism spectroscopy and molecular modeling study.Using the Job's plot, a single class of binding sites for theD-IB on DNA was put in evidence. The Stern–Volmer analysis of fluorescence quenching data shows the presence of both the static and dynamic quenching mechanisms. The binding constants, K_b were calculated at different temperatures, and the thermodynamic parameters ?G^°, ?H^° and ?S^° were given. The experimental results showed that D-IB molecules could bind with DNA via groove binding mode as evidenced by: I. DNA binding constant from spectrophotometric studies of the interaction of D-IB with DNA is comparable to groove binding drugs. II. Competitive fluorimetric studies with Hoechst 33258 have shown that D-IB exhibits the ability of this complex to displace with DNA-bounded Hoechst, indicating that it binds to DNA in strong competition with Hoechst for the groove binding. III. There is no significantly change in the absorption of the MB-DNA system upon adding the D-IB, indicates that MB molecules are not released from the DNA helix after addition of the D-IB and are indicative of a non-intercalative mode of binding. IV. Small changes in DNA viscosity in the presence of D-IB, indicating weak link to DNA, which is consistent with DNA groove binding. As well as, induced CD spectral changes, and the docking results revealed that groove mechanism is followed by D-IB to bind with DNA.     

Methylene blue binding to DNA with alternating GC base sequence: continuum treatment of salt effects

Methylene blue (MB), an efficient singlet oxygen generating photoactive dye, binds to DNA and allows photosensitized reactions to be used for sequence-specific cleavage of the DNA backbone. Intercalation and groove binding are possible binding modes of the dye, depending on base sequences and environmental conditions. In a recent modeling study of methylene blue binding to a double stranded DNA decamer with an alternating GC sequence, six structural models for intercalation structures and for minor and major groove binding have been obtained. By estimating the binding energies (including electrostatic reaction field contributions of a salt-free aqueous solvent), symmetric intercalation at the 5'-CpG-3' and 5'-GpC-3' steps was found as the predominant binding mode, followed by a slightly weaker binding of the dye in the minor groove. In this study, the stability of the modeled structures has been analysed as a function of salt concentration. The results of finite difference numerical solutions of the non-linear Poisson-Boltzmann equation show that the stabilizing effect of salt is larger for free DNA than for the modeled MB-DNA complexes. Accordingly, the estimated binding energies decrease with increasing ionic strength. A slightly higher stabilization of the groove binding complexes results in comparable binding energies for symmetric intercalation and minor groove binding at high salt concentration. Both results are in qualitative agreement with experimental data.     

1H NMR studies of selective interactions of norfloxacin with double-stranded DNA

The interaction of the antibiotic drug norfloxacin with double-stranded DNA containing interior 5'-CpG-3', 5'-GpC-3', and 5'-GpG-3' steps was studied by 1H NMR. The drug is in fast exchange on the NMR timescale. A highly selective broadening of the imino proton resonances assigned to central CpG steps was observed after addition of drug, indicating an intercalation-like interaction. DNA sequences with central CpG steps also displayed broadening of non-hydrogen-bonded cytosine amino protons in the major groove upon addition of norfloxacin. Furthermore, a sequence-independent selective broadening of the adenine H2 resonance and an upfield shift of the guanine amino proton resonance, both protons located in the minor groove, was observed. Two-dimensional-NOESY spectra showed that no significant structural changes were induced in the DNA by the drug. The results suggest that the planar two-ring system of norfloxacin partially intercalates into CpG steps and that the drug also exhibits non-specific groove binding.     

How dopamine transporter interacts with dopamine: insights from molecular modeling and simulation

By performing homology modeling, molecular docking, and molecular dynamics simulations, we have developed three-dimensional (3D) structural models of both dopamine transporter and dopamine transporter-dopamine complex in the environment of lipid bilayer and solvent water. According to the simulated structure of dopamine transporter-dopamine complex, dopamine was orientated in a hydrophobic pocket at the midpoint of the membrane. The modeled 3D structures provide some detailed structural and mechanistic insights concerning how dopamine transporter (DAT) interacts with dopamine at atomic level, extending our mechanistic understanding of the dopamine reuptake with the help of Na(+) ions. The general features of the modeled 3D structures are consistent with available experimental data. Based on the modeled structures, our calculated binding free energy (DeltaG(bind) = -6.4 kcal/mol) for dopamine binding with DAT is also reasonably close to the experimentally derived DeltaG(bind) value of -7.4 kcal/mol. Finally, a possible dopamine-entry pathway, which involves formation and breaking of the salt bridge between side chains of Arg(85) and Asp(476), is proposed based on the results obtained from the modeling and molecular dynamics simulation. The new structural and mechanistic insights obtained from this computational study are expected to stimulate future, further biochemical and pharmacological studies on the detailed structures and mechanisms of DAT and other homologous transporters.     

Binding site identification of anticancer drug gefitinib to HSA and DNA in the presence of five different probes

This study was carried out to evaluate the binding interaction of gefitinib (GEF) with human serum albumin (HSA) and calf thymus DNA (ct-DNA) using fluorescence, UV–Visible, zeta potential measurements and molecular docking methods in order to understand its pharmacokinetic mechanism. By increasing the temperature, a steady decrease in Stern–Volmer quenching constants was observed for HSA binding properties; this indicates a static type of fluorescence quenching. Negative values were calculated for Gibbs free energy (ΔG), enthalpy (ΔH), and entropy (ΔS) changes, indicating that the reaction is spontaneous and enthalpy-driven. Probe competitive experimental results showed that GEF contains the same binding site as warfarin and are consistent with modeling results. The zeta potential of the HSA increased with increasing GEF, which represents the presence of electrostatic interactions in the system. DNA binding properties were investigated in the presence of three probes. The experimental results showed that by increasing GEF to DNA-AO (acridine-orange) and DNA-MB (methylene-blue) system, the fluorescence intensity and absorbance spectra had no considerable change. Furthermore, with the addition of GEF to DNA, the zeta potential decreased gradually, indicating that the hydrophobic interaction between the GEF and the bases of DNA is the major factor. Thus, GEF can bind to DNA via a groove binding mode. It was also found that GEF entered into the minor groove in the A–T rich region of DNA fragment and bind via van der-Waals forces and three H-bond with double strands of DNA. This is in good agreement with experimental results.     

Binding properties of palmatine to DNA: spectroscopic and molecular modeling investigations

Palmatine, an isoquinoline alkaloid, is an important medicinal herbal extract with diverse pharmacological and biological properties. In this work, spectroscopic and molecular modeling approaches were employed to reveal the interaction between palmatine and DNA isolated from herring sperm. The absorption spectra and iodide quenching results indicated that groove binding was the main binding mode of palmatine to DNA. Fluorescence studies indicated that the binding constant (K) of palmatine and DNA was ~ 10⁴ L·mol^?1. The associated thermodynamic parameters, ΔG, ΔH, and ΔS, indicated that hydrogen bonds and van der Waals forces played major roles in the interaction. The effects of chemical denaturant, thermal denaturation and pH on the interaction were investigated and provided further support for the groove binding mode. In addition to experimental approaches, molecular modeling was conducted to verify binding pattern of palmatine–DNA. Copyright © 2015 John Wiley & Sons, Ltd.     

Control of DNA minor groove width and Fis protein binding by the purine 2-amino group

The width of the DNA minor groove varies with sequence and can be a major determinant of DNA shape recognition by proteins. For example, the minor groove within the center of the Fis–DNA complex narrows to about half the mean minor groove width of canonical B-form DNA to fit onto the protein surface. G/C base pairs within this segment, which is not contacted by the Fis protein, reduce binding affinities up to 2000-fold over A/T-rich sequences. We show here through multiple X-ray structures and binding properties of Fis–DNA complexes containing base analogs that the 2-amino group on guanine is the primary molecular determinant controlling minor groove widths. Molecular dynamics simulations of free-DNA targets with canonical and modified bases further demonstrate that sequence-dependent narrowing of minor groove widths is modulated almost entirely by the presence of purine 2-amino groups. We also provide evidence that protein-mediated phosphate neutralization facilitates minor groove compression and is particularly important for binding to non-optimally shaped DNA duplexes.     

Molecular dynamics simulations and free energy calculations of netropsin and distamycin binding to an AAAAA DNA binding site

Molecular dynamics simulations have been performed on netropsin in two different charge states and on distamycin binding to the minor groove of the DNA duplex d(CGCGAAAAACGCG)·d(CGCGTTTTTCGCG). The relative free energy of binding of the two non-covalently interacting ligands was calculated using the thermodynamic integration method and reflects the experimental result. From 2 ns simulations of the ligands free in solution and when bound to DNA, the mobility and the hydrogen-bonding patterns of the ligands were studied, as well as their hydration. It is shown that even though distamycin is less hydrated than netropsin, the loss of ligand–solvent interactions is very similar for both ligands. The relative mobilities of the ligands in their bound and free forms indicate a larger entropic penalty for distamycin when binding to the minor groove compared with netropsin, partially explaining the lower binding affinity of the distamycin molecule. The detailed structural and energetic insights obtained from the molecular dynamics simulations allow for a better understanding of the factors determining ligand–DNA binding.