1. Introduction

Computational Chemistry

2332-5968

Scientific Research Publishing

10.4236/cc.2017.52007

CC-76024

Articles

Chemistry&Materials Science

Theoretical Investigation of the Newly Designed Benzimidazole Based Metal Mediated DNA Base Couples with DFT Method

Fatma

Sevin

¹^*Mostafa

Asghari Dilmani

¹Kübra

Sarıkavak

¹Samira

Farshbaf Jeddi

Chemistry Department, Faculty of Science, Hacettepe University, Ankara, Turkey

* E-mail:sevin@hacettepe.edu.tr(FS);

28042017

05027490February 7, 2017Accepted: April 27, April 30, 2017

2014

This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

DNA has the genetic information storage and transmission capacity according to the sequential order of the monomer and creates a central role in the chemical evolution by copying itself with the proliferation feature. Watson-Crick base pairs define two base pairs with hydrogen bonds. If metal coordination bonds replace hydrogen bonds more stable alternative metallo-DNA sequence can be established. If the replication feature can be obtained for the metallo-DNA, this will greatly benefit the creation of DNA computer keys. In this study, a new type of benzimidazole based metallo-DNA sensors consisting of a connector unit that unsaturated azinil bridge linked to Watson-Crick base pairs with Ni ²⁺, Hg ²⁺, Zn ²⁺, Ag ⁺, Pt ²⁺, Pd ²⁺ metal cations and a benzimidazole has been designed. Absorption and emission spectrum of the newly designed aqua medium based fluorophore and their metallo-DNA sensors with selected cations have been theoretically investigated by using DFT method. The logic gates of selected possible sensors which response in the visible region have also been examined in detail in acidic and water phase. As a result of calculated absorption-emission spectrum data show that T-Hg-A-Bnz, A-Ni-T-Bnz, C-Pt-G-Bnz, C-Ni-C-Bnz complexes produce OR gate. T-Zn-T-Bnz and G-Pt-C-Bnz results demonstrated XOR and AND logic gate, respectively.

Flourescence Metallo-DNA Sensor Benzimidazole Logic Gates DFT

1. Introduction

Today, development of biomolecular structures is generally based on supra- molecules that include non-covalent interactions, such as hydrogen bonds, hydrophobic effects and metal coordination bonds [1] [2] [3] [4] . Without a doubt, the most basic structures in origins of the chemical evolution are the nucleic acids. The most basic indicators that allow storing, transferring and copying of the genetic information within these nuclear acids are hydrogen bonded Watson-Crick base pairs. Having metal coordination bonds instead of hydrogen bonds presents alternative and completely different base pairs [5] [6] [7] . As known, the event of binding a metal ion to a molecule affects the characteristics of UV-visible spectrum of molecule. Fluorescence technics are also good tool to identify biomolecular interactions and uses widely in many research. For these technics, some fluorophores have been used to analyze cations and anions that can be found in the nature. Among various fluorophores, benzimidazole has drawn attention due to its optical properties and high stability. Benzimidazole ring is a very good type of fluorophore as it can produce “scorpion type” complexes [8] .

Theoretically suggested one dimensional TMn (benzimidazole)_n+1 (TM = Sc, Ti, V, Cr, Mn) system’s electronic and magnetic characteristics were investigated with DFT method [9] . It was detected that benzimidazole can be stable while also retaining the helix and one-dimensional characteristics of the DNA with the Ti, V and Cr transition atoms. Experimentally [10] and theoretically [11] studies were established about structural-electronic characteristics and UV absorption spectrums of the Thymine-Hg²⁺-Thymine base pair (T-Hg-T), which is a mercury (II) linked metallo-DNA complex. In the experimental study, the stability of the thymine pair with the mercury ion was determined in various temperatures and pH’s and it was seen that its maximum absorption was at 260 nm. In the theoretic study, the base pairs and dimers that thymine and its derivatives of cis and trans form with the mercury ion were calculated with the TD-B3LYP and TD-PCM-B3LYP methods, it was determined that maximum absorption occurs in values that are closer to red, which is to say 263 nm in solvent phase and 276 nm in gas phase.

Takezawa and Shionoya presented an abstract about the chemistry of metal-linked base pairs which includes primary approaches to the DNA based molecular systems, molecular designs,structures, characteristics and their applications in their research [12] . One of the examples that were provided Cu²⁺ linked hydroxipyron base pair. The H-Cu²⁺-H base pair has a square planar characteristic while H defines a nucleoacid that carries a hydroxipyron. The most distinctive characteristic of the metallo-base pairs added to the DNA is their potential to increase the duplex stability. As expected, metal coordination bond energies are two or three times greater than the hydrogen bonds. For this reason, their stability effects are worked out from the melting temperatures of the DNA base pairs. An example for this is as the following: when there is a Cu²⁺ ion, this temperature has increased from 37˚ Celsius to 50˚ Celsius. Similar stabilities are also achieved in Salen-type. Structures like similar to these are synthesized and characterized as monomer structures at first and then they were added to, for example, 15-mer DNA duplexes. Also, Li et al. have examined the aggregation behavior of the silver molecules with the DFT method by using polymorphic DNA models that contain Watson-Crick base pairs, i-motif and G-quadrulex [13] . Leutwyler group examined double hydrogen bonded complexes of 2-pyridones cytosine and 1-methyl cytosine with mass, UV and IR spectroscopic techniques and observed five different cytosine 2-pyridone isomers in their theoretical and experimental study [14] .

Usage of sensors as logic gates in biochemical researches has begun only recently, but it has started developing [15] [16] . In chemical logic systems, binding a guest molecule to the host molecule corresponds to a logical input and the result includes physical changes as an output which corresponds to absorption or fluorescence spectrum. Whenever multiple chemical inputs provide a single output independent of each other, the system is defined as OR logic gates, which is to say that this is a very weak chemical selection system. On the other hand, AND logic gate identifies multiple chemical inputs based on luck and it provides an input that requires high chemical selectivity.

In this study, metallo-DNA sensors that have new fluorescent characteristics and the capability to work in aqueous mediums have been designed by binding metals (Hg²⁺, Ag⁺, Ni²⁺, Pb²⁺, Pt²⁺, Zn²⁺) that can provide them coordination, especially planar coordination, to Watson-Crick base pairs and their reversible and changeable optical characteristics in acidic mediums have been investigated and possible logic gates have been suggested. As it can be seen on Figure 1 and Figure 2, a new type of benzimidazole based metallo-DNA base pair sensors have been designed which consisting of a connector unit and unsaturated azinil bridge linked to Watson-Crick (T = Thymine, A = Adenine, C = Cytosine and G = Guanine) base pairs bonded with Ni²⁺, Hg²⁺, Zn²⁺, Ag⁺, Pt²⁺, Pd²⁺ metal cations and a fluorophore. These newly designed metallo-sensors have the characteristic that allows molecular identification with visible changes in color (colorimeter) and increases and decreases in emission wavelength (fluorescence) due to their coordination with various cations.

To explain the structural and electronic characteristics of these sensors in different media, their energies, absorption and emission spectrums, energy

