Next Article in Journal
Modelling of Octahedral Manganese II Complexes with Inorganic Ligands: A Problem with Spin-States
Previous Article in Journal
Studies on Electronic Charge of the Hydrogen Bond Proton in Model Molecular Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Structure, Stability and Interaction Studies on Nucleotide Analogue Systems

1
Department of Physics, Bharathiar University, Coimbatore-641046, India
2
Department of Physics, Nallamuthu Gounder Mahalingam College, Pollachi-642 001, India
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2003, 4(8), 486-502; https://0-doi-org.brum.beds.ac.uk/10.3390/i4080486
Submission received: 27 January 2003 / Accepted: 25 August 2003 / Published: 30 August 2003

Abstract

:
Most of the biological molecules have a good interaction with water molecules. The hydrogen bonding interactions with the structural analogues of nucleic acid phosphate group namely dimethyl phosphate anion (DMP) and diethyl phosphate anion (DEP) are studied employing the ab initio and density functional theory methods. Inspections have been made to locate the reactive sites for the interactions of isomeric forms of mono, di and tri hydrates of alkyl phosphate anion using the above theories. It reveals, water molecules have a very strong interaction with the phosphate group in both the molecules and their interactions induce the changes in the structural parameters of the PO4 group for both the DMP and DEP anions. The optimized structural parameters, total energy, dipole moment and rotational constants are calculated and are compared with the available experimental values. The chemical hardness and chemical potential for these complexes have been calculated at HF/6-31G* level of theory and discussed the conformational stability of these complexes.

Introduction

The functional diversity of DNA or RNA involves specialized base sequences, but the flexibility of the sugar-phosphate backbone is expected to play a key role in nucleic acid polymorphism. Although x-ray crystallography and NMR spectroscopy can provide detailed information about the structural parameters of the nucleic acid phosphodiester backbone, the precise information on conformations of phosphodiester linkages in many biologically important nucleic acids, remain beyond the scope of the above methods. On the basis of quantum chemical methods, extensive studies have been carried out on phosphodiester model compounds. The dimethyl phosphate anion [(CH3O)2PO2-, DMP] serves as the simplest structural analogue for the nucleic acid phosphate group. Raman and IR spectra of DMP were employed to develop a comprehensive molecular force field for the phosphate group [1] and to explore the conformational dependence of C-O-P-O-C skeletal vibrations. Changes in the phosphodiester torsions of DMP, which correspond to the nucleic acid backbone torsions α and ζ, generate the most significant frequency shifts in Raman bands and have thus provided an initial basis for correlating the empirical conformation marker bands of nucleic acids with specific internal coordinates of the phosphodiester moiety. The another molecule, diethylphospate anion [(CH3CH2O)2PO2-, DEP] is also considered here as an improved structural analogue for the nucleic acid phosphate group. The back bone molecular skeleton, C-C-O-P(O2-)-O-C-C, incorporates two C-C bonds capable of vibrational interaction with modes of the phosphodiester linkages and two C-C bonds with torsion angles analogous to the nucleic acid backbone torsions β(O5’-C5’) and ε(C3’-O3’) [2].
It is always interesting and important to study the water interactions with real building blocks of DNA. Since dimethylphosphate and diethylphosphate anions are the analogue systems of nucleotides, the water interaction with these systems will help to model the real building blocks of DNA molecules. To make the study more effective, we have considered the interaction of single, double, and triple water molecules with the DMP and DEP anions. We have not imposed any constraints for these complex molecules and optimized them using the supra molecular approach. The interaction of single water molecule can be possible either O4 or O5 atom of PO4 group which is the most powerful proton acceptor, where the water molecules moves towards the PO4 group and makes very strong interaction. Even di and tri hydrates are having strong interactions with PO4 group. So, all water molecules have a very strong interaction with the PO4 group of DMP and DEP anions. The present study helps to understand the most possible sites for the interaction of water molecules with the DMP and DEP anions. The interaction of mono-, di- and tri-hydrates with PO4 group will elucidate the structural deformation of DMP and DEP anions. The HF, Post Hartree-Fock and density functional theory methods have been used for the above study. Even though the appreciable changes have been observed in the structural parameters of the PO4 group, the water interactions occurred only with the PO2- anion.

2. Computational methods

The ab initio and density functional theory methods have been used to optimize the geometries of these complexes. The restricted Hartree-Fock (RHF), Møller-Plesset perturbation theory (MP2) [3] of ab initio method and Becke’s three parameter exact exchange functional (B3) [4] combined with gradient corrected correlation functional of Lee- Yang- Parr (LYP) [5] and Perdew and Wang’s 1991 (PW91) [6] functional of DFT have been employed to optimize the isolated and complex molecules by implementing 6-31G* basis set. Boys and Bernardi’s counter poise procedure [7] have been applied to correct the basis set superposition error for the interaction energy:
Eint(corr)=E(AB)-[EA(AB)+EB(AB)]
where E(AB), EA(AB), and EB(AB) are the energies of complex and monomer structures in complex geometry. The chemical hardness (η) and chemical potential (μ) have been calculated using the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) energies determined by the HF/6-31G* method. The chemical hardness and chemical potential can be expressed within orbital basis as
η = I A 2       μ = I + A 2
where I=-EHOMO and A=-ELUMO: I and A are the ionization potential and electron affinity of a molecule. All these calculations were performed using the Gaussian 98W program [8].

