Biophysical Characterization and Molecular Docking Studies of


Biophysical Characterization and Molecular Docking Studies of...

0 downloads 103 Views 1MB Size

Subscriber access provided by NEW YORK UNIV

Article

Biophysical Characterization and Molecular Docking Studies of Imidazolium based Polyelectrolytes–DNA complexes: Role of Hydrophobicity Manojkumar Kasina, Prabhu Charan K.T, Sivaramakrishna Akella, Prakash C. Jha, Vijay M Khedkar, Siva Ramamoorthy, Jayaraman Gurunathan, and Vijayakrishna Kari Biomacromolecules, Just Accepted Manuscript • DOI: 10.1021/bm5018029 • Publication Date (Web): 11 Feb 2015 Downloaded from http://pubs.acs.org on February 16, 2015

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Biomacromolecules is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

Biophysical Characterization and Molecular Docking Studies of Imidazolium based Polyelectrolytes–DNA complexes: Role of Hydrophobicity Kasina Manojkumar1, K. T. Prabhu Charan1, Akella Sivaramakrishna1, Prakash C. Jha2, Vijay M. Khedkar3, Ramamoorthy Siva4, Gurunathan Jayaraman,5 and Kari Vijayakrishna*,1 1

Organic Chemistry Division, School of Advanced Sciences, VIT University, Vellore-632014,

Tamil Nadu, India. 2

School of Chemical Sciences, Central University of Gujarat, Sector-30, Gandhinagar-38200,

Gujarat, India. 3

Combi-Chem Resource Centre, CSIR-National Chemical Laboratory, Pune, 411008, India.

4

Plant Biotechnology Division, School of Biosciences and Technology, VIT University, Vellore-

632014, Tamil Nadu, India.

ACS Paragon Plus Environment

1

Biomacromolecules

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

5

Page 2 of 37

Bioinformatics Division, School of Biosciences and Technology, VIT University, Vellore-

632014, Tamil Nadu, India. KEYWORDS Polyelectrolytes; Poly(ionic liquids); DNA; binding interactions; electrostatic interactions; hydrophobic interactions; molecular docking; fluorescence quenching

ABSTRACT

Non-viral gene delivery vectors are acquiring greater attention in the field of gene therapy by replacing the biological viral vectors. DNA-cationic polymer complexes are one of the most promising systems to find application in gene therapy. Hence, a complete insight of their biophysical characterization and binding energy profile is important in understanding the mechanism involved in non viral gene therapy. In this investigation, the interaction between calf thymus DNA (ctDNA) and imidazolium-based poly(ionic liquids) (PILs) also known as polyelectrolytes with three different alkyl side chains (ethyl, butyl and hexyl) in physiological conditions using various spectroscopic experiments with constant DNA concentration and varying polyelectrolyte concentrations is reported. UV-Visible absorption, fluorescence quenching studies, gel electrophoresis, circular dichroism (CD) and Fourier transform infrared spectroscopy (FTIR) have confirmed the binding of polyelectrolytes with DNA. UV-Vis absorption measurements and fluorescence quenching revealed that the binding between DNA and polyelectrolyte is dominated by electrostatic interactions. Additionally, CD and FTIR results indicated that the DNA retained its B-form with minor perturbation in the phosphate backbone without significant change in the conformation of its base pairs. Preference for alkyl side chains

ACS Paragon Plus Environment

2

Page 3 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

(KPIL-Ethyl Br < KPIL -Butyl Br < KPIL -Hexyl Br) towards efficient binding between polyelectrolyte and DNA was inferred from the binding and quenching constants calculated from the absorption and emission spectra, respectively. Further, in silico molecular docking studies not only validated the observed binding trend but also provided insight into the binding mode of the polyelectrolyteDNA complex.

INTRODUCTION Gene therapy, a pivotal method involving the transfer of genetic material into target cells in order to rectify the expression of a particular gene; has attracted significant interest throughout the research community since the past few decades for curing genetic diseases from inherited disorders to cancers.1, 2 The development of safe, efficient and biocompatible gene delivery vectors has become a crucial challenge in the field of gene therapy and clinical applications. Viral vectors are found to be potential gene transfer vehicles but acquired severe drawbacks causing immune reactions, mutations and cancer in patients.3 These safety concerns have encouraged the researchers towards the development of non viral vectors possessing minimal mutagenesis.4 These synthetic vectors which mostly include cationic polymers and cationic lipids have gained much popularity over their counter parts and emerged as a recent topic of interest in the field of gene therapy.5-9 The DNA-cationic polymer complexes involving electrostatic interactions between the cationic polymeric chains and negatively charged phosphate backbone of DNA are being developed as promising non viral gene delivery systems due to their salient features such as aqueous solubility, charge neutralization, pH sensitivity and tunable physiochemical properties.10-14 Several cationic polymers such as linear poly(amidoamine)s (PAA),15 poly[2-(dimethylamino)ethyl methacrylate],16, 17 poly(3-guanidinopropyl

ACS Paragon Plus Environment

3

Biomacromolecules

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 37

methacrylate),18 ester-functionalized linear poly-(ethylenimine),19 polyhistidine,20 poly(ethylenimine) and poly(L-lysine),21 and triazole based click polymers22 have been employed so far to study their interaction with DNA. Their transfection efficiency and cell viability were also evaluated and comprehended with their DNA binding characteristics. On the other hand, PILs containing a polymerized vinylic backbone and an ionic moiety quarternized at the heteroatom with hydrophobic alkyl groups possess unique properties of ionic liquids (ILs) as well as polymers.23-25 The potential and innovative applications of polyelectrolytes include, biosensors, catalytic membranes, absorbing membranes, support for catalysts, polymeric surfactants, electrochemical devices, solid phase microextraction and phase transfer medium.23, 26-28 These polyelectrolytes are advancing as an interesting class of cationic polymers having significant binding affinity to DNA. So far in the literature, the interaction of phosphonium,29 sulfonium,30 ammonium,31, 32 and pyridinium33, 34 based polyelectrolytes with DNA has been reported and their binding mode and transfection efficiency has been discussed. To date, very little data are available in literature regarding the interaction between imidazoliumbased polyelectrolytes and DNA. Recently, Long et al. reported the DNA binding studies of imidazolium-based copolymer quaternized with hydroxyl-functionalized alkyl chains.35 The authors have discussed the contribution of charge density and hydroxyl concentration on DNA binding, cytotoxicity, and in vitro transfection efficiency. It was found that DNA binding increased with the increase in charge density and hydroxy levels, whereas, a systematic change in both the parameters had a significant impact on transfection efficiency. However, a detailed biophysical characterization, structure activity relationship and role of alkyl side chain’s hydrophobicity on the interactions between imidazolium-based polyelectrolytes and DNA are still unclear. We hypothesize that the tuning of polyelectrolyte’s

ACS Paragon Plus Environment

4

Page 5 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

