Quantum Mechanical Investigation of 2 , 3-dihydro-1-methyl-2-thioxopyrimidin-4 ( 1 H )-one Tautomers and their Complexes with Iodine

DFT and AB initio theoretical methods were used to calculate the relative stability of tautomers and complexes with iodine in the 2,3-dihydro-1-methyl-2-thioxopyrimidin-4(1H)-one. This compound can be used to treat hyper thyroidism due to their ability to make complexes with iodine.All tautomers and complexes are optimized using the B3LYP Method with two different energies, the relative energies shows that in all tautomers and complexes. Thione forms are more stable than thieolforms.The NBO calculation is carried out for tautomers and complexes to obtain atomic charges and accept donor interactions. These result confirm the ability of T tautomers to form complexes and show that the planar complexes has more effective interaction than the perpendicular complex.The atom in molecule (AIM) analysis show that the charge density and its laplacian at the S-I bond critical point of the planar complex is greater than the perpendicular complex.

ability to complex with iodine, which can be used as a new drug for hyper thyroidism.
In the present work, we used computational method to investigate tautomerismand the stability of 2,3-dihydro-1methyl-2-thioxopyrimidin-4(1H)-one complexes in the gas phase.Another important aspect of this molecule is its diversity in tautomerism 20,21 .
g.1.possibletautomers and numbering schem for T and structuers of complexes.
A study of tautomerism, especially in biologically active molecules, has been one of the most interesting in computational researches because tautomerism plays an important role in the determination of compound application [22][23][24][25][26][27][28][29][30][31][32][33][34] , and different tautomers of each molecule behave differently in both chemical and biological systems.Thus, since the biological, chemical, and complexation properties of T are different in different tautomers, we have attempted a thorough analysis of this molecule in both the tautomery scheme and complexation properties in our study.
In this research, molecular parameters (Table 1), relative energie (Table 2), of T tautomers were calculated usingB3LYP/6-311++G ** and B3LYP/6-31+G ** levels of theory to study its tautomery scheme and properties of all tautomers.In addition, NBO calculations were made to obtain natural atomic charges (Table 4).Occupation numbers (Table 5), and acceptor_donorinteraction (Table 6,7)of all tautomers and complexes using the B3LYP/6-311++g ** level of theory then, AIM analyses (Table 8) of the complexes were performed to determine the complexation properties of different tautomers in the interaction with iodine.Computation details and the result obtained in this work are presented below.

Methods
Density functional theory (DFT) has been widely applied by chemists to study the electronic structure of molecule in the past year 35,36 .In this work, all calculation were carried out using Becke'sthree parameter density functional 37 and Lee, Yang, and Parr functional [Table 13] to describe gradient-corrected correlation effects, which leads to the well-known B3LYPmethod.The B3lyp method has been validatedto give result similar to those of the more computationally expensive MP2 theory for molecular geometry and frequency calculation 39- 41 .The geometry optimizations were performed for all tautomers by the B3LYP method with 6-31++G ** and 6-31++G ** basis sets.And the NBO analysis 42 was carried out at the B3LYP/6-31++G ** level of theory.All optimizations and NBO analysis was carried out using the Gaussion03 program package 43 .AIM analyses were performed using the AIM 2000 program 44 .This method has presented useful information about intermolecular interactions and the characterization of bonds through the analysis of the electron density 45 .