Figure 1 Designed benzimidazole based metallo-DNA based pair sensorsFigure 2 Some azinil benzimidazole mediated tymine base pair examples

differences between frontier molecular orbitals (HOMO: highest occupied molecular orbital and LUMO: lowest unoccupied molecular orbital) which means HOMO-LUMO gap has been calculated as well as color and emission changes. Same calculations have been made with protonation of sp2 hybrid nitrogen atom of benzimidazole portions and logic gates have been presented for acidic medium and aqueous phase. As a result, since these designed sensors are expected to be potential keys for the DNA computer technology, we expect them to contribute greatly to science and technology applications.

2. Method

In this study, all the calculations have been performed with Gaussian 09W [17] and GaussView 5.0.8 molecular modeling software [18] . The methods used in the calculation of the organometallic compounds with mercury in the literature have been chosen [19] [20] [21] . It is especially known, M06 is one of the best functional for the study of organometallic and inorganomatellic thermochemistry and noncovalent interactions. In the frame of these informationsM06 (Hybrid meta exchange-correlation functionals) [22] , B3LYP (Becke, three-parameter, Lee-Yang-Parr) [23] , PBE0 (Perdew, Burke and Ernzerhof) [24] methods and the double-zeta pseudo-potential LandL2DZ [25] basis set has been selected for geometry optimization and evaluate the absorption and emission spectroscopy. The selection of appropriate method among the aforecited methods has been based on the experimental results of T-Hg²⁺-T base pair which has a known maximum absorption value in the literature as 260 nm. These calculated absorption results have been compared with the corresponding experimental values. The calculation results show that M06 functional has been found better than the other methods as in Table S1 (it is given in Supporting info). After determining the appropriate method, all subsequent calculations performed with M06 functional LandL2DZ basis set. Vibrational frequency analyses have been carried out to confirm local minima of the structures. In order to compute the solvation effect self-consistent reaction field (SCRF) theory with polarizable continuum method (PCM) used in the water phase calcultaions. The dielectric constant was chosen as the standard value for water, (ε = 78.39). Calculations corresponding to acidic medium implemented with nitrogen protonation of benzimizadole fragment. Absorption and emission spectra calculated and logic gates determined through these constructions.

The total energy of the structures (E), Gibbs free energy (G) and enthalpy (H) has been calculated. Along with gas phase calculation, Conductor-Like Screening Model [26] has been used to calculate theoretical absorption wavelength in aqueous media. Also, 7-digit TD-DFT method has been used for emission calculation [27] .

3. Results and Discussion

Investigation of the DNA bases and their structures with metal and benzimidazole has been divided into two categories as T-A, A-T, C-G and G-C base combinations. All calculations have been conducted in the water phase. Formation energies, frontier molecular orbital band gaps have been calculated and studies particularly have focused on the photophysical properties.

3.1. Properties of Designed T-T and T-A and T-A Base Pairs and Their Benzimidazole Base Pairs with and without Metal Cations

In the first stage, the complexation energies of T-T, T-A and A-T base pairs with metal cations have been calculated and obtained results given in Table 1. Also, total energies and entaphy have been demonstrated in Tables S2-S4. Complexation energies of benzimidazole (Bnz) based T-T,T-A and A-T base pairs with metal cations have been calculated and results given in the same table. As can be seen from the table, complexation of T-T occurs easily with Ni²⁺. This sequence followed by Hg²⁺. The complexation abilities of metal cation in this sequence for T-T is Ni²⁺ > Hg²⁺ > Zn²⁺ > Ag⁺ > Pt²⁺ > Pd²⁺. This sequence is the

Table 1 The complexation energies of T-T, T-A and A-T base pairs with metal cations and Bnz (Kcal/mol)

Base Pairs	Energy Values
Base Pairs	ΔG	Base Pairs		ΔG	Base Pairs
T-T → T-Hg-T T-T → T-Zn-T T-T → T-Ag-T T-T → T-Ni-T T-T → T-Pd-T T-T → T-Pt-T T-A-Bnz → T-Hg-A-Bnz T-A-Bnz → T-Zn-A-Bnz T-A-Bnz → T-Ag-A-Bnz T-A-Bnz → T-Ni-A-Bnz T-A-Bnz → T-Pd-A-Bnz T-A-Bnz → T-Pt-A-Bnz	−122.66 −98.06 −78.26 −183.55 −59.08 −71.98 96.05 −77.00 −64.22 −122.41 −45.89 −56.00		T-A → T-Hg-A T-A → T-Zn-A T-A → T-Ag-A T-A → T-Ni-A T-A → T-Pd-A T-A → T-Pt-A A-T-Bnz → A-Hg-T-Bnz A-T-Bnz → A-Zn-T-Bnz A-T-Bnz → A-Ag-T-Bnz A-T-Bnz → A-Ni-T-Bnz A-T-Bnz → A-Pd-T-Bnz A-T-Bnz → A-Pt-T-Bnz	-101.32 −81.46 −72.19 −134.85 −50.00 −61.35 −105.22 −84.92 −76.00 −104.55 −54.02 −67.26	T-T-Bnz → T-Hg-T-Bnz T-T-Bnz → T-Zn-T-Bnz T-T-Bnz → T-Ag-T-Bnz T-T-Bnz → T-Ni-T-Bnz T-T-Bnz → T-Pd-T-Bnz T-T-Bnz → T-Pt-T-Bnz	−112.08 −99.02 −79.37 −184.62 −60.81 −73.26

same with T-T for T-A base pair, too.

Complexation of benzimidazole with T-T base pair and its complexation energies with metal cations prefers the Ni²⁺ cation. Its sequence is as following: Ni²⁺ > Hg²⁺ > Zn²⁺ > Ag⁺ > Pt²⁺ > Pd²⁺. T-A-Bnz base pair shows the same trend. This sequence is the same as T-T and benzimidazole did not have an impact on adenine.

It is clearly observed from the results, benzimidazole based A-T designed base pair results show that when benzimidazole bonds with the thymine part of this base pair, the first two rows of the sequence slightly changes: Hg²⁺ > Ni²⁺ > Zn²⁺ > Ag⁺ > Pt²⁺ > Pd²⁺. It can be considered that benzimidazole has more effects on thymine. The result shows that, with the exception of the A-T-Bnz base pair which prefers Hg²⁺ cation, other benzimidazole based DNA base pairs primarily prefer Ni²⁺ cation.