3. Results and discussion

3.1 DMP – Water

The optimized geometrical parameters, rotational constants and dipole moment of the DMP-anion are given in Table 1 and shown in Figure 1. The optimized structures are compared with the crystal structure of NH4DMP molecule [9]. Though the optimized C-O bond length is shorter than the experimental value of NH4DMP in all the four levels of theory, the value calculated at MP2/6-31G* level of theory is in good agreement with the crystal data. The calculated P-O2 and P-O6 bond lengths are longer by 0.1 Å in all the four levels of theory and the P-O4 and P-O5 bond lengths are in excellent agreement with the experimental values. Single water molecules can interact with PO4 group in three different ways. The optimized structural parameters for the single water isomeric form with DMP anion are given in Table 2. The O4 or O5 atoms of the PO4 group can interact with H atom of the water molecule and H atoms of water molecules simultaneously interact with O4 and O5 atoms. The interaction induces the change in the P-O2 and P-O4 bond distances. There is a decrease in the bond distance of P-O2 atoms and increase in the P-O6 bond distance, when the water molecule interacts with O4 atom. Suppose, if water molecule interacts with O5 atom, there is not much change in the P-O2 distance, where as P-O6 bond distance is decreasing. Another observation is that the bond distance P-O4 is increasing while the water interaction is with O4 atom and P-O5 bond distance is decreasing. The same trends have been found while the water interaction is with O5 atom. All the P-O bond distances are same when the water molecule interacts symmetrically with the O4 and O5 atoms of PO4. Recently many authors [10,11,12] have explained the charge transfer process during the hydrogen bond formation. PO4 group is the proton acceptor and O-H bond of water molecule is proton donor group. So the charge is transferred from O4 atom to H atom of water molecule. This induces the electronic
Table 1. Optimized geometrical parameters (bond length in Å, bond angle in degrees), rotational constants RA, RB, RC(in GHz), dipole moment µM (in Debye), chemical potential μ (in eV) and chemical hardness η (in eV) for DMP and DEP anions.
Table 1. Optimized geometrical parameters (bond length in Å, bond angle in degrees), rotational constants RA, RB, RC(in GHz), dipole moment µM (in Debye), chemical potential μ (in eV) and chemical hardness η (in eV) for DMP and DEP anions.
ParametersDMP DEP
HFMP2B3LYPB3PW91Crystal data (NH4DMP) HFMP2B3LYPB3PW91Crystal data (BaDEP)
R(P3-O4)1.4701.5011.4981.4961.498 1.4721.5031.5001.4981.520
R(P3-O6)1.6421.6841.6841.6801.582 1.6421.6841.6841.6791.620
R(C1-O2)1.3931.4211.4121.4061.425 1.3981.4251.4181.4121.440
R(C1-C9) 1.5211.5221.5301.5251.510
θ(O4-P3-O5)124.8126.1125.7125.7117.2 124.7126.1125.6125.7121.6
θ(O2-P3-O6)99.398.498.698.6104.8 108.399.599.7108.8103.5
θ(P3-O2-C1)118.5114.4116.2115.8118.3 121.2116.9118.7118.3116.0
θ(P3-C1-C9) 112.6110.7111.5112.2114.0
μM5.76025.74565.03695.0128 6.20066.15855.62165.5709
RA2.96852.81052.82792.8588 2.13452.03391.99722.0300
RB2.19822.22512.17382.1795 1.06951.10791.08311.0820
RC1.76321.75791.73291.7460 0.87050.88790.85850.8638
η8.35 8.16
μ2.28 1.97
Figure 1. Structures of the DMP and DEP anions.
Figure 1. Structures of the DMP and DEP anions.
Ijms 04 00486 g001
delocalization in the PO2- group, so to compensate this, the charge is transferred from O2 and O6 to PO2- group. Due to this process, the structural reorganization has occurred in the PO4 and thus the P-O distances are found to be different where the water interaction takes place. The other bond length changes are found to be minimum. The various possible interactions of two water molecules with the DMP anion have been studied. The three isomeric forms of interaction have been identified for the two water molecule from the optimization of the complex using the supramolecular approach. The optimized structural parameters are tabulated in the Table 3 and the various interactions of water with DMP anion are shown in Fig. 2. The DMP-2W1 figure shows that both the water molecules point towards the same side. Here, the P-O4 bond distance is found to be more than the P-O5 bond distance by 0.001 Å and again P-O6 bond distance is more by 0.02 Å than the bond distance between P-O2 atoms. The reason for the above differences may be due to the possible dispersion force of second water molecule with O4 and O6 atoms. In the DMP-2W2, the interactions of water molecules are
Table 2. Optimized geometrical parameters (bond length in Å, bond angle in degrees), rotational constants RA, RB, RC(in GHz), dipole moment µM (in Debye), chemical potential μ (in eV) and chemical hardness η (in eV) for DMP-1water complexes.
Table 2. Optimized geometrical parameters (bond length in Å, bond angle in degrees), rotational constants RA, RB, RC(in GHz), dipole moment µM (in Debye), chemical potential μ (in eV) and chemical hardness η (in eV) for DMP-1water complexes.
ParametersDMP-1W1DMP-1W2DMP-1W3
HFMP2B3LYPB3PW91HFMP2B3LYPB3PW91HFMP2B3LYPB3PW91
R(P-O2)1.6311.6691.6681.6631.6441.6841.691.691.6321.6711.6711.666
R(P-O4)1.4781.5091.5071.5051.4871.4971.4941.4921.4751.5061.5031.501
R(P-O5)1.4671.4971.4941.4921.4781.5051.5071.5091.4751.5061.5031.501
R(P-O6)1.6441.691.691.6841.6311.6691.6681.6631.6321.6711.6711.666
R(C1-O2)1.3991.4291.4221.4161.3971.411.4161.4241.3961.4241.4161.41
R(C7-O6)1.3971.4241.4161.411.3991.4291.4221.4161.3961.4241.4161.41
θ(Ο4−Π3−Ο5)108.5125.7125.2125.2124.6125.2125.2125.7106.7106106105.9
θ(Ο2−Π3−Ο6)99.89999.499.599.899.599.499100.299.7100.1100.1
θ(Π3−Ο6−Η7)118.6114.6116.5116.2118.7114.9117116.4118.8114.7116.6116.2
θ(Π3−Ο2−Χ1)118.8114.9111116.5118.5116.2116.5114.6118.8114.7116.7116.3
θ(Ο14−Η15−Ο4)154.9153.4155.5155.7
θ(Η16−Ο14−Η15)100.497.997.897.8
θ(Η16−Ο14−Η15)
θ(Ο5−Η16−Ο14)100.497.897.897.9
θ(Ο15−Η14−Ο4)154.7155.7155.4153.9
θ(Η16−Ο15−Η14)145.1147.6147.8147.5
R(H15…O4)100.198.197.897.8
R(H16…O5)1.951.931.8761.8631.9531.8621.8771.926
R(H14…O4)
Eint(1W1)2.1572.0812.062.052
Eint(1W2)17.0921.220.6819.35
Eint(1W3)17.0921.220.6819.35
µM17.9820.619.9618.95
RA5.37575.1114.49914.51565.37155.12094.49354.5137.11577.19736.61916.6145
RB1.81131.83251.81781.83181.81091.83261.81781.83152.17532.17522.1222.1308
RC1.3931.40481.41251.42141.39341.40421.41281.42181.25461.24611.25861.2716
η1.20161.23521.22841.23591.20191.23521.22811.23570.96990.97640.97330.9822
μ8.398.398.5
1.71.71.76
Table 3. Optimized geometrical parameters (bond length in Å, bond angle in degrees), rotational constants RA, RB, RC(in GHz), dipole moment µM (in Debye), chemical potential μ (in eV) and chemical hardness η (in eV) for DMP-2 water complexes.
Table 3. Optimized geometrical parameters (bond length in Å, bond angle in degrees), rotational constants RA, RB, RC(in GHz), dipole moment µM (in Debye), chemical potential μ (in eV) and chemical hardness η (in eV) for DMP-2 water complexes.
ParametersDMP-2W1DMP-2W2DMP-2W3
HFMP2B3LYPB3PW91HFMP2B3LYPB3PW91HFMP2B3LYPB3PW91
R(P-O4)1.4811.5131.5121.5091.4711.5021.5001.4971.4741.5051.5031.501
R(P-O5)1.4711.5021.5001.4971.4811.5131.5121.5091.4741.5051.5031.501
R(P-O6)1.6351.6781.6771.6711.6221.6581.6561.6521.6341.6751.6731.668
R(P-O2)1.6221.6581.6561.6521.6351.6781.6771.6721.6341.6751.6731.668
R(C1-O2)1.4021.4331.4271.4201.4281.4281.4201.4141.4031.4331.4261.420
R(C7-O6)1.4001.4281.4201.4141.4021.4331.4261.4201.4031.4331.4261.420
θ(O4-P3-O5)122.6123.8123.2123.3122.6123.8123.2123.3124.4125.5124.9124.9
θ(O2-P3-O6)100.6100.3100.8100.8100.6100.3100.8100.8100.299.6100.0100.0
θ(P3-O2-C1)119.1115.2117.4116.9118.8114.8116.9116.5118.8115.0117.1116.6
θ(P3-O6-C7)116.0114.8116.9116.5119.1115.3117.5116.9118.8115.0117.1116.6
θ(O15-H14-O4)143.0145.9145.3145.2146.3148.0149.0149.1
θ(O17-H14-O4)153.7152.0154.0154.7
θ(H18-O17-H19)100.898.498.398.4100.898.498.398.3100.898.498.398.2
θ(H14-O15-H16)100.698.698.498.4100.698.698.498.4100.898.498.398.2
θ(O17-H19-O5)153.6152.0153.9154.2
θ(O4-H18-O17)154.2152.6154.3154.9
θ(O14-H16-O5)154.2152.5154.3154.9
R(O4…H18)1.9801.9581.9071.8891.9711.9541.9061.887
R(H16…O5)2.1412.0842.0482.0342.2162.1182.1152.1051.9711.9541.9061.887
R(O4…H14)2.2162.1182.1152.1082.1422.0852.0482.036
R(O5…H19)1.9821.9581.9071.894
Eint(H14O15H16)16.4718.9218.1217.09
Eint(O17 H18H19)15.6119.5018.8717.47
Eint(H14O15H16)16.4718.9218.1217.10
Eint(O17 H18H19)15.6119.5018.1617.48
Eint(O14H15H16)15.9419.8119.2918.00
Eint(O17 H18H19)15.9419.8119.2918.00
µM6.70306.64566.08846.14016.70006.64476.09176.11754.59644.01063.41103.4750
RA1.26451.32111.31791.31931.26491.32251.31541.32201.74101.69461.66991.6933
RB0.92220.92520.92760.93530.92210.92530.92700.93520.84710.88490.89770.8998
RC0.82110.82600.82640.83450.82140.82490.82810.83480.76970.81600.81930.8228
η8.538.538.47
μ1.221.221.21
Figure 2. Structures of the DMP anions interacting with one, two, and three water molecules.
Figure 2. Structures of the DMP anions interacting with one, two, and three water molecules.
Ijms 04 00486 g002
symmetrically opposite to each other. But one water molecule is near to the methyl group of DMP anion. This may be the reason that P-O5 bond distance is more by 0.01 Å than the P-O4 bond distance and P-O2 bond distance is more by 0.02 Å than the P-O6 bond distance. The one more isomeric form of DMP-2W3 is that two water molecules are symmetrically opposite to each other and interacts with O4 and O5 atoms of PO4 group. It is noted that P-O4 and P-O5 bond distances are found to be same and other two bonds P-O2 and P-O6 are also found to be the same. It clearly shows that the two water molecules do not have any Van der Waals interactions with the nearby groups of the DMP anions.
Finally, we have considered the interaction of three water molecules with PO4 group. Two hydrogen atoms of one water molecule symmetrically interact with the O4 and O5 atoms and other two water molecules, again interact with O4 and O5 atoms. It is interesting to note that P-O4, P-O5 distances are same and P-O2 and P-O6 bond distances are also found to be the same. Since all the water molecular interactions are symmetrical, the charge transfers from PO4 group to water molecules are uniform, so there is no change in the structural parameters of the DMP anion. As the bonding in the hydrogen atoms could not be accurately located in the x-ray structure, their optimized parameters are compared with standard experimental values of similar model compounds. The bond angles O-P-O and P-O-C are smaller compared to the corresponding experimental values in all the levels of theory except for P-O-C bond angle at HF level, which is very close to the experimental value.
The interaction energies of water molecule with DMP anion are given in Table 2, Table 3, Table 4. The interaction energy is corrected for BSSE using full counterpoise method of Boys and Bernardi [7]. The interaction energies of single water molecule is same when the interaction is O4 or O5 atoms. The hydrogen bond lengths are also found to be the same. When the two water molecules interact with the O4 and O5 atoms, the interaction energies are found to be slightly smaller than the single water interaction. The hydrogen bond lengths are found to be slightly higher than the single water interaction. If the three water molecules are having interaction with PO4 group, again the interaction energies are found to be slightly less and the hydrogen bond lengths are found to be slightly higher. The reason for the above variation may be due to the charge transfer from PO2- to OH of water molecule. If, more water molecules are interacting, then the charge transfer becomes less, due to the equal distribution of charge in the PO2- anion. The HF/6-31G* level of theory underestimated the interaction energy. But MP2/6-31G*, B3LYP/6-31G* and B3PW91/6-31G* levels of theory predicted the interaction energies and hydrogen bond lengths very reasonably.