hydrophobicity may be of help in studying the role of hydrophobic interaction between the polyelectrolyte’s alkyl side chains and DNA, apart from the electrostatic interaction between cationic polymer moieties and phosphate anions of DNA. Moreover, there are several reports36-39 in the literature which has provided the evidence for the important role of cationic polymer’s hydrophobicity in gene delivery where it was found to enhance the DNA condensation,40, 41 promote the polymer-DNA complex interaction with cell membranes40 and facilitate the release of nucleic acid from the polymeric carriers.42 Also, the incorporation or increase in hydrophobic alkyl chains in the cationic polymer was found to increase the binding with DNA.43 In this study, imidazolium-based polyelectrolytes with different alkyl side chain length (ethyl, butyl, and hexyl) were synthesized and their electrostatic and hydrophobic interactions with DNA were examined. For the first time to our knowledge, the binding characteristics, mode of interaction and the effect of imidazolium-based polyelectrolyte alkyl side chains hydrophobicity on the binding affinity to DNA were investigated in detail using different experiments such as UVVisible absorption, fluorescence, CD, agarose gel electrophoresis and FTIR spectroscopy. In addition, molecular docking studies were carried out to visualize the binding mode and preferential docking position of polyelectrolytes in DNA. Furthermore, the structure-activity relationship was explained by comparing the binding affinity of polyelectrolytes with DNA. EXPERIMENTAL SECTION Materials Alkyl halides, such as, n-bromoethane, n-bromobutane, and n-bromohexane; organic solvents; ethidium bromide; TRIS-HCl; NaOH; N-vinylimidazole; AIBN; and sodium salt of calf thymus DNA (ctDNA) were purchased from Sigma Aldrich. AIBN was used after recrystallization. All

ACS Paragon Plus Environment

5

Biomacromolecules

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 37

other purchased chemicals were of analytical grade and used as such without further purification. The imidazolium-based ionic liquids (IL-Ethyl Br, IL-Butyl Br and IL-Hexyl Br) were synthesized according to the literature report.44 10 mM TRIS-HCl buffer was prepared using double distilled deionized water and the pH was adjusted to 7.4 using 0.1 M NaOH solution. 1% w/w of ctDNA in 10 mM TRIS-HCl buffer (pH 7.4) was prepared with occasional stirring and stored in refrigerator at 5 °C. The final concentration of ctDNA was determined by UV-Visible absorption spectroscopy at 260 nm using known molar extinction coefficient, ε260 = 6600 M-1 cm-1 (expressed as molarity of phosphate groups). The absorbance of a dilute solution of ctDNA in buffer solution at 260 nm was recorded to be 0.42091, the concentration, calculated using ε260 was found to be 63.8 µM and the final concentration of ctDNA stock solution was estimated to be 12.7 mM. The ctDNA used for all the experiments was sufficiently free of protein with the ratio of UV absorbance of ctDNA in buffer solution at 260 nm and 280 nm (A260/A280) ᵙ 1.8.45 General procedure for the preparation of polyelectrolytes (PIL-Alkyl Br) Imidazolium-based polyelectrolytes of the type poly(N-alkyl-3-vinylimidazolium bromide)s were synthesized according to the literature report.28 N-vinylimidazole was homopolymerized using free radical polymerization with the degree of polymerization targeted as 50 repeating units. The obtained poly(N-vinylimidazole) was purified by precipitation in cold chloroform and then, quaternized with excess of three different alkyl bromides and used for the DNA binding studies after purification. The quaternization of poly(N-vinylimidazole) with different alkyl halides were followed by 1H NMR and the relative integration confirm the quantitative conversions as shown in Figure S1 (see Supplementary Information). The hydrophilic polyelectrolytes PIL-Alkyl Br were converted to hydrophobic PIL-Alkyl Tf2 by simple anionic metathesis using LiNTf2. Molecular weights of these PIL-Alkyl NTf2 polymers were calculated

ACS Paragon Plus Environment

6

Page 7 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

by Waters gel permission chromatography instrument with Styragel columns equipped with RI detector in THF solvent with a flow rate of 1 mL/minute using polystyrene internal standards. GPC in THF. The molecular weights, Mn (polydispersity index, PDI) of [PIL-Ethyl NTf2], [PIL-Butyl NTf2], and [PIL-Hexyl NTf2] were found to be 20,800(1.40), 24,400 (1.56) and 29,950(1.64) respectively as shown in Figure S2 (see Supplementary Information). DNA-BINDING STUDIES UV-Visible absorption The absorption spectra were recorded on HITACHI U-2910 UV-Visible spectrophotometer with 1 cm quartz cuvette at 298.15 ± 0.15 K. The absorbance was measured at pH 7.4 by keeping the DNA concentration constant (63.5 µM) and varying the polyelectrolyte concentration (0, 0.375, 0.75, 1.25, and 1.6 µM). The solutions were incubated for 5 minutes and then scanned at a wavelength between 200 and 400 nm. The interference due to absorbance of polyelectrolyte was subtracted during the calculation of absorbance of DNA- polyelectrolyte complex. Fluorescence quenching studies The fluorescence quenching studies of ethidium bromide (EtBr)-bound ctDNA were carried out using the HITACHI F-700 fluorescence instrument. The ctDNA at 254 µM concentration and the EtBr at 25.3 µM concentration were prepared in 10 mM TRIS-HCl buffer (pH 7.4) and used in the experiment. The EtBr-DNA complex was titrated using different concentrations (0, 0.25, 0.75, 1.5, 2.5, 4, and 6 µM) of polyelectrolyte. The solutions were gently mixed, incubated for 5 minutes and then the readings recorded at 298.15 ± 0.15 K. Salt effect

ACS Paragon Plus Environment

7

Biomacromolecules

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 37

The effect of ionic strength on the fluorescence quenching of EtBr-DNA complex by polyelectrolyte was studied using various concentrations of NaCl (0.2 M, 0.6 M and 0.8 M). The experiment was performed using 254 µM concentration of ctDNA, 25.3 µM concentration of EtBr and different concentrations (0.25, 0.45, 0.75, 1, 1.25, 1.45, 1.75, 2, 2.25, and 2.45 µM) of PIL-Ethyl Br in the presence of NaCl. Circular dichroism Circular dichroism experiments were performed on JASCO J-715 spectropolarimeter and the spectra were recorded using 1 cm path length rectangular quartz cuvette at 298.15 ± 0.15 K. The circular dichroism titrations were performed using fixed concentration of ctDNA (635 µM) and various concentrations (0, 0.16, 0.33, 0.5, 0.66, 1.0 and 1.3 µM) of polyelectrolyte prepared in 10 mM Tris-HCl buffer and scanned at a wavelength between 220 and 320 nm, bandwidth of 1 nm, response time of 1 sec and an average of three scans to reduce the signal to noise ratio. The baseline of the ctDNA- polyelectrolyte complexes CD spectra was corrected by subtracting with that of 10 mM TRIS-HCl buffer and the experiment was carried out in an atmosphere of nitrogen. Agarose gel electrophoresis The electrophoretic mobility of ctDNA- polyelectrolyte complexes at varying concentration of polyelectrolytes was determined by gel electrophoresis using 1% agarose gel in a buffer containing 45 mM TRIS-Borate buffer and 1 mM EDTA buffer (TBE). The ctDNApolyelectrolyte complexes were prepared by mixing appropriate amount of 12.7 mM DNA and various concentrations of polyelectrolytes (0, 0.1, 0.3, 0.5, 0.9 and 1.1 mM) and allowed to incubate for 30 minutes at 25 °C. These solutions were placed carefully inside the wells of 1%

ACS Paragon Plus Environment

8

Page 9 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

