Biosolvation Nature of Ionic Liquids: Molecular Dynamics Simulation of Methylated Nucleobases in Hydrated 1-Ethyl-3-methylimidazolium Acetate

Toward Greater DNA Stability by Leveraging the Proton-Donating Ability of Protic Ionic Liquids

Deoxyribonucleic acid (DNA) stability is a prerequisite in many applications, ranging from DNA-based vaccines and data storage to gene therapy. However, the strategies to enhance DNA stability are limited, and the underlying mechanisms are poorly understood. Ionic liquids (ILs), molten salts of organic cations and organic/inorganic anions, are showing tremendous prospects in myriads of applications. With a judicious choice of constituent ions, the protic nature of ILs can be tuned. In this work, we investigate the relative stability of full-length genomic DNA in aqueous IL solutions of increasing protic nature. Our experimental measurements show that the protic ionic liquids (PILs) enhance the DNA melting temperature significantly while unaltering its native B-conformation. Molecular dynamics simulations and quantum mechanical calculation results suggest that the intramolecular Watson-Crick H-bonding in DNA remains unaffected and, in addition, the PILs induce stronger H-bonding networks in solution through their ability to make multiple intermolecular H-bonds with the nucleobases and among its constituent ions, thus aiding greater DNA stability. The detailed understanding obtained from this study could bring about the much-awaited breakthrough in improved DNA stability for its sustained use in the aforesaid applications!

Ionic liquids: prospects for nucleic acid handling and delivery

Operations with nucleic acids are among the main means of studying the mechanisms of gene function and developing novel methods of molecular medicine and gene therapy. These endeavours usually imply the necessity of nucleic acid storage and delivery into eukaryotic cells. In spite of diversity of the existing dedicated techniques, all of them have their limitations. Thus, a recent notion of using ionic liquids in manipulations of nucleic acids has been attracting significant attention lately. Due to their unique physicochemical properties, in particular, their micro-structuring impact and tunability, ionic liquids are currently applied as solvents and stabilizing media in chemical synthesis, electrochemistry, biotechnology, and other areas. Here, we review the current knowledge on interactions between nucleic acids and ionic liquids and discuss potential advantages of applying the latter in delivery of the former into eukaryotic cells.

Conformational Free-Energy Landscapes of Alanine Dipeptide in Hydrated Ionic Liquids from Enhanced Sampling Methods

Understanding the interaction of the ionic liquid (IL) with protein is vital to find the origin of the conformational changes of proteins in these alternative solvents. Here, we performed biased molecular dynamics simulations of alanine dipeptide (ADP), a widely used model for protein backbone structure, in water and two hydrated ionic liquids (ILs): 80% (w/w) 1-ethyl-3-methylimidazolium acetate ([EMIm][Ac]) and 80% (w/w) choline dihydrogen phosphate ([Cho][DHP]). We employed three different biasing methods, metadynamics (metaD), well-tempered metadynamics (WT-metaD), and adaptive biasing force (ABF), to construct the free-energy landscapes of the ADP conformations using the backbone dihedral angles (ϕ and ψ) as the collective variables. The calculations were also performed in water; the free-energy landscapes of ADP in water obtained from three methods are similar and agree well with the previously reported results. In hydrated [EMIm][Ac], α-planar conformation emerges as a minimum, which is comparable to that of α and β conformations corresponding to α-helix and β-sheet-like conformations of proteins. Investigation of corresponding conformations suggests that the imidazolium ring of [EMIm] cation is stacked with the amide bonds of ADP. Acetate anion makes hydrogen bonds with the amide hydrogens of the ADP. The amide-π stacking interaction is the driving force for α-planar conformation to become one of the minimum energy conformations in this IL, which destabilizes the protein conformation. However, α and β conformations are more stable in hydrated [Cho][DHP] compared to α-planar and β-planar conformations; therefore, this IL stabilizes the protein conformation. These findings are in good correlation with the previous study of proteins in these ILs. Our study helps to understand the interaction of proteins with the ionic entities and their stability in ILs.

Conformational dynamics of amyloid-β (16–22) peptide in aqueous ionic liquids

Molecular dynamics simulations of amyloid-β (16–22) peptide dimer in water as well as at two different experimentally studied concentrations of hydrated ionic liquids (ILs), ethylammonium mesylate (EAM), ethylammonium nitrate (EAN), and triethylammonium mesylate (TEAM), were carried out employing an umbrella sampling method. We used the average Ψ angle of the peptide backbone as the reaction coordinate to observe the conformational changes of a peptide dimer. Secondary structural element values were calculated for the peptide dimer along the reaction coordinate to see the transition of the peptide dimer between β-sheet and α-helix conformations. We observe the β-sheet conformation as the global minimum on the free energy surfaces in both EAM and EAN ILs at both the concentrations and at a low concentration of TEAM. However, we observe α-helix conformation as the global minimum at a high concentration of TEAM. Our results are in good correlation with the experimental findings. We calculated the average number of intramolecular and intermolecular hydrogen bonds of α-helix and β-sheet conformations in all solutions, and they are in correlation with the secondary structure element values. To understand the peptide–IL interactions, atom–atom radial distribution functions of cation, anion, and water around amide oxygen and hydrogen atoms were calculated. The solvent-accessible surface area of the peptide dimer was calculated to understand the exposure of the peptide towards the solvent during conformational changes. Finally, van der Waals (vdW) and Coulomb interaction energies were calculated between peptide–cation, peptide–anion, and peptide–water to understand the stability of conformations in different concentrations. We find that the TEA cation has more vdW interaction energy compared to Coulomb interaction energy with peptide in 70% (w/w) TEAM, which mimics a membrane-like environment to induce α-helix conformation rather than β-sheet conformation.

Molecular dynamics simulations of amyloid-β (16–22) peptide dimer at two different experimentally studied concentrations of hydrated ethylammonium mesylate, ethylammonium nitrate, and triethylammonium mesylate were carried out employing an umbrella sampling method.

Association of Nucleobases in Hydrated Ionic Liquid from Biased Molecular Dynamics Simulations

We employed metadynamics-based classical molecular dynamics simulations to methylated adenine-thymine (mA-mT) and guanine-cytosine (mG-mC) base pairs to see favorable conformations in various concentrations of hydrated 1-ethyl, 3-methyl imidazolium acetate. We investigated various stacked and hydrogen-bonded conformations of association of base pairs through appropriately chosen collective variables. Stacked conformations more favored in water for both base pairs, whereas Watson-Crick (WC) hydrogen-bonding conformations are favored in pure and hydrated ionic liquids (ILs) except for 0.75 mol fraction IL. We observe that EMIm cations surround the base pairs in WC conformations creating a kind of hydrophobic cavity and protect the hydrogen bonds between base pairs. However, the five-membered heteroaromatic rings of cations stack with the nucleobases in the cation-base-cation (π-π-π) model, which resembles the base-base-base stacking in a DNA duplex. Interestingly, from additional simulations of 0.5 mol fraction hydrated choline dihydrogen phosphate IL, we observe that the stacked conformations become more favored than the WC conformation due to the absence of π-bonds in cations. The calculated values of relative solubility of base pairs in pure and hydrated ionic liquids compared to those in pure water correlate well with the free energy values of WC and stacked conformations.