Journal of Chemistry

Journal of Chemistry / 2013 / Article

Research Article | Open Access

Volume 2013 |Article ID 126502 |

G. R. Ramkumaar, S. Srinivasan, T. J. Bhoopathy, S. Gunasekaran, "Vibrational Spectroscopic Studies of Tenofovir Using Density Functional Theory Method", Journal of Chemistry, vol. 2013, Article ID 126502, 12 pages, 2013.

Vibrational Spectroscopic Studies of Tenofovir Using Density Functional Theory Method

Academic Editor: Isabel Seiquer
Received25 Jun 2012
Revised23 Aug 2012
Accepted03 Sep 2012
Published24 Oct 2012


A systematic vibrational spectroscopic assignment and analysis of tenofovir has been carried out by using FTIR and FT-Raman spectral data. The vibrational analysis was aided by electronic structure calculations—hybrid density functional methods (B3LYP/6-311++G(d,p), B3LYP/6-31G(d,p), and B3PW91/6-31G(d,p). Molecular equilibrium geometries, electronic energies, IR intensities, and harmonic vibrational frequencies have been computed. The assignments proposed based on the experimental IR and Raman spectra have been reviewed and complete assignment of the observed spectra have been proposed. UV-visible spectrum of the compound was also recorded and the electronic properties such as HOMO and LUMO energies and were determined by time-dependent DFT (TD-DFT) method. The geometrical, thermodynamical parameters, and absorption wavelengths were compared with the experimental data. The B3LYP/6-311++G(d,p)-, B3LYP/6-31G(d,p)-, and B3PW91/6-31G(d,p)-based NMR calculation procedure was also done. It was used to assign the 13C and 1H NMR chemical shift of tenofovir.

1. Introduction

A recent addition to the antiretroviral armamentarium is the nucleotide analogue tenofovir disoproxil fumarate (Viread), approved for use in the USA and the European Union. Tenofovir is unique among the NRTIs in that it is an acyclic nucleoside phosphonate, analogous to the monophosphate form of the other NRTIs [1]. Tenofovir disoproxil fumarate is an oral prodrug of tenofovir that is rapidly converted into tenofovir upon absorption [2, 3]. Tenofovir has activity in vitro against both HIV-1 and HIV-2 [4, 5], and in resting and activated T cells, monocytes and macrophages [5, 6]. Cross-resistance within the NRTI class of drugs has important clinical consequences for patients who are highly treatment experienced, or for those patients in whom primary HIV infection is associated with the transmission of a resistant virus [7].

Literature survey reveals that to the best of our knowledge, neither the complete IR, Raman spectra, nor the quantum mechanical calculations for molecular structure of tenofovir have been reported so far. In the present communication, we report detailed quantum chemical studies of the molecular structure and IR, Raman spectra of tenofovir using basis set. Information about the geometry and structure of the molecule, together with analysis of the IR, Raman spectra, based on frequency and intensity should help in understanding the structural and spectral characteristics. The band assignments have made by assuming C1 point group symmetry. Density Functional Theory (DFT) calculations have been performed to support our wave number assignments.

2. Experimental Details

The compound tenofovir in powder form was procured from Reputed pharmaceutical company, Chennai, India with more than 98% purity and was used as such without further purification to record FTIR and FT-Raman spectra. The FTIR spectrum of the compound was recorded in the region 4000–450 cm−1 in evacuation mode on Bruker IFS 66V spectrophotometer using KBr pellet technique (solid phase) with 4.0 cm−1 resolution. The FT-Raman spectrum was recorded using 1064 nm line of Nd: YAG laser as excitation wavelength in the region 5000–10 cm−1 on Bruker IFS 66V spectrometer equipped with FRA 106 Raman module was used as an accessory. The UV-vis spectral measurements were carried out using a varian cary 5E-UV-NIR spectrophotometer. The spectral measurements were carried out at sophisticated instrumentation Analysis Facility, IIT, Chennai, India.

3. Computational Details

To provide complete information regarding the structural characteristics and the fundamental vibrational modes of tenofovir of B3LYP/6-311++G(d,p), B3LYP/6-31G(d,p), and B3PW91/6-31G(d,p) correlation functional calculations have been carried out. The calculations of geometrical parameters in the ground state were performed using the Gaussian 03 [9] programs, invoking gradient geometry optimization [10] on Intel core i4/2.93 GHZ processor. The geometry optimization was carried out using the initial geometry generated from standard geometrical parameters at B3LYP and B3PW91 methods adopting 6-31G(d,p) and 6-311++G(d,p) basis sets to characterize all stationary points as minima. The optimized structural parameters of the compound tenofovir were used for harmonic vibrational frequency calculations resulting in IR and Raman frequencies together with intensities. In DFT methods, Becke’s three parameter exchange-functional (B3) [11, 12] combined with gradient-corrected correlation functional of Johnson and Frisch [13] and Burke et al. [1419] by implementing the split-valence polarized 6-31G(d,p) basis set [20, 21] have been utilized for the computation of molecular structure optimization, vibrational frequencies, thermodynamic properties, and energies of the optimized structures. The time-dependent density functional theory (TD-DFT) methods were used for the calculation of the UV-vis spectra. The IR and UV-vis spectra were calculated and visualized using the SWizard program [22]. The 13C and 1H NMR chemical shifts of the title compound were calculated using the keyword NMR in the DFT-B3LYP and B3PW91 levels.

4. Results and Discussion

4.1. Molecular Geometry

To compare the structural parameters of tenofovir and with the available experimental data, they are subjected to geometry optimization in their ground state. For each method, geometry optimizations were performed. The requested convergence on the maximum density matrix was 1026 a.u. the threshold value of the maximum displacement was 0.0018 A° and that of the maximum force was 0.00045 Hartree/Bohr using the Berny analytical gradient optimization routine [10, 23]. The nature of stationary points is checked by diagonalising the Hessian matrix to determine the number of imaginary wave numbers (zero for local minimum). The optimized bond lengths and bond angles of tenofovir calculated by B3LYP/6-311++G(d,p), B3LYP/6-31G(d,p), and B3PW91/6-31G(d,p) levels are listed in Table 1, in accordance with the atom numbering scheme shown in Figure 1. Table 1 shows the comparison of the calculated optimized structure parameters for tenofovir with those experimentally available data [8]. It is observed that most of the optimized structure parameters are slightly larger than the available data values. The bond lengths and bond angles computed by the DFT-B3LYP and B3PW91 levels show excellent agreement with the available computed values.

Structural parametersTenofovir
Experimental [8]B3LYP/6-311++G(d,p)B3LYP/6-31G(d,p)B3PW91/6-31G(d,p)

Internuclear distance (A°)

R(1-2) 1.4021.4121.4271.421
R(1-16) 1.8561.8371.8411.838
R(1-20) 1.1131.0961.0961.096
R(1-21) 1.1131.0931.0951.096
R(2-4) 1.4021.4431.4441.437
R(3-4) 1.5231.5321.5311.527
R(3-8) 1.4701.4561.4541.448
R(3-22) 1.1131.0911.0931.093
R(3-23) 1.1131.0921.0951.095
R(4-5) 1.5231.5241.5261.521
R(4-24) 1.1131.0951.0981.099
R(5-25) 1.1131.0911.0931.093
R(5-26) 1.1131.0941.0951.095
R(5-27) 1.1131.0941.0961.096
R(6-7) 1.3691.3101.3121.310
R(6-14) 1.3661.3821.3851.380
R(7-8) 1.3151.3811.3831.378
R(7-28) 1.1001.0801.0811.082
R(8-9) 1.3691.3791.3811.377
R(9-10) 1.3451.3371.3391.337
R(9-14) 1.3841.3961.3981.396
R(10-11) 1.3401.3351.3371.335
R(11-12) 1.3521.3411.3431.340
R(11-29) 1.1001.0861.0881.089
R(12-13) 1.3491.3431.3461.343
R(13-14) 1.3761.4081.4111.410
R(13-15) 1.2661.3541.3561.350
R(15-30) 1.0501.0061.0071.006
R(15-31) 1.0501.0071.0081.006
R(16-17) 1.6151.6171.6191.616
R(16-18) 1.6151.6261.6191.614
R(16-19) 1.4801.4841.4861.484
R(17-32) 0.9420.9650.9690.967
R(18-33) 0.9420.9650.9720.972

Bond angle (°)

A(2-1-16) 109.5111.9109.1109.2
A(2-1-20) 109.4112.5112.1112.2
A(2-1-21) 109.5108.3108.1108.0
A(1-2-4) 120.0117.6116.9116.5
A(16-1-20) 109.4108.0110.1110.2
A(16-1-21) 109.5107.7108.9108.7
A(1-16-17) 109.5102.7103.7103.8
A(1-16-18) 109.4105.2101.1100.5
A(1-16-19) 109.5117.7116.3116.6
A(20-1-21) 109.5108.3108.5108.4
A(2-4-3) 109.5106.5105.9105.9
A(2-4-5) 109.4111.1112.0112.2
A(2-4-24) 109.5109.0108.5108.5
A(4-3-8) 109.5113.0113.2113.2
A(4-3-22) 109.4109.6109.5109.4
A(4-3-23) 109.5109.4109.4109.5
A(3-4-5) 109.4111.5111.9111.8
A(3-4-24) 109.5108.6108.3108.2
A(8-3-22) 109.4107.7107.0107.0
A(8-3-23) 109.5108.4109.1109.2
A(3-8-7) 127.5128.0128.5128.6
A(3-8-9) 127.5126.1125.6125.5
A(22-3-23) 109.5108.6108.5108.4
A(5-4-24) 109.5110.0110.1110.1
A(4-5-25) 109.5110.3110.9111.0
A(4-5-26) 109.4110.3110.3110.3
A(4-5-27) 109.5110.4110.6110.7
A(25-5-26) 109.4108.5108.5108.4
A(25-5-27) 109.5108.8108.2108.2
A(26-5-27) 109.5108.4108.2108.2
A(7-6-14) 107.5104.0103.8103.7
A(6-7-8) 111.8113.9114.0114.0
A(6-7-28) 124.1125.3125.4125.3
A(6-14-9) 104.8111.0111.3111.3
A(6-14-13) 135.0133.0133.0133.0
A(8-7-28) 124.1120.8120.7120.7
A(7-8-9) 105.0105.8105.8105.9
A(8-9-10) 125.7128.1127.9127.8
A(8-9-14) 110.9105.2105.1105.0
A(10-9-14) 104.8126.6127.0127.2
A(9-10-11) 112.6111.6111.2111.0
A(9-14-13) 120.3116.0115.8115.6
A(10-11-12) 128.1128.4128.8128.9
A(10-11-29) 115.9116.1115.9115.9
A(12-11-29) 115.9115.4115.3115.2
A(11-12-13) 118.1118.7118.4118.4
A(12-13-14) 117.6118.7118.9118.9
A(12-13-15) 121.2118.9118.8118.8
A(14-13-15) 121.2122.4122.3122.3
A(13-15-30) 120.0119.9119.0119.5
A(13-15-31) 120.0118.7117.9118.3
A(30-15-31) 120.0119.9119.4120.0
A(17-16-18) 109.4102.8101.8102.2
A(17-16-19) 109.5116.9114.2113.9
A(16-17-32) 120.0114.2109.7109.3
A(18-16-19) 109.5113.1117.5117.6
A(16-18-33) 120.0111.9108.6107.9

The bond distance of carbon-carbon in the benzene ring in the range of 1.524 Å confirms the double-bond character and in the range of 1.370 Å supports the single-bond character. The bond distances of C–N and C–O were in the expected range and agree with the experimental values. Although there is a slight difference between the calculated and experimental structure parameters, the calculated geometric parameters represent a good approximation, and they are the bases for calculating other parameters such as atomic charges, NMR chemical shift values, vibrational modes, and thermodynamic properties.

4.2. Bond-Order Analysis

The bond order of tenofovir is presented in Table 2. Bond order is related to bond strength. The bonds with the higher bond order values have short bond length and vice versa. The bond-order analysis may predict that the weakest bonds may be cleaved preferentially and they may possess a relatively low pi-bond character. From Table 2, it is noted that bond between O2–C4 possess relatively low pi-bond character with low-bond-order value of 0.848 obtained DFT-B3LYP method, respectively. The P16–O19 bond order values are in the range 1.947, which depict the double bond character, while the C4–C5 and N8–C9 bond-order value is approximately unity, which shows the single-bond character. The optimized geometrical values are in support of the bond-order analysis.

Bond orderB3LYP/6-31G(d,p)

C1–O2 0.913
C1–P16 0.845
C1–H20 0.921
C1–H21 0.909
O2–C4 0.848
C3–C4 0.923
C3–N8 0.862
C3–H22 0.94
C3–H23 0.934
C4–C5 1.001
C4–H24 0.815
C5–H25 0.96
C5–H26 0.954
C5–H27 0.956
N6–C7 1.517
N6–C14 1.138
C7–N8 1.185
C7–H28 0.949
N8–C9 1.04
C9–N10 1.123
C9–C14 1.311
N10–C11 1.464
C11–N12 1.396
C11–H29 0.948
N12–C13 1.3
C13–C14 1.187
C13–N15 1.109
N15–H30 0.876
N15–H31 0.878
P16–O17 1.021
P16–O18 1.038
P16–O19 1.947
O17–H32 0.856
O18–H33 0.853

4.3. Electronic Properties

The energies of four important molecular orbitals of tenofovir: the highest and second highest occupied MOs (HOMO and HOMO−1), the lowest and the second lowest unoccupied MOs (LUMO and LUMO+1) were calculated and are presented in Table 3. The lowest singlet→singlet spin-allowed excited states of tenofovir were taken into account for the TD-DFT calculation in order to investigate the properties of electronic absorption. The experimental values are obtained from the UV/visible spectra recorded in methanol as reported earlier [24]. The calculations were also performed with methanol solvent effect. The calculated absorption wavelengths (), oscillator strength, excitation energies, and the experimental wavelengths are also given in Table 3. The energy gap between HOMO and LUMO is a critical parameter in determining molecular electrical transport properties [25]. In the electronic absorption spectrum of tenofovir, there are three absorption bands with a maximum 268.67, 250, and 230.69 nm. The strong absorption band 268.67 nm is caused by the and the other two moderately intense bands are due to transitions. The transitions are expected to occur relatively at lower wavelength, due to the consequence of the extended aromaticity of the benzene ring. The HOMO and LUMO of tenofovir are represented in Figure 2.

(Exp.; nm)
(Cal.; nm)





268.67286.264.33110.0004 π*
250252.034.91940.0019ππ −*−5.95−0.6631−5.98820.0117

4.4. Natural Population Analysis

The calculation of effective atomic charges plays an important role in the application of quantum mechanical calculations to molecular systems. Our interest here is in the comparison of different methods (B3LYP and B3PW91) to describe the electron distribution in tenofovir as broadly as possible, and to assess the sensitivity of the calculated charges to changes in the choice of the quantum chemical method. The calculated natural atomic charge values from the natural population analysis (NPA) and Mulliken population analysis (MPA) procedures using the B3LYP/6-311++G(d,p), B3LYP/6-31G(d,p), and B3PW91/6-31G(d,p) methods are listed in Table 4. The NPA from the natural-bonding orbital (NBO) method is better than the MPA scheme. Table 4 compares the atomic charge site of tenofovir from both MPA and NPA methods. The NPA of tenofovir shows that the presence of three oxygen atoms in the nitrate moiety [O17 = −0.979 (B3LYP/6-311++G(d,p)), −1.023 (B3LYP/6-31G(d,p)) and −1.023 (B3PW91/6-31G(d,p)); O18 = −0.983 (B3LYP/6-311++G(d,p)), −1.027 (B3LYP/6-31G(d,p)) and −1.027 (B3PW91/6-31G(d,p)); O19 = −1.070 (B3LYP/6-311++G(d,p)), −1.092 (B3LYP/6-31G(d,p)) and −1.087 (B3PW91/6-31G(d,p))] imposes large positive charges on the Phosphorus atom [P16 = 2.265 (B3LYP/6-311++G(d,p)), 2.388 (B3LYP/6-31G(d,p) and 2.373 (B3PW91/6-31G(d,p))]. However, the nitrogen atoms N6, N8, N10, and N12 possess large negative charges, resulting in the positive charges on the carbon atoms C7, C9, C11, C13, and C14. Moreover, there is no difference in charge distribution observed on all hydrogen atoms except the H30, H31, H32, and H33 hydrogen atoms. The positive charge on H30, H31, H32, and H33 hydrogen atoms is due to the negative charge accumulated on the N15, O17, and O18 atoms.

Atom with numberingMPANPA

O2 −0.527−0.508−0.506−0.617−0.599−0.596
C3 −0.059−0.077−0.122−0.184−0.287−0.305
C4 0.1630.1330.0900.0500.0440.035
C5 −0.347−0.333−0.404−0.577−0.703−0.729
N6 −0.539−0.587−0.602−0.516−0.510−0.510
C7 0.2850.2470.2420.2480.2050.198
N8 −0.513−0.479−0.516−0.412−0.393−0.391
C9 0.5120.5110.5420.3830.3690.366
N10 −0.524−0.564−0.596−0.583−0.587−0.586
C11 0.2270.2260.2170.3130.2640.257
N12 −0.521−0.548−0.568−0.558−0.557−0.556
C13 0.4800.5610.5940.4010.4180.412
C14 0.1950.2480.223−0.0020.0070.004
N15 −0.639−0.662−0.709−0.727−0.788−0.796
P16 1.0721.0541.0812.2652.3882.373
O17 −0.556−0.548−0.564−0.979−1.023−1.023
O18 −0.555−0.552−0.567−0.983−1.027−1.027
O19 −0.550−0.549−0.554−1.070−1.092−1.087
H20 0.1370.1290.1560.2040.2380.247
H21 0.1560.0890.1120.1700.1990.207
H22 0.1530.1280.1570.2080.2480.257
H23 0.1230.1280.1550.2080.2440.253
H24 0.1130.2230.2600.2590.2950.301
H25 0.1650.1240.1510.2190.2540.263
H26 0.1160.1040.1300.2040.2390.249
H27 0.1060.0950.1200.1940.2250.234
H28 0.1220.1270.1590.2030.2380.246
H29 0.1000.1090.1410.1790.2250.232
H30 0.2880.3040.3260.4190.4450.452
H31 0.2800.2890.3110.4110.4380.445
H32 0.3530.3490.3620.5140.5370.541
H33 0.3400.3530.3660.5200.5410.545

4.5. Vibrational Spectra

The experimental FTIR, calculated (B3LYP, B3PW91) and FT-Raman and vibrational spectra were shown in Figures 3 and 4. Since the calculated vibrational wave numbers were known to be higher than the experimental ones, they were scaled down by the wavenumber linear scaling procedure of Yoshida et al. [26] by using the expression: Comparison of the frequencies calculated at B3LYP and B3PW91 with experimental values (Table 5) reveals the overestimation of the calculated vibrational modes due to the neglect anharmonicity in real system. According to the theoretical calculations, tenofovir has a nonplanar structure of C1 point group symmetry. The molecule has 33 atoms and 93 normal modes of vibration active in both IR and Raman. Since the title molecule possess C1 point group symmetry, all the modes of vibration belong to A species only. The Chemcraft program was used to display the vibrational modes. Vibrational wave number assignments were made on the basis of combining the results of Chemcraft program with the symmetry, and taking the atomic displacements into consideration based on the frequency calculation and also made in analogy with the structurally related molecules.

ExperimentalTheoretical frequency
FTIRFT-RAMANB3LYP/6-311++G(d,p)B3LYP/6-31G(d,p)B3PW91/6-31G(d,p)Vibrational assignment

383536291738253620203859365021O17-H32 Stretching
38243619193761356315378335838O18-H33 stretching
3735354083750355373769357018N15-H30H31 asymmetric stretching
3227360534251536133432143640345615N15-H30H31 symmetric stretching
3106324130980326831220327731300C7-H28 stretching
30513052316630303317830414318630484C11-H29 stretching
312429920314230080315330180C5-H25H26 asymmetric stretching
311729861313730040314930150C3-H23H22 asymmetric stretching
29862991309229633311629852313630032C5-H26H27 asymmetric stretching
309129620311129800312629940C1-H20H21 asymmetric stretching
29392942306429384307229455308229544C3-H22H23 symmetric stretching
303429102305129261306129351C1-H20H21 symmetric stretching
302829050304529201305329281C5-H26H27 symmetric stretching
2836302729043303729131304529201C4-H24 stretching
16791656165516251001673164210016871655100N15-H30H31 scissoring
1574160615783161715881162415950N15-H30H31 scissoring
15041505152515005154315186155615308CH asymmetric stretching in methylene C5
149414712151014860150414800CH symmetric stretching in methylene C5
1472146914471148014570147814552C3-H22H23 + C1-H20H21 scissoring
1447145914372147814551147214492C3-H22H23 + C1-H20H21 scissoring
14211424141913993143514141142814072C5-H25H26H27 wagging
1330133013134133913212133413173C1-H20H21 wagging

The CH and OH stretching modes are expected to be observed at the high wave number region. The CH stretching bands are observed between 3106 cm−1 and 3051 cm−1 in FTIR spectrum. The CH stretching modes in the ring and methylene groups were identified and their effect in tenofovir was examined. Let us start considering the contribution of theoretical methods B3LYP/6-311++G(d,p), B3LYP/6-31G(d,p) and B3PW91/6-31G(d,p), the OH stretching mode is predicted at 3650 cm−1, 3583 cm−1 in B3PW91 and 3620 cm−1, 3563 cm−1 in B3LYP and 3629 cm−1, 3619 cm−1 in B3LYP/6-311++G(d,p) after scaling. This OH stretching mode is affected in the presence of water molecules. The CH2 scissoring and wagging mode was also identified. A major coincidence of theoretical values with that of experimental evaluations is found in the asymmetric and symmetric vibrations of the methylene (–CH2–) moiety. The –CH2– wagging mode at 1330 cm−1 in FT-Raman deviates positively by 20 cm−1 from the reported value of 1350 cm−1 [25] and theoretically coincide with 1321 cm−1 in B3LYP and 1334 cm−1 in B3PW91. The vibrational frequency 1455 cm−1 in B3LYP and 1449 cm−1 in B3PW91 is found neater to FTIR–CH2– scissoring mode of 1447 cm−1.

The scaled NH2 asymmetric and symmetric stretching vibrations of theoretical values are 3570 cm−1 and 3456 cm−1. The computed NH2 scissoring vibration at 1588 cm−1 in B3LYP and 1595 cm−1 in B3PW91 is in agreement with the expected characteristic value, 1600 cm−1 [27, 28]. In the presence of water molecules in the tenofovir, the NH stretching vibrations were red shifted by 10 cm−1. But however, the NH2 scissoring mode does not get affected. Notice that the most evident discrepancies between the experimental and calculated spectra are associated with the stretching modes showing the effect of water molecules. List of selected observed and calculated bands of tenofovir was presented in Table 5.

5. 13C and 1H-NMR Chemical Shift Assignment

The 13C and 1H NMR simulated theoretically with the aid of ChemDraw Ultra 10.0 is shown in Figures 5(a) and 5(b). Table 6 present the predicted chemical shift values of tenofovir obtained by the B3LYP/6-311++G(d,p), B3LYP/6-31G(d,p), B3PW91/6-31G(d,p) and ChemDraw Ultra 10.0 software package and its assignment along with the shielding values. In general, highly shielded electrons appear downfield and vice versa. The predicted chemical shift values by the theoretical methods both DFT values slightly deviates from the computed values of ChemDraw Ultra. The carbon atom C11 appearing at very higher chemical shift value (177.0 in B3LYP/6-311++G(d,p), 157.1 ppm in B3LYP/6-31G(d,p) and 152.6 ppm in B3PW91/6-31G(d,p)) due to negative charge of two nitrogen atoms (N10 and N12). Similarly C7, C9, and C13 appearing at higher chemical shift values (163.5, 174.7 and 166.4 ppm in B3LYP/6-311++G(d,p), 143.3, 155.1 and 144.9 ppm in B3LYP/6-31G(d,p), 141.1, 152.6 and 143.1 ppm in B3PW91/6-31G(d,p)) are due to nitrogen atoms N6, N8, and N15, respectively.

Atom positionB3LYP/6-311++B3LYP/6-31B3PW91/6-31ChemDraw UltraAssignment
Absolute shieldingChemical shiftAbsolute shieldingChemical shiftAbsolute shieldingChemical shift

1114.785.3126.673.3129.770.373.4C1 in Aliphatic
3124.775.3136.663.4139.360.758.6C3 in Aliphatic
4122.577.4134.765.3137.562.476C4 in Aliphatic
5153.946.1163.336.7166.433.620C5 in Aliphatic
736.4163.556.7143.358.8141.1143C7 in Purine
925.3174.744.9155.147.3152.6149.8C9 in Purine
1123.0177.042.9157.145.1154.9152.4C11 in Purine
1333.6166.455.0144.956.9143.1156.1C13 in Purine
1462.8137.281.3118.683.6116.4119.4C14 in Purine
2027.84.828.04.627.94.73.4H20 in Methylene
2127.94.728.04.628.04.63.4H21 in Methylene
2227. in Methylene
2326.85.826.75.926.66.03.7H23 in Methyene
2416.416.216.316.316.216.412.0H24 in Methine
2529.72.929.72.929.63.01.2H25 in Methyl
2630.32.330.22.430.22.41.2H26 in Methyl
2730.32.330.32.330.22.41.2H27 in Methyl
2823.88.823.88.823.69.08.1H28 in Purine
2923.39.323.19.522.99.78.2H29 in Purine
3025.07.625.27.425.07.67.0H30 in Aromatic
3125. in Aromatic
3229.23.429.03.628.93.712.0H32 in Alcohol
3328.34.328.24.428.14.512.0H33 in Alcohol

The carbon atoms C7, C9, C11, C13, and C14 are electropositive and possess more positive charges than the other carbon atoms, and hence the shielding is very small and appears upfield (see Table 6). In both molecules, the DFT-calculated atomic charges revealed that the more electron-rich atoms are C1, C3, C4, and C5; they are highly shielded atoms and hence appear at downfield (lower chemical shift). The carbon atoms in the benzene ring are deshielded than the carbon atoms in Aliphatic, so that the benzene carbon atoms in purine appear at higher chemical shift values than the aliphatic carbon atoms that were made by DFT methods. In this study, a good correlation between atomic charges and chemical shift was made. It is to be noted that 13C NMR chemical shifts for tenofovir which agree with the ChemDraw Ultra.

Table 6 gives the 1H NMR predicted chemical shift values obtained by the DFT methods and ChemDraw Ultra 10.0 software program [8] along with assignments. The predicted shielding values for each atom in the tenofovir molecule by B3LYP/6-311++G(d,p), B3LYP/6-31G(d,p), and B3PW91/6-31G(d,p) are given in Table 6. The predicted chemical shift values by the ChemDraw Ultra software program are in good agreement with the two DFT methods.

The spectrum of tenofovir showed a singlet at 7.0 ppm for the proton of the aromatic (H31) group, which is in good agreement with the chemDraw ultra value. Triplet is predicted at 2.9, 2.3, and 2.3 ppm for the methyl group of hydrogen atoms (H25, H26, and H27). This higher chemical shift for methyl group hydrogen is mainly due to the carbon atom (C5). The predicted value of singlet peak at higher chemical shift at 16.2 ppm for methine group of the indole proton (H24). Doublet is predicted at 7.0 ppm for the aromatic group of hydrogen atoms (H30 and H31). This higher chemical shift for methyl group hydrogen is mainly due to the nitrogen atom (N15). The hydrogen atoms of methylene group attached with C1 and C3 atoms show a multiplet at 3.4–3.9 ppm, which is due to the presence of oxygen and nitrogen atoms, respectively. In all the DFT methods predict the chemical shifts value of the hydrogen in hydroxyl group(H32, H33) is contradicting to chemDraw ultra value, but other hydrogen atoms fairly agrees.

6. Thermodynamic Properties

Several calculated thermodynamic parameters at room temperature are presented in Table 7. According to Koopmanns’ theorem, ionization potential (I) is the negative of the highest occupied molecular orbital (HOMO) energy, that is, I = −EHOMO, and affinity potential (A) is the negative of the lowest unoccupied molecular orbital (LUMO) energy, that is, A = −ELUMO, which are summarized in Table 7.


Total energy (a.u)−1267.50−1267.25−1266.851702
Zero-point energy (Kcal/Mol)158.540160.017160.67362
Rotational constants (GHz)0.16250.14630.1466
 Dipole moment (Debye)2.31521.59731.6778
 HOMO (eV)−6.27−5.94−6.04
 LUMO (eV)−0.95−0.51−0.57
 Energy gap (eV)5.315.435.46

Knowledge of permanent dipole moment of a molecule provides a wealth of information to determine the exact molecular conformation. The total dipole moment value of tenofovir in both DFTs and B3PW91 methods was observed. The variations in the entropy and zero-point vibrational energies seem to be insignificant.

7. Conclusion

The geometry of tenofovir was optimized with both DFT and B3PW91 methods using 6-31G(d,p) and 6-311++G(d,p) basis sets. The complete molecular structural parameters and thermodynamic properties of the optimized geometry of the compound have been obtained from DFT calculations. The vibrational frequencies of the fundamental modes of the compound have been precisely assigned and analyzed and the theoretical results were compared with the experimental vibrations. The bond order and atomic charges of the title molecule have been studied by DFT and B3PW91 methods. The energies of important MOs, absorption wavelength (), oscillator strength, and excitation energies of the compound were also determined from TD-DFT method and compared with the experimental values. This study predicted that the molecular geometry, vibrational wave numbers and 13C and 1H NMR chemical shifts for tenofovir could be successfully elucidated by the B3LYP/6-311++G(d,p), B3LYP/6-31G(d,p) and B3PW91/6-31G(d,p) methods using Gaussian program. Thus, the present investigation provides complete vibrational assignments, structural informations and electronic properties of the compound which may be useful to upgrade the knowledge on tenofovir.


  1. L. Naesens, R. Snoeck, G. Andrei, J. Balzarini, J. Neyts, and E. De Clercq, “HPMPC (cidofovir), PMEA (adefovir) and related acyclic nucleoside phosphonate analogues: a review of their pharmacology and clinical potential in the treatment of viral infections,” Antiviral Chemistry & Chemotherapy, vol. 8, no. 1, pp. 1–23, 1997. View at: Google Scholar
  2. M. Arimilli, C. Kim, and N. Bischofberger, “Synthesis, in vitro biological evaluation and oral bioavailability of 9-[2-(phosphonomethoxy)propyl]adenine (PMPA) prodrugs,” Antiviral Chemistry & Chemotherapy, vol. 8, no. 6, pp. 557–564, 1997. View at: Google Scholar
  3. L. Naesens, N. Bischofberger, P. Augustijns et al., “Antiretroviral efficacy and pharmacokinetics of oral bis(isopropyloxycarbonyloxymethyl)-9-(2-phosphonylmethoxypropyl)adenine in mice,” Antimicrobial Agents and Chemotherapy, vol. 42, no. 7, pp. 1568–1573, 1998. View at: Google Scholar
  4. J. Balzarini, A. Holy, J. Jindrich et al., “Differential antiherpesvirus and antiretrovirus effects of the (S) and (R) enantiomers of acyclic nucleoside phosphonates: potent and selective in vitro and in vivo antiretrovirus activities of (R)-9-(2-phosphonomethoxypropyl)- 2,6-diaminopurine,” Antimicrobial Agents and Chemotherapy, vol. 37, no. 2, pp. 332–338, 1993. View at: Google Scholar
  5. B. L. Robbins, R. V. Srinivas, C. Kim, N. Bischofberger, and A. Fridland, “Anti-human immunodeficiency virus activity and cellular metabolism of a potential prodrug of the acyclic nucleoside phosphonate 9-R-(2- phosphonomethoxypropyl)adenine (PMPA), bis(isopropyloxymethylcarbonyl)PMPA,” Antimicrobial Agents and Chemotherapy, vol. 42, no. 3, pp. 612–617, 1998. View at: Google Scholar
  6. B. L. Robbins, J. J. Greenhaw, M. C. Connelly, and A. Fridland, “Metabolic pathways for activation of the antiviral agent 9-(2- phosphonylmethoxyethyl)adenine in human lymphoid cells,” Antimicrobial Agents and Chemotherapy, vol. 39, no. 10, pp. 2304–2308, 1995. View at: Google Scholar
  7. V. Miller and B. A. Larder, “Mutational patterns in the HIV genome and cross-resistance following nucleoside and nucleotide analogue drug exposure,” Antiviral Therapy, vol. 6, supplement 3, pp. 25–44, 2001. View at: Google Scholar
  8. Carl Kemnitz 2002 Chemoffice ultra 10, Trial version.
  9. M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., Gaussian Inc., Woodbridge, Conn, USA, 2004.
  10. H. B. Schlegel, “Optimization of equilibrium geometries and transition structures,” Journal of Computational Chemistry, vol. 3, no. 2, pp. 214–218, 1982. View at: Publisher Site | Google Scholar
  11. A. D. Becke, “Density-functional exchange-energy approximation with correct asymptotic behavior,” Physical Review A, vol. 38, no. 6, pp. 3098–3100, 1988. View at: Publisher Site | Google Scholar
  12. C. Lee, W. Yang, and R. G. Parr, “Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density,” Physical Review B, vol. 37, no. 2, pp. 785–789, 1988. View at: Publisher Site | Google Scholar
  13. B. G. Johnson and M. J. Frisch, “Analytic second derivatives of the gradient-corrected density functional energy. Effect of quadrature weight derivatives,” Chemical Physics Letters, vol. 216, no. 1-2, pp. 133–140, 1993. View at: Publisher Site | Google Scholar
  14. K. Burke, J. P. Perdew, and Y. Wang, “Derivation of a generalized gradient approximation: the PW91 density functional,” in Electronic Density Functional Theory; Recent Progress and New Directions, J. F. Dobson, G. Vignale, and M. P. Das, Eds., Plenum Press, 1998. View at: Google Scholar
  15. J. P. Perdew, “Unified theory of exchange and correlation beyond the local density approximation,” in Electronic Structure of Solids '91, P. Ziesche and H. Eschrig, Eds., pp. 11–20, Academic Verlag, Berlin, Germany, 1991. View at: Google Scholar
  16. J. P. Perdew, J. A. Chevary, S. H. Vosko et al., “Atoms, molecules, solids, and surfaces: applications of the generalized gradient approximation for exchange and correlation,” Physical Review B, vol. 46, no. 11, pp. 6671–6687, 1992. View at: Publisher Site | Google Scholar
  17. J. P. Perdew, J. A. Chevary, S. H. Vosko et al., “Erratum: atoms, molecules, solids, and surfaces: applications of the generalized gradient approximation for exchange and correlation,” Physical Review B, vol. 48, no. 7, p. 4978, 1993. View at: Publisher Site | Google Scholar
  18. J. P. Perdew, K. Burke, and Y. Wang, “Generalized gradient approximation for the exchange-correlation hole of a many-electron system,” Physical Review B, vol. 54, no. 23, pp. 16533–16539, 1996. View at: Google Scholar
  19. J. P. Perdew and Y. Wang, “Accurate and simple analytic representation of the electron-gas correlation energy,” Physical Review B, vol. 45, no. 23, pp. 13244–13249, 1992. View at: Publisher Site | Google Scholar
  20. W. J. Hehre, L. Radom, P. V. R. Schleyer, and A. J. Pople, Ab Initio Molecular Orbital Theory, John Wiley & Sons, New York, NY, USA, 1989.
  21. Z. Dega-Szafran, A. Katrusiak, and M. Szafran, “Molecular structure of the complex of N-methylmorpholine betaine with 2,4-dinitrophenol,” Journal of Molecular Structure, vol. 741, no. 1–3, pp. 1–9, 2005. View at: Publisher Site | Google Scholar
  22. S. I. Gorelsky, SWizard program, 2005,
  23. H. B. Schlegel, “An efficient algorithm for calculating ab initio energy gradients using s, p Cartesian Gaussians,” The Journal of Chemical Physics, vol. 77, no. 7, pp. 3676–3681, 1982. View at: Google Scholar
  24. S. A. Ghorpade, M. S. Sali, A. H. Kategaonkar, D. M. Patel, V. P. Choudhari, and B. S. Kuchekar, “Simultaneous determination of emtricitabine and tenofovir by area under curve and dual wavelength spectrophotometric method,” Journal of the Chilean Chemical Society, vol. 55, no. 1, pp. 115–117, 2010. View at: Google Scholar
  25. R. M. Silverstein, G. C. Bassler, and T. C. Morrill, Spectrometric Identification of Organic Compounds, John Wiley & Sons, New York, NY, USA, 5th edition, 1981.
  26. H. Yoshida, K. Takeda, J. Okamura, A. Ehara, and H. Matsuura, “A new approach to vibrational analysis of large molecules by density functional theory: wavenumber-linear scaling method,” Journal of Physical Chemistry A, vol. 106, no. 14, pp. 3580–3586, 2002. View at: Publisher Site | Google Scholar
  27. R. A. Yadav and I. S. Sing, “Hydrogen bonding in o- and m-ethyl phenols,” Indian Journal of Pure and Applied Physics, vol. 23, pp. 626–627, 1985. View at: Google Scholar
  28. K. B. Wibery and A. Shrake, “A vibrational study of cyclohexane and some of its isotopic derivatives-III. A vibrational analysis of cyclohexane, cyclohexane-d12, cyclohexane-1,1,4,4-d4 and cyclohexane-1,1,2,2,4,4,5,5-d8,” Spectrochimica Acta A, vol. 29, no. 3, pp. 583–594, 1973. View at: Publisher Site | Google Scholar

Copyright © 2013 G. R. Ramkumaar et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.