agarose gel along with the gel loading dye to visualize the gel run. The gel was run at 100 V for 45 minutes and then the DNA was visualized under the UV transilluminator after staining the gels with ethidium bromide overnight at 25 °C. The images of DNA were obtained using gel documentation instrument. FTIR spectroscopy FTIR spectra of ctDNA and its complexes with polyelectrolyte were recorded with SHIMADZUIR AFFINITY-1 spectrophotometer equipped with DLATGS (deuterated L-alanine triglycene sulfate) detector and KBr beam splitter. The ctDNA-polyelectrolyte complexes were made using 10 mM TRIS-HCl buffer solution to give the desired polyelectrolyte/ctDNA ratios of 1/25, 1/10 and 1/5 at the final ctDNA concentration of 12.7 mM. The prepared solutions were incubated for 30 minutes at 25 °C and then recorded using KBr pellets. Molecular docking Molecular docking studies were carried out using the Glide (Grid-Based Ligand Docking With Energetics)46, 47 program incorporated in the Schrodinger molecular modeling package (Schrodinger, Inc., USA) which is an interactive molecular graphics program for the interaction, docking calculations and identification of the possible binding site of the biomolecules.48, 49 The receptor, DNA duplex (PDB ID 425D) retrieved from Protein Data Bank (www.rcsb.org) contained 12 base pairs.50 The B-DNA structure was optimized using the Protein Preparation Wizard. The B-DNA crystal structure was preprocessed by deleting the crystallographically observed water molecules (water without H bonds) as no water molecule was found to be conserved, and optimizing the hydrogen bonds. The hydrogen atoms were added to the DNA structure corresponding to pH 7.0. After assigning the charge and protonation state, energy

ACS Paragon Plus Environment

9

Biomacromolecules

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 37

minimization with root mean square deviation (RMSD) value of 0.30Å was carried out using Optimized Potentials for Liquid Simulations (OPLS-2005) force field. On the other hand, the 3D geometries of polyelectrolytes (pentamer) were optimized using LigPrep module and the partial charges were ascribed using OPLS-2005 force-field as performed in Glide. The binding pocket of B-DNA is defined by a 14× 14 × 14Å box that is centered on the geometric centroid of BDNA structure. Sufficiently large grids were chosen to include a significant part of B-DNA. Extra-precision scoring function in Glide was used to rank the docking poses and measure the binding affinity of polyelectrolytes (pentamer) to B-DNA.46, 47 RESULTS AND DISCUSSION Poly(N-alkyl-3-vinylimidazolium bromide)s with three different alkyl side chains (ethyl, butyl and hexyl) were synthesized by quaternizing the radically-synthesized poly(N-vinylimidazole) according to the literature report (Scheme 1).28 Since, the back bone of all the employed polymers (PIL-Ethyl Br, PIL -Butyl Br, and PIL -Hexyl Br) of this study are from the same unit i,e poly(N-vinylimidazole), it facilitates a comparative study with the respective alkyl side chain length. (Scheme 1 here) DNA BINDING STUDIES The thermodynamic parameters and electrostatic interactions between ionic liquids of the type Nalkyl-3-methylimidazoliums and DNA were already reported in the literature.51-54 As a preliminary investigation, we have studied the fluorescence quenching of EtBr-DNA complex using N-hexyl-3-methylimidazolium bromide (IL-Hexyl Br) (Figure S3). Surprisingly, IL-

ACS Paragon Plus Environment

10

Page 11 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

Hexyl Br found to quench the fluorescence only at high concentrations (0.1 to 1.7 mM) whereas PIL-Hexyl Br showed significant fluorescence quenching at minute concentrations (0.25 to 6 µM) and the fluorescence quenching constants are as follows: KIL-Hexyl Br = 0.0742 × 103 M-1 and KPIL-Hexyl Br = 0.5926 × 106 M-1 (Figure S4). Moreover, the relative binding energy from the molecular docking studies between ILs and DNA was found to be less than that of polyelectrolytes (The detailed docking modes (Figure S5) and docking studies for ILs (IL-Ethyl Br, IL-Butyl Br and IL-Hexyl Br) and DNA are given in electronic supporting information). Hence, these findings have encouraged us to continue the DNA binding studies with imidazolium-based polyelectrolytes. UV-Visible absorption The interaction between polyelectrolytes and ctDNA was examined by UV-Vis absorption spectra of ctDNA with varying concentrations of polyelectrolytes (0.375 to 1.6 µM), which has been illustrated in Figure 1. The absorption spectra of ctDNA showed increase in absorbance at 260 nm with the gradual addition of polyelectrolytes. This hyperchromic shift is ascribed to the electrostatic mode of interaction between the cationic imidazolium groups of polyelectrolytes and phosphate anions of ctDNA.55 The absorbance of ctDNA increased dramatically with 0.375 and 0.75 µM of polyelectrolyte and negligible or minimal increase was noted with 1.25 and 1.6 µM of polyelectrolyte. These observations concurred well with the DNA binding studies of cationic polymers available in the literature.56, 57 It is assumed that the interaction between DNA and ligand (polyelectrolyte) is 1:1 and the equilibrium between them could be established according to the following relationship.58, 59  +  ↔  ∶ 

ACS Paragon Plus Environment

11

Biomacromolecules

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

=

Page 12 of 37

[ ∶  ] [][ ]

The binding constants of polyelectrolytes were calculated according to the literature reports.58-61 A double reciprocal plot of 1/A-A0 and 1/CPolyelectrolyte is linear and the binding constant, K was calculated from the ratio of the intercept to the slope. A0 and A correspond to the absorbance of free ctDNA at 260 nm in the absence of polyelectrolyte and absorbance recorded with the addition of different concentrations of polyelectrolyte respectively. The graphs (Figure 1) followed linear regression and the binding constants increased with the increase in the alkyl chain length of the polyelectrolytes thereby suggesting the contribution of alkyl chain’s hydrophobicity in the interaction between polyelectrolyte and DNA. The estimated binding constants of polyelectrolytes are KPIL-Ethyl Br = 1.6008 × 106 M-1, KPIL-Butyl Br = 1.6493 × 106 M-1 and KPIL-Hexyl Br = 3.2087 × 106 M-1. (Figure 1 here) Fluorescence quenching studies EtBr is a dye that binds to ctDNA by intercalation with DNA bases and generally used as a fluorescent probe in DNA-binding studies. In this experiment, EtBr was used to study the interactions between ctDNA and polyelectrolyte through fluorescent quenching studies. The fluorescence intensity of EtBr-ctDNA complex decreases when ctDNA gets condensed or compacted by the quencher (competing molecule) thereby restricting the intercalation with EtBr.21 In this process, EtBr is leached out freely into the bulk solvent and therefore, decreases the fluorescence intensity of the complex.

ACS Paragon Plus Environment

12

Page 13 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

It is expected that polyelectrolyte would condense the ctDNA and form compact macromolecular structures through electrostatic interaction between polymeric imidazolium cations and phosphate anions of ctDNA, resulting in insufficient space for intercalation with EtBr. The emission spectra of EtBr-ctDNA complex with varying concentrations of polyelectrolytes are represented in Figure 2. The fluorescence intensity of EtBr-ctDNA complex decreased concomitantly with the addition of increasing concentrations (0.25 to 6 µM) of polyelectrolytes indicative of fluorescence quenching by the polyelectrolytes employed. The fluorescence intensity of EtBr-DNA complex decreased steadily from 0.25 to 2.5 µM of polyelectrolyte and then, decreased abruptly from 4 to 6 µM of polyelectrolyte. The fluorescence quenching constant, Ksv was calculated using the following Stern-Volmer equation:62  = 1 +  []  Here, F0 is the steady-state fluorescence intensity of EtBr-DNA complex in the absence of polyelectrolyte, F is the intensity of EtBr-DNA complex in the presence of polyelectrolyte (quencher), [] is the concentration of quencher, and Ksv was obtained from the slope of [] vs. F0/F plot. The graphs (Figure 2.d) follow almost linear regression and the quenching constant increases with the increase in the polyelectrolyte’s alkyl side chain length, which is in consistency with the binding trend observed in UV-Vis absorption spectral results. KPIL-Ethyl Br = 0.2676 × 106 M-1, KPIL--Butyl Br = 0.3167 × 106 M-1 and KPIL--Hexyl Br = 0.5926 × 106 M-1 are the calculated values of polyelectrolyte quenching constants. This effect of alkyl side chain length on the quenching of fluorescence intensity may be attributed to the increase in hydrophobic interactions between polyelectrolyte alkyl side chains and ctDNA along with the electrostatic interactions.53 This trend is in agreement with the reported literature regarding binding studies between ctDNA and ILs with different alkyl side chains.54

