Fluorescence and Molecular Simulation Studies on the Interaction between Imidazolium-Based Ionic Liquids and Calf Thymus DNA

: This work presents a molecular level investigation on the nature and mode of binding between imidazolium-based ionic liquids (ILs) ([C n bim]Br where n = 2, 4, 6) with calf thymus DNA. This investigation o ﬀ ers valuable insight into the mechanisms of interactions that can a ﬀ ect the structural features of DNA and possibly cause the alteration or inhibition of DNA function. To expedite analysis, the study resorted to using molecular docking and COnductor like Screening MOdel for Real Solvents (COSMO-RS) in conjunction with ﬂuorescence spectroscopic data for conﬁrmation and validation of computational results. Both the ﬂuorescence and docking studies consistently revealed a weak interaction between the two molecules, which corresponded to the binding energy of a stable docking conformation in the range of − 5.19 to − 7.75 kcal mol − 1 . As predicted, the rod-like structure of imidazolium-based ILs prefers to bind to the double-helix DNA through a minor groove. Interestingly, the occurrence of T-shape π - π stacking was observed between the amine group in adenine that faces the aromatic ring of imidazole. In addition, data of COSMO-RS for the interaction of individual nucleic acid bases to imidazolium-based ILs a ﬃ rmed that ILs showed a propensity to bind to di ﬀ erent bases, the highest being guanine followed by cytosine, thymine, uracil, and adenine.


Introduction
In recent years, the study on ionic liquids (ILs) has been a fascinating area of chemical research due to their unique properties such as low toxicity, high thermal and chemical stability, high boiling point, negligible vapor pressure, and excellent conductivity [1]. The applicability and versatility of ILs extend to fields such as organic synthesis [2,3], electrochemistry [4], catalysis and biocatalyst [2,5], separations [6], pharmaceuticals, and medicine [7,8], etc. As a matter of fact, the past decade has seen ILs being utilized in life science applications concerning biomolecules viz. proteins, enzymes, and nucleic acids [9,10].
It is undeniable to say that the utilization of ILs in DNA innovation has been broadly reported in recent years. However, most studies were focused on comprehending the progression and properties of DNA in ILs. Be that as it may, there exist a few conflicting reports with regards to the binding mode

Steady-State Fluorescence
Steady state fluorescence was carried out using a Cary Eclipse Fluorescence Spectrophotometer (Agilent, Santa Clara, CA, USA). The measurements were carried out using a 5.0 nm fixed band slits through a 1.0 cm path length quartz cuvette. All spectra were corrected for background intensities using a solution without DNA. The experiments were performed by continuous addition of 1 µL of 0.3 mM calf thymus DNA solution into a fixed 0.5 M [C n bim]Br (where n = 2, 4, 6) ILs solution for each titration, and the emission intensity was observed. All triplicate titrations were recorded at 25 • C. The excitation wavelength and maximum emission for each [C n bim]Br ILs were detected at 320 and 358 nm, respectively.