The energy gap reflects the reactivity or stability of the molecule. HOMO- LUMO energy gap of molecules is considered as a measure of charge transfer and is regarded as an important parameter in determining the properties such as electrical conductivity. Table 2 shows that T-T-Bnz has made it stable than T-T about 0.35 eV. However, in case of complexation with metal cations, Pd²⁺ and Pt²⁺ is more thermodynamically stable while the gap energies of other metal cations decrease. This situation shows that conjugation is also increased on the structure. The highest amount conjugation occurs with Ni²⁺ and Ag⁺among the selected cations. Hence, HUMO-LUMOs of T-Ag-T seen that the contribution of metal orbitals is greater in linked benzimidazole and conjugation shifts to benzimidazole ring. The significant effect has not been observed by linking benzimidazole for T-A base pairs. Complexation energies with Pd²⁺and Pt²⁺ metal cations are higher than the cations for T-A and also with Hg²⁺ and Zn²⁺ cations are decreased their gap energies. In terms of A-T base pairs, bonding of benzimidazole to structure increases molecular stability about 0.20 eV. Although Ag⁺ and Ni²⁺ cations more stable in case of complexation with metal cations, the energy gap values have decreased only for Hg²⁺ cation. Thus, base pair orbitals in T-Hg-T HOMO-LUMOs show benzimidazole orbitals in LUMO when benzimi-

Table 2 The HOMO-LUMO band gap energies of T-T-, T-A and their complexes with Bnz metal cations (eV)

Base Pairs	HOMO-LUMO gap energies
Base Pairs	-	Hg	Zn	Ag	Ni	Pd	Pt
T-T T-T-Bnz T-T/T-T-Bnz T-A T-A-Bnz T-A/T-A-Bnz A-T A-T-Bnz A-T/A-T-Bnz	1.178 1.552 −0.347 0.940 0.980 −0.040 0.940 1.139 −0.199	4.602 3.050 1.552 3.432 1.413 2.019 3.432 1.337 2.095	4.124 3.107 1.017 3.785 1.392 2.393 3.785 4.372 −0.587	4.736 3.009 1.727 3.230 4.565 −1.335 3.230 4.566 −1.336	4.649 3.013 1.636 3.263 4.763 −1.500 3.263 4.120 −0.857	0.824 0.912 −0.088 3.395 4.680 −1.285 3.395 3.517 −0.122	0.926 2.354 −1.428 3.671 4.884 −1.213 3.671 4.227 −0.556

Table 3 Maximum absorption/emission wavelength (nm) and differences between them of the designed T-T, T-A, A-T pairs and their Bnz complexes

Pairs	T-T- Based
Pairs	λ Abs		Osc		λ Ems	Osc	∆λ	Pairs λ Abs			Osc		λ Ems			Osc		∆λ
T-T T-Hg-T T-Zn-T T-Ag-T T-Ni-T T-Pd-T T-Pt-T	264 335 460 586 283 866 603	0.060 0.230 0.100 0.010 0.004 0.003 0.025			293 397 471 709 314 901 738	0.020 0.110 0.080 0.009 0.001 0.001 0.013	28 43 10 122 31 34 13	T-T-Bnz T-Hg-T-Bnz T-Zn-T-Bnz T-Ag-T-Bnz T-Ni-T-Bnz T-Pd-T-Bnz T-Pt-T-Bnz	319 358 468 741 331 967 807			0.620 0.560 0.009 0.090 0.046 0.001 0.078		337 407 482 808 359 1109 921	0.430 0.380 0.005 0.070 0.019 0.008 0.047		18 49 13 67 105 141 113
		T-A- Based								T-A-Bnz Based
T-A T-Hg-A T-Zn-A T-Ag-A T-Ni-A T-Pd-A T-Pt-A	301 473 602 645 343 864 767	0.030 0.047 0.051 0.088 0.027 0.002 0.700			339 509 668 861 457 985 789	0.028 0.044 0.040 0.082 0.023 0.001 0.650	38 35 65 216 114 121 21	T-A-Bnz T-Hg-A-Bnz T-Zn-A-Bnz T-Ag-A-Bnz T-Ni-A-Bnz T-Pd-A-Bnz T-Pt-A-Bnz	334 468 692 781 351 986 801			0.060 0.019 0.029 0.067 0.004 0.058 0.061		357 509 715 898 477 1108 827	0.051 0.016 0.017 0.056 0.003 0.054 0.540		22 41 23 116 126 121 26
		A-T- Based							A-T-Bnz Based
A-T A-Hg-T A-Zn-T A-Ag-T A-Ni-T A-Pd-T A-Pt-T	301 473 602 645 343 864 767	0.030 0.047 0.051 0.088 0.027 0.002 0.700		339. 509 668 861 457 985 789		0.028 0.044 0.040 0.082 0.023 0.001 0.650	38 35 65 216 114 121 21	A-T-Bnz A-Hg-T-Bnz A-Zn-T-Bnz A-Ag-T-Bnz A-Ni-T-Bnz A-Pd-T-Bnz A-Pt-T-Bnz	365 439 658 772 443 822 789			0.013 0.019 0.020 0.060 0.002 0.001 0.670		404 488 670 893 661 1043 839	0.012 0.014 0.018 0.056 0.002 0.001 0.590		38 48 12 121 217 221 40

dazole enters. According to calculations; Hg²⁺, Ag⁺ and Zn²⁺ cations can increase resonance stability of structure and Pt²⁺ and Pd²⁺ cations provides thermody- namic stability.

In the scope of this section, absorption and emission spectrum has been studied and results of the complexes formed by T-T, T-A, A-T bases and their benzimidazole based pairs with metals has been given in Table 3. This table shows that T-T-Bnz and T-M-T complexes have maximum absorption and emission wavelength (λ_max) with Zn²⁺ and Ag⁺ cations, respectively. For Ag⁺ cations, Stokes shift value are greater than Zn²⁺ cations. When the results are evaluated for T-A and T-A-Bnz, it has seen that maximum absorption and emission wavelength λ_max values belong to T-Hg-A, T-Zn-A and T-Zn-A?Bnz complexes. Δλ value of the T-A based complex with Ni²⁺ cations is higher than Hg²⁺ cations one. The results obtained by A-T based structures are the same with the obtained for T-A structure. Maximum absorption and emission wavelength for A-T-Bnz belong to the A-Zn-T-Bnz complex

In accordance with spectral data of T-T, T-A, A-T base pairs with selected metals and benzimidazole their colors have been determined before and after radiation. The color difference caused by the binding of metal cations to T-T and T-T-Bnz structures is given in the Table S5. The color change from blue to green has been observed in T-Hg-A, T-Zn-T-Bnz, A-Hg-T and T-Hg-A-Bnz complexes. The color change from violet to green has been seen in A-Hg-T-Bnz and A-Ni-T-Bnz pairs. Yellow coloured T-Ag-T and orange coloured T-Pt-T base pairs has losts their colours after radiation like red coloured T-Ag-A, T-Zn-A-Bnz and A-Ag-T pairs. T-Ni-A, T-Ni-A-Bnz and A-Ni-T complexes that are colourless initially,have gained blue colour after radiation. Colorless A-T-Bnz and T-Hg-T-Bnz pairs have gained violet colour, too. The color change from yellow to red has been observed for T-Ag-T and A-Zn-T pairs

3.2. Properties of Designed C-G and G-C Base Pairs and Their Benzimidazole Base Pairs with and without Metal Cations

As in the previous section, complexation energies, band gaps and spectral properties of the targeted structures has been studied in this part, too. Table 4 displays complexation energies of C-G and C-C base pairs and their combinations with benzimidazole and selected cations.