RESULT AND DISCUSSION
All geometries of T were optimized with the B3LYP methods as a part of the Gaussian 03 3 .
A theorical rational method cannot presently be reached for complexes with so many electrons.Therefore, for T-I 2 complex we used DFT methods with two basis sets.
As it can be seen from the calculated Gibbs free energy values listed in Table (1 ), the B3LYP methods all suggest that thionetautomer T is more stable than the thiol tautomer in the gas phase.By observing Table 1, important aspect of the molecular structure can be followed.The C 3 -S bond lengths are listed in the first row.The value of this bond length in T is 1.672A R" and1.693-1.696AR" in their complexes.
The difference between the C = S bond lengthin complexes and monomer can result interactions between anti-bonding orbital Lone pair of sulfur and iodine.In the presence interaction ,the charge-transfer of the ring nitrogen atom to sulfur increases and the C = S bond length decreases relative to the monomer.
In next rows, S-I and I-I bond length of the complexes are listed.The S-I bond is slightly shorter in the planar complexes than that in the perpendicular complexes, but any important difference does not exist between different complexes in I-I bond lengths.
Next 4 rows of table 1consist of bond angles.Observing bond angle variation, we can follow hybridization changes in the central action of each angle.For example, when the central atom is SP 2 , its angle is near 120 degree (as C 3 in all tautomers and N 4 in T 2 and T 4 ).
But for C 3 , the exocyclic angle (N 2 -C 3 -N 4 ) is smaller than the exocyclic angle (S-C 3 -N 4 ) because of a higher mobility of the external sulfur atom.
Dihedral angles are listed in the last columns of table 1.The value of the N 2 -C 1 -C 5 -N 4 dihedral angle can show the degree of planarity of the molecules.This value is near zero in all tautomers and complexes, so that these values confirm the planarity of all structures.
Moreover, the I-I-S-C 3 and I-S-C 3 -N 4

Energies
The relative zero point energies, enthalpies, and Gibbs free energies for the most stable tautomers and isomers.T has the lowest enthalpy and Gibbs free energy at both levels of theory.T has both C=S and C=O strong double bonds, so these strong bonds can be responsible for the most stability of this tautomer.The some observation about the relative stability of cyclic tautomers was observed in various pyrimidine and purine bases [47][48][49] .
The relative stability of other tautomers and geometric isomers obtained at the B3LYP/6-31++G ** Level of theory is found to be as follows: T> T 1 >T 2 >T 3 >T 4 .
Moreover, the relative stability of the complexes is found to be as follows: T-I 2 (planar)>T-I 2 (perp).The data in the table ( 2) show that "G for formation complex (T-I2) planer by 10.7 kJ/ mol lower than another T-I2complex .The equilibrium constant of the complex formed T-I2planar is more than T-I2 perpendicular.Enthalpy and entropy reduced with formation complex.The reduction in the enthalpy of the complex plane is more the complexperpendicular.
In order to support this hypothesis, we have calculated the NHstretching wave number of the T and T-I2 complexes.The NH stretching wave numbers of the T and T-I2 complexes as computed at B3LYP/6-31++G** and B3LYP/6-31++G* (3-21G** for iodine atom) levelof theory, are reported in Table 3.The comparison of the calculated vibrational wave  In tautomers (T,T 1 ,T 2 ,T 3 ,T 4 ), all nitrogen and oxygen atoms have negative charges and carbon and hydrogen atoms have positive charges.These charges are usual, however, interesting charges were observed in the sulfur atom.In T and complex sulfur has a negative charge, while in the other tautomers it has a positive charge because in T, the sulfur atom is bound only to C 3 with a double bond while in the others, sulfur is bond to two atoms by single bonds.Other atomic charges of the atoms in tautomers have usual values with a little variation between the tautomers.In the complexes, the sum of net charges of all atoms in larger in the tautomer than those in simple tautomers, because negative charges are placed on the iodine atoms.These charges provide another proof for the existence of a real complex between the tautomers and iodine.Also, in the planar complex, the absolute values of charges placed on the iodine atoms are larger than perpendicular complex.This shows that a more effective interaction (between the tautomer and iodine) seems to exist in the planar complex.
One of the important results obtained from NBO calculations is orbital occupancies and another is acceptor-donor interactions.In table 5, selected occupation numbers for all tautomers and complexes are listed.The occupancies given show that most bounding orbitals and lone pair consist of more than 1.9 electrons, but some lone pairs have occupation numbers between 1.6 and 1.98.these occupancies display that important donor properties can exist on these lone pairs.Also, a significant change in the occupancies of the ´ or À bond between the tautomers and complexes has not been observed.This observation confirms that the  Or  bond do not contribute in electron donation to iodine.Onbehalf , a significant decrease was observed in the occupation number of lone pair 2 of sulfur in the complexes versus tautomer T, which shows that this lone pair is donated to the iodine atom.The occupation number of  * C 5 -O 8 in the tautomers show that resonance is between the N 4 lone pair and the C 5 =O 8 double bond.Otherwise, the N 4 lone pair has resonance with  * C 5 -O 8 The same result were observed in the  * C 3 -S 6 Occupation number in T and another tautomers.
Because of a higher capacity of the sulfur atom in electron acceptance, these occupancies are higher in the previous example.
The list of important donor-acceptor interactions is given in table 6.As shown in this table, the lone pairs of N 2 , N 4 , and S6 are important donors, and C-N, C-S, and C-O bonds are the main acceptors.The most powerful interactions are LP 1 N 2 to  * C 3 -S 6 and LP 1 N 4 to  * C 4 -O 8 , and powerful resonance exists between the nitrogen lone pair and the C=S or C=O double bond.Another important interaction is electron donation by LP 2 S 6 to the  * C 3 -N 4 .The table data listed shows that although both sulfur lone pairs have been in electron donation to iodine, however, one of them (LP 2 ) is more effective.Moreover in the planar complexes the acceptor-donor energy is significantly higher than in the perpendicular complex in both interactions (LP 1 and LP 2 to  * I-I).This observation confirms our previous estimate about higher efficiency .AIM results.In the final part of our study, AIM analyses performed for the complexes to obtain important complexation properties, and the results are listed in table 8.

