Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Binding of Phenazinium Dye Safranin T to Polyriboadenylic Acid: Spectroscopic and Thermodynamic Study

  • Ankur Bikash Pradhan ,

    Contributed equally to this work with: Ankur Bikash Pradhan, Lucy Haque

    Affiliation Department of Chemistry, Jadavpur University, Jadavpur, Kolkata, India

  • Lucy Haque ,

    Contributed equally to this work with: Ankur Bikash Pradhan, Lucy Haque

    Affiliation Department of Chemistry, Jadavpur University, Jadavpur, Kolkata, India

  • Snigdha Roy,

    Affiliation Department of Chemistry, Jadavpur University, Jadavpur, Kolkata, India

  • Suman Das

    sumandas10@yahoo.com

    Affiliation Department of Chemistry, Jadavpur University, Jadavpur, Kolkata, India

Abstract

Here, we report results from experiments designed to explore the association of the phenazinium dye safranin T (ST, 3,7-diamino-2,8-dimethyl-5-phenylphenazinium chloride) with single and double stranded form of polyriboadenylic acid (hereafter poly-A) using several spectroscopic techniques. We demonstrate that the dye binds to single stranded polyriboadenylic acid (hereafter ss poly-A) with high affinity while it does not interact at all with the double stranded (ds) form of the polynucleotide. Fluorescence and absorption spectral studies reveal the molecular aspects of binding of ST to single stranded form of the polynucleotide. This observation is also supported by the circular dichroism study. Thermodynamic data obtained from temperature dependence of binding constant reveals that association is driven by negative enthalpy change and opposed by negative entropy change. Ferrocyanide quenching studies have shown intercalative binding of ST to ss poly-A. Experiments on viscosity measurements confirm the binding mode of the dye to be intercalative. The effect of [Na+] ion concentration on the binding process suggests the role of electrostatic forces in the complexation. Present studies reveal the utility of the dye in probing nucleic acid structure.

Introduction

The use of small molecules to specifically control important cellular functions through binding to nucleic acids is an area of active interest at the interface of chemical biology and medicinal chemistry. The characterization of interaction of small molecules with nucleic acid not only provides insights in biology but also gives the opportunity for developing effective therapeutic agents [1], [2]. Some small molecules have been proven to be useful as sensitive probes for nucleic acid structures [3]. Since structural polymorphism of nucleic acids in vivo presents numerous opportunities for the development of new therapeutic agents, the design of small molecules that bind to noncanonical nucleic acid structures represents an active area of rational drug design [3][5]. The role of RNA in the progression of many diseases particularly in viral infections like HIV, AIDS and hepatitis B has led to growing interest in RNA as a potential target for therapeutic intervention [6][8]. Hence, the study of interaction between ligands with different polymorphic forms of RNA is not only of immense importance to understand the basic fundamentals of RNA-ligand interaction, but it is a necessity for the development of RNA targeted therapeutic agents.

Among the single stranded nucleic acids, polyriboadenylic acid [poly-A] is of particular biological relevance due to its role in mRNA functioning and gene expression [9], [10]. Single stranded poly-A has critical roles in cell biology [11], [12]. So it is important to understand how nucleic acid-targeting ligands affect such structures. Virtually all mRNAs in eukaryotic cells have a poly-A tail at the 3′ end that is an important determinant in the maturation and stability of mRNA and also the initiation of translation as well as the production of alternate forms of proteins [13][15]. Molecules that can interact with the polyadenylate tail may inhibit mRNA function and impair protein production in the cell. This switching off of the protein production by targeting poly-A tail may be an avenue for the development of RNA based therapeutic agents. Since the discovery that Neo-PAP (a recently identified human poly-A polymerase) is significantly over expressed in some human cancer cells [16], [17], it has been suggested that the poly-A tails of mRNA might interfere with the full processing of mRNA by PAP and switch off protein synthesis. This suggests that poly-A tail is a potential tumor-specific target [18]. Further, poly-A has been reported to exist in single stranded helical structure and parallel double stranded helix depending on a narrow pH variation and temperature [18][22]. Such a polymorphic conversion of poly-A makes it a prospective target for investigation to understand the binding and structural aspects through ligand interaction [23][27]. Polyadenylation process also plays a significant regulatory role in the production of alternative forms of proteins [28], [29]. A possible biological role for ds poly-A structure has been proposed by Zarudnaya et al. [30] and other group [31]. They have suggested the involvement of such structure in intracellular process as termination of mRNA–poly-A synthesis and auto regulation of poly-A binding protein synthesis.

Phenazine derivatives are a kind of antibiotic. Studies have shown that some phenazinium dyes have antimalarial potency and selectivity and they also inhibit many bacteria from growing [32]. Safranin T is a water soluble phenazinium dye which features a planar phenazine ring and a positive charge (Figure 1). This compound is mainly used as food dye in flavoring and colouring candies and cookies. It is also used for dyeing tannin, cotton, bast fibers, wool, silk, leather and paper [33]. A very wide range of biological application of ST has been reported [34], [35].

DNA binding aspects of ST has been reported and it has been shown that this ligand binds to DNA by the mechanism of intercalation [36][38]. Binding studies of this compound with RNA are scarce [39]. Due to the relevance of poly-A to mRNA stability, protein synthesis, virology, as well as the potential importance of poly-A to cancer biology, we choose to assess the ability of ST to target poly-A. We investigated in details the binding of ST to poly-A and our studies reveal its preferential binding to ss poly-A structure compared to the ds poly-A. This present study would be a good addition to the literature to understand the fundamentals of small molecule interactions with various RNA structures in order to develop RNA targeted therapeutic agents.

Materials and Methods

Biochemicals

Poly-A was a product from Sigma-Aldrich Corporation (St. Louis, MO, USA). Concentration of poly-A in terms of nucleotide was determined by UV absorbance measurements at 257 nm using molar extinction coefficient values of 10,000 M−1cm−1 [27]. Average molar mass of the poly-A used in our study was 2–2.7×105 Da as determined from viscometric measurements. ST was purchased from Fluka (USA). Purity of ST was verified by thin layer chromatography and melting point determination [40]. Since no impurities were seen, no further purification was done.

Preparation of stock solutions