ACS Paragon Plus Environment

13

Biomacromolecules

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 37

(Figure 2 here) Salt Effect Monitoring the spectral changes with varying ionic strength is an efficient method to identify the binding mode between molecules and DNA.63 The addition of NaCl would weaken the electrostatic interactions between phosphate backbone of DNA and cationic molecules due to their competition with sodium ions. The emission spectra of EtBr-DNA complex with the addition of polyelectrolyte in the presence of varying concentrations of NaCl (0.2, 0.6 and 0.8 M) were monitored in order to evaluate the contribution of electrostatic interaction in the binding between ctDNA and polyelectrolyte. PIL-Ethyl Br was used as an example and the emission spectra were represented in Figure S6 (see Supplementary Information). As shown in Figure 3, the fluorescence quenching of EtBr-DNA complex by polyelectrolyte decreased drastically with the increase in NaCl concentration, which is ascribed to the competition between cationic imidazolium moiety of polyelectrolyte and sodium ions. The fluorescence quenching constants of the polyelectrolyte also decreased with the increase in the ionic strength suggesting the strong dependence of fluorescence intensity on ionic strength and the sequence is as follows: K0 M = 0.1455 × 106 M-1> K0.2 M NaCl = 0.1010 × 106 M-1> K0.6 M NaCl = 0.04828 × 106 M-1> K0.8 M NaCl = 0.03510 × 106 M-1. The results clearly confirm the mode of interaction between DNA and polyelectrolytes as mainly electrostatic.56 (Figure 3 here) Circular dichroism

ACS Paragon Plus Environment

14

Page 15 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

The conformational changes in the secondary structure of DNA before and after the addition of polyelectrolytes were studied using CD spectropolarimeter. The CD spectra of ctDNA after the addition of varying concentrations (0.16 to 1.3 µM) of polyelectrolytes are shown in Figure 4. The characteristic elliptic peaks at 247 nm (negative) and 275 nm (positive) are caused due to ππ base stacking and helicity of ctDNA respectively. This further suggests that the ctDNA retained the B-form after binding with polyelectrolytes.64 On the other hand, the employed polyelectrolytes do not have any optical activity in the present investigation. The molar ellipticity progressively decreased at both positive (247 nm) and negative (275 nm) bands with increase in the concentration of polyelectrolytes. This shift at both the positive and negative bands indicated an electrostatic mode of interaction between the phosphate backbone of DNA and cationic imidazolium groups of polyelectrolytes that is in accordance with UV-Vis absorption and salt effect results. Moreover, a slight red shift was observed at the cross over region in all the CD spectra of ctDNA bound to polyelectrolytes. Also, the decrease in the intensity of CD signals with increase in the concentration of polyelectrolyte may be due to the condensation of ctDNA by the polyelectrolyte.65 From the CD results, it is evident that there are no significant alterations in the secondary structure of DNA after binding to polyelectrolytes and the DNA’s helicity was well retained. Though the change in molar ellipticity at 247 nm did not follow a regular trend with respect to alkyl chains; the ellipticity at 275 nm decreased appreciably with the increase in polyelectrolyte’s alkyl side chain length as shown in Figure S7.a and b (see Supplementary Information). This trend further supported the results discussed previously in the UV-Vis absorption and fluorescence studies. (Figure 4 here) Agarose gel electrophoresis

ACS Paragon Plus Environment

15

Biomacromolecules

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 37

Gel electrophoresis allows the visualization of the interaction between ctDNA and polyelectrolyte. The representative gel images of ctDNA-polyelectrolyte complexes with different concentrations and alkyl side chains of polyelectrolytes are depicted in Figure 5. Lane 0 represents the illumination of free ctDNA and lanes 1-6 in Figure 5.a, b and lanes 1-4 in Figure 5.c represent the illumination of incompletely neutralized ctDNA migrating in the electric field toward the anode. The illumination decreased with the increase in the concentration and alkyl side chain length of polyelectrolyte, thus indicating the formation of ctDNA-polyelectrolyte complexes. No illumination was observed in lanes 5 and 6 in Figure 5.c showing complete neutralization of anionic phosphate groups of ctDNA with cationic imidazolium groups of polyelectrolyte. From the gel electrophoresis, it could be summarized that ctDNA-PIL-Hexyl Br complex showed complete neutralization with no illumination in 5 & 6 lanes as seen in Figure 5c; whereas, ctDNA-PIL-Ethyl Br and ctDNA-PIL-Butyl Br complexes showed partial neutralization. (Figure 5 here) FTIR spectroscopy FTIR spectra of free ctDNA and ctDNA-polyelectrolyte complexes with different molar ratios of polyelectrolyte are represented in Figure 6. The vibrational bands of free ctDNA at 1261, 1402, and 1639 cm-1 are assigned to thymidine, cytosine, and adenine bases, respectively, according to the literature reports.66-69 Bands at 1087 and 991 cm-1 denote symmetric stretching and bending frequencies of phosphate groups, respectively.69 After the addition of polyelectrolyte to ctDNA, no characteristic change was observed in DNA base-pair bands (thymidine,1261 cm-1 and adenine, 1639 cm-1), except for a slight shift in cytosine band from 1402 to 1400 cm-1 at 1/10

ACS Paragon Plus Environment

16

Page 17 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

molar ratio of polyelectrolyte/ctDNA implying no considerable intercalative mode of interaction between polyelectrolyte and ctDNA bases. A prominent shift in phosphate asymmetric stretching from 1083 to 1087 cm-1 was seen at 1/10 and 1/5 molar ratios of polyelectrolyte (PIL-Ethyl Br, PIL -Butyl Br and PIL -Hexyl Br); whereas, shift from 1083 to 1085 cm-1 was noted at 1/25 molar ratio of PIL-Butyl Br and PIL-Hexyl Br. Also, phosphate bending at 990 cm-1 shifted slightly to 989 cm-1 at all the molar ratios of all the three polyelectrolytes employed and the shifts of bands are close to the experimental error of the equipment. It may be concluded that there was no considerable change in the bands of DNA base pairs whereas the contrary was observed in the phosphate asymmetric bands implying an electrostatic interaction between the oppositely charged groups of polyelectrolytes and DNA. (Figure 6 here) Molecular Docking In an effort to understand and interpret the molecular mechanism of the interaction between polyelectrolytes with DNA, molecular docking was performed to simulate the modes of interactions between the polyelectrolytes with five repeating monomeric units and B-DNA. Structures of polyelectrolyte and DNA were kept flexible to attain different conformations in order to predict the best docked complex. The study revealed that groove binding makes intimate contacts with the walls of the groove, and as a result of this interaction, numerous electrostatic and van der Waals interactions occur between the units of polyelectrolyte and DNA bases and its phosphate backbone. From the ensuing docked structures, it is clear that all the polyelectrolytes fit snugly into the curved contour of the targeted DNA in the minor groove, with the walls of the groove in close contact with imidazolium scaffold and associated alkyl side chains stabilized by

ACS Paragon Plus Environment

17

Biomacromolecules

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 37