Calculations show that the complexation of C-G and C-C occurs easily with Ni²⁺ like T-T and A-T ba-ses. The metal cation sequence for C-G is Ni²⁺ > Hg²⁺ > Pt²⁺ > Pd²⁺ > Ag⁺ > Zn²⁺. The metal cation sequence for C-C is Ni²⁺ > Ag⁺ > Hg²⁺ > Zn²⁺ > Pd²⁺ > Pt²⁺. Unlike the others, this sequence is followed by Ag⁺. When the complexation energies analyzed for created by complexation of Bnz based C-G pairs with metal cations, it has been seen that the bonding of Bnz with guanine is the same as its bonding with G-C base pair and this has not changed the preference of it. The same results have been obtained for G-C-Bnz complex and C-C-Bnz complex. The results of all the complexation calculations

Table 4 The complexation energies of C-G -Bnz and G-C-Bnz base pairs with metal cations (Kcal/mol)

Base Pairs	Energy Values
Base Pairs	ΔG		Base Pairs	ΔG	Base Pairs
C-G → C-Hg-G C-G → C-Zn-G C-G → C-Ag-G C-G → C-Ni-G C-G → C-Pd-G C-G → C-Pt-G G-C-Bnz → G-Hg-CBnz G-C-Bnz → G-Zn-C-Bnz G-C-Bnz → G-Ag-C-Bnz G-C-Bnz → G-Ni-C-Bnz G-C-Bnz → G-Pd-C-Bnz G-C-Bnz → G-Pt-C-Bnz	−140.33 −63.55 −77.77 −154.64 −101.04 −113.36 −132.45 −69.01 −81.34 −150.09 −99.58 −110.78	C-C → C-Hg-C C-C → C-Zn-C C-C → C-Ag-C C-C → C-Ni-C C-C → C-Pd-C C-C → C-Pt-C C-C-Bnz → C-Hg-CBnz C-C-Bnz → C-Zn-C-Bnz C-C-Bnz → C-Ag-C-Bnz C-C-Bnz → C-Ni-C-Bnz C-C-Bnz → C-Pd-C-Bnz C-C-Bnz → C-Pt-C-Bnz		−91.33 −82.22 −94.27 −110.01 −80.00 −79.88 −132.45 −69.01 −81.34 −150.09 −99.58 −110.78	C-G-Bnz → C-Hg-G-Bnz C-G-Bnz → C-Zn-G-Bnz C-G-Bnz → C-Ag-G-Bnz C-G-Bnz → C-Ni-G-Bnz C-G-Bnz → C-Pd-G-Bnz C-G-Bnz → C-Pt-G-Bnz C-G-Bnz → C-Hg-G-Bnz C-G-Bnz → C-Zn-G-Bnz C-G-Bnz → C-Ag-G-Bnz C-G-Bnz → C-Ni-G-Bnz C-G-Bnz → C-Pd-G-Bnz C-G-Bnz → C-Pt-G-Bnz	−125.05 −61.04 −73.24 −147.69 −94.26 −104.97

Table 5 The HOMO-LUMO band gap energies of C-G, C-C and their complexes with Bnz and metal cations (eV)

Base Pairs	HOMO-LUMO gap energies
Base Pairs	-	Hg	Zn	Ag	Ni	Pd	Pt
C-G C-G-Bnz C-G/C-G-Bnz G-C G-C-Bnz G-C/G-C-Bnz C-C C-C-Bnz C-C/C-C-Bnz	1.573 0.973 0.600 1.573 0.580 0.993 0.829 0.914 −0.085	1.629 0.630 0.999 1.629 2.319 −0.690 2.377 2.051 0.326	3.198 3.489 −0.291 3.198 3.452 −0.254 3.800 3.652 0.148	1.201 2.708 −1.507 1.201 3.078 −1.877 2.977 2.029 0.948	1.814 2.046 −0.232 1.814 3.342 −1.528 2.993 2.916 0.077	2.329 3.470 −1.141 2.329 3.211 −0.882 2.979 3.274 −0.295	0.825 3.476 −2.651 0.825 3.188 −2.363 2.998 2.706 0.292

shows that except for A-T-Bnz base pair which prefers Hg²⁺ cation, other benzimidazole based DNA base pairs primarily prefer Ni²⁺ cation. Also, the HOMO-LUMO band gaps of the molecules can be seen in Table 5. It is observed from Table 5 binding benzimidazole to C-G has created 0,60 eV of conjugation andthis effect has been found for the Hg²⁺ cation while the band gap of other cations has in-creased. The greatest difference has been observed with Pt²⁺ cation. Binding of Bnz to G-C provides 0.993 eV of resonance stability. However, when complexation occurs with other metal cations, gap energy values have increased and this value is greater than what is for Pt²⁺ cation. Binding benzimidazole to C-C does not bring about a large effect and it is complexed with any metal cation other than Pd²⁺, gap energy value decreases and this decrease is greater than the decrease for Ag⁺. All the results related to band gaps show that resonance stability is increased with Hg²⁺, Ag⁺ and Zn²⁺ while thermodynamic stability is increased with Pt²⁺ and Pd²⁺ cations.

The absorption and emission spectrum has been studied and results of the complexes have been formed by C-G, C-G and C-C bases and benzimidazole based pairs with metals have been given in Table 6.

According to Table 6, maximum absorption and emission wavelength values have been calculated for the Pt²⁺ cation, which is the most probable cation for both C-M-G and C-M-G-Benzimidazole complexes. The results gathered for G-C based structures are the same with the data which obtained by C-G pair. Examination of the maximum absorption and emission values in the visible area shows that the biggest Δλ value belongs to the C-Ni-C and C-Ni-C-Bnz complexes.

Base pairs that show whatsoever colour change has been determined as in the previous section. Colour changes of the base pairs depending on before and after radiation has been presented in Table S9. As can be seen, colour changes of C-M-G and C-M-G-Bnz pairs are same with G-M-C and G-M-C-Bnz pairs because of absorption and emission wavelength of these pairs have equal values. Red coloured C-Hg-C, C-Ag-C and Bnz pairs have disappeared colour. A similar situation has occurred in blue coloured C-Pt-G and violet coloured C-Ni-C base

Table 6 Maximum absorption/emission wavelength (nm) and differences between them of the designed C-G, G-C, C-Cpairs and their Bnz complexes

Pairs	C-G- Based									C-G-Bnz Based
Pairs	λ Abs	Osc		λ Ems	Osc	∆λ	Pairsλ		Abs		Osc		λ Ems		Osc	∆λ
C-G C-Hg-G C-Zn-G C-Ag-G C-Ni-G C-Pd-G C-Pt-G	263 324 830 794 297 473 454		0.037 0.024 0.020 0.045 0.026 0.026 0.029	294 450 1047 815 514 491 789	0.020 0.170 0.015 0.030 0.025 0.024 0.026	30 126 216 21 217 18 334	C-G-Bnz C-Hg-G-Bnz C-Zn-G-Bnz C-Ag-G-Bnz C-Ni-G-Bnz C-Pd-G-Bnz C-Pt-G-Bnz		297 349 799 743 311 479 439			0.005 0.020 0.020 0.038 0.005 0.007 0.006	310 462 921 761 735 491 663	0.005 0.016 0.017 0.039 0.004 0.006 0.006			13 112 122 18 423 11 194
	G-C- Based									G-C-Bnz Based