Molecular Docking Studies
The crystal structure of calf thymus DNA (PDB ID: 425D) was downloaded from the Research Collaboratory for Structural Bioinformatics (RCSB) Protein Data Bank (https://www.rcsb.org/pdb) [22]. The TmoleX package (Version 4.1.1) (COSMOlogic, Leverkusen, Germany) was employed for geometry optimization of 3 imidazolium-based ILs' cations using the Hartree-Fock theory with 6-31G* basic set [23]. A geometry optimization calculation using the Hartree-Fock level provided more meaningful, accurate values. The * accounts for the polarization effect of the species and was further optimized using the Density Functional Theory (DFT) level with the BP function that consisted of resolution of identity approximation and a triple-ξ valence polarized basis set. A similar optimization method has been employed previously [24,25]. The required pdbqt file for the DNA and cations were prepared using Auto Dock Tools (ADT) 1.5.4 [26]. Partial atomic charges for each atom in DNA were adopted from the Optimized Potential for Liquid Simulations (OPLS) force field [27], while the atomic charges for ILs were obtained from a previous study [18]. Parameters of the optimized ILs structures were successfully used in previous MD simulations and provided reliable data against the experimental evidence. The grid size for the x-, y-and z-direction was set to 100, 56, and 52, respectively, while the center of the search space was set to 29.58, 29.82, and 29.71 for the x-, y-and z-direction. AutoDock4 using Lamarckian Genetic Algorithm (LGA) was used for the docking studies between the DNA and the imidazolium-cation [28], and default values were used for other setting parameters. In this study, a total of 10 different conformations with different binding energies were generated. The conformation that exhibited the lowest energy value and most stable conformation (obtained from the clustering analysis) was identified as the best binding mode in each docking study. The 3D model of the docking results was visualized using PyMOL [29].

COSMO-RS Modelling
COnductor like Screening MOdel for Real Solvents (COSMO-RS) is the integration of the COSMO model used to predict thermos-physico-chemical properties of fluids in its pure and mixture state. Details on the calculation and procedure to estimate the excess enthalpy were obtained from the literature [30]. Initially, the nucleic acid bases (A, C, G, T, and U) were geometry optimized using the TmoleX package (Version 4.1.1). Then, COSMO-RS calculated the thermodynamic properties of fluid mixtures based on the 3D polarized charge distribution (σ, sigma) on the molecular surface of the individual bases, obtained from the quantum chemical calculations [31]. A 3D screening of the charge distribution on the molecular surface can be easily visualized on the histogram σ-profile. The information was useful for a qualitative description of the molecule and anticipated possible interactions of the compounds in a fluid mixture. Conversely, ILs were modeled as an equimolar mixture of cation and anion, considering that they were composed of 2 distinct ions, rather than of a single ionic molecule that conceptually resembled the actual situation. The thermodynamic calculation of the molecular interaction in this work was conducted using COSMOtherm software using the BP_TZVP_C30_1301 parameter file [24,25,31].

Steady State Fluorescence
fSteady state fluorescence spectroscopy is an effective technique to study the interaction of macromolecules and small compounds. In this work, the fluorescence spectroscopic method was applied to the imidazolium-based ILs as the compounds are known for non-negligible absorption, as described by Paul et al. [32] who observed the emission intensities of imidazolium-based ILs. They discovered fascinating fluorescence behaviors by the ILs at wavelengths beyond 300 nm, in which the [C n bim]Br was excitation wavelength-dependent, exhibiting the highest peak at an excitation wavelength of 320 nm.
The type of binding can be identified by observing changes in the spectrum. According to Nakamoto et al. [33], the formation of a covalent bond tended to cause a red shift and hyperchromism, whereas intercalative binding contributed to bathochromism and hypochromism [34]. Moreover, electrostatic interaction was characterized by the lower hypochromicity without a bathochromic shift. Figure 1 illustrates the emission intensity of fixed concentration of 0.5 M [C 2 bim]Br during continuous addition of ctDNA. Most importantly, elevation in fluorescence intensity suggested that the ILs were capable of binding with the DNA. calculation of the molecular interaction in this work was conducted using COSMOtherm software using the BP_TZVP_C30_1301 parameter file [24,25,31].

Steady State Fluorescence
fSteady state fluorescence spectroscopy is an effective technique to study the interaction of macromolecules and small compounds. In this work, the fluorescence spectroscopic method was applied to the imidazolium-based ILs as the compounds are known for non-negligible absorption, as described by Paul et al. [32] who observed the emission intensities of imidazolium-based ILs. They discovered fascinating fluorescence behaviors by the ILs at wavelengths beyond 300 nm, in which the [Cnbim]Br was excitation wavelength-dependent, exhibiting the highest peak at an excitation wavelength of 320 nm.
The type of binding can be identified by observing changes in the spectrum. According to Nakamoto et al. [33], the formation of a covalent bond tended to cause a red shift and hyperchromism, whereas intercalative binding contributed to bathochromism and hypochromism [34]. Moreover, electrostatic interaction was characterized by the lower hypochromicity without a bathochromic shift.  The study predicted that the smaller sized cation molecules with shorter hydrocarbon chain lengths could enter the DNA grooves to create a more hydrophobic microenvironment between the base pairs. Figure 2 represents the DNA bases showing numerous recognition sites of C-G and A-T pairs. Adenine of the C-6 amine acts as a hydrogen bond donor, while Adenine N-7 and thymine C-4 carbonyl in the A-T base pair act as hydrogen bond acceptors. In contrast, the thymine methyl group acts as a stabilizer for van der Waals interaction. In the C-G base pair, cytosine C-4 amine can function as a hydrogen bond donor, whereas both guanine C-6 carbonyl and guanine N-7 can act as hydrogen bond acceptors. Non-covalent interactions are usually formed between a small molecule, for instance, the rod-shaped imidazolium cation, and DNA molecules at the DNA bases boundary and the grooves. Since the fluorescence emission is very sensitive, any changes in the local environment of ILs can invoke changes in the intensity of the signal. In this perspective, the increasing intensity of The study predicted that the smaller sized cation molecules with shorter hydrocarbon chain lengths could enter the DNA grooves to create a more hydrophobic microenvironment between the base pairs. Figure 2 represents the DNA bases showing numerous recognition sites of C-G and A-T pairs. Adenine of the C-6 amine acts as a hydrogen bond donor, while Adenine N-7 and thymine C-4 carbonyl in the A-T base pair act as hydrogen bond acceptors. In contrast, the thymine methyl group acts as a stabilizer for van der Waals interaction. In the C-G base pair, cytosine C-4 amine can function as a hydrogen bond donor, whereas both guanine C-6 carbonyl and guanine N-7 can act as hydrogen bond acceptors. Non-covalent interactions are usually formed between a small molecule, for instance, the rod-shaped imidazolium cation, and DNA molecules at the DNA bases boundary and the grooves. Since the fluorescence emission is very sensitive, any changes in the local environment of ILs can invoke changes in the intensity of the signal. In this perspective, the increasing intensity of the fluorescence signifies changes in the local environment of the imidazolium cation from being well-solvated in aqueous solution to a hydrophobic surrounding.

Docking Studies
Docking is one of the well-founded methods to observe molecular level interactions between biomolecules and ligands. In this work, the binding energy was predicted to result from the interaction between the different alkyl chain length of imidazolium cation ([Cnbim] + where n = 2, 4, and 6) and ctDNA. The searching grid space was set to cover the whole duplex DNA in order to get reliable data. By referring to Figure 3, the docking study analysis revealed that the cation of ILs preferably binds to the minor groove of DNA. Pertinently, the results of the docking study and fluorescence spectra agreed well with the observed increase in emission intensity of ILs.   Table 1 shows the tabulated data for the binding energy of docked structures between cations and DNA. Based on the predicted results, the change in alkyl chain lengths clearly influenced the

Docking Studies
Docking is one of the well-founded methods to observe molecular level interactions between biomolecules and ligands. In this work, the binding energy was predicted to result from the interaction between the different alkyl chain length of imidazolium cation ([C n bim] + where n = 2, 4, and 6) and ctDNA. The searching grid space was set to cover the whole duplex DNA in order to get reliable data. By referring to Figure 3, the docking study analysis revealed that the cation of ILs preferably binds to the minor groove of DNA. Pertinently, the results of the docking study and fluorescence spectra agreed well with the observed increase in emission intensity of ILs.

Docking Studies
Docking is one of the well-founded methods to observe molecular level interactions between biomolecules and ligands. In this work, the binding energy was predicted to result from the interaction between the different alkyl chain length of imidazolium cation ([Cnbim] + where n = 2, 4, and 6) and ctDNA. The searching grid space was set to cover the whole duplex DNA in order to get reliable data. By referring to Figure 3, the docking study analysis revealed that the cation of ILs preferably binds to the minor groove of DNA. Pertinently, the results of the docking study and fluorescence spectra agreed well with the observed increase in emission intensity of ILs.   Table 1 shows the tabulated data for the binding energy of docked structures between cations and DNA. Based on the predicted results, the change in alkyl chain lengths clearly influenced the  Table 1 shows the tabulated data for the binding energy of docked structures between cations and DNA. Based on the predicted results, the change in alkyl chain lengths clearly influenced the binding energy (∆G) between the ILs and ctDNA. The value of binding energy (∆G) decreased as the length of the alkyl chain was increased. [C 6 bim] + had the lowest energy (−7.75 kcal mol −1 ), which corresponded to a binding constant of 3.44 × 10 3 L mol −1 , with van der Waals and desolvation energies being the main sources of energy. In this scenario, the electrostatic energy was insignificant for the interaction. It is important to note that the desolvation energy in this study referred to the hydrophobic energy that resulted from the hydrophobic interaction between the hydrocarbon length of the cation and non-polar sites of the bases inside the groove. An earlier work that assessed the [Br]) to DNA also described a similar observation [35]. They reported that total energies of both van der Waals and desolvation were more negative than the electrostatic energy.  [33]. By comparing to Haque et al. [36], their estimated binding energy of the [Ch][Cl]-DNA complex was −3.86 kcal mol −1 , which was two times lower as compared to ours. Possible reasoning for this discrepancy was the existence of a planar geometry in the imidazolium ring that offered multiple sites of interaction in the minor grooves. In addition, the non-bonded cation-π interaction can be established from the charge delocalization among polar atoms in the planar group [37]. In our case, such an interaction occurred between the oxygen atom in cytosine that faced the aromatic imidazole cation (Figure 3b). Conversely, the binding between the cholinium-based ILs and DNA was found to be weak-largely due to the steric hindrance of the bulky choline cation. Similarly, morpholinium-based ILs exhibited slightly lower binding energy than imidazolium-based ILs as the ring of morpholinium cation was non-planar and large.
Based on the stable docked complex in Figure 3, the cation had a higher tendency to bind to A-T rich regions of residues DA19-DT6 and DT18-DA7. The binding occurred between the polar hydrogen atoms (C-H) in the imidazolium ring and oxygen atom (O4) at the backbone of DNA. Based on previous work [37], such an interaction was the consequence of the electrostatic attraction between the ring of imidazole and the backbone of DNA. Data of the docking study demonstrated that ILs with similar structures to the [C n bim] + cations tended to have a similar binding pattern. Notably, the addition of a hydrocarbon chain in the imidazolium cation could favorably improve the binding of ILs to DNA. On the contrary, the employment of increasingly longer alky chains reduced the binding energy, thereby inferring the formation of stronger hydrophobic interactions. The findings in this study are consistent with an earlier MD simulation that described a duplex DNA in hydrated imidazolium-based ILs that was stabilized by an electrostatic attraction and groove binding mechanism [21].