ST solution was prepared freshly each day; the concentration of the compound was measured using molar extinction coefficients of 29,000 M−1 cm−1 at 520 nm. The compound obeyed Beer's law in the concentration range employed in the study. All the studies with ss poly-A were performed in Citrate-Phosphate (CP) buffer, pH 7.0 containing 5 mM of disodium hydrogen phosphate and 0.75 mM of citric acid. Salt dependent studies were performed in the same buffer containing different concentrations of Na+ ion.

The double stranded form of poly-A was prepared by very slow addition of ss poly-A solution into CP buffer of pH 4.5 under constant stirring and allowing three hour incubation at room temperature for the transition to be completed [27]. The formation of ds poly-A was verified by UV spectroscopy and circular dichroism (CD) study. Deionized and triple distilled water were used throughout our studies. All the buffer solutions were filtered through Millipore membrane filter of 0.45 µm before use.

Absorption spectrophotometric measurements

All the UV-VIS absorption studies were made on a Shimadzu model UV-1800 spectrophotometer (Shimadzu Corporation, Japan) in matched quartz cells of 1 cm path length. A thermo programmer was attached to it to maintain the temperature of the spectrometer by Peltier effect. Spectrophotometric titrations were performed by keeping fixed concentration of compound and varying the concentration of polynucleotide. The change in the absorption at the λmax of the dye was noted at each P/D ratio (RNA nucleotide/dye molar ratio) till saturation was obtained. This spectrophotometric data were then converted to titration curve of absorbance versus P/D ratio. The concentrations of free and bound dye were calculated from the derived points of the graph.

Spectrofluorimetric measurements

Steady state fluorescence measurements were done in Horiba Jobin Yvon Fluoromax-4 spectrofluorimeter which was attached to highly sensitive temperature controller. Measurements were done in fluorescence free quartz cell of 1 cm path length. A fixed concentration of the dye was titrated by increasing concentration of poly-A under constant stirring condition.

Evaluation of binding parameters from spectroscopic results

The results of absorption and fluorescence titration of ST with ss poly-A were expressed in terms of Scatchard plots as r/Cf versus r [41]. All the Scatchard plots were nonlinear and showed negative slopes at low r values as observed in non-cooperative binding isotherms and hence were analyzed by excluded site model for non-linear non-cooperative ligand binding phenomenon using McGhee and von Hippel equation [42]. (1)

Here, K′ is the intrinsic binding constant to an isolated site; n is the number of nucleotides occluded after the binding of one single ligand molecule; r is the number of moles of ligand bound per mole of nucleotide and Cf is the molar concentration of free ligand. The binding data were analyzed using the Origin 7.0 software to determine the best-fit parameters of K′ and n.

Determination of fluorescence polarization anisotropy

The steady state fluorescence anisotropy was also measured using the same spectrofluorimeter. Steady state anisotropy (r′), was defined by [43], (2)

Where, G is the ratio IHV/IHH used for instrumental correction. IVV, IVH, IHV and IHH represent the fluorescence signal for excitation and emission with the polarizer positions set at (0°,0°), (0°,90°), (90°, 0°) and (90°,90°), respectively.

Fluorescence lifetime measurements

Time correlated single photon counting (TCSPC) measurements were carried out in 10 mM CP buffer of pH 7.0 at 20°C for the fluorescence decay of ST in the absence and in the presence of increasing concentration of poly-A. For the TCSPC measurements the photoexcitation was made at 450 nm using a picosecond diode laser (IBH Nanoled-07) in an IBH fluorocube apparatus. The fluorescence decay data were collected on a Hamamatsu MCP photo multiplier (R3809) and were analyzed by using IBH DAS6 software using the equation, (3)where, αi is the ith preexponential factor and τ is the decay time. The decay time is life time of the excited species. The quality of fits was judged from χ2 criterion and visual inspection of the residuals of the fitted function to the data.

Mode of binding: fluorescence quenching studies

Quenching studies were carried out with the anionic quencher potassium ferrocyanide. Solutions of KCl and K4[Fe(CN)6] were mixed in different ratios to give a fixed total ionic strength. Fluorescence quenching experiments were performed at a constant P/D ratio and fluorescence intensity was measured as a function of changing concentration of the ferrocyanide [44]. The data were plotted as Stern-Volmer plot of relative fluorescence intensity (Fo/F) versus [Fe(CN)6]4− using the equation [45], (4)where, Fo and F are the fluorescence intensities in the absence and in the presence of the quencher (Q) K4[Fe(CN)6]; KSV is the Stern Volmer constant.

Circular dichroic titration

Circular dichroism (CD) measurements were carried out on a PC-driven JASCO J815 spectropolarimeter (Jasco International Co.) attached with a temperature controller and a thermal programmer (model PFD-425L/15) interfaced in a rectangular quartz cuvette of 1 cm path length. All CD spectra were recorded in the wavelength range of 200–600 nm with a scan speed of 100 nm min−1. Each spectrum was averaged from five readings for each sample. Final CD spectra were expressed in terms of molar ellipticity ([θ]) in unit of deg cm2 dmol−1 by using the software provided with the spectropolarimeter. The molar ellipticity is based on nucleic acid concentration for the intrinsic CD and dye concentration for the extrinsic CD. All measurements were done at 20°C.

Influence of increasing ionic strength on the perturbation of ss poly-A CD spectrum by ST was investigated at four different salt concentrations, viz. 5, 10, 50 and 100 mM [Na+].

Denaturation and renaturation studies

The denaturation and renaturation profiles of ss poly-A in the absence and in the presence of ST were taken at a rate of 0.5°C min−1 in the temperature range of 15°C to 60°C. Thermal melting profile of ds poly-A was taken in the temperature range 20°C to 110°C at the same rate of heating.

Solution viscosity measurements

Viscometric measurements were carried out using a Cannon-Manning semi micro dilution viscometer type 75 (Cannon Instruments Co., State College, PA, USA) submerged vertically in a constant temperature bath maintained at 20±0.5°C. The molecular weight of the ss poly-A sample was estimated to be in the order 2–2.7×105 Da with an intrinsic viscosity of 2.8 dL/g. 700 µL of RNA solution (500 µM for ss poly-A and 400 µM for ds poly-A) was placed in the viscometer and aliquots of stock solution of ST were directly added into the viscometer to obtain increasing D/P (ST/nucleotide molar ratio) values. Flow times of ss poly-A and ds poly-A in absence and in presence of increasing concentration of ST were measured in triplicate with an accuracy of ±0.01 s and the relative specific viscosity was calculated using the equation: (5)where η'sp and ηsp were the specific viscosity of ss poly-A or ds poly-A in presence and in absence of ST; tcomplex and tcontrol were the time of flow of complex and control solution and t0 is the same for buffer solution as described previously (44).