G-C G-Hg-C G-Zn-C G-Ag-C G-Ni-C G-Pd-C G-Pt-C	263 324 830 794 297 473 454		0.037 0.024 0.020 0.045 0.026 0.026 0.029	294 450 1047 815 514 491 789	0.020 0.170 0.015 0.030 0.025 0.024 0.026	30 126 216 21 217 18 334	G-C-Bnz G-Hg-C-Bnz G-Zn-C-Bnz G-Ag-C-Bnz G-Ni-C-Bnz G-Pd-C-Bnz G-Pt-C-Bnz		279 376 784 718 293 468 403			0.015 0.015 0.012 0.081 0.098 0.076 0.006	290 597 891 774 417 489 615	0.014 0.011 0.011 0.079 0.086 0.075 0.005			10 220 107 56 178 21 212
	C-C Based							C-C-Bnz Based
C-C C-Hg-C C-Zn-C C-Ag-C C-Ni-C C-Pd-C C-Pt-C	398 699 778 670 403 775 786		0.001 0.019 0.002 0.009 0.014 0.017 0.086	419 742 797 985 729 892 802	0.001 0.018 0.002 0.008 0.012 0.011 0.070	21 42 18 314 326 117 16	C-C-Bnz C-Hg-C-Bnz C-Zn-C-Bnz C-Ag-C-Bnz C-Ni-C-Bnz C-Pd-C-Bnz C-Pt-C-Bnz		343 693 722 691 414 751 769			0.001 0.014 0.003 0.014 0.016 0.008 0.023	374 711 775 813 635 860 781	0.001 0.012 0.002 0.012 0.012 0.005 0.021			31 18 52 122 220 108 11

pairs. Blue-green change exhibited in C-Pd-G and its benzimidazole pair. Violet coloured C-Pt-G-Bnz, G-Pt-C-Bnz and C-Ni-C-Bnz pairs have changed their colours to red, orange and again orange, respectively. Colourless pairs C-Hg-G and its Bnz pair has turned their colour to blue while C-Ni-G green.

3.3. Logic Gates

The most probable DNA base pairs with and without metal cations that could be demonstrate Stokes shift has been selected for logic gate calculations. Acidic media effect is included by the protonation of the nitrogen atom on benzimidazole fragment. For this purpose, a proton has been linked to the nitrogen atom of benzimidazole and calculated their absorption and emission spectrums aqueous phase. The results of the selected T-Zn-T, T-Hg-A, A-Ni-T, C-Pt-G, G-Pt-C and C-Ni-C base pairs have been given in the Table 7 and Table 8.

Examination of the tables show that T-Hg-A, A-Ni-T, C-Pt-G, C-Ni-C base pairs represents an OR logic gate, while T-Zn-T, G-Pt-C produces AND gate and XOR gate, respectively. In the AND gate for T-Zn-T pair, fluorescence takes

Table 7 Logic Gates for T-T, T-A, A-T base pairs and their complexes with Bnz

Pairs	λ Abs (nm)	λ Ems (nm)		Δλ (nm)		Input H⁺	Input Zn⁺	Input Hg²⁺	Input Ni⁺	Output AND	Output OR	λ Max (nm)
T-T-Bnz T-H⁺-T-Bnz T-Zn-T-Bnz T-Zn,H⁺-T-Bnz T-A-Bnz T-H⁺-A-Bnz T-Hg-A-Bnz T-Hg,H⁺-A-Bnz A-T-Bnz A-H⁺-T-Bnz A-Ni-T-Bnz A-Ni,H⁺-T-Bnz	319 517 468 843 334 412 468 574 365 487 443 364		337 531 482 912 357 431 509 581 404 356 661 381		18.23 24.81 13.60 91.02 22.87 54.18 41.68 73.20 38.69 43.26 217.40 57.05	0 1 0 1 0 1 0 1 0 1 0 1	0 0 1 1	0 0 1 1	0 0 1 1	0 0 0 1	0 1 1 1 0 1 1 1	337.73 531.34 482.10 912.80 357.04 431.13 509.08 581.34 404.40 536.38 661.20 381.09

Table 8 Logic Gates for C-G, G-C, C-C base pairs and their complexes with Bnz

Pairs	λ Abs (nm)	λ Ems (nm)	Δλ (nm)	Input H⁺	Input Pt⁺	Input Ni⁺	Output OR	Output XOR		λ Max (nm)
C-G-Bnz C-H⁺-G-Bnz C-Pt-G-Bnz C-Pt,H⁺-G-Bnz G-C-Bnz G-H⁺-C-Bnz G-Pt-C-Bnz G-Pt,H⁺-C-Bnz C-C-Bnz C-H⁺-C-Bnz C-Ni-C-Bnz C-Ni,H⁺-C-Bnz	297 363 439 812 279 964 403 321 343 363 414 276	310 401 663 940 290 1012 615 340 374 382 635 616	13.77 74.43 194.76 217.60 10.38 148.90 212.16 29.03 31.53 37.66 220.34 89.47	0 1 0 1 0 1 0 1 0 1 0 1	0 0 1 1 0 0 1 1	0 0 1 1	0 1 1 1		0 1 1 0	310.80 401.06 663.76 940.62 290.03 1012.1 615.2 340.6 374.60 382.17 635.24 316.43

place when the proton and the Zn²⁺ cation available in the system at the same time. XOR gate for G-Pt-C occurs in the event of being proton or Pt²⁺ cation. For the afore-mentioned other base pairs have shown OR logic gateand it means that fluorescence can be seen when the metal cation, proton and both of them are in the system.

Also, bonding of metal cation and protonation of selected pairs has been caused to red shift in absorption spectrum as well as emission spectrum for all the interested base pairs. C-Ni-C-Bnz pair shows the most largest Stokes shift among the molecules that have been selected for logic gates calculations. The C and G base pairs have higher Stokes shift values than the thymine and adenine pairs.

4. Conclusion