strong electrostatic interaction (Figure 7). Relative binding energy docking scores of the docked polymeric structures was found to be -74.66 kcal/mol (-10.56) , -83.08 kcal/mole (-11.67) and 90.91 kcal/mol (-13.32) for pentamers (PIL-Ethyl Br)5, (PIL-Butyl Br)5, and (PIL-Hexyl Br)5 respectively. The binding energy of PIL-Ethyl Br and PIL-Butyl Br was found to be lower than that of PIL-Hexyl Br. Higher negative binding energy indicates a more stable combination with DNA. While a consistent interaction was observed between the imidazolium nucleus and the DNA base pairs across the three polyelectrolytes, it was the varying alkyl chain length that contributed to the variation in the binding strengths of these polyelectrolytes. Thus, it is concluded that the binding ability of PIL-Hexyl Br with DNA is stronger than that of PIL-Ethyl Br and PIL-Butyl Br. The docking scores and the corresponding binding energy values revealed an increase in polyelectrolytes binding affinity with increase in their alkyl side chain length (KPIL-Ethyl Br > KPIL-Butyl Br > KPIL-Hexyl Br). Furthermore, analysis of the binding modes of these polyelectrolytes revealed that the B-form was well retained by the DNA with minor perturbation in its phosphate backbone and no significant change in the base pairs after complexation with polyelectrolytes, which is in good agreement with the experimental observations. (Figure 7 here)

CONCLUSION Imidazolium-based polyelectrolytes with three different alkyl side chains (ethyl, butyl, and hexyl) derived radically, were successfully employed to study their interaction with ctDNA. The absorption spectra of ctDNA showed hyperchromic shift with the increase in polyelectrolyte concentration that may be ascribed to the electrostatic interaction between the cationic

ACS Paragon Plus Environment

18

Page 19 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

imidazolium groups of polyelectrolyte and phosphate anions of ctDNA. This electrostatic mode of interaction was further confirmed by the effect of ionic strength on the binding of polyelectrolyte with ctDNA. The infrared and CD spectroscopic results showed evidence that the DNA remained in the B-form without much perturbation to the secondary structure after binding to polyelectrolytes. The binding constants of polyelectrolytes obtained from the absorption spectra followed the same trend of KPIL-Ethyl Br < KPIL -Butyl Br < KPIL -Hexyl Br as that of the quenching constants calculated from the emission spectra. This binding trend demonstrates an increase in binding between polyelectrolytes and DNA with an increase in polyelectrolyte’s alkyl side chain length suggesting the contribution of hydrophobic interaction between polyelectrolyte’s alkyl side chains and DNA apart from the electrostatic interaction. The above binding trend was further confirmed by molecular docking studies with respect to the binding energy of polyelectrolytes towards DNA. Furthermore, it was visualized that polyelectrolytes (pentamer) showed groove binding with DNA along with predominant electrostatic interaction between the oppositely charged groups.

ACS Paragon Plus Environment

19

Biomacromolecules

Figure 1. UV-Vis absorption spectra of free ctDNA (63.5 µM) and ctDNA- polyelectrolyte complexes with different concentrations (0, 0.375, 0.75, 1.25 and 1.6 µM) of polyelectrolyte (i) PIL-Ethyl Br (ii) PIL-Butyl Br and (iii) PIL-Hexyl Br. PIL-Butyl Br - DNA complex

PIL-Ethyl Br - DNA Complex

0.18

K = 1.6008 x 10 6 M-1 b

0.12

a

1/(A-A0)

0.15

Normalized Absorbance

40 d

30 a

20

d b

c

10

0 0.0 0.5 1.0 1.5 2.0 2.5 3.0

0.09

6

1/[Polyelectrolyte] x 10 M

0.06

-1

c DNA

0.03

K = 1.6493 x 106 M-1

d b

0.12

a

30 20

a

d c

b

10 0 0.0 0.5 1.0 1.5 2.0 2.5 3.0

0.09

1/[Polyelectrolyte] x 10 6 M -1

0.06 c

DNA

0.03 0.00

240 260 280 300 320 340 360 380 400 Wavelength (nm)

240 260 280 300 320 340 360 380 400 Wavelength (nm)

(i)

(ii) PIL-Hexyl Br - DNA complex 0.18 35

0.15

d

c b

0.12

K = 3.2087 x 106 M-1 1/(A-A0)

Normalized Absorbance

0.00

40

0.15

1/(A-A0)

0.18

Normalized Absorbance

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 37

a

30 25 d

0.09

c

b 20 0.0 0.5 1.0 1.5 2.0 2.5 3.0

a

1/[Polyelectrolyte] x 10 6 M-1

0.06 DNA

0.03 0.00

240 260 280 300 320 340 360 380 400 Wavelength (nm)

(iii)

ACS Paragon Plus Environment

20

Page 21 of 37

Figure 2. Fluorescence emission spectra of EtBr-ctDNA complex with the addition of different concentrations (0, 0.25, 0.75, 1.5, 2.5, 4, and 6 µM) of (i) PIL-Ethyl Br, (ii) PIL-Butyl Br and (iii) PIL-Hexyl Br. The comparative graph (iv) represents the decrease in relative fluorescence intensity of EtBr-ctDNA complex with increase in alkyl side chain length of polyelectrolytes.

PIL-Ethyl Br - DNA Complex

PIL-Butyl Br - DNA Complex

35

30

b

c

20 d e

15

3

f

d 2 b e 1 c a 0 0 1 2 3 4 5 6 7 8 [Polyelectrolyte] x 10 -6 M

f

10 5 0

5

DNA

25

6 -1 4 KSV = 0.3167 x 10 M

a

F0/F

a

25

Fluorescence Intensity

6 -1 4 KSV = 0.2676 x 10 M

DNA

30

F0/F

Fluorescence Intensity

5

b

20

c

15

d e

10

d

2 b e 1 a c 0 0 1 2 3 4 5 6 7 8 [Polyelecctrolyte] x 10 -6 M

f

5 0

550 575 600 625 650 675 700 725 750 Wavelength (nm)

f

3

550 575 600 625 650 675 700 725 750 Wavelength (nm)

(i)

(ii)

PIL-Hexyl Br - DNA complex 30

5

5

20

a

b c

15 10

d e

4 e

3 2

a b

c

f

PIL-Hexyl Br PIL-Butyl Br PIL-Ethyl Br

d

1

KSV = 0.5926 x 106 M-1 0 0 1 2 3 4 5 6 7 8 [Polyelectrolyte] x 10 -6 M

f

3 2 1

5 0

4 F0/F

DNA

25

F0/F

Fluorescence Intensity

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

550 575 600 625 650 675 700 725 750 Wavelength (nm)

0 0

1

2

3

4

5

6

7

[Polyelectrolyte] x 10-6 M

(iii)

(iv)

ACS Paragon Plus Environment

21

Biomacromolecules

Figure 3. Graph representing the decrease in relative fluorescence intensity of EtBr bound ctDNA with increase in the ionic strength.

Salt Effect 1.5 PIL-Ethyl Br 0.2 M NaCl 0.6 M NaCl 0.8 M NaCl

1.4 1.3

F0/F

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 37

1.2 1.1 1.0 0.9 0.0

0.5

1.0

1.5

2.0

2.5

3.0

[Polyelectrolyte] x 10-6 M

ACS Paragon Plus Environment

22

Page 23 of 37

Figure 4. CD spectra of DNA (635 µM) with different concentrations (0, 0.16, 0.33, 0.5, 0.66, 1.0, and 1.3 µM) of polyelectrolytes (a) PIL-Ethyl Br (b) PIL-Butyl Br and (c) PIL-Hexyl Br.