3.2 DEP-Water Complex

In accordance with the x-ray crystal structure of BaDEP [13], the DEP anion and DEP-water have been fully optimized at HF, MP2, B3LYP, and B3PW91 levels of theory using 6-31G* basis set and are shown in Table 1, Table 5, Table 6, Table 7. Figure 3 shows various interactions of water with DEP anion. The calculated structural parameters are compared with those obtained by x-ray crystallography of BaDEP. In the DEP anion, the optimized P-O bond lengths are slightly larger than the corresponding values in the x-ray structure of BaDEP. Conversely, the optimized P-O4 and C-O bonds are shorter than the values obtained from the x-ray methods. The C-C bond length is in good agreement with the experimental values. The optimized O2-P-O6 bond angle is smaller and P-O-C bond angles on the other hand, are generally close to the experimental values. Similar differences can be attributed to the effect of metal cation in the crystal structure [1,14,15]. In the present study the interaction of single,
Table 4. Optimized geometrical parameters (bond length in Å, bond angle in degrees), rotational constants RA, RB, RC(in GHz), dipole moment µM (in Debye), chemical potential μ (in eV) and chemical hardness η (in eV) for DMP-3 water complex.
Table 4. Optimized geometrical parameters (bond length in Å, bond angle in degrees), rotational constants RA, RB, RC(in GHz), dipole moment µM (in Debye), chemical potential μ (in eV) and chemical hardness η (in eV) for DMP-3 water complex.
ParametersHFMP2B3LYPB3PW91
R(P-O4)1.4781.5101.5081.505
R(P-O5)1.4781.5101.5081.505
R(P-O2)1.6251.6641.6621.657
R(P-O6)1.6251.6641.6621.657
R(C1-O2)1.4061.4371.4301.423
R(C7-O6)1.4061.4371.4301.423
θ(O5-P3-O4)122.5123.6123.0123.1
θ(O4-P3-O2)101.1100.9101.3101.3
θ(P3-O2-C1)119.1115.4117.5117.0
θ(P3-O6-C7)119.1115.4117.5117.0
θ(O15-H14-O4)144.3146.6146.7146.8
θ(O20-H21-O4)153.1151.2152.9153.8
θ(H18-O17-H19)101.298.898.898.8
θ(H14-O15-H16)101.099.198.998.9
θ(O17-H19-O5)153.1151.2152.9153.8
θ(H22-O20-H21)101.298.898.898.8
R(O4…H21)1.9991.9811.9351.915
R(O5…H16)2.1972.1172.0992.091
R(O4…H14)2.1972.1172.0992.091
R(O5…H19)1.9991.9811.9351.915
Eint(H14O15H16)15.1217.4816.5915.51
Eint(O17 H18H19)14.6118.3417.6816.27
Eint(O20 H21H22)14.6118.3417.6816.27
µM5.99855.63885.12855.2152
RA0.92090.90920.90510.9170
RB0.77040.81280.81560.8165
RC0.59350.61270.62040.6218
η8.60------
μ0.75------
double and triple water molecules with the PO4 group have been considered. The interaction of single water molecule either O4 or O5 atom induces the change in the structural parameters of PO4 group. When the O4 atom of PO4 group forms the hydrogen bond with H atom of water molecule, the P-O4 bond length is found to be increasing while P-O5 is decreasing. The bond distance P-O2 is decreasing while the other bond P-O6 is increasing. When the water molecule interacts with the O5 atom, the reverse trend has been observed. Two water molecules can interact with the PO2- in three different ways which have been shown in fig. In the first case, the P-O4 distance is larger than the P-O5 bond distance. The corresponding hydrogen bond lengths are 1.917 Å and 1.912 Å respectively at MP2/6-31G* level of theory. Both the water molecules pointed towards the P-O4 bond, so the latter bond distance is more. In the second case, both water molecules are in the opposite direction, so the P-O4
Table 5. Optimized geometrical parameters (bond length in Å, bond angle in degrees), rotational constants RA, RB, RC(in GHz), dipole moment µM (in Debye), chemical potential μ (in eV) and chemical hardness η (in eV) for DEP-1water complexes.
Table 5. Optimized geometrical parameters (bond length in Å, bond angle in degrees), rotational constants RA, RB, RC(in GHz), dipole moment µM (in Debye), chemical potential μ (in eV) and chemical hardness η (in eV) for DEP-1water complexes.
ParametersDEP-1W1DEP-1W2DEP-1W3
HFMP2B3LYPB3PW91HFMP2B3LYPB3PW91HFMP2B3LYPB3PW91
R(P-O2)1.631.6691.6681.6631.6431.6891.6881.6811.6311.6721.6721.666
R(P-O4)1.4791.5111.5091.5071.4691.4991.4961.4941.4761.5081.5051.503
R(P-O5)1.4681.4981.4951.4931.4781.511.5081.5061.4751.5071.5041.502
R(P-O6)1.6441.6911.691.6841.631.6691.6671.6631.6321.6721.671.665
R(C1-O2)1.4051.4361.4311.4241.4041.4311.4251.4191.4031.4321.4251.419
R(C7-O6)1.4031.431.4231.4171.4031.4311.4251.4191.4021.4291.4231.416
R(C1-C9)1.521.5181.5251.5211.521.5181.5261.5211.521.5181.5261.521
R(C7-C8)1.5231.5211.5281.5241.5231.5211.5281.5231.5231.5211.5291.524
θ(Ο4−Π3−Ο5)124.4125.6125125124.3125.5124.9124.9106.7105.7105.8105.8
θ(Ο6−Π3−Ο2)100.6100.1100.4100.5100.599.9100.1100.3100.9100.7101101.6
θ(Π3−Ο6−Χ7)121.1117118.8118.4121.9117.7119.7119.2121.5117.1119118.6
θ(Π3−Ο2−Χ1)120.2116118.2117.7119.7115.7117.7117.3120.1115.8118.1117.6
θ(Ο6−Χ7−Χ8)112.4111.8112.2112.1113112.9113.4113.4112.4112.1112.1112.2
θ(Ο2−Χ1−Χ9)110.8110.1110.9110.9111.1110.6111.4111.4111110.5111.4111.3
θ(Ο20−Η22&Ο4)
θ(Ο20−Η21&Ο5)155.4153.7155.8156.4
θ(Ο21−Η20&Ο4)155154.2156156.8
θ(Η21−Ο20−Η22)145.2147.9147.6147.7
θ(Η21−Ο20−Η22)100.597.997.997.9
θ(Η22−Ο21−Η20)100.497.997.897.8
R(H20…O4)100.298.297.997.9
R(H22…O4)2.1542.0782.0622.05
R(H21…O5)1.9371.9171.8661.849
1.9451.9121.8691.85
Eint(DEP-1W1)
Eint(DEP-1W2)17.0721.1824.119.05
Eint(DEP-1W3)17.2321.2720.7219.27
µM17.7220.219.4318.45
RA5.95675.69335.12985.14495.88015.70674.94114.96027.81657.95977.49137.476
RB1.18681.21381.21561.22381.21491.22871.26641.2461.1191.13281.13931.1387
RC0.85660.86040.83660.84160.86440.87820.85070.8580.95030.95750.9290.9408
η0.75030.79170.78020.78060.74670.77010.74690.74960.59450.60670.59290.5971
μ8.258.278.33
1.481.51.5
Table 6. Optimized geometrical parameters (bond length in Å, bond angle in degrees), rotational constants RA, RB, RC(in GHz), dipole moment µM (in Debye), chemical potential μ (in eV) and chemical hardness η (in eV) for DEP-2 water complexes.
Table 6. Optimized geometrical parameters (bond length in Å, bond angle in degrees), rotational constants RA, RB, RC(in GHz), dipole moment µM (in Debye), chemical potential μ (in eV) and chemical hardness η (in eV) for DEP-2 water complexes.
ParametersDEP-2W1DEP-2W2DEP-2W3
HFMP2B3LYPB3PW91HFMP2B3LYPB3PW91HFMP2B3LYPB3PW91
R(P-O2)1.6201.6581.6551.6511.6361.6391.6761.6311.6331.6751.6721.666
R(P-O4)1.4831.5151.5131.5111.4731.5051.5021.5001.4751.5061.5051.502
R(P-O5)1.4721.5031.5001.4981.4831.5221.5201.5171.4741.5051.5031.501
R(P-O6)1.6341.6791.6331.6711.6131.6471.6461.6421.6341.6761.6731.667
R(C1-O2)1.4091.4401.4351.4281.4091.4391.4331.4261.4091.4401.4341.427
R(C7-O6)1.4061.4341.4271.4211.4071.4371.4321.4251.4071.4371.4321.425
R(C1-C9)1.5191.5171.5241.5201.5231.5211.5291.5231.5191.5171.5241.520
R(C7-C8)1.5221.5201.5271.5231.5181.5151.5221.5171.5181.5161.5231.518
θ(O4-P3-O5)122.3123.6123.0123.0120.9121.3121.0121.1124.2125.2124.5124.5
θ(O4-H20-O21)143.2146.2145.7145.9145.8147.0147.3147.6
θ(O5-H22-O21)146.0147.7148.4148.4142.5146.1146.1146.3
θ(O5-H24-O23)154.6152.8154.2154.4
θ(H20-O21-H22)100.698.798.498.4100.598.398.198.0
θ(H25-O23-H24)100.998.498.498.5101.198.398.498.3101.098.498.598.5
θ(O2-P3-O6)101.4101.3101.5101.699.397.797.797.8100.399.5100.0100.1
θ(O23-H25-O4)154.4152.3154.3155.0
θ(O4-H25-O23)155.4153.2155.6156.5
θ(O5-H21-O20)154.3152.6154.7155.5
θ(H21-O20-H22)100.898.498.398.3
R(O4…H25)1.9641.9441.8951.833
R(H22…O5)2.1462.0912.0582.0432.1902.0652.0422.029
R(H24…O5)1.9841.9511.9161.899
R(H20…O4)2.2142.1162.1102.0992.1142.0622.0482.028
R(H25…O4)1.9491.9351.8801.860
R(H21…O5)1.9701.9521.9001.881
Eint(H20O21H22)17.6620.4519.0118.04
Eint(O23H24H25)15.8419.8618.0518.04
Eint(H20O21H22)15.9419.3019.8117.60
Eint(O23H24H25)15.2920.2419.8718.08
Eint(O20H21H22)15.8619.8119.2117.85
Eint(O23H24H25)15.9419.9119.1017.74
µM7.55677.46756.98376.99375.72305.08744.75554.79745.43295.00084.31574.3933
RA0.77200.79830.79610.79511.02431.08281.10001.10300.81200.85200.85730.8593
RB0.69820.72330.71810.72450.60440.63000.61460.61890.75310.75260.74290.7483
RC0.57850.58230.57090.57630.53480.53980.53620.54130.52840.56510.55430.5566
η8.248.358.29
μ1.031.131.07
Table 7. Optimized geometrical parameters (bond length in Å, bond angle in degrees), rotational constants RA, RB, RC(in GHz), dipole moment µM (in Debye), chemical potential μ (in eV) and chemical hardness η (in eV) for DEP-3 water complex.
Table 7. Optimized geometrical parameters (bond length in Å, bond angle in degrees), rotational constants RA, RB, RC(in GHz), dipole moment µM (in Debye), chemical potential μ (in eV) and chemical hardness η (in eV) for DEP-3 water complex.
ParametersHFMP2B3LYPB3PW91
R(P-O2)1.6241.6641.6601.655
R(P-O4)1.4791.5111.5091.507
R(P-O5)1.4781.5091.5081.505
R(P-O6)1.6251.6651.6611.656
R(C1-O2)1.4131.4441.4391.432
R(C7-O6)1.4111.4411.4361.430
R(C1-C9)1.5181.5161.5231.519
R(C7-C8)1.5171.5151.5211.517
θ(O4-H20-O21)144.2146.6146.7146.7
θ(H20-O21-H22)101.099.098.998.9
θ(O4-H28-O26)154.2152.0154.4155.3
θ(H28-O26-H27)101.398.998.999.0
θ(O5-H22-O21)144.3146.6146.7146.7
θ(H24-O23-H25)101.298.798.798.8
θ(O5-H24-O23)153.1151.0152.9153.5
R(O4…H20)2.2022.1182.1002.091
R(O4…H28)1.9771.9591.9021.884
R(O5…H24)1.9971.9801.9301.913
R(O5…H22)2.1942.1172.0952.086
Eint(H20O21H22)15.0117.2916.3215.25
Eint(O23H24 H25)14.5818.4017.5616.17
Eint(O26H27 H28)14.6718.5117.5816.20
µM7.70356.83026.22276.3122
RA0.57220.60530.60430.6061
RB0.55150.56260.55280.5573
RC0.49220.50050.50100.5043
η8.44
μ0.59
and P-O5 bond distances are found to be same. Moreover the hydrogen bond lengths are also almost the same. In the third isomeric form, one water molecule pointed towards the ethyl form. Here, again the P-O4 and P-O5 bond distance are the same, but the hydrogen bond lengths are different. In the last case, three water molecules are interacting with the PO2-, the bond lengths are almost same. Both the ab initio and DFT methods were produced the similar trends of result. The B3LYP/6-31G* and B3PW91/6-31G* levels of theory produced the structural parameters, which are very well coincide with the MP2/6-31G* level of theory. The bond angle O4-P3-O5 in DEP- water molecule is decreased by 2o compared to that of DEP anion due to the interaction of hydrogen atoms in the molecule and the oxygen atoms of phosphodioxy group. The same amount of decrease in the bond angle is also experienced in the DMP-water complex. Similarly all the P-O bond distances of DEP-water are slightly
Figure 3. Structures of the DEP anions interacting with one, two, and three water molecules
Figure 3. Structures of the DEP anions interacting with one, two, and three water molecules
Ijms 04 00486 g003
shorter than those of DEP anion. The interaction energies and hydrogen bond lengths of DEP-water are given in Table 5, Table 6 and Table 7. The interaction energies of mono-, di-, tri-hydrates with DEP anion is almost found to be the same. The hydrogen bond lengths are also approximately equal to 2 Å. For all the DEP-water interactions the HF/6-31G* level of theory underestimated the interaction energies, whereas B3LYP/6-31G*, B3PW91/6-31G* and MP2/6-31G* levels of theory predicted the interaction energies and hydrogen bond lengths reasonably.
The chemical hardness and chemical potential have been calculated for these complexes at HF/6-31G* level of theory and are tabulated in Table 1, Table 2, Table 3, Table 4, Table 5, Table 6 and Table 7. In both the anions, the chemical hardness values are found to be increased after the formation of the complexes. This result indicates that the system tends to stabilized after the formation of the complex. For the past few years, the chemical hardness is one of the important parameters to study the stability of the molecules. In the exothermic exchange reactions of the AB+CD→AC+BD type would produce the hardest possible species in one of the products [16]. Chattaraj and Schleyer [17] have concluded that in general this is not true for acid base reactions leading to binary complexes. Nath et el., [18] have shown the hardness profiles of hydrogen bond complexes, that the formation of H3N…HF complex is exothermic to the extent of 11.82 kcal/mol and is accompanied by an increase in η of 0.28 eV for proton-transfer reaction (F-H…Cl)- → (F…H-Cl)-. Arulmozhiraja and Kolandaivel, who studied [19] hardness profiles for hydrogen-bonded system, indicate that hydrogen-bonding interactions do not create higher hardness complex species in the product side. But in the present case for both the complexes DMP-water and DEP-water, the chemical hardness values are found to be higher than the monomer DMP and DEP anions. The chemical hardness value is found to be more for tri hydrates interaction with DEP and DMP anions. While forming the hydrogen bond in complexes, the change of bond length and bond angles of PO4 and the more charge transfer from O atoms induce the change in the electronic chemical potential, which in turn increases the chemical hardness value. So the proton transfer reactions, the minimum energy structure of the complex will have the maximum hardness, which is in agreement with the other studies also. In a similar way the dipole moment of the complexes is increased compared to that of the corresponding isolated anion. The earlier studies also show that dipole moment of the dimer is enhanced during the formation of hydrogen bonds in complexes [20]. Both the DMP and DEP water interactions, PO2- anion act as proton acceptor group. Moreover PO2- group is more polar. The charge from the σ orbital of the P-O bond is transferred to OH group of water molecule. Due to the charge transfer, the electron delocalization has been occurred in the PO2- anion, to compensate this, the other charges are reoriented and structural reorganization has been occurred in the PO4 group of DMP and DEP anions. So the water interaction induces the reorientation in the PO4 group rather than the other parts of the molecules.