In this study, we have theoretically studied complexation energies, their band gap and electronic absorption-emission spectral behaviors of targeted DNA base pairs in the water phase. In order to determine the appropriate method for UV-vis absorption wavelength maxima of pairs, it has been calculated with three different functionals and compared with the experimental data and each other. Results show that M06/Lanl2dz level is good agreement with experimental absorption wavelength. The energy calculation results give that all the base pairs primarily prefer Ni²⁺ cation for complexation, whereas A-T base pairs prefer Hg²⁺ cation. The same trend has been observed in case of being benzimidazole in pairs. HOMO-LUMO band gap results have shown that the resonance stability increased with Hg²⁺, Ag⁺ and Zn²⁺ cations, while the Pt²⁺ and Pd²⁺ provide thermodynamically stable. Our calculated electronic spectrum presents that designed DNA base pairs can be use as a probe to detected selected cations. The fluorescence of C-G, G-C pairs answers for Pt^2+, Pd²⁺ and Zn²⁺ cations, while T-T, T-A and A-T are given for Hg²⁺, Ag⁺ and Zn²⁺ cations. The calculated results for their logic gates have been given that addition of protons to designed DNA pairs causes a red shift for all pairs in the water phase. Also, the presence of metal cations causes a red shift like protonation, except for Ni²⁺ cation. As a result of calculated absorption-emission spectrum data show that T-Hg-A-Bnz, A-Ni-T-Bnz, C-Pt-G-Bnz, C-Ni-C-Bnz complexes produce OR gate. T-Zn-T- Bnz and G-Pt-C-Bnz results demonstrated XOR and AND logic gate, respectively. In brief, this theoretical study manifestes important electronic and photophysical behaviors of designed metallo-DNA pairs which can be use to determining of selected cations.

Acknowledgements

We acknowledge the support of TUBITAK (Scientific and Technical Research Council of the Turkish Republic) under grant no. 214Z022.

Cite this paper

Sevin, F. Dilmani, M.A., Sarıkavak, K. and Jeddi, S.F. (2017) Theoretical Investigation of the Newly Designed Benzimidazole Based Metal Mediated DNA Base Couples with DFT Method. Computational Chemistry, 5, 74-90. https://doi.org/10.4236/cc.2017.52007

Supplementary Tables

Table S1 HOMO-LUMO molecular orbital energies and complexation energies in eV for T-Hg-T, calculated with MO6, B3LYP and PBE0 method (Kcal/mol)

Methods	T-Hg-T
Methods	HOMO	LUMO	ΔGap(L-H)	ΔG	ΔH	ΔE	Absorbance
B3LYP PBE0 M06	−6.46 −5.56 −6.76	−1.15 −1.78 −0.95	5.30 5.78 5.81	−595906.58 −595417.20 −595556.71	−595862.38 −595373.64 −595513.17	−595874.80 −595386.05 −595525.35	269.34 331.99 264.48

Table S2 The complexation energies of thymine-thymine (T-T), thymine-adenine (T-A) base pairs with metal cations (Kcal/mol)

Base Pairs	Energy Values
Base Pairs	ΔE	ΔG	ΔH	Base Pairs	ΔE	ΔG	ΔH
T-T→ T-Hg-T T-T→ T-Zn-T T-T→ T-Ag-T T-T→ T-Ni-T T-T→ T-Pd-T T-T→ T-Pt-T	−113.11 −89.06 −69.98 −147.64 −52.12 −63.41	−122.66 −98.06 −78.26 −183.55 −59.08 −71.98	−122.64 −88.49 −69.69 −143.50 −51.23 −61.54	T-A → T-Hg-A T-A → T-Zn-A T-A → T-Ag-A T-A → T-Ni-A T-A → T-Pd-A T-A → T-Pt-A	−97.02 −77.75 −64.99 −123.26 −47.78 −56.36	−101.32 −81.46 −72.19 −134.85 −50.00 −61.35	95.98 −77.02 −61.76 −119.63 −47.22 −54.49

Table S3 The complexation energies of T-T-Bnz and T-A-Bnz base pairs with metal cations (Kcal/mol)

	Energy Values
Base Pairs	ΔE	ΔG	ΔH	Base Pairs	ΔE	ΔG	ΔH
T-T-Bnz → T-Hg-T-Bnz T-T-Bnz → T-Zn-T-Bnz T-T-Bnz → T-Ag-T-Bnz T-T-Bnz → T-Ni-T-Bnz T-T-Bnz → T-Pd-T-Bnz T-T-Bnz → T-Pt-T-Bnz	−112.91 −89.37 −70.01 −147.98 −52.18 −63.74	−112.08 −99.02 −79.37 −184.62 −60.81 −73.26	−112.53 −88.62 −69.91 −143.72 −51.39 −62.35	T-A-Bnz → T-Hg-A-Bnz T-A-Bnz → T-Zn-A-Bnz T-A-Bnz → T-Ag-A-Bnz T-A-Bnz → T-Ni-A-Bnz T-A-Bnz → T-Pd-A-Bnz T-A-Bnz → T-Pt-A-Bnz	−97.02 −77.75 −64.99 −123.26 −47.78 −56.36	−96.05 −77.00 −64.22 −122.41 −45.89 −56.00	−94.01 −75.67 −60.80 −121.33 −45.46 −55.03

Table S4 The complexation energies of A-T-Bnz base pairs with metal cations (Kcal/mol)

Base Pairs	Energy Values
Base Pairs	ΔE	ΔG	ΔH
A-T-Bnz → A-Hg-T-Bnz A-T-Bnz → A-Zn-T-Bnz A-T-Bnz → A-Ag-T-Bnz A-T-Bnz → A-Ni-T-Bnz A-T-Bnz → A-Pd-T-Bnz A-T-Bnz → A-Pt-T-Bnz	−101.11 −83.10 −71.26 −102.95 −51.01 −63.53	−105.22 −84.92 −76.00 −104.55 −54.02 −67.26	−98.88 −83.00 −69.77 −100.00 −50.13 −63.06

Table S5 The colour changes of of T-T, T-A, A-T and their complexes with Bnz and metal cations

T-M-T	B.L	A.L	T-M-T-Bnz	B.L	A.L
T-T T-Hg-T T-Zn-T T-Ag-T T-Ni-T T-Pd-T T-Pt-T	- - Blue Yellow - - Orange	- - Blue - - - -	T-T-Bnz T-Hg-T-Bnz T-Zn-T-Bnz T-Ag-T-Bnz T-Ni-T-Bnz T-Pd-T-Bnz T-Pt-T-Bnz	- - Blue - - - -	- Violet Green - - - -
T-M-A	B.L	A.L	T-M-A-Bnz	B.L	A.L
T-A T-Hg-A T-Zn-A T-Ag-A T-Ni-A T-Pd-A T-Pt-A	- Blue Yellow Red - - -	- Green Red - Blue - -	T-A-Bnz T-Hg-A-Bnz T-Zn-A-Bnz T-Ag-A-Bnz T-Ni-A-Bnz T-Pd-A-Bnz T-Pt-A-Bnz	- Blue Red - - - -	- Green - - Blue - -
A-M-T	B.L	A.L	A-M-T-Bnz	B.L	A.L
A-T A-Hg-T A-Zn-T A-Ag-T A-Ni-T A-Pd-T A-Pt-T	- Blue Yellow Red - - -	- Green Red - Blue - -	A-T-Bnz A-Hg-T-Bnz A-Zn-T-Bnz A-Ag-T-Bnz A-Ni-T-Bnz A-Pd-T-Bnz A-Pt-T-Bnz	- Violet Red - Violet - -	Violet Green Red - Red - -

B.L: Before luminescence; A.L: After luminescence

Table S6 The complexation energies of C-G and C-C base pairs with metal cations (Kcal/mol)