Sigma Profiles of the Individual Compounds (Nucleic Acid Bases and ILs Molecules)
COSMO-RS was a reliable tool due to its ability to produces results with good agreement to the experimental data. COSMO-RS was utilized in this study to estimate the interaction energy in the ternary blend of water + IL with nucleic acid bases in a mixture state. The COSMO-RS theory prescribed that the σ-profile of molecules ranged from −0.03 to 0.03 e/Å2, which was the probability of the relative amount of surface segments with polarity, σ, on the surface of the molecule. The COSMO-RS histogram of σ-profile can be qualitatively divided into three main regions according to the cut-off values: The hydrogen bond donor (HBD) (σ < −0.008 e/Å2), the hydrogen bond acceptor (HBA) (σ > +0.008 e/Å2), and the nonpolar regions (−0.008 < σ < +0.008 e/Å2). of the relative amount of surface segments with polarity, σ, on the surface of the molecule. The COSMO-RS histogram of σ-profile can be qualitatively divided into three main regions according to the cut-off values: The hydrogen bond donor (HBD) (σ < −0.008 e/Å2), the hydrogen bond acceptor (HBA) (σ > +0.008 e/Å2), and the nonpolar regions (−0.008 < σ < +0.008 e/Å2).   Based on the σ-profiles (Figure 4), all nucleic acid bases exhibited peaks in both HBD and HBA regions. These peaks can be interpreted as the locations where H-bonds may be formed between the base pairs, which resulted in the double helix structure. The peaks in the HBD region originated from polar H atoms attached to the N or O atom in the DNA structure (light blue colored sigma surface). Whereas, peaks in the HBA region identify polar N or O atoms in the structure of the bases (red colored sigma surface). The aromatic structure of the bases produced broad peaks in the nonpolar region (green colored sigma surface).
The main peak in the nonpolar region corresponding to charge distribution of −0.003 e/Å2 represented the non-polar attribute that originated from the aliphatic alkyl chain and the aromatic head group (green colored sigma surface) ( Figure 5). Peak with charge distribution of −0.013 e/Å2 in the HBD region represented the acidic hydrogen atoms in the aromatic ring (light blue colored sigma COSMO-RS histogram of σ-profile can be qualitatively divided into three main regions according to the cut-off values: The hydrogen bond donor (HBD) (σ < −0.008 e/Å2), the hydrogen bond acceptor (HBA) (σ > +0.008 e/Å2), and the nonpolar regions (−0.008 < σ < +0.008 e/Å2).   Based on the σ-profiles (Figure 4), all nucleic acid bases exhibited peaks in both HBD and HBA regions. These peaks can be interpreted as the locations where H-bonds may be formed between the base pairs, which resulted in the double helix structure. The peaks in the HBD region originated from polar H atoms attached to the N or O atom in the DNA structure (light blue colored sigma surface). Whereas, peaks in the HBA region identify polar N or O atoms in the structure of the bases (red colored sigma surface). The aromatic structure of the bases produced broad peaks in the nonpolar region (green colored sigma surface).
The main peak in the nonpolar region corresponding to charge distribution of −0.003 e/Å2 represented the non-polar attribute that originated from the aliphatic alkyl chain and the aromatic head group (green colored sigma surface) ( Figure 5). Peak with charge distribution of −0.013 e/Å2 in the HBD region represented the acidic hydrogen atoms in the aromatic ring (light blue colored sigma Based on the σ-profiles (Figure 4), all nucleic acid bases exhibited peaks in both HBD and HBA regions. These peaks can be interpreted as the locations where H-bonds may be formed between the base pairs, which resulted in the double helix structure. The peaks in the HBD region originated from polar H atoms attached to the N or O atom in the DNA structure (light blue colored sigma surface). Whereas, peaks in the HBA region identify polar N or O atoms in the structure of the bases (red colored sigma surface). The aromatic structure of the bases produced broad peaks in the nonpolar region (green colored sigma surface).
The main peak in the nonpolar region corresponding to charge distribution of −0.003 e/Å2 represented the non-polar attribute that originated from the aliphatic alkyl chain and the aromatic head group (green colored sigma surface) ( Figure 5). Peak with charge distribution of −0.013 e/Å2 in the HBD region represented the acidic hydrogen atoms in the aromatic ring (light blue colored sigma surface). The increasing peak intensity in the nonpolar region denoted the increasing lengths of the aliphatic alkyl chain. In contrast, the single dominant peak in the HBA region signified the [Br] − anion at +0.018 e/Å2 (red shaded sigma surface) ( Figure 5).

Prediction of Excess Enthalpies and Evaluation of Molecular Interactions in Ternary Mixture of Water + ILs with DNA Bases Systems Using COSMO-RS
Excess enthalpy, H E m refers to the change of enthalpy following the introduction of DNA into the mixture of water and ILs. This value is useful to predict the strength of the interaction between nucleic acid bases and ILs in the mixture. The predicted H E m whole composition for the water and IL mixtures is illustrated in Figures 6 and A1 (Appendix A), respectively. Results showed that ILs bind significantly to various bases, the negative H E m value indicated a favorable interaction between aqueous IL with free nucleic acid bases. A similar trend was observed in the aqueous mixture of [Ch]Cl + DNA bases [36]. The order of interaction for all types of ILs and the bases were ranked as follows: Guanine > cytosine > thymine > uracil > adenine, thus consistent with findings of another similar study [37]. A similar order of interaction was also obtained based on the excess Gibbs energy as shown in Figures 7 and A2 in Appendix A. surface). The increasing peak intensity in the nonpolar region denoted the increasing lengths of the aliphatic alkyl chain. In contrast, the single dominant peak in the HBA region signified the [Br] -anion at +0.018 e/Å2 (red shaded sigma surface) ( Figure 5).

Prediction of Excess Enthalpies and Evaluation of Molecular Interactions in Ternary Mixture of Water + ILs with DNA Bases Systems Using COSMO-RS
Excess enthalpy, refers to the change of enthalpy following the introduction of DNA into the mixture of water and ILs. This value is useful to predict the strength of the interaction between nucleic acid bases and ILs in the mixture. The predicted whole composition for the water and IL mixtures is illustrated in Figure 6 and Figure A1 (Appendix A), respectively. Results showed that ILs bind significantly to various bases, the negative value indicated a favorable interaction between aqueous IL with free nucleic acid bases. A similar trend was observed in the aqueous mixture of [Ch]Cl + DNA bases [36]. The order of interaction for all types of ILs and the bases were ranked as follows: Guanine > cytosine > thymine > uracil > adenine, thus consistent with findings of another similar study [37]. A similar order of interaction was also obtained based on the excess Gibbs energy as shown in Figure 7 and Figure A2 in Appendix A.  surface). The increasing peak intensity in the nonpolar region denoted the increasing lengths of the aliphatic alkyl chain. In contrast, the single dominant peak in the HBA region signified the [Br] -anion at +0.018 e/Å2 (red shaded sigma surface) ( Figure 5).