The relative increase in length L/Lo were obtained from the corresponding increase in relative viscosity by using the following equation [46]. (6)

Here, L and Lo are the contour lengths of RNA in presence and absence of the dye, η and η0 are the corresponding values of intrinsic viscosity (approximated by the reduced viscosity η = ηsp/C, where C is the RNA concentration and β is the slope when L/Lo is plotted against r.

Temperature dependent spectrophotometric study

Temperature dependent absorption spectra were recorded by using Shimadzu UV-1800 of spectrophotometer attached with thermometric cell temperature programmer and temperature controller. These measurements were performed at 20, 25, 30 and 35°C allowing an equilibrium period of 15 minutes for each addition.

Evaluation of thermodynamic parameters

The values of K′ and n were determined at different temperatures. Thermodynamic parameters were estimated by analysis of van't Hoff plot (ln K′ vs. 1/T) obtained over the temperature range of the study. The slope of the plot gives the binding enthalpy change (ΔH0) as (7)

The Gibbs free energy change (ΔG0) was determined from the binding constant at a particular temperature according to this relation (8)

The entropy change (ΔS0) was estimated from the following relation (9)

Salt effect on binding

ST carries a positive charge and poly-A has a negative phosphate skeleton. Therefore in the binding process electrostatic contribution is thought to be one of the driving forces. Cations are condensed around the polyanionic poly-A helix and charged dye compete to expel the cations for phosphate neutralization; these are thermodynamically linked processes. To provide insight into such molecular details absorption and fluorimetric titrations were carried out at four different salt concentrations, viz. 5, 10, 25 and 50 mM [Na+], and the association constants (K′) were evaluated. The following relationship between K′ and [Na+] has been derived previously linking the charge to the variation of binding affinity with [Na+] [47]. (10)

Here, Z is the apparent charge of the bound ligand per phosphate binding. Ψ is the fraction of the [Na+] bound per RNA phosphate and SK is equivalent to the number of counterions released upon binding of a dye. From the salt dependence of binding constant the electrostatic (ΔGpe) and nonelectrostatic (ΔGt) contributions to the overall Gibbs free energy was calculated. The polyelectrolytic contribution to the overall observed free energy can be quantitatively estimated from the relationship, (11)

At a given salt concentration, the non electrostatic contribution can be calculated as the difference between ΔG and ΔGpe.

Results and discussion

Spectral characteristics of poly-A structure

Characteristic circular dichroism (CD) spectra of poly-A at pH 7.0 and 4.5 for single stranded and double stranded structures are presented in Figure 2. Corresponding absorption spectra are shown in the inset. The CD spectrum of the ss poly-A structure at pH 7.0 has a positive band around 264 nm and a strong negative band around 248 nm followed by a positive band 220 nm. When pH was lowered to 4.5 the positive band was enhanced in ellipticity and was blue shifted while the negative band was blue shifted and became weaker with a peak at around 242 nm. The positive band around 220 nm in the ss poly-A structure was disappeared and a negative shoulder appeared around 234 nm. The absorption spectrum (inset of Figure 2) of ss poly-A at pH 7.0 had a maximum at 257 nm and at pH 4.5 this band maximum was blue shifted to 252 nm with a hypochromic effect. All these changes characterized the transition of poly-A from a single stranded helix to a double stranded helix [48].

thumbnail
Figure 2. Circular Dichorism spectral titration of poly-A with ST.

CD spectra of ss poly-A (17.0 µM, curve 1) and ds poly-A (17.0 µM, curve 2) at 20°C in 10 mM CP buffer, pH 7.0 and 4.5 respectively. Inset: UV spectra ss poly-A (17.0 µM, curve 1) and ds poly-A (17.0 µM, curve 2) at 20°C in 10 mM CP buffer, pH 7.0 and 4.5 respectively.

https://doi.org/10.1371/journal.pone.0087992.g002

Binding characteristics of ST–poly-A interaction: Spectrophotometric titration

The effect of increasing concentration of ss poly-A on the absorption spectrum of ST is presented in Figure 3. The interaction of the ST with ss poly-A resulted in marked changes in the visible region (370–650 nm) of the dye absorption spectrum with hypochromic and bathochromic effects that essentially indicated strong intermolecular interaction involving effective overlap of the π electron cloud of ST with the RNA bases. Pronounced hypochromic effect in the dye spectrum could be the indication of intercalative binding of ST between the bases of single stranded RNA. Presence of a sharp isosbestic point at 537 nm clearly indicates the equilibrium between bound and the free dye molecules [44]. Similar type of spectral features on the interaction of ST with calf thymus (CT) DNA have also been reported by Cao and He [36]. In their study two isosbestic points were observed at 432 nm and 538 nm respectively. In the present study on the interaction of ss poly-A we have observed only one isosbestic point.

thumbnail
Figure 3. Absorption spectral change of ST in presence of poly-A.

Representative absorption spectrum of ST (3.03 µM) treated with increasing concentration of ss poly-A in 10 mM CP buffer pH 7.0 at 20°C: Spectra 1–6 denote ss poly-A concentrations of 0, 12.02, 16.40, 26.02, 40 and 72.50 µM. Inset: Scatchard plot for the binding of ST to ss Poly-A. The solid line is the non linear least squares best fit of the experimental points to the Von Hippel equation.

https://doi.org/10.1371/journal.pone.0087992.g003

The absorption titration data were expressed in the form of Scatchard plot (inset of Figure 3) which was observed to be nonlinear and noncooperative in nature. Therefore the Scatchard plot was fitted to a theoretical curve drawn according to the excluded site model developed by McGhee and von Hippel [42] for a nonlinear non-cooperative ligand binding. The best fit values of the binding parameters calculated from Equation 1 are presented in Table 1. The striking result that emerged from this experiment and analysis was that ST binds to ss poly-A in a noncooperative manner with an intrinsic binding affinity (K′) of 5.98×105 M−1 and an exclusion site (n) of 7.4 bases (at 20°C, Table 1). Similar stoichiometry (∼7) was also obtained from continuous variation method (Job's plot) using both spectrofluorimetry and UV-vis spectrophotometry (data not shown). Binding of ST with CT DNA, fish sperm (FS) DNA and yeast RNA has been reported by Huang et al [39]. In our study on the interaction of ST with ss poly-A, binding constant has been found to be of the same order (∼105 M−1) as that of the data reported by them (K′ = 1.25×105, 2.40×105 and 2.41×105 M−1 with CT DNA, FS DNA and yeast RNA respectively). In comparison with the data reported by Huang et al on the binding of ST with CT DNA (n = 4.2) FS DNA (n = 4.3) and yeast RNA (n = 2.9), in our case the binding stoichiometry (n∼7) has been found to be large for ST-ss poly-A interaction. This may be attributed to the differences in the sequences of polymers used and their structures in the solutions. In an another study [36] the binding constant for the association of ST with CT DNA has been reported to less (∼104 M−1). This may be due to the difference in the components, pH and ionic strength of the buffer solution used. In our study no effect of ds poly-A on the absorption spectrum of ST was noticed up to a P/D ratio of 100 (Figure S1, panel A). This confirmed the absence of any binding of ST with ds poly-A.

thumbnail
Table 1. Binding parameters for the interaction of ST with ss poly-A in 10 mM CP buffer, pH 7.0 at 20°C obtained from spectrophotometry and spectrofluorimetrya.

https://doi.org/10.1371/journal.pone.0087992.t001

Fluorescence spectral study

ST is a fluorescent molecule. The emission spectrum of ST was recorded in the wavelength range 530–750 nm with maximum around 580 nm when excited at 520 nm. Complex formation was monitored by titration studies keeping a constant concentration of the dye and increasing the concentration of ss poly-A. With increasing concentration of ss poly-A a progressive enhancement in the fluorescence intensity of the dye was observed and eventually reaching a saturation point with approximately 3 nm blue shift in the wavelength maximum (Figure 4).The results of the spectrofluorimetric titrations were analyzed by constructing Scatchard plot. The Scatchard plot (inset of Figure 4) exhibited non-cooperative behavior as revealed by negative slope and hence were analyzed further by the McGheee and von Hippel methodology [42] for noncooperative binding using Equation 1 for the evaluation of the binding parameters. Analysis yielded a binding constant of 6.15×105 M−1 and n value of 7.60 (at 20°C), which are in excellent agreement with the spectrophotometric results (table 1). The values of binding parameters are presented in the Table 1. Here also no effect of ds poly-A on the fluorescence spectrum of ST up to a P/D ratio of 100 confirmed the absence of binding of the dye to ds poly-A (Figure S1, panel B).

thumbnail
Figure 4. Fluorescence spectra of ST in presence of poly-A.

Fluorescence spectra of ST in presence of ss poly-A in 10 mM CP buffer pH 7.0 at 20°C: Spectra 1–8 denote the fluorescence spectrum of ST (5.0 µM) treated with 0, 15.0, 27.0, 37.5, 50.0, 75.0, 100.0 and 125.0 µM of the ss poly-A. Inset: Scatchard plot for the binding of ST to ss poly-A. The solid line is the non linear least squares best fit of the experimental points to the Von Hippel equation.

https://doi.org/10.1371/journal.pone.0087992.g004

Mode of binding: fluorescence quenching studies

Fluorescence quenching experiments provide an effective method to address the mode of binding of small ligands to nucleic acids [44]. In principle, molecules that are either free or bound on the surface of nucleic acid are easily accessible for the quencher while those are inserted between bases of the polynucleotide may not be accessible to the quencher. Negative charge on the phosphate backbone of the polynucleotide causes an electrostatic barrier at the helix surface and limits the penetration of an anionic quencher into the interior core of the helix. As a result very little or no quenching may be observed in the presence of such quencher if the binding involves strong stacking or intercalation and consequently the magnitude of the Stern-Volmer quenching constant (KSV) of the ligands that are bound inside will be lower than that of the free molecules. It is observed that binding to the poly-A resulted in an increase of the fluorescence intensity of ST (Figure 4). Representative Stern-Volmer plots for free and RNA bound ST are shown in Figure 5. KSV values for free ST and its complex with ss poly-A were 32.5 and 8.0 L mol−1 respectively. This indicated that bound ST was less accessible to the quencher or in other words were considerably protected and sequestered away from the solvent suggesting intercalative binding with ss poly-A. Similar kind of results have also been reported for the binding of methylene blue, ethidium bromide, propidium iodide with ss poly-A [48], [49].

thumbnail
Figure 5. Stern–Volmer plots.

Quenching of ST fluorescence by [Fe(CN)6]4− at 20°C in the absence (•) and in presence (○) of ss poly-A in 10 mM CP buffer, pH 7.0.

https://doi.org/10.1371/journal.pone.0087992.g005

Steady state fluorescence polarization anisotropy study

Fluorescence anisotropy measurements provide effective information about the nature of the environment of probes. Any factor that affects the shape, size and flexibility of a molecule alters the observed anisotropy [45]. Increase in the rigidity of the environment surrounding the fluorophore causes an increase in the fluorescence anisotropy. Change in anisotropy can thus help in finding the probable location of a probe in environments like nucleic acids [50]. Figure 6 shows the variation of fluorescence anisotropy with increasing concentration of ss poly-A for the fluorophore ST. A significant increase in the fluorescence anisotropy on binding with ss poly-A suggests that the dye was trapped in a motionally restricted region within the polynucleotide. It has been found that the fluorescence polarization anisotropy of ST upon binding to the ss poly-A showed a value of 0.18 at saturation against a value of 0.03 for free ST under identical condition. This increase in fluorescence polarization anisotropy provided insight towards strong interaction, probably intercalation, of ST to ss poly-A but further studies like quenching and viscosity measurements were done to illustrate the mode of binding.

thumbnail
Figure 6. Variation of the anisotropy.

Variation of the anisotropy of ST fluorescence as a function of concentration of ss poly-A. λex and λem for ST was 520 nm and 580 nm respectively.

https://doi.org/10.1371/journal.pone.0087992.g006

Lifetime measurements

The fluorescence decay profiles of ST in presence and in absence of ss poly-A are shown in Figure 7. The lifetimes obtained from the best fittings to the decay traces are presented in Table 2. Goodness of the fits has been evaluated from the χ2 criterion (χ2 within 1.0–1.1). It was observed that ST had monoexponential decays in absence as well as in presence of the added polynucleotide. With increase in concentration of ss poly-A lifetime increased gradually up to the saturating condition. Increase in values of lifetime showed strong binding of the compound to the polymer.

thumbnail
Figure 7. Time correlated single photon counting measurements.

Time-resolved fluorescence decay curves (logarithm of normalised intensity versus time in ns) for ST (8.0 µM) in absence (red) and in presence of 40.0 (blue), 80.0 (green), 160.0 (purple) µM of ss poly-A at 20°C in 10 mM CP buffer of pH 7.0. (black) indicates decay curve for the scatterer.

https://doi.org/10.1371/journal.pone.0087992.g007

thumbnail
Table 2. Fluorescence lifetime of ST in absence and in presence of ss poly-A in 10 mM CP buffer, pH 7.0 at 20°C.

https://doi.org/10.1371/journal.pone.0087992.t002

Circular dichroism spectral study

The CD spectrum of ss poly-A was characterized by a large positive band around 264 nm and a negative band around 248 nm. With increasing concentration of ST the CD spectral nature was strongly perturbed (Figure 8). The intensity of 264 nm band with an initial molar ellipticity of around 44,000 deg cm2 dmol−1 decreased when titrated with ST. In presence of the dye the peak was red shifted to around 269 nm and ellipticity became around 12,500 deg cm2 dmol−1. The characteristic negative CD band around 248 nm increased from an initial molar ellipticity of −35,000 deg cm2 dmol−1. Finally at saturation this peak was red shifted to 253 nm and ellipticity went upto −9,000 deg cm2 dmol−1. A clear isoelliptic point at 254 nm was observed for the ss poly-A-ligand interaction. ST molecule itself does not show any characteristic CD in its absorption range due to its molecular symmetry. Figure 8 shows that there is appearance of induced CD band in the absorption region of ST. This is due to the molecular asymmetry induced by the RNA polynucleotide on ST.

thumbnail
Figure 8. Circular dichroism spectral titration of poly-A with ST.

Circular dichroic spectra of ss Poly-A (100.0 µM) with varying concentrations of ST in 10 mM CP buffer pH 7.0 at 20°C. The curves 1–7 represent the ST concentrations of 0, 4.28, 10.68, 21.32, 31.91, 42.44 and 52.92 µM respectively. Inset represents the induced CD spectra of (21.0 µM) of ST with increasing concentrations of ss poly-A. Curves 1–7 represent ss poly-A concentrations of 0, 38.3, 76.6, 152.8, 228.6, 305.2 and 380.2 µM respectively.

https://doi.org/10.1371/journal.pone.0087992.g008

Inset of Figure 8 represents the induced CD of ST in presence of added single stranded polynucleotide. Here, to a fixed concentration of ST, ss poly-A was added until saturation was reached. In presence of added single stranded form of the polymer there was formation of induced CD band with enhanced ellipticity. These induced CD spectra occurred in a wavelength region where the RNA had no absorption band and thus exclusively monitored changes in the environment of the dye molecules. It is pertinent to note here that ST being achiral and planar is not CD active by itself. Therefore the induced CD spectrum was generated exclusively from the asymmetric arrangement of the dye that had been intercalated into the RNA bases. The presence of isoelliptic point at 466 nm in the series of spectra suggested that the bound and free molecules were in equilibrium as revealed from absorption spectral studies also. Appearance of such strong CD band in presence of the polynucleotide confirmed the intercalative mode of binding of ST with ss poly-A. Similar observation has been reported for ST binding with CT DNA [38]. No change in the CD spectrum of ds poly-A in presence of ST clearly indicated the absence of any interaction of ST with double stranded form of the polynucleotide (Figure S1, panel C). This supports our earlier observation in spectrophotometric and spectrofluorimetric data.

CD melting study

The effect of temperature on ss poly-A structure in absence and in presence of the dye was monitored using CD spectroscopy (Figure S2). It was observed that ss poly-A went denaturation (structural transition from helix to random coil) with rise of temperature from 20 to 60°C and a complete renaturation occurred with lowering in temperature to 20°C. This clearly indicated a complete transition of the helical structure of ss poly-A to random coil structure and back. No change in molar ellipticity value of the bound form was observed on heating or cooling at D/P∼0.53 in the same temperature range as mentioned earlier. This indicates that the dye binds to ss poly-A. In our experimental conditions there was absence of any self-structure formation of poly-A in presence of ST. Similar kind of observation has also been reported for ethidium bromide and propidium iodide [48]. Thermal melting profiles of ds poly-A in absence and in presence of the dye are shown in panel C of Figure S2. Double helical to single helical transition was observed around 74°C. No change in the profile in presence of the dye again confirmed the absence of any binding of ST with the double stranded form of the polynucleotide.

Viscometric study

The mode of binding by which ST interacts with poly-A was investigated by viscometric techniques. Hydrodynamic measurements are very sensitive to length changes and stiffening in rod like nucleic acid by insertion of ligands between adjacent bases [51] and are regarded as the most decisive and reliable test for elucidating the binding mode of small molecules to nucleic acids in solution. Ligand insertion between the bases is associated with the lengthening of the nucleic acid and this increase in length manifests itself as an increase in the viscosity of the nucleic acid solution. Thus ligand induced enhancement in the viscosity in the solution of nucleic acid is consistent with an intercalative mode of binding. On the other hand, the lack of a ligand induced increase in viscosity is consistent with a non-intercalative mode of binding (e.g. groove binding). In Figure 9, a plot of (η/η0)1/3 versus r for the complexation of ST to ss and ds poly-A is presented. Increase in viscosity of ss poly-A solution was observed in presence of ST while no change was noted in case of solution of ds poly-A. The value of the slope obtained for ss poly-A-ST interaction was ∼1.4 and this in turn lead to a helix length extension (ΔL) of 0.24 nm at the binding site calculated from a standard value of β = 2 corresponding to a length enhancement of 0.34 nm. This result clearly established intercalation type of binding of ST to ss poly-A helical structure. It is to be noted that, since ss poly-A does not have any base pairing and has only stacked helical structure, a true intercalation model where planar ligand molecules are fully sandwiched between hydrogen-bonded base pairs of double stranded DNA cannot be ideated.

thumbnail
Figure 9. Viscometric measurements.

Plot of change of relative specific viscosity of ss poly-A (blue) and ds poly-A (red) with increasing concentration of ST in 10 mM CP buffer pH 7.0 or 4.5 at 20°C. The concentrations of ss poly-A and ds poly-A were 500 and 400 µM, respectively.

https://doi.org/10.1371/journal.pone.0087992.g009

Thermodynamics of the interaction

Thermodynamic parameters for the association of ST with ss poly-A were calculated from the temperature dependence of binding constants using UV-Vis absorption spectroscopic studies at four different temperatures namely 20, 25, 30 and 35°C respectively at constant salt molarity of 10 mM [Na+].The binding parameters are presented in Table 3. It can be seen that with increase in temperature the binding constant decreased while the number of nucleotide occluded sites (n) changed marginally. The van't Hoff plot for binding is shown in Figure 10. Linear fit of the data indicates a very small value of heat capacity change (ΔCP≈0). The values of the thermodynamic parameters are presented in Table 3. It can be seen that the binding of ST with ss poly-A was characterized by both negative enthalpy and entropy changes.

thumbnail
Figure 10. Temperature dependence of binding constant.

Van't Hoff plot for the complexation of ST with ss poly-A in 10 mM CP buffer of pH 7.0. The data were fitted to straight line with correlation coefficient of 0.97.

https://doi.org/10.1371/journal.pone.0087992.g010

thumbnail
Table 3. Binding and thermodynamic parameters for the interaction of ST with ss poly-A in 10 mM CP buffer, pH 7.0 obtained from absorption spectrophotometrya.

https://doi.org/10.1371/journal.pone.0087992.t003

Conceptually the thermodynamic parameters for the binding process may be divided into three contributions: [i] contributions from hydrogen bonding and hydrophobic part due to interactions between the bound dye and nucleic acid binding site; [ii] contribution from the conformational changes upon binding in either the nucleic acid or the ligand and [iii] contributions from coupled processes like ion release, proton transfer or changes in the water of hydration. Our data revealed that the binding process was favoured by negative enthalpy change and was opposed by negative entropy change. The possible contribution of negative enthalpy change for the binding process of ST to ss poly-A may be explained in terms of the van der Waal's stacking interactions, hydrophobic as well as weak electrostatic interactions. Negative entropy change for the interaction process may be attributed to the increase in order due to the formation of rigid polynucleotide-dye complex.

Dependence of the binding on ionic strength

It has been observed that aromatic systems comprised of two condensed rings can enforce destacking of the nucleobases leading to intercalation only if their association to nucleic acid is facilitated by positive charges in the ring system. RNA condenses counter ions on its surfaces to screen the polyionic charge on the phosphodiester backbone. The phenazinium dye ST is a cationic molecule that may compete with sodium ions for binding to RNA phosphates. To assess the salt dependence of binding of ST with ss poly-A we have done binding studies using spectrophotometry at varying concentration of Na+. Variation of K′ versus [Na+] is plotted in Figure 11. Data on salt dependence of binding are presented in Table 4. The binding constants were found to decrease with increase of salt concentration. Here the plot was found to be linear with slope (−ZΨ) being −0.73. The physical meaning of this value is that the thermodynamic extent of counterion released from each phosphate on binding of single ST molecule is 0.73. Similar data have been reported for other cationic molecules binding to double stranded DNAs and RNAs [48], [52]. The salt dependence of binding clearly emphasized the role of the electrostatic component in the binding process. It has been observed that as the salt concentration was increased, the ΔGpe contribution decreased while nonelectrolytic contribution remained almost unchanged. Data presented in Table 4 clearly revealed that there was a remarkably large contribution from the nonelectrostatic forces to the overall binding free energy. This clearly suggested the larger role of hydrophobic forces in the binding process of ST to ss poly-A. This is in good agreement with the report on cationic molecule binding to double stranded DNAs and RNAs [48], [52][54].

thumbnail
Figure 11. Salt dependence of binding constant.

ln–ln plot of intrinsic binding constant (K′) as a function of ionic strength [Na+] for ST binding to ss poly-A in CP buffer of pH 7.0 at 20°C.

https://doi.org/10.1371/journal.pone.0087992.g011

thumbnail
Table 4. Ionic strength dependence of binding constant for the interaction of ST with ss poly-A obtained from absorption spectrophotometrya.

https://doi.org/10.1371/journal.pone.0087992.t004

Salt dependent CD studies

To understand the role of electrostatic interaction in the binding process, salt dependence of the binding of ST with ss poly-A was performed by CD experiment. The CD spectral changes of the interaction were followed at four Na+ concentrations viz. 05, 10, 20 and 50 mM to complement the observed decrease in binding as the salt increased. The data is presented in Figure S3. It can be seen that the conformational changes in poly-A were less pronounced as the salt concentration enhanced. This was also complemented by the lower intensity for the induced CD bands of ST in the complex. This result confirmed that the binding was favoured at lower salt concentration probably due to the fact that below 50 mM of [Na+] no self-assembled structure was formed in poly-A in presence of ST.

Conclusion

The phenazinium dye has been shown to interact strongly with B-form DNA by the mechanism of intercalation [36][38]. A few studies have been reported on the interaction of phenazinium dye with RNA structures [39]. The strong interaction of ST with ss poly-A was evident from the bathochromic and hypochromic effect in the absorption spectra; significant enhancement of fluorescence intensity, fluorescence polarization anisotropy and life lime; remarkable perturbation of CD spectra; enhancement of viscosity of the polynucleotide solution and magnitude of thermodynamic parameters. From the data presented here it can be concluded that the compound binds to the single stranded form of the polynucleotide by the mechanism of intercalation while it does not bind at all with the double stranded form polyriboadenylic acid. CD melting study shows ST could not induce self structure formation in poly-A. In the binding process there was a remarkably large magnitude from the nonelectrostatic forces to the binding free energy which clearly suggested the larger role of hydrophobic forces in the binding process. The specific ss poly-A binding features of ST would turn into further opportunity for RNA based therapeutic agents.

Supporting Information

Figure S1.

Characterization of interaction of ST with double stranded poly-A. [A].Absorption spectrum of ST (3.05 µM) in absence (curve 1, black) and in presence of 301.09 µM ds poly-A (curve 2, red) in 10 mM CP buffer, pH 4.5 at 20°C. [B] Fluorescence spectrum of ST (5.0 µM) in absence (curve 1, black) and in presence of 510.50 µM ds poly-A (curve 2, red) in 10 mM CP buffer, pH 4.5 at 20°C. [C] Circular dichroic spectra of ds poly-A (100.0 µM) in absence (curve 1, black) and in presence of 55.60 µM ST (curve 2, red) in 10 mM CP buffer, pH 4.5 at 20°C.

https://doi.org/10.1371/journal.pone.0087992.s001

(TIF)

Figure S2.

Denaturation and renaturation study. Spectropolarimetric measurements on heating (red) and cooling (blue) of 100.0 µM of ss poly-A in absence (A) and in presence of ST (B, D/P = 0.53) in 10 mM CP buffer, pH 7.0. The arrows indicate the direction of heating (red) and cooling (blue) process. (C) Melting profile of 100.0 µM ds poly-A in absence (red) and in presence of ST (blue, D/P = 0.53) in10 mM CP buffer, pH 4.5.

https://doi.org/10.1371/journal.pone.0087992.s002

(TIF)

Figure S3.

Salt dependence of CD spectral Changes. Representative CD spectra resulting from the interaction of ss poly-A (100.0 µM) in (A) 5 mM [Na+] treated with 0, 4.28, 10.68, 21.32, 31.91, 42.44, 52.92 and 63.32 µM of ST (curves 1–8); (B) 10 mM [Na+] treated with 0, 4.28, 10.68, 21.32, 31.91, 42.44 and 52.92 µM of ST (curves 1–7); (C) 50 mM [Na+] treated with 0, 10.68, 31.91, 52.92 and 63.32 µM of ST (curves 1–5); and (D) 100 mM [Na+] treated with 0, 10.68, 31.91, 52.92 and 63.32 µM of ST (curves 1–5) in CP buffer, pH 7.0 at 20°C.

https://doi.org/10.1371/journal.pone.0087992.s003

(TIF)

Acknowledgments

The authors thank Professor Nitin Chattopadhyay, Department of Chemistry, Jadavpur University for providing laboratory facilities.

Author Contributions

Conceived and designed the experiments: SD ABP LH. Performed the experiments: ABP LH SR. Analyzed the data: SD ABP LH SR. Contributed reagents/materials/analysis tools: SD. Wrote the paper: SD ABP LH.

References

  1. 1. Dervan PB (2001) Molecular Recognition of DNA by Small Molecules. Biorg Med Chem 9: 2215–2235.
  2. 2. Maiti M, Kumar GS (2007) Molecular aspects on the interaction of protoberberine, benzophenanthridine, and aristolochia group of alkaloids with nucleic acid structures and biological perspectives. Med Res Rev 27: 649–695.
  3. 3. Ihmels H, Otto D (2005) Topics in Current Chemistry: Supramolecular Dye Chemistry. Berlin: Springer. pp 161–204.
  4. 4. Chaires JB (2005) DNA Binders and Related Subjects. Berlin: Springer. pp 33–53.
  5. 5. Palchaudhuri R, Hergenrother PJ (2007) DNA as a target for anticancer compounds: Methods to determine the mode of binding and the mechanism of action. Curr Opin Biotechnol 18: 497–503.
  6. 6. Foloppe N, Matassova N, Aboul-ela F (2006) Towards the discovery of drug-like RNA ligands? Drug Discov Today 11: 1019–1027.
  7. 7. Liu Z, Sall A, Yang D (2008) MicroRNA: An emerging therapeutic target and intervention tool. Int J Mol Sci 9: 978–999.
  8. 8. Fulle S, Gohlke H (2010) Molecular recognition of RNA: challenges for modelling interactions and plasticity. J Mol Recognit 23: 220–231.
  9. 9. Wickens M, Anderson P, Jackson RJ (1997) Life and death in the cytoplasm: messages from the 3′-end. Curr Opin Genet Dev 7: 220–232.
  10. 10. Dower K, Kuperwasser N, Merrikh H, Rosbash M (2004) A synthetic A tail rescues yeast nuclear accumulation of a ribozyme-terminated transcript. RNA 10: 1888–1899.
  11. 11. Svitkin YV, Sonenberg N (2004) An efficient system for cap- and poly(A)-dependent translation in vitro. Meth Mol Biol 257: 155–170.
  12. 12. Alvarez E, Castello A, Menendez-Arias L, Carrasco L (2006) HIV protease cleaves poly(A)-binding protein. Biochem J 396: 219–226.
  13. 13. Munroe D, Jacobson A (1990) mRNA poly(A) tail, a 3′ enhancer of translation initiation. Mol Cell Biol 10: 3441–3455.
  14. 14. Wang Z, Day N, Trifillis P, Kiledjian M (1999) An mRNA stability complex functions with poly(A)-binding protein to stabilize mRNA in vitro. Mol Cell Biol 19: 4552–4560.
  15. 15. Tian B, Hu J, Zhang H, Lutz CS (2005) A large-scale analysis of mRNA polyadenylation of human and mouse genes. Nucleic Acids Res 33: 201–212.
  16. 16. Topalian SL, Kaneko S, Gonzales MI, Bond GL, Ward Y, et al. (2001) Identification and functional characterization of neo-poly(A) polymerase. Mol Cell Biol 21: 5614–5623.
  17. 17. Topalian SL, Gonzales MI, Ward Y, Wang X, Wang RF (2002) Revelation of a cryptic major histocompatibility complex class II-restricted tumor epitope in a novel RNA processing enzyme. Cancer Res 62: 5505–5509.
  18. 18. Saenger W (1984) Principles of nucleic acid structure. New York: Springer-Verlag.
  19. 19. Fresco JR, Doty P (1957) Polynucleotides. I. Molecular properties and configurations of polyriboadenylic acid in solution. J Am Chem Soc 79: 3928–3929.
  20. 20. Rich A, Davies DR, Crick FHC, Watson JD (1961) The molecular structure of polyadenylic acid. J Mol Biol 3: 71–86.
  21. 21. Holcomb DN, Tinoco, Jr I (1965) Conformation of polyriboadenylic acid: pH and temperature dependence. Biopolymers 3: 121–133.
  22. 22. Petrovic AG, Polavarapu PL (2005) Structural transitions in polyriboadenylic acid induced by the changes in pH and temperature: vibrational circular dichroism study in solution and film states. J Phys Chem B 109: 23698–23705.
  23. 23. Yadav RC, Kumar GS, Bhadra K, Giri P, Sinha R, et al. (2005) Berberine, a strong polyriboadenylic acid binding plant alkaloid: spectroscopic, viscometric and thermodynamic study. Bioorg Med Chem 13: 165–174.
  24. 24. Xing F, Song G, Ren J, Chaires JB, Qu X (2005) Molecular recognition of nucleic acids: coralyne binds strongly to poly(A). FEBS Lett 579: 5035–5039.
  25. 25. Giri P, Kumar GS (2007) Specific binding and self-structure induction to poly(A) by the cytotoxic plant alkaloid sanguinarine. Biochimica et Biophysica Acta 1770: 1419–1426.
  26. 26. Cetinkol ÖP, Hud NV (2009) Molecular recognition of poly(A) by small ligands: an alternative method of analysis reveals nanomolar, cooperative and shape-selective binding. Nucleic Acids Research 37: 611–621.
  27. 27. Giri P, Kumar GS (2008) Self-structure induction in single stranded poly(A) by small molecules: studies on DNA intercalators, partial intercalators and groove binding molecules. Arch Biochem Biophys 474: 183–192.
  28. 28. Alt FW, Bothwell ALM, Knapp M, Siden E, Mather E, et al. (1980) Synthesis of secreted and membrane-bound immunoglobulin mu heavy chains is directed by mRNA's that differ at their 3′ ends. Cell 20: 293–301.
  29. 29. McDevitt MA, Hart RP, Wong WW, Nevins JR (1986) Sequences capable of restoring poly(A) site function define two distinct downstream elements. EMBO J 5: 2907–2913.
  30. 30. Zarudnaya MI, Hovorun DM (1999) Hypoterical double-helical poly(A) formation in a cell and its possible biological significance. IUBMB Life 48: 581–584.
  31. 31. Xi H, Gray D, Kumar S, Arya DP (2009) Molecular Recognition of Single-Stranded RNA: Neomycin Binding to Poly(A). FEBS Lett 583: 2269–2275.
  32. 32. Vennerstrom JL, Makler MT (1995) Antimalarial dyes revisited: xanthenes, azines, oxazines, and thiazines. Antimicrob Agents Chemother 39: 2671–2677.
  33. 33. Zaghbani N, Hafiane A, Dhahbi M (2008) Removal of Safranin T from wastewater using micellar enhanced UF. Desalination 222: 348–356.
  34. 34. Krings M (2000) Technical note: The use of biological stains in the analysis of Late Palaeozoic pteridosperm cuticles. Rev Palaeobot Palynol 108: 143–150.
  35. 35. Penney DP, Powers JM, Franck M, Willis C, Churukian C (2002) Analysis and testing of biological stains– The Biological Stain Commission Procedures. Biotech Histochem 77: 237–275.
  36. 36. Cao Y, He XW (1998) Studies of interaction between safranine T and double helix DNA by spectral methods. Spectrochim Acta, Part A 54: 883–892.
  37. 37. Saha I, Kumar GS (2011) Spectroscopic characterization of the interaction of Phenosafranin and Safranin O with double stranded, heat denatured and single stranded Calf Thymus DNA. J Fluoresc 21: 247–255.
  38. 38. Sarkar D, Das P, Basak S, Chattopadhyay N (2008) Binding interaction of cationic phenazinium dyes with calf thymus DNA:a comparative study. J Phys Chem B 112: 9243–9249.
  39. 39. Huang CZ, Li YF, Liu XD (1998) Determination of nucleic acids at nanogram levels with safranine T by a resonance light-scattering technique. Analytica Chimica Acta 375: 89–97.
  40. 40. Proevska LI, Pojarlieff IG (1998) 1H NMR spectra and structure of safranines. Hindered rotation of the 3-dialkylamino group in 7-azo derivatives. Dyes and Pigments 36: 177–190.
  41. 41. Scatchard G (1949) The attraction of proteins of small molecules and ions. Ann. N Y Acad Sci 51: 660–762.
  42. 42. McGhee JD, von Hippel PH (1974) Theoretical aspects of DNA-protein interactions: Co-operative and non-co-operative binding of large ligands to a one-dimensional homogeneous lattice. J Mol Biol 86: 469–489.
  43. 43. Bose D, Ghosh D, Das P, Girigoswami A, Sarkar D, et al. (2010) Binding of a cationic phenazinium dye in anionic liposomal membrane: A spectacular modification in the photophysics. Chem Phys Lipids 163: 94–101.
  44. 44. Das S, Kumar GS (2008) Molecular aspects on the interaction of phenosafranine to deoxyribonucleic acid: model for intercalative drug–DNA binding. J Mol Struct 872: 56–63.
  45. 45. Lakowicz JR (1983) Principles of Fluorescence Spectroscopy. New York: Plenum Press. pp 257–295.
  46. 46. Müller W, Crothers DM (1968) Studies of the binding of actinomycin and related compounds to DNA. J Mol Biol 35: 251–290.
  47. 47. Record Jr MT, Anderson CF, Lohman TM (1978) Thermodynamic analysis of ion effects on the binding and conformational equilibria of proteins and nucleic acids: the roles of ion association or release, screening, and ion effects on water activity. Q Rev Biophys 11: 103–178.
  48. 48. Das S, Parveen S, Pradhan AB (2014) An insight into the interaction of phenanthridine dyes with polyriboadenylic acid: Spectroscopic and thermodynamic approach. Spectrochimica Acta Part A 118: 356–366.
  49. 49. Hossain M, Kabir A, Kumar GS (2012) Binding of the phenothiazinium dye methylene blue with single stranded polyriboadenylic acid. Dyes and Pigments 92: 1376–1383.
  50. 50. Das P, Chakrabarty A, Haldar B, Mallick A, Chattopadhyay N (2007) Effect of Cyclodextrin Nanocavity Confinement on the Photophysics of a β-Carboline Analogue: A Spectroscopic Study. J Phys Chem B 111: 7401–7408.
  51. 51. Lerman LS (1961) Structural considerations in the interaction of DNA and acridines. Journal of Molecular Biology 3: 18–30.
  52. 52. Wilson WD, Lopp IG (1979) Analysis of cooperativity and ion effects in the interaction of quinacrine with DNA. Biopolymers 18: 3025–3041.
  53. 53. Nandi R, Chaudhuri K, Maiti M (1985) Effects of ionic strength and pH on the binding of sanguinarine to deoxyribonucleic acid. Photochem Photobiol 42: 497–503.
  54. 54. Kumar GS, Debnath D, Sen A, Maiti M (1993) Thermodynamics of the interaction of berberine with DNA. Biochem Pharmacol 46: 1665–1667.