Base Pairs	Energy Values
Base Pairs	ΔE	ΔG	ΔH	Base Pairs	ΔE	ΔG	ΔH
C-G → C-Hg-G C-G → C-Zn-G C-G → C-Ag-G C-G → C-Ni-G C-G → C-Pd-G C-G → C-Pt-G	−137.74 −62.16 −75.82 −151.04 −97.76 −110.26	−140.33 −63.55 −77.77 −154.64 −101.04 −113.36	−135.80 −60.06 −75.10 −149.86 −95.55 −109.55	C-C → C-Hg-C C-C → C-Zn-C C-C → C-Ag-C C-C → C-Ni-C C-C → C-Pd-C C-C → C-Pt-C	−88.10 −81.07 −93.00 −108.75 −77.42 −79.55	−91.33 −82.22 −94.27 −110.01 −80.00 −79.88	−85.90 −79.56 −89.90 −107.00 −76.20 −76.59

Table S7 The complexation energies of C-G -Bnz and G-C-Bnz base pairs with metal cations (Kcal/mol)

Base Pairs	Energy Values
Base Pairs	ΔE	ΔG	ΔH	Base Pairs	ΔE	ΔG	ΔH
C-G-Bnz → C-Hg-G-Bnz C-G-Bnz → C-Zn-G-Bnz C-G-Bnz → C-Ag-G-Bnz C-G-Bnz → C-Ni-G-Bnz C-G-Bnz → C-Pd-G-Bnz C-G-Bnz → C-Pt-G-Bnz	−124.05 −61.04 −73.24 −147.69 −94.26 −104.97	−125.05 −61.04 −73.24 −147.69 −94.26 −104.97	−123.86 −58.97 −73.00 −146.11 −92.50 −104.00	G-C-Bnz → G-Hg-CBnz G-C-Bnz → G-Zn-C-Bnz G-C-Bnz → G-Ag-C-Bnz G-C-Bnz → G-Ni-C-Bnz G-C-Bnz → G-Pd-C-Bnz G-C-Bnz → G-Pt-C-Bnz	−130.99 −66.11 −79.05 −148.27 −96.13 −108.60	−132.45 −69.01 −81.34 −150.09 −99.58 −110.78	−128.94 −64.89 −78.22 −148.00 −94.99 −108.01

Table S8 The complexation energies of C-C-Bnz base pairs with metal cations (Kcal/mol)

Base Pairs	Energy Values
Base Pairs	ΔE	ΔG	ΔH
C-C-Bnz→ C-Hg-CBnz C-C-Bnz→ C-Zn-C-Bnz C-C-Bnz→ C-Ag-C-Bnz C-C-Bnz→ C-Ni-C-Bnz C-C-Bnz→ C-Pd-C-Bnz C-C-Bnz→ C-Pt-C-Bnz	−130.99 −66.11 −79.05 −148.27 −96.13 −108.60	−132.45 −69.01 −81.34 −150.09 −99.58 −110.78	−128.94 −64.89 −78.22 −148.00 −94.99 −108.01

Table S9 The colour changes of of C-G, G-C, C-C and their complexes with Bnz and metal cations

C-M-G	B.L	A.L	C-M-G-Bnz	B.L	A.L
C-G C-Hg-G C-Zn-G C-Ag-G C-Ni-G C-Pd-G C-Pt-G	- - Green - - Blue Blue	- Blue Green - Green Green -	C-G-Bnz C-Hg-G-Bnz C-Zn-G-Bnz C-Ag-G-Bnz C-Ni-G-Bnz C-Pd-G-Bnz C-Pt-G-Bnz	- - - - - Blue Violet	- Blue - - - Green Red-
G-M-C	B.L	A.L	C-M-G-Bnz	B.L	A.L
G-C G-Hg-C G-Zn-C G-Ag-C G-Ni-C G-Pd-C G-Pt-C	- - Green - - Blue Blue-	- Blue Green - Green Green --	G-C-Bnz G-Hg-C-Bnz G-Zn-C-Bnz G-Ag-C-Bnz G-Ni-C-Bnz G-Pd-C-Bnz G-Pt-C-Bnz	- - - - - Blue Violet	- Yellow - - Violet Green Orange-
C-M-C	B.L	A.L	C-M-C-Bnz	B.L	A.L
C-C C-Hg-C C-Zn-C C-Ag-C C-Ni-C C-Pd-C C-Pt-C	- Red - Red Violet - -	Violet - - - - - -	C-C-Bnz C-Hg-C-Bnz C-Zn-C-Bnz C-Ag-C-Bnz C-Ni-C-Bnz C-Pd-C-Bnz C-Pt-C-Bnz	- Red - Red Violet - --	- - - - Orange - -

References1

Holliday, B.J. and Mirkin, C.A. (2001) Strategies for the Construction of Supramolecular Compounds through Coordination Chemistry. Angewandte Chemi, 40, 2022-2043. http://dx.doi.org/10.1002/1521-3773(20010601)40:11<2022::AID-ANIE2022>3.0.CO;2-D

Panja, K.S., Dwivedi, N. and Saha, S. (2016) Highly Stable Naphthalene Core Based Novel Cleftshaped Strain Molecule: Influence of Intermolecular H-Bonding Architectures. The Royal Society of Chemistry, 6, 59574-59581. http://dx.doi.org/10.1039/c6ra06855c

Mahmudov, K.T., Kopylovich, M.N., Silva, F.C. and Pombeiro, A.J.L. (2016) Non-Covalent Interactions in the Synthesis of Coordination Compounds: Recent Advances. Coordination Chemistry Reviews. http://dx.doi.org/10.1016/j.ccr.2016.09.002

Hisamatsu, Y., Miyazawa, Y., Yoneda, K., Miyauchi, M., Zulkefeli, M. and Aoki, S. (2016) Supramolecular Complexes Formed by the Self-Assembly of Hydrophobic Bis(Zn2+-cyclen) Complexes, Copper, and Dior Triimide Units for the Hydrolysis of Phosphate Mono- and Diesters in Two-Phase Solvent Systems (Cyclen = 1,4,7,10-Tetraazacyclododecane). Chemical and Pharmaceutical Bulletin, 64, 451-464. http://dx.doi.org/10.1248/cpb.c15-01014

Smith, S.J., Radford, R.J., Subramanian, R.H., Barnett, B.R., Figueroa, J.S. and Tezcan, F.A. (2016) Tunable Helicity, Stability and DNA-Binding Properties of Short Peptides with Hybrid Metal Coordination Motifs. Chemical Science, 7, 5453-5461. http://dx.doi.org/10.1039/c6sc00826g

Santamaria-Diaz, N., Mendez-Arriaga, M.J., Salas, J.M. and Galindo, M.A. (2016) Highly Stable Double-Stranded DNA Containing Sequential Silver (I)-Mediated 7-Deazaadenine/Thymine Watson-Crick Base Pairs. Angewandte Chemie, 55, 6170-6174. http://dx.doi.org/10.1002/anie.201600924