Prediction of Excess Enthalpies and Evaluation of Molecular Interactions in Ternary Mixture of Water + ILs with DNA Bases Systems Using COSMO-RS
Excess enthalpy, refers to the change of enthalpy following the introduction of DNA into the mixture of water and ILs. This value is useful to predict the strength of the interaction between nucleic acid bases and ILs in the mixture. The predicted whole composition for the water and IL mixtures is illustrated in Figure 6 and Figure A1 (Appendix A), respectively. Results showed that ILs bind significantly to various bases, the negative value indicated a favorable interaction between aqueous IL with free nucleic acid bases. A similar trend was observed in the aqueous mixture of [Ch]Cl + DNA bases [36]. The order of interaction for all types of ILs and the bases were ranked as follows: Guanine > cytosine > thymine > uracil > adenine, thus consistent with findings of another similar study [37]. A similar order of interaction was also obtained based on the excess Gibbs energy as shown in Figure 7 and Figure A2 in Appendix A.    (1) Figure 7. Calculated total excess Gibb's Free energy, ,, of water + [C2bim]Br + nucleic acid bases at 298.15 K by COSMO-RS. Figure 8 shows the special contribution for aqueous [C2bim]Br with the guanine at 298.15 K at = 0.8 (for other studied ILs see Appendix A Figures A3-A5). Using COSMO-RS model, specific contribution to the can be estimated. is the sum of specific molecular interaction, which represents the electrostatic-misfit ( , ), hydrogen bonds ( , ) and van der Waals forces ( , ). The total contribution of each specific interaction can be summarized by the following equation:   Figure 10. Based on the COSMO-RS model, the most favorable interaction between aqueous ILs + nucleic acid bases mixtures is aqueous [C2bim]Br + guanine. The more negative the excess Gibbs energy, the more attractive the interaction between the ILs and the nucleic acid bases.