4. Conclusion

The anionic forms of dimethyl phosphate and diethyl phosphate are optimized at HF, MP2, B3LYP and B3PW91 levels of theory using 6-31G* basis set. The interaction energies of isomeric forms of mono-, di-, and tri-hydrates with the above anions are further studied. The optimized structural parameters have been compared with the available experimental results and they are in good agreement with each other. The effect of charge transfer within these complexes makes the structural reorganization in the PO4 group of DMP and DEP anions. The chemical hardness values have been calculated for these complexes and it is found that interaction of more number of water molecules increase the stability of the complex.

Acknowledgement

One of the authors (R.K.) expresses his sincere thanks to the UGC, New Delhi, for the award of FIP fellowship during the IX plan period and he is thankful to the Management and the Principal, Nallamuthu Gounder Mahalingam College, Pollachi, for allowing him to undergo this program. P.K. is thankful to DST, Government of India for the financial support in the form of Research Project.

References

  1. Guan, Y.; Choy, G.; Glaser, R.; Thomas, G. J., Jr. J. Phys. Chem. 1995, 99, 12054. [CrossRef]
  2. Guan, Y.; Thomas, G. J., Jr. Biopolymers. 1996, 39, 813.
  3. Møller, C.; Plesset, M.S. Phys. Rev. 1934, 46, 618. [CrossRef]
  4. Becke, A.D. Phys. Rev. A. 1988, 38, 3098. [CrossRef]
  5. Lee, C.; Yang, W.; Parr, R.G. Phys. Rev. B. 1988, 37, 785. [CrossRef]
  6. Perdew, J.P.; Wang, Y. Phys. Rev. B. 1992, 45, 244. [CrossRef]
  7. Boys, S.F.; Bernardi, F. Mol. Phys. 1970, 19, 553.
  8. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Zakrzewski, V.G.; Montgomery, J.A., Jr.; Stratmann, R.E.; Burant, J.C.; S. Dapprich, J. M.; Millam, A. D.; Daniels, K. N.; Kudin, M. C.; Strain, Farkas, O.; Tomasi, J.; Barone, V.; Cossi, M.; Cammi, R.; Mennucci, B.; Pomelli, C.; Adamo, C.; Clifford, S.; Ochterski, J.; Petersson, G.A.; Ayala, P.Y.; Cui, Q.; Morokuma, K.; Rega, N.; Salvador, P.; Dannenberg, J.J.; Malick, D.K.; Rabuck, A.D.; Raghavachari, K.; Foresman, J.B.; Cioslowski, J.; Ortiz, J.V.; Baboul, A.G.; Stefanov, B.B.; Liu, G.; Liashenko, A.; Piskorz, P.; Komaromi, I.; Gomperts, R.L.; Martin, R.L.; Fox, D.J.; Keith, T.; Al-Laham, M.A.; Peng, C.Y.; Nanayakkara, A.; Challacombe, M.; Gill, P.M.W.; Johnson, B.; Chen, W.; Wong, M.W.; Andres, J.L.; Gonzalez, C.; Head-Gordon, M.; Replogle, E.S.; Pople, J.A. Gaussian98, Revision A.11.2; Gaussian, Inc.: Pittsburgh, PA, 2001. [Google Scholar]
  9. Giarda, L.; Garbassi, F.; Calcaterra, M. Acta Crystallogr. 1973, B29, 1826.
  10. Hobza, P; Havlas, Z. Chem. Rev. 2000, 100, 4253.
  11. Hobza, P. Phys. Chem. Chem. Phys. 2001, 3, 2555. [CrossRef]
  12. Hermansson, K. J. Phys. Chem. A 2002, 106, 4695. [CrossRef]
  13. Kyogoku, Y.; Iitaka, Y. Acta Cryst. 1966, 21, 49.
  14. Hadzi, D.; Hodoscek, M.; Grdadolnik, J.; Avbelj, F. J.Mol.Struct. 1992, 266, 9. [CrossRef]
  15. Liang, C.; Ewig, C.S.; Stouch, T.R.; Hagler, A.T. J. Am. Chem. Soc. 1993, 115, 1537. [CrossRef]
  16. Datta, D. J. Phys. Chem. 1992, 96, 2400. [CrossRef]
  17. Chattaraj, P.K.; Schleyer, P.V.R. J. Am. Chem. Soc. 1994, 116, 1067. [CrossRef]
  18. Nath, S.; Sannigrahi, A.B.; Chattaraj, P.K. J. Mol. Struct. (Theochem). 1994, 309, 65. [CrossRef]
  19. Arulmozhiraja, S.; Kolandaivel, P. Int. J. Quant. Chem. 1996, 64, 221.
  20. Kanakaraju, R.; Kolandaivel, P.; Gowenlock, B.G. J. Mol. Struct. (Theochem) 2002, 577, 121. [CrossRef]

Share and Cite

MDPI and ACS Style

Kolandaivel, P.; Kanakaraju, R. Structure, Stability and Interaction Studies on Nucleotide Analogue Systems. Int. J. Mol. Sci. 2003, 4, 486-502. https://0-doi-org.brum.beds.ac.uk/10.3390/i4080486

AMA Style

Kolandaivel P, Kanakaraju R. Structure, Stability and Interaction Studies on Nucleotide Analogue Systems. International Journal of Molecular Sciences. 2003; 4(8):486-502. https://0-doi-org.brum.beds.ac.uk/10.3390/i4080486

Chicago/Turabian Style

Kolandaivel, P., and R. Kanakaraju. 2003. "Structure, Stability and Interaction Studies on Nucleotide Analogue Systems" International Journal of Molecular Sciences 4, no. 8: 486-502. https://0-doi-org.brum.beds.ac.uk/10.3390/i4080486

Article Metrics

Back to TopTop