PIL-Ethyl Br - DNA complex

PIL-Butyl Br - DNA complex

40 30

CD (mdeg)

20

0 0.5 µM

-20

-40 220

240

0 µM

0.66 µM

0.16 µM

1.0 µM

0.33 µM

1.3 µM

260 280 Wavelength (nm)

300

15 0 0.5 µM

-15 -30

320

220

240

0 µM

0.66 µM

0.16 µM

1.0 µM

0.33 µM

1.3 µM

260 280 Wavelength (nm)

(a)

300

320

(b) PIL-Hexyl Br - DNA complex

20

CD (mdeg)

CD (mdeg)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

0 0.5 µM

-20

220

240

0 µM

0.66 µM

0.16 µM

1.0 µM

0.33 µM

1.3 µM

260 280 Wavelength (nm)

300

320

(c)

ACS Paragon Plus Environment

23

Biomacromolecules

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 37

Figure 5. Electrophoretic mobility of ctDNA- polyelectrolyte complexes on 1% agarose gel (a) PIL-Ethyl Br (b) PIL-Butyl Br (c) PIL-Hexyl Br. Lanes 0, 1, 2, 3, 4, 5 and 6 correspond to ctDNA- polyelectrolyte complexes with 12.7 mM DNA and 0, 0.1, 0.3, 0.5, 0.9 and 1.1 mM concentrations of polyelectrolyte respectively.

(a)

(b)

(c)

ACS Paragon Plus Environment

24

Page 25 of 37

Figure 6. Stacked view of FTIR spectra of polyelectrolyte/ctDNA complexes in the molar ratios of 1/25, 1/10 and 1/5 in the region of 1800 to 800 cm-1.

1400

1200

1000

800

1800

Wavelength (cm-1)

991

1083

1402

1639

989

1085

1261

1402

1639

989

1087

1261

1402

Arbitrary Unit

989

1085

989

1087

1261

989

1087

1261

1400

1639

989

1087

1261

1400

1/5

1600

1/25

1/10

1639

989

1087

1261

1400

1639

1/10

1/5

1261

991

1083

1261 1261

1402

1639 1/25

1402

989

1087

1261

1402

1639 1/10

DNA

1639

Arbitrary Unit

989

1083

1261

1402

1639

DNA

1/25

1800

PIL-Hexyl Br

1639

DNA

1402

1639

991

PIL-Butyl Br

1083

1261

1402

1639

PIL-Ethyl Br

Arbitrary Unit

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

1/5

1600

1400

1200

1000

800

Wavelength (cm-1)

(a)

1800

1600

1400

1200

1000

800

Wavelength (cm-1)

(b)

(c)

Figure 7. Binding mode of pentamer of polyelectrolytes (a) (PIL-Ethyl Br)5, (b) (PIL-Butyl Br)5, and (c) (PIL-Hexyl Br)5 with B-DNA.

(a)

(b)

(c)

ACS Paragon Plus Environment

25

Biomacromolecules

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 37

Scheme 1. Imidazolium-based poly(ionic liquids) used in this study.

ASSOCIATED CONTENT Supporting Information In the electronic supplementary information we have provided results of fluorescence quenching of EtBr-DNA complex using IL-Hexyl Br and the detailed docking studies for ILs (IL-Ethyl Br, IL-Butyl Br and IL-Hexyl Br) and DNA. This material is available free of charge via the Internet at http://pubs.acs.org. AUTHOR INFORMATION Corresponding Author *Kari Vijayakrishna. Mailing address: Organic Chemistry Division, School of Advanced Sciences, VIT University, Vellore 632014, Tamil Nadu, India. Phone: +91 416 224 2334. Fax: +91 416224 3092. E-mail: [email protected]. Author Contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript.

ACS Paragon Plus Environment

26

Page 27 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

ACKNOWLEDGMENT Kari thank DST-SERB, India (Project NO: SR/S1/OC-22/2012) for the financial support. Authors also thank DST-VIT-FIST for NMR, VIT-SIF for other instrumentation facilities. REFERENCES (1)

Partridge, K. A.; Oreffo, R. O. C. Gene Delivery in Bone Tissue Engineering: Progress and Prospects Using Viral and Nonviral Strategies. Tissue Eng. 2004, 10, 295-307.

(2)

Haider, M.; Megeed, Z.; Ghandehari, H. Genetically engineered polymers: status and prospects for controlled release. J. Control. Release 2004, 95, 1-26.

(3)

Thomas, C. E.; Ehrhardt, A.; Kay, M. A. Progress and problems with the use of viral vectors for gene therapy. Nat. Rev. Genet. 2003, 4, 346-358.

(4)

Pathak, A.; Patnaik, S.; Gupta, K. C. Recent trends in non-viral vector-mediated gene delivery. Biotechnol. J. 2009, 4, 1559-1572.

(5)

Park, T. G.; Jeong, J. H.; Kim, S. W. Current status of polymeric gene delivery systems. Adv. Drug Deliv. Rev. 2006, 58, 467-486.

(6)

Yue, Y.; Wu, C. Progress and perspectives in developing polymeric vectors for in vitro gene delivery. Biomaterials Science 2013, 1, 152-170.

(7)

Schaffert, D.; Wagner, E. Gene therapy progress and prospects: synthetic polymer-based systems. Gene Ther. 2008, 15, 1131-1138.

(8)

Mintzer, M. A.; Simanek, E. E. Nonviral Vectors for Gene Delivery. Chem. Rev. 2008, 109, 259-302.

ACS Paragon Plus Environment

27

Biomacromolecules

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(9)

Page 28 of 37

Samal, S. K.; Dash, M.; Van Vlierberghe, S.; Kaplan, D. L.; Chiellini, E.; van Blitterswijk, C.; Moroni, L.; Dubruel, P. Cationic polymers and their therapeutic potential. Chem. Soc. Rev. 2012, 41, 7147-7194.

(10) Boussif, O.; Lezoualc'h, F.; Zanta, M. A.; Mergny, M. D.; Scherman, D.; Demeneix, B.; Behr, J. P. A versatile vector for gene and oligonucleotide transfer into cells in culture and in vivo: polyethylenimine. Proc. Natl. Acad. Sci. 1995, 92, 7297-7301. (11) Kabanov, A. V.; Kabanov, V. A. DNA Complexes with Polycations for the Delivery of Genetic Material into Cells. Bioconjugate Chem. 1995, 6, 7-20. (12) Newland, B.; Zheng, Y.; Jin, Y.; Abu-Rub, M.; Cao, H.; Wang, W.; Pandit, A. Single Cyclized Molecule Versus Single Branched Molecule: A Simple and Efficient 3D “Knot” Polymer Structure for Nonviral Gene Delivery. J. Am. Chem. Soc. 2012, 134, 4782-4789. (13) Zhou, D.; Li, C.; Hu, Y.; Zhou, H.; Chen, J.; Zhang, Z.; Guo, T. PLL/pDNA/P(His-coDMAEL) ternary complexes: assembly, stability and gene delivery. J. Mater. Chem. 2012, 22, 10743-10751. (14) Zhao, T.; Zhang, H.; Newland, B.; Aied, A.; Zhou, D.; Wang, W. Significance of Branching for Transfection: Synthesis of Highly Branched Degradable Functional Poly(dimethylaminoethyl methacrylate) by Vinyl Oligomer Combination. Angew. Chem., Int. Ed. 2014, 53, 6095-6100. (15) Jones, N. A.; Hill, I. R. C.; Stolnik, S.; Bignotti, F.; Davis, S. S.; Garnett, M. C. Polymer chemical structure is a key determinant of physicochemical and colloidal properties of polymer–DNA complexes for gene delivery. Biochim. Biophys. Acta 2000, 1517, 1-18.

