The optimized electronic structure of Phosphorylation of Nucleosides and Nucleotides of deoxy ribonucleic acid (DNA) (a to c) are shown in Scheme- 1 to 4. The calculated values of Heat of formation (∆HfO in kcal/mol), dipole moment (µ in Debye), frontier molecular orbital energy (HOMO & LUMO) (Table-I), quantum chemical descriptors (Table-II), Calculation of Heat of Reaction (∆Hr) (Table-III), Induction effect(µind),Cosmo area (Sq.Ang) and Cosmo volume (Cu.Ang) (Table-IV) are included.
Phosphorylation of deoxyadenosine (I), as per Scheme – 1, and calculated values are incorporated in Table – I. It is observed the stability in the order of 3’5’-ADP > 3’AMP > 5’AMP > deoxyadenosine as per Heat of formation data. The dipole moment depends on the nature of the atoms and bonds comprising the molecules and on their arrangement. It is also found that the dipole moment has increased in the order of deoxyadenosine < 3’AMP < 5’AMP < 3’5’-ADP.
Table-I - IHeat of formation (∆HfO in kcal/mol), dipole moment (µ in Debye), frontier orbital energy (HOMO & LUMO), Cosmo area and Cosmo volume for nucleosides (I to IV) and nucleotides (a to c) of DNA from AM1 calculation
Comp
No.
|
DNA Nucleotides
|
ΔHfO (Kcal/mol)
|
Dipole moment (Debye)
|
HOMO (ev)
|
LUMO (ev)
|
Cosmo Area (Sq.Ang)
|
Cosmo volume
(Cu.Ang)
|
Mol wt
|
|
H3PO4
|
-238.1240
|
5.8959
|
-11.716
|
+ 0.598
|
100.39
|
82.58
|
97.9951
|
|
H2O
|
-59.2499
|
1.8602
|
-12.463
|
+ 4.412
|
42.61
|
26.03
|
18.0152
|
I
|
deoxyadenosine
|
+ 60.6205
|
4.1752
|
-7.176
|
-1.807
|
280.64
|
297.10
|
251.2444
|
Ia
|
3’-AMP
|
-170.8813
|
5.0619
|
-8.616
|
-0.492
|
351.76
|
362.93
|
331.2243
|
Ib
|
5’-AMP
|
-131.0727
|
5.3107
|
-8.184
|
-1.037
|
352.60
|
361.92
|
331.2243
|
Ic
|
3’,5’-ADP
|
-384.3956
|
9.3849
|
-6.338
|
-1.424
|
383.30
|
416.70
|
411.2041
|
II
|
deoxyguanosine
|
-24.7427
|
3.8273
|
-8.585
|
-0.742
|
305.74
|
319.85
|
267.2438
|
IIa
|
3’-GMP
|
-176.1358
|
6.0821
|
-8.341
|
-0.634
|
359.47
|
377.44
|
347.2237
|
IIb
|
5’-GMP
|
-284.8335
|
6.1391
|
-8.754
|
-0.744
|
363.77
|
382.68
|
347.2237
|
IIc
|
3’,5’-GDP
|
-480.0820
|
2.5209
|
-8.975
|
-0.808
|
401.50
|
435.94
|
427.2035
|
III
|
deoxycytidine
|
-52.6407
|
4.6162
|
-7.524
|
-2.021
|
272.69
|
280.78
|
227.2194
|
IIIa
|
3’-CMP
|
-272.7677
|
8.3389
|
-8.758
|
-0.748
|
309.41
|
335.30
|
307.1993
|
IIIb
|
5’-CMP
|
-257.1856
|
3.4799
|
-6.591
|
-1.091
|
313.53
|
331.79
|
307.1993
|
IIIc
|
3’,5’-CDP
|
-498.8033
|
7.1174
|
-9.135
|
-1.154
|
384.63
|
403.02
|
387.1791
|
IV
|
deoxythymidine
|
-112.9364
|
6.6805
|
-7.714
|
-2.058
|
263.41
|
291.83
|
242.2310
|
IVa
|
3’-TMP
|
-288.6573
|
8.4049
|
-8.731
|
-0.683
|
314.81
|
330.36
|
307.1993
|
IVb
|
5’-TMP
|
-257.1517
|
3.4015
|
-6.605
|
-1.097
|
313.51
|
333.14
|
307.1993
|
IVc
|
3’,5’-TDP
|
-498.3798
|
7.4077
|
-9.035
|
-1.051
|
387.14
|
401.64
|
387.1791
|
The heats of reactions are calculated and values incorporated in Table – III. It is investigated that Heat of Reaction (∆Hr) is in the order of 4 > 1 > 3 > 2 as per Scheme – 1. A reaction can occur spontaneously only if the change in free energy is negative. Free energy is a measure of capacity of a system to do useful work at the constant temperature and pressure.
Table- II - Quantum Chemical descriptors i.e. Energy gap or Electron excitation energy (ΔE, in eV), Ionization Potential (IP), Electron Affinity (EA), Electro Negativity(EN) (χ) , Global hardness (ƞ), Softness (s) , Chemical potential (μ0) and Electrophilicity Index (w) for nucleosides ( I to IV) and nucleotides (a to c) of DNA are calculated from frontier molecular orbital energy (HOMO & LUMO)
Comp. No.
|
DNA
Nucleotides
|
HOMO
eV
|
LUMO
eV
|
Energy gap (ΔE) eV
|
IP
eV
|
EA
eV
|
EN
eV
|
Global Hardness
|
Softness
|
Chemical Potential
|
Eletro philicity Index
|
|
H3PO4
|
-11.716
|
0.598
|
12.314
|
11.716
|
-0.598
|
5.559
|
6.157
|
1.903
|
-5.559
|
2.510
|
|
H2O
|
-12.463
|
4.412
|
16.875
|
12.463
|
-4.412
|
4.026
|
8.438
|
1.477
|
-4.026
|
0.960
|
I
|
deoxyadenosine
|
-7.176
|
-1.807
|
5.369
|
7.176
|
1.807
|
4.492
|
2.685
|
2.673
|
-4.492
|
3.757
|
a
|
3’-AMP
|
-8.616
|
-0.492
|
8.124
|
8.616
|
0.492
|
4.554
|
4.062
|
2.121
|
-4.554
|
2.553
|
b
|
5’-AMP
|
-8.184
|
-1.037
|
7.147
|
8.184
|
1.037
|
4.611
|
3.574
|
2.290
|
-4.611
|
2.974
|
c
|
3’,5’-ADP
|
-6.338
|
-1.424
|
4.914
|
6.338
|
1.424
|
3.881
|
2.457
|
2.580
|
-3.881
|
3.065
|
II
|
deoxyguanosine
|
-8.585
|
-0.742
|
7.843
|
8.585
|
0.742
|
4.664
|
3.922
|
2.189
|
-4.664
|
2.773
|
a
|
3’-GMP
|
-8.341
|
-0.634
|
7.707
|
8.341
|
0.634
|
4.488
|
3.854
|
2.165
|
-4.488
|
2.613
|
b
|
5’-GMP
|
-8.754
|
-0.744
|
8.010
|
8.754
|
0.744
|
4.749
|
4.005
|
2.186
|
-4.749
|
2.816
|
c
|
3’,5’-GDP
|
-8.975
|
-0.808
|
8.167
|
8.975
|
0.808
|
4.892
|
4.084
|
2.198
|
-4.892
|
2.930
|
III
|
deoxycitidine
|
-7.524
|
-2.021
|
5.503
|
7.524
|
2.021
|
4.773
|
2.752
|
2.735
|
-4.773
|
4.139
|
a
|
3’-CMP
|
-8.758
|
-0.748
|
8.010
|
8.758
|
0.748
|
4.753
|
4.005
|
2.187
|
-4.753
|
2.820
|
b
|
5’-CMP
|
-6.591
|
-1.091
|
5.500
|
6.591
|
1.091
|
3.841
|
2.750
|
2.397
|
-3.841
|
2.682
|
c
|
3’,5’-CDP
|
-9.135
|
-1.154
|
7.981
|
9.135
|
1.154
|
5.145
|
3.991
|
2.289
|
-5.145
|
3.316
|
IV
|
deoxythymidine
|
-7.714
|
-2.058
|
5.656
|
7.714
|
2.058
|
4.886
|
2.828
|
2.728
|
-4.886
|
4.221
|
a
|
3’-TMP
|
-8.731
|
-0.683
|
8.048
|
8.731
|
0.683
|
4.707
|
4.024
|
2.170
|
-4.707
|
2.753
|
b
|
5’-TMP
|
-6.605
|
-1.097
|
5.508
|
6.605
|
1.097
|
3.851
|
2.754
|
2.398
|
-3.851
|
2.692
|
c
|
3’,5’-TDP
|
-9.035
|
-1.051
|
7.984
|
9.035
|
1.051
|
5.043
|
3.992
|
2.263
|
-5.043
|
3.185
|
The study of heat is produced from the system to surroundings are termed as exothermic process and Heat of Reaction (∆Hr) is negative (-ve). When heat is observed by the system from the surroundings is termed as endothermic process and Heat of Reaction (∆Hr) is positive (+ ve). All reactions are exothermic and their Heat of Reactions (∆Hr) are negative (-ve).
Table - III - Calculation of Heat of Reaction (∆Hr) as per Schemes (1 to 4) from AM1 Method.
Reactions as per Schemes (1 to 4)
|
Heat of Reaction (∆Hr)
(k.cal/mol)
|
Scheme − 1
|
1
|
deoxyadenosine + H3PO4 → 3’-AMP + H2O
|
-52.6277
|
2
|
deoxyadenosine + H3PO4 → 5’-AMP + H2O
|
-12.8191
|
3
|
3’-AMP + H3PO4 → 3’5’-ADP + H2O
|
-34.6402
|
4
|
5’-AMP + H3PO4 → 3’5’-ADP + H2O
|
-74.4488
|
Scheme − 2
|
1
|
deoxyguanosine + H3PO4 → 3’-GMP + H2O
|
+ 27.4810
|
2
|
deoxyguanosine + H3PO4 → 5’-GMP + H2O
|
-81.2167
|
3
|
3’-GMP + H3PO4 → 3’5’-GDP + H2O
|
-125.0721
|
4
|
5’-GMP + H3PO4 → 3’5’-GDP + H2O
|
-16.3744
|
Scheme − 3
|
1
|
deoxycytidine + H3PO4 → 3’-CMP + H2O
|
-41.2529
|
2
|
deoxycytidine + H3PO4 → 5’-CMP + H2O
|
-25.6708
|
3
|
3’-CMP + H3PO4 → 3’5’-CDP + H2O
|
-47.1615
|
4
|
5’-CMP + H3PO4 → 3’5’-CDP + H2O
|
-62.7436
|
Scheme − 4
|
1
|
deoxythymidine + H3PO4 → 3’-TMP + H2O
|
+ 3.1532
|
2
|
deoxythymidine + H3PO4 → 5’-TMP + H2O
|
+ 34.6588
|
3
|
3’-TMP + H3PO4 → 3’5’-TDP + H2O
|
-30.8484
|
4
|
5’-TMP + H3PO4 → 3’5’-TDP + H2O
|
-62.3540
|
The strong kind of dipole-dipole attractions make hydrogen bonding in the formation of the double helix of DNA and these are important in determining the nature of the hydrogen bonding in nucleic acids4. It is expected that the dipole moment has increased in the conversion of deoxyadenosine (a to c) due to the induction effect. So that the magnitude of the induction effect11 (µind) can be estimated as per Eq. (1). The electronegative nitrogen atoms cause displacement of elections in the π-frame work to induce an additional dipole moment in the conversion of one to another. The order of induction effect is found and incorporated values in Table – IV. It is also observed in the case of Cosmo area has increased and Cosmo volume has decreased as 3’-AMP > 5’-AMP
Table-IV - Change of Induction effect, Cosmo area (Sq.Ang) and Cosmo volume (Cu.Ang) in the formation of nucleotides.
Conversions
|
Induction effect
(µind)
|
Cosmo area (Sq.Ang)
|
Cosmo volume (Cu.Ang)
|
Deoxyadenosine → 3’-AMP
|
0.8867
|
71.12
|
65.83
|
Deoxyadenosine → 5’-AMP
|
1.1355
|
71.96
|
64.82
|
Deoxyadenosine → 3’5’-ADP
|
5.2097
|
102.64
|
119.60
|
Deoxyguanosine → 3’-GMP
|
2.2548
|
53.73
|
57.59
|
Deoxyguanosine → 5’-GMP
|
2.3118
|
58.03
|
62.83
|
Deoxyguanosine → 3’5’-GDP
|
-1.3064
|
95.76
|
116.09
|
Deoxycytidine → 3’-CMP
|
3.7227
|
36.72
|
54.52
|
Deoxycytidine → 5’-CMP
|
-1.1363
|
40.84
|
51.01
|
Deoxycytidine → 3’5’-CDP
|
2.5012
|
111.94
|
122.24
|
Deoxythymidine → 3’-TMP
|
1.7245
|
51.10
|
38.53
|
Deoxythymidine → 5’-TMP
|
-3.2789
|
50.10
|
41.31
|
Deoxythymidine → 3’5’-TDP
|
0.7273
|
123.73
|
109.81
|
Phosphorylation of deoxyguanosine (II) as per Scheme – 2, and calculated values are incorporated in Table – I. It is observed the stability in the order of 3’5’-GDP > 5’-GMP > 3’-GMP > deoxyguanosine as per Heat of formation data. The dipole moment depends on the nature of the atoms and bonds comprising the molecules and on their arrangement. It is also found that the dipole moment has increased in the order of 3’5’-GDP < deoxyguanosine < 3’-GMP < 5’-GMP. The heats of reactions are calculated and values incorporated in Table – III. It is investigated that Heat of Reaction (∆Hr) is in the order of 3 > 2 > 4 > 1 > as per Scheme – 2.
The dipole moment has increased in the conversion of deoxyguanosine (a to c) due to the induction effect. So that the magnitude of the induction effect11 (µind) can be estimated as per Eq. (1). The order of induction effect is found and incorporated values in Table – IV. It is also observed that Cosmo area and Cosmo volume have increased as 3’-GMP < 5’-GMP. The study of heat is produced from the system to surroundings are termed as exothermic process and Heat of Reaction (∆Hr) is negative (-ve). When heat is observed by the system from the surroundings is termed as endothermic process and Heat of Reaction (∆Hr) is positive (+ ve). All reactions are exothermic and their Heat of Reactions (∆Hr) are negative (-ve) except Scheme – 2, reaction − 1 is endothermic process as in Table – III.
Phosphorylation of deoxycytidine (III) as per Scheme – 3, and calculated values are incorporated in Table – I. It is observed the stability in the order of 3’5’-CDP > 3’-CMP > 5’-CMP > deoxycytidine as per Heat of formation data. The dipole moment depends on the nature of the atoms and bonds comprising the molecules and on their arrangement. It is also found that the dipole moment has in the order of 5’-CMP < 3’5’-CDP < deoxycytidine < 3’-CMP. The heats of reactions are calculated and values incorporated in Table – III. It is investigated that Heat of Reaction (∆Hr) is in the order of 4 > 3 > 1 > 2 as per Scheme – 3. All reactions are exothermic and their Heat of Reactions (∆Hr) are negative (-ve).
The dipole moment has increased in the conversion of deoxycytidine (a to c) due to the induction effect. So that the magnitude of the induction effect (µind) can be estimated as per Eq. (1). The order of induction effect is found and incorporated values in Table – IV. It is also observed that Cosmo area has increased (3’-CMP < 5’-GMP) and Cosmo volume have decreased as 3’-CMP > 5’-GMP.
Phosphorylation of deoxythymidine (IV) as per Scheme – 4, and calculated values are incorporated in Table – I. According to the heat of formation (∆Hfo) data, the stability is found in the increasing order of 3’5’-TDP > 3’-TMP > 5’-TMP > deoxythymidine as per Heat of formation data. The dipole moment depends on the nature of the atoms and bonds comprising the molecules and on their arrangement. It is also found that the dipole moment has increased in the order of 5’-TMP < deoxythymidine < 3’5’-TDP < 3’-TMP.
The heats of reactions are calculated and values incorporated in Table – III. It is investigated that Heat of Reaction (∆Hr) is in the order of 4 > 3 > 1 > 2 as per Scheme – 4. The dipole moment has increased in the conversion of deoxythymidine (a to c) due to the induction effect. So that the magnitude of the induction effect (µind) can be estimated as per Eq. (1). The order of induction effect is found and incorporated values in Table – IV. It is also observed that Cosmo area has decreased (3’-TMP > 5’-TMP) and Cosmo volume increased (3’-TMP < 5’-TMP). The study of heat is produced from the system to surroundings are termed as exothermic process and Heat of Reaction (∆Hr) is negative (-ve). When heat is observed by the system from the surroundings is termed as endothermic process and Heat of Reaction (∆Hr) is positive (+ ve). As per Scheme – 2, reactions 1 and 2 are endothermic process and reactions 3 and 4 are exothermic process as in Table – III.