Léon, J.C., Sinha, I. and Müller, J. (2016) 6-Pyrazolylpurine as an Artificial Nucleobase for Metal-Mediated Base Pairing in DNA Duplexes. International Journal of Molecular Sciences, 17, 554-564. http://dx.doi.org/10.3390/ijms17040554

Camacho, C.C., Ruiz, A.E., Hueso P.A., Mijangos, E., Garcia, I.R., Contreras, R. and Parra, A.F. (2013) Organometallic Tin Compounds Derived from 2-Benzimidazole Propionic Acid. Zeitschrift für anorganische und allgemeine Chemie, 7, 1122-1128. http://dx.doi.org/10.1002/zaac.201300108

Mallajosyula, S.S. and Pati, S.K. (2008) Benzimidazole-Modified Single-Stranded DNA: Stable Scaffolds for 1-Dimensional Spintronics Constructs. Journal of Physical. Chemistry B, 112, 16982-16989. http://dx.doi.org/10.1021/jp8080782.

Miyake, Y., Togashi, M., Tashiro, M., Yamaguchi, H., Oda, S., Kudo, M., Tanaka, Y., Kondo, Y., Sawa, R., Fujimoto, T., Machinami, T. and Ono, A. (2005) MercuryII-Mediated Formation of Thymine-HgII-Thymine Base Pairs in DNA Duplexes. Journal of the American Chemical Society, 128, 2172-2173. http://dx.doi.org/10.1021/ja056354d

Miyachi, H., Matsui, T., Shigetac, Y. and Hirao, K. (2009) Effects of Mercury (II) on Structural Properties, Electronic Structure and UV Absorption Spectra of a Duplex Containing Thymine-Mercury(II)-Thymine Nucleobase Pairs. Physical Chemistry Chemical Physics, 12, 909-917. https://doi.org/10.1039/B912807G

Takezawa, Y. and Shionoya, M. (2012) Metal-Mediated DNA Base Pairing: Alternatives to Hydrogen-Bonded Watson Crick Base Pairs. Accounts of Chemical Research, 45, 2066-2076. https://doi.org/10.1021/ar200313h

Wu, J., Fu, Y., He, Z., Han, Y., Zheng, L., Zhang, J. and Li, W. (2012) Growth Mechanisms of Fluorescent Silver Clusters Regulated by Polymorphic DNA Templates: A DFT Study. The Journal of Physical Chemistry B, 116, 1655-1665. https://doi.org/10.1021/jp206251v

Müller, A., Talbot, F. and Leutwyler, S. (2002) Hydrogen Bond Vibrations of 2-Aminopyridinea2-Pyridone, Watson-Crick Analogue of Adenine Uracil. Journal of the American Chemical Society, 124, 14486-14494. https://doi.org/10.1021/ja0209969

Brown, G.J. (2004) Encyclopedia of Supramolecular Chemistry. In: Atwood, J.L. and Steed, J.W., Eds., Encyclopedia of Supramolecular Chemistry, CRC Press, Taylor and Francis Group, New York, 893.

Kubo, K. (2005) PET Sensors. In: Geddes, C.D. and Lakowicz, J.R., Eds., Topics in Fluorescence Spectroscopy: Advanced Concepts in Fluorescence Sensing, Springer Science Business Media, New York, 241.

Frisch, M.J., Trucks, G.W., Schlegel, H.B., Scuseria, G.E., Robb, M.A., Cheeseman, J.R., Scalmani, G., Barone, V., Mennucci, B., Petersson, G.A., Nakatsuji, H., Caricato, M., Li, X., Hratchian, H.P., Izmaylov, A.F., Bloino, J., Zheng, G., Sonnenberg, J.L., Hada, M., Ehara, M., Toyota, K., Fukuda, R., Hasegawa, J., Ishida, M., Nakajima, T., Honda, Y., Kitao, O., Nakai, H., Vreven, T., Montgomery, J.A., Peralta, J.E., Ogliaro, F., Bearpark, M., Heyd, J.J., Brothers, E., Kudin, K.N., Staroverov, V.N., Kobayashi, R., Normand, J., Raghavachari, K., Rendell, A., Burant, J.C., Iyengar, S.S., Tomasi, J., Cossi, M., Rega, N., Millam, J.M., Klene, M., Knox, J.E., Cross, J.B., Bakken, V., Adamo, C., Jaramillo, J., Gomperts, R., Stratmann, R.E., Yazyev, O., Austin, A.J., Cammi, R., Pomelli, C., Ochterski, J.W., Martin, R.L., Morokuma, K., Zakrzewski, V.G., Voth, G.A., Salvador, P., Dannenberg, J.J., Dapprich, S., Daniels, A.D., Farkas, O., Foresman, J.B., Ortiz, J.V., Cioslowski, J. and Fox, D.J. (2009) Gaussian 09, Revision E.01. Gaussian Inc., Wallingford.

Dennington, R., Keith, T. and Millam, J. (2009) Gauss View.

Zhang, B., Ma, P., Gao, D., Wang, Z., Sun, Y., Song, D. and Li, X. (2016) A FRET-Based Fluorescent Probe Formercury Ions in Water and Living Cells. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 165, 99-105.

Sicilia, V., Borja, P., Bayab, M. and Casasb, J. (2015) Selective Turn-Off Phosphorescent and Colorimetric Detection of Mercury (II) in Water ByhalfLantern Platinum(II) Complexes. Dalton Transactions, 44, 6936-6943. https://doi.org/10.1039/C5DT00087D

Karhu, A., Rautiainen, M., Oilunkaniemi, R., Chivers, T. and Laitinen, R. (2015) Mercury- and Cadmium-Assisted [2 + 2] Cyclodimerization of Tert-Butylselenium Diimide. Inorganic Chemistry, 54, 9499-9508. https://doi.org/10.1021/acs.inorgchem.5b01434

Zhao, Y. and Truhlar, D.G. (2006) The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States and Transition Elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 Other Functionals. Theoretical Chemistry Accounts, 120, 215-241. https://doi.org/10.1007/s00214-007-0310-x

Rives, J.T. and Jorgensen, W.L. (2008) Performance of B3LYP Density Functional Methods for a Large Set of Organic Molecules. Journal of Chemical Theory and Computation, 4, 297-306. https://doi.org/10.1021/ct700248k

Adamo, C. and Barone, V. (1999) Toward Reliable Density Functional Methods without Adjustable Parameters: The PBE0 Model. The Journal of Chemical Physics, 113, 6158-6170. https://doi.org/10.1063/1.478522

Chiodo, S., Russo, N. and Sicilia, E. (2006) LANL2DZ Basis Sets Recontracted in the Framework Of Density Functional Theory. The Journal of Chemical Physics, 125, 104-107. https://doi.org/10.1063/1.2345197

Klamt, A. (1994) Conductor-Like Screening Model for Real Solvents: A New Approach to the Quantitative Calculation of Solvation Phenomena. The Journal of Physical Chemistry, 99, 2224-2235. https://doi.org/10.1021/j100007a062

http://www.gaussian.com/g_tech/g_ur/k_scrf.htm