ACS Paragon Plus Environment

28

Page 29 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

(16) Rungsardthong, U.; Ehtezazi, T.; Bailey, L.; Armes, S. P.; Garnett, M. C.; Stolnik, S. Effect of Polymer Ionization on the Interaction with DNA in Nonviral Gene Delivery Systems. Biomacromolecules 2003, 4, 683-690. (17) Synatschke, C. V.; Schallon, A.; Jérôme, V.; Freitag, R.; Müller, A. H. E. Influence of Polymer

Architecture

and

Molecular

Weight

of

Poly(2-(dimethylamino)ethyl

methacrylate) Polycations on Transfection Efficiency and Cell Viability in Gene Delivery. Biomacromolecules 2011, 12, 4247-4255. (18) Funhoff, A. M.; van Nostrum, C. F.; Lok, M. C.; Fretz, M. M.; Crommelin, D. J. A.; Hennink, W. E. Poly(3-guanidinopropyl methacrylate):  A Novel Cationic Polymer for Gene Delivery. Bioconjugate Chem. 2004, 15, 1212-1220. (19) Liu, X.; Yang, J. W.; Miller, A. D.; Nack, E. A.; Lynn, D. M. Charge-Shifting Cationic Polymers

That

Promote

Self-Assembly

and

Self-Disassembly

with

DNA.

Macromolecules 2005, 38, 7907-7914. (20) Pack, D. W.; Putnam, D.; Langer, R. Design of imidazole-containing endosomolytic biopolymers for gene delivery. Biotechnol. Bioeng. 2000, 67, 217-223. (21) Vuorimaa, E.; Urtti, A.; Seppänen, R.; Lemmetyinen, H.; Yliperttula, M. Time-Resolved Fluorescence Spectroscopy Reveals Functional Differences of Cationic Polymer−DNA Complexes. J. Am. Chem. Soc. 2008, 130, 11695-11700. (22) Gao, Y.; Yin, Q.; Chen, L.; Zhang, Z.; Li, Y. Linear Cationic Click Polymers/DNA Nanoparticles: In Vitro Structure–Activity Relationship and In Vivo Evaluation for Gene Delivery. Bioconjugate Chem. 2011, 22, 1153-1161.

ACS Paragon Plus Environment

29

Biomacromolecules

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 37

(23) Mecerreyes, D. Polymeric ionic liquids: Broadening the properties and applications of polyelectrolytes. Prog. Polym. Sci. 2011, 36, 1629-1648. (24) Yuan, J.; Mecerreyes, D.; Antonietti, M. Poly(ionic liquid)s: An update. Prog. Polym. Sci. 2013, 38, 1009-1036. (25) Yuan, J.; Antonietti, M. Poly(ionic liquid)s: Polymers expanding classical property profiles. Polymer 2011, 52, 1469-1482. (26) Vijayakrishna, K.; Jewrajka, S. K.; Ruiz, A.; Marcilla, R.; Pomposo, J. A.; Mecerreyes, D.; Taton, D.; Gnanou, Y. Synthesis by RAFT and Ionic Responsiveness of Double Hydrophilic Block Copolymers Based on Ionic Liquid Monomer Units. Macromolecules 2008, 41, 6299-6308. (27) Vijayakrishna, K.; Mecerreyes, D.; Gnanou, Y.; Taton, D. Polymeric Vesicles and Micelles Obtained by Self-Assembly of Ionic Liquid-Based Block Copolymers Triggered by Anion or Solvent Exchange. Macromolecules 2009, 42, 5167-5174. (28) Prabhu Charan, K. T.; Pothanagandhi, N.; Vijayakrishna, K.; Sivaramakrishna, A.; Mecerreyes, D.; Sreedhar, B. Poly(ionic liquids) as “smart” stabilizers for metal nanoparticles. Eur. Polym. J. 2014, 60, 114-122. (29) Hemp, S. T.; Allen, M. H.; Green, M. D.; Long, T. E. Phosphonium-Containing Polyelectrolytes for Nonviral Gene Delivery. Biomacromolecules 2011, 13, 231-238.

ACS Paragon Plus Environment

30

Page 31 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

(30) Hemp, S. T.; Allen, M. H.; Smith, A. E.; Long, T. E. Synthesis and Properties of Sulfonium Polyelectrolytes for Biological Applications. ACS Macro Lett. 2013, 2, 731735. (31) Ren, X.; Xu, Q.-H. Label-Free DNA Sequence Detection with Enhanced Sensitivity and Selectivity Using Cationic Conjugated Polymers and PicoGreen. Langmuir 2008, 25, 4347. (32) Davies, M. L.; Burrows, H. D.; Cheng, S.; Morán, M. C.; Miguel, M. d. G.; Douglas, P. Cationic Fluorene-Based Conjugated Polyelectrolytes Induce Compaction and Bridging in DNA. Biomacromolecules 2009, 10, 2987-2997. (33) San Juan, A.; Letourneur, D.; Izumrudov, V. A. Quaternized Poly(4-vinylpyridine)s as Model Gene Delivery Polycations:  Structure−Function Study by Modification of Side Chain Hydrophobicity and Degree of Alkylation. Bioconjugate Chem. 2007, 18, 922-928. (34) Han, F.; Lu, Y.; Zhang, Q.; Sun, J.; Zeng, X.; Li, C. Homogeneous and sensitive DNA detection based on polyelectrolyte complexes of cationic conjugated poly(pyridinium salt)s and DNA. J. Mater. Chem. 2012, 22, 4106-4112. (35) Allen, M. H.; Green, M. D.; Getaneh, H. K.; Miller, K. M.; Long, T. E. Tailoring Charge Density and Hydrogen Bonding of Imidazolium Copolymers for Efficient Gene Delivery. Biomacromolecules 2011, 12, 2243-2250. (36) Zhou, J.; Liu, J.; Cheng, C. J.; Patel, T. R.; Weller, C. E.; Piepmeier, J. M.; Jiang, Z.; Saltzman, W. M. Biodegradable poly(amine-co-ester) terpolymers for targeted gene delivery. Nat. Mater. 2012, 11, 82-90.

ACS Paragon Plus Environment

31

Biomacromolecules

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 37

(37) Kurisawa, M.; Yokoyama, M.; Okano, T. Transfection efficiency increases by incorporating hydrophobic monomer units into polymeric gene carriers. J. Controlled Release 2000, 68, 1-8. (38) Yu, S.; Chen, J.; Dong, R.; Su, Y.; Ji, B.; Zhou, Y.; Zhu, X.; Yan, D. Enhanced gene transfection efficiency of PDMAEMA by incorporating hydrophobic hyperbranched polymer cores: effect of degree of branching. Polym. Chem. 2012, 3, 3324-3329. (39) Liu, Z.; Zhang, Z.; Zhou, C.; Jiao, Y. Hydrophobic modifications of cationic polymers for gene delivery. Prog. Polym. Sci. 2010, 35, 1144-1162. (40) Kuhn, P. S.; Levin, Y.; Barbosa, M. C. Charge inversion in DNA–amphiphile complexes: possible application to gene therapy. Physica A 1999, 274, 8-18. (41) Alvarez-Lorenzo, C.; Barreiro-Iglesias, R.; Concheiro, A.; Iourtchenko, L.; Alakhov, V.; Bromberg,

L.; Temchenko, M.; Deshmukh,

S.; Hatton, T. A. Biophysical