This interaction has not been observed in
The first column consists of S-I interaction data.In this part, P (electron destiny) and this Laplacian may be very useful parameters to estimate the strengths of the S-I interaction.
The low values for the change destiny reflect the weak character of this bond in our complexes.
It is noticeable that the planar complex has the larger value than perp complex of the electron density, confirming the higher efficiency of the planar complex.The negative values for the laplacian of the C-S bond show the covalent bond.

CONCLUSIONS
In this report2,3-dihydro-1-methyl-2thioxopyrimidin-4(1H)-one and their tautomers have been studies to find new structures with the effective ability to make a complex with iodine (for treatment of hyperthyroidism).
All structures and their complex with iodine have been optimized using the B3LYP method with 6-31+G ** and 6-31++G ** Basis sets (for Iodine 3-21G**)to obtain their molecular parameters and relative energies.The relative stability of the complexes is found to be:T>T1> T2> T3> T4 The optimized structures show that T menomer make planar complex and perpendicular complex with iodine.The relative stability of the complexes is found to be: T-I 2 (planar)>T-I 2 (perp).Then, NBO calculations were performed for the tautomers and complexes to obtain atomic charges, occupation numbers, and acceptor-donor interactions.These results confirm the ability of (T) tautomers to form complexes and show that the planar complexes have more effetely interaction than the perp complex.Finally the AIM analyses performed to prove them and obtain complexation properties.This calculation of the T with iodine is more powerful in the planar complex versus the perpendicular complex. dihydro-1-methyl-2-thioxopyrimidin-4(1H)-one(T)

Table 6 : Acceptor-Donor Interaction Energies (in kcal/mol) in the tautomers and complexes
The numbering scheme of this table is the same as the scheme presented in fig .1.

Table 8 .Calculated critical point properties of the complexes
complexes because of the donation of this lone pair to iodine.It is obvious that in T 1 and T 2 , the C 3 =N 4 double bond exists.So this interaction is replaced by LP 2 S 6 donation to  * C 3 -N 4 .In addition, special donor-acceptor interactions were listed only for the complexes in table 7. most