Conclusions
The study performed a systematic assessment of the mode and strength of binding between calf thymus DNA and imidazolium-based ILs. Data on molecular docking revealed that ILs and calf thymus DNA formed a weak interaction with ∆G in the range of −5.19 to −7.75 kcal mol −1 . Aside from electrostatic attraction, the occurrence of T-shape π-π stacking between the bases and the aromatic ring of imidazole was also observed in the stable docking conformation. Despite that ILs interacted weakly with calf thymus DNA, the preparations exhibited potential affinity towards individual nuclei acid bases, as predicted by COSMO-RS. Most importantly, data of the molecular docking study corroborated the results of the fluorescence experiments. The study effectively demonstrated the mode of binding of imidazolium-based ILs to bases in which the binding preference increases from adenine, uracil, thymine, cytosine, and guanine. The findings of this study also emphasize the role of [Cnbim]Br ILs cation absorption for long term/high-temperature stabilization of DNA material. The present work demonstrated the role of the developed [Cnbim]Br ILs as solvents to stabilize DNA stored at room temperature. It was made apparent by the preferable binding of smaller cations of [Cnbim]Br ILs to the minor groove of ct-DNA, causing a formidable

Conclusions
The study performed a systematic assessment of the mode and strength of binding between calf thymus DNA and imidazolium-based ILs. Data on molecular docking revealed that ILs and calf thymus DNA formed a weak interaction with ∆G in the range of −5.19 to −7.75 kcal mol −1 . Aside from electrostatic attraction, the occurrence of T-shape π-π stacking between the bases and the aromatic ring of imidazole was also observed in the stable docking conformation. Despite that ILs interacted weakly with calf thymus DNA, the preparations exhibited potential affinity towards individual nuclei acid bases, as predicted by COSMO-RS. Most importantly, data of the molecular docking study corroborated the results of the fluorescence experiments. The study effectively demonstrated the mode of binding of imidazolium-based ILs to bases in which the binding preference increases from adenine, uracil, thymine, cytosine, and guanine. The findings of this study also emphasize the role of [C n bim]Br ILs cation absorption for long term/high-temperature stabilization of DNA material. The present work demonstrated the role of the developed [C n bim]Br ILs as solvents to stabilize DNA stored at room temperature. It was made apparent by the preferable binding of smaller cations of [C n bim]Br ILs to the minor groove of ct-DNA, causing a formidable change in the microenvironment that to be more hydrophobic to expel water from immediate contact with DNA. This dehydration effect can extend the lifetime of ct-DNA by circumventing a hydrolytic attack that instigates rapid deamination and depurination while increasing the thermal stability of the base pairs. Notwithstanding safe storage and extending the lifetime of genetic blueprints, it is envisioned that the applicabilities of [C n bim]Br ILs can be extended to the manipulation of nucleic acids and the fine-tuning of DNA microenvironment for compaction and possible applications in nanotechnology.