Characterization of Complexation of DNA with Block Copolymers of Poly(2dimethylaminoethyl) Methacrylate, Poly(ethylene oxide), and Poly(propylene oxide). Langmuir 2005, 21, 5142-5148. (42) Gabrielson, N. P.; Pack, D. W. Acetylation of Polyethylenimine Enhances Gene Delivery via Weakened Polymer/DNA Interactions. Biomacromolecules 2006, 7, 2427-2435. (43) Matulis, D.; Rouzina, I.; Bloomfield, V. A. Thermodynamics of Cationic Lipid Binding to DNA and DNA Condensation:  Roles of Electrostatics and Hydrophobicity. J. Am. Chem. Soc. 2002, 124, 7331-7342.

ACS Paragon Plus Environment

32

Page 33 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

(44) Dzyuba,

S.

V.;

Bartsch,

R.

A.

Efficient

synthesis

of

1-alkyl(aralkyl)-3-

methyl(ethyl)imidazolium halides: Precursors for room-temperature ionic liquids. J. Heterocyclic Chem. 2001, 38, 265-268. (45) Marmur, J. A procedure for the isolation of deoxyribonucleic acid from micro-organisms. J. Mol. Biol. 1961, 3, 208-218. (46) Friesner, R. A.; Banks, J. L.; Murphy, R. B.; Halgren, T. A.; Klicic, J. J.; Mainz, D. T.; Repasky, M. P.; Knoll, E. H.; Shelley, M.; Perry, J. K.; Shaw, D. E.; Francis, P.; Shenkin, P. S. Glide:  A New Approach for Rapid, Accurate Docking and Scoring. 1. Method and Assessment of Docking Accuracy. J. Med. Chem. 2004, 47, 1739-1749. (47) Halgren, T. A.; Murphy, R. B.; Friesner, R. A.; Beard, H. S.; Frye, L. L.; Pollard, W. T.; Banks, J. L. Glide:  A New Approach for Rapid, Accurate Docking and Scoring. 2. Enrichment Factors in Database Screening. J. Med. Chem. 2004, 47, 1750-1759. (48) Bhakta, D.; Siva, R. Morindone, an Anthraquinone, Intercalates DNA Sans Toxicity: a Spectroscopic and Molecular Modeling Perspective. Appl. Biochem. Biotechnol. 2012, 167, 885-896. (49) Ghosh, P.; Devi, G. P.; Priya, R.; Amrita, A.; Sivaramakrishna, A.; Babu, S.; Siva, R. Spectroscopic and In Silico Evaluation of Interaction of DNA with Six Anthraquinone Derivatives. Appl. Biochem. Biotechnol. 2013, 170, 1127-1137. (50) Rozenberg, H.; Rabinovich, D.; Frolow, F.; Hegde, R. S.; Shakked, Z. Structural code for DNA recognition revealed in crystal structures of papillomavirus E2-DNA targets. P. Natl. Acad. Sci. 1998, 95, 15194-15199.

ACS Paragon Plus Environment

33

Biomacromolecules

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 34 of 37

(51) Nishimura, N.; Nomura, Y.; Nakamura, N.; Ohno, H. DNA strands robed with ionic liquid moiety. Biomaterials 2005, 26, 5558-5563. (52) Xie, Y.-N.; Wang, S.-F.; Zhang, Z.-L.; Pang, D.-W. Interaction between Room Temperature Ionic Liquid [bmim]BF4 and DNA Investigated by Electrochemical Micromethod. J. Phys. Chem. B 2008, 112, 9864-9868. (53) Ding, Y.; Zhang, L.; Xie, J.; Guo, R. Binding Characteristics and Molecular Mechanism of Interaction between Ionic Liquid and DNA. J. Phys. Chem. B 2010, 114, 2033-2043. (54) Singh, P. K.; Sujana, J.; Mora, A. K.; Nath, S. Probing the DNA–ionic liquid interaction using an ultrafast molecular rotor. J. Photochem. Photobio. B 2012, 246, 16-22. (55) Sirajuddin, M.; Ali, S.; Badshah, A. Drug–DNA interactions and their study by UV– Visible, fluorescence spectroscopies and cyclic voltametry. J. Photochem. Photobio. B 2013, 124, 1-19. (56) Dey, D.; Kumar, S.; Banerjee, R.; Maiti, S.; Dhara, D. Polyplex Formation between PEGylated Linear Cationic Block Copolymers and DNA: Equilibrium and Kinetic Studies. J. Phys. Chem. B 2014, 118, 7012-7025. (57) Ottaviani, M. F.; Furini, F.; Casini, A.; Turro, N. J.; Jockusch, S.; Tomalia, D. A.; Messori, L. Formation of Supramolecular Structures between DNA and Starburst Dendrimers Studied by EPR, CD, UV, and Melting Profiles. Macromolecules 2000, 33, 7842-7851.

ACS Paragon Plus Environment

34

Page 35 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

(58) Saito, S. T.; Silva, G.; Pungartnik, C.; Brendel, M. Study of DNA–emodin interaction by FTIR and UV–vis spectroscopy. J. Photochem. Photobio. B 2012, 111, 59-63. (59) Nafisi, S.; Bonsaii, M.; Maali, P.; Khalilzadeh, M. A.; Manouchehri, F. β-Carboline alkaloids bind DNA. J. Photochem. Photobio. B. 2010, 100, 84-91. (60) Klotz, I. M.; Hunston, D. L. Properties of graphical representations of multiple classes of binding sites. Biochemistry 1971, 10, 3065-3069. (61) Connors, K. A. Binding constants: the measurement of molecular complex stability. J. Wiley & Sons. 1987. (62) Li, N.; Ma, Y.; Yang, C.; Guo, L.; Yang, X. Interaction of anticancer drug mitoxantrone with DNA analyzed by electrochemical and spectroscopic methods. Biophys. Chem. 2005, 116, 199-205. (63) Misra, V. K.; Sharp, K. A.; Friedman, R. A.; Honig, B. Salt Effects on Ligand-DNA Binding: Minor Groove Binding Antibiotics. J. Mol. Biol. 1994, 238, 245-263. (64) Kypr, J.; Kejnovská, I.; Renčiuk, D.; Vorlíčková, M. Circular dichroism and conformational polymorphism of DNA. Nucleic Acids Res. 2009, 37, 1713-1725. (65) Mercado, C. M.; Tomasz, M. Circular dichroism of mitomycin-DNA complexes. Evidence for a conformational change in DNA. Biochemistry 1977, 16, 2040-2046. (66) Mello, M. L. S.; Vidal, B. C. Changes in the Infrared Microspectroscopic Characteristics of DNA Caused by Cationic Elements, Different Base Richness and Single-Stranded Form. PLOS ONE 2012, 7, 43169.

ACS Paragon Plus Environment

35

Biomacromolecules

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 36 of 37

(67) Marty, R.; N'Soukpoé-Kossi, C. N.; Charbonneau, D.; Weinert, C. M.; Kreplak, L.; Tajmir-Riahi, H. -A. Structural analysis of DNA complexation with cationic lipids. Nucleic Acids Res. 2009, 37, 849-857. (68) Alex, S.; Dupuis, P. FT-IR and Raman investigation of cadmium binding by DNA. Inorg. Chim. Acta 1989, 157, 271-281. (69) Hembram, K. P. S. S.; Rao, G. M. Studies on CNTs/DNA composite. Mat. Sci. Eng. C 2009, 29, 1093-1097.

ACS Paragon Plus Environment

36

Page 37 of 37

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

TOC

ACS Paragon Plus Environment

37