Hydration Analysis of Antiviral Agent AZT by Means of DFT and MP2 Calculations

Effect of the sulfur atom on S2 and S4 positions of the uracil ring in different DNA:RNA hybrid microhelixes with three nucleotide base pairs

The effect of the sulphur atom on the uracil ring was analyzed in different DNA:RNA microhelixes with three nucleotide base-pairs, including uridine, 2-thiouridine, 4-thiouridine, 2,4-dithiouridine, cytidine, adenosine and guanosine. Distinct backbone and helical parameters were optimized at different density functional (DFT) levels. The Watson-Crick pair with 2-thiouridine appears weaker than with uridine, but its interaction with water molecules appears easier. Two types of microhelixes were found, depending on the H-bond of H2' hydroxyl atom: A-type appears with the ribose ring in ³ E-envelope C_3' -endo, and B-type in ² E-envelope C_2' -endo. B-type is less common but it is more stable and with higher dipole-moment. The sulphur atoms significantly increase the dipole-moment of the microhelix, as well as the rise and propeller twist parameters. Simulations with four Na atoms H-bonded to the phosphate groups, and further hydration with explicit water molecules were carried out. A re-definition of the numerical value calculation of several base-pair and base-stacking parameters is suggested.

FT-IR and FT-Raman spectra of 5-chlorocytosine: Solid state simulation and tautomerism. Effect of the chlorine substitution in the Watson-Crick base pair 5-chlorodeoxycytidine-deoxyguanosine

The laser Raman and IR spectra of 5-chlorocytosine have been recorded and accurately assigned in the solid state using Density functional calculations (DFT) together with the linear scaling equation procedure (LSE) and the solid state simulation of the crystal unit cell through a tetramer form. These results remarkably improve those reported previously by other authors. Several new scaling equations were proposed to be used in related molecules. The six main tautomers of the biomolecule 5-chlorocytosine were determined and optimized at the MP2 and CCSD levels, using different basis sets. The relative stabilities were compared with those obtained in cytosine and their 5-halo derivatives. Several relationships between energies, geometric parameters and NBO atomic charges were established. The effect of the chlorine substitution in the fifth position was evaluated through the stability of the Watson-Crick (WC) base pair of 5-chlorodeoxycytidine with deoxyguanosine, and through their vibrational spectra.

Structure and conformational analysis of the anti-HIV reverse transcriptase inhibitor AZT using MP2 and DFT methods. Differences with the natural nucleoside thymidine. Simulation of the 1st phosphorylation step with ATP

A comprehensive quantum-chemical investigation of the conformational landscape of the HIV-1 reverse transcriptase inhibitor AZT (3'-azido-3'-deoxythymidine) nucleoside analogue was carried out. The whole conformational parameters (χ, γ, β, δ, ϕ, P, νmax) were analysed as well as the NBO charges. The search located at least 55 stable structures, 9 of which were by MP2 within a 1 kcal mol(-1) electronic energy range of the global minimum. Most conformers were anti or high-anti around the glycoside bond and with North sugar ring puckering angles. The distribution of all the conformers according to the ranges of stability of the characteristic torsional angles was established. The results obtained were in accordance with those found in related anti-HIV nucleoside analogues. The best conformer in the anti form corresponded to the calculated values by MP2 of χ = -126.9°, β = 176.4° and γ = 49.1°. An analysis of the lowest vibrations in conformer C1 was carried out. The first hydration shell was simulated and the structural differences with the natural nucleoside deoxythymidine (dT) were determined. The first phosphorylation step was simulated by interacting ATP with the best hydrated clusters of AZT and dT. The Na cations act as a bridge between the phosphate moieties of ATP making it easy for -P3O3 to receive the H5' proton from AZT or dT. A proton-transfer mechanism is proposed through the water molecules. When the number of the water molecules surrounding AZT is lower than 8, the first phosphorylation step of AZT can be carried out. However, the appropriate orientation of the O5'-H in dT avoids this limitation and it can be performed with large numbers of water molecules.

Conformational analysis, molecular structure and solid state simulation of the antiviral drug acyclovir (zovirax) using density functional theory methods

The five tautomers of the drug acyclovir (ACV) were determined and optimised at the MP2 and B3LYP quantum chemical levels of theory. The stability of the tautomers was correlated with different parameters. On the most stable tautomer N1 was carried out a comprehensive conformational analysis, and the whole conformational parameters (R, β, Φ, φ1, φ2, φ3, φ4, φ5) were studied as well as the NBO Natural atomic charges. The calculations were carried out with full relaxation of all geometrical parameters. The search located at least 78 stable structures within 8.5 kcal/mol electronic energy range of the global minimum, and classified in two groups according to the positive or negative value of the torsional angle j1. In the nitrogen atoms and in the O2' and O5' oxygen atoms of the most stable conformer appear a higher reactivity than in the natural nucleoside deoxyguanosine. The solid state was simulated through a dimer and tetramer forms and the structural parameters were compared with the X-ray crystal data available. Several general conclusions were emphasized.

The biomolecule of 5-bromocytosine: FT-IR and FT-Raman spectra and DFT calculations. Identification of the tautomers in the isolated state and simulation the spectra in the solid state

An accurate assignment of the IR spectrum in Ar matrix of 5-bromocytosine and of the IR and Raman spectra in the solid state was carried out. For this purpose Density functional calculations (DFTs) were performed to clarify wavenumber assignments of the experimental observed bands. The calculated values were scaled using scaling equations and they were compared with IR and Raman experimental data. Good reproduction of the experimental wavenumbers is obtained and the% error is very small in the majority of cases. In the isolated state all the tautomer forms of 5-bromocytosine were determined and optimized. The wavenumbers corresponding to C1 and C2b tautomers were identified and assigned in the IR experimental spectrum reported in Ar matrix. Our study confirms the existence of at least two tautomers, the amino-oxo and the amino-hydroxy in the isolated state. In the solid state the FT-IR and FT-Raman spectra of 5-bromocytosine in the powder form were recorded in the region 400-4000 cm(-1) and 50-3500 cm(-1), respectively. The unit cell found in the crystal was simulated as a tetramer form in three tautomers. Thus, it has been possible to assign all the 33 normal modes of vibration. The study indicates that the features, that are the characteristic of the vibrational spectra of cytosine, are retained by the spectra of 5-bromocytosine and it exists in the solid phase in the amino-oxo form.

Solid state simulation of tetramer form of 5-aminoorotic acid: the vibrational spectra and molecular structure study by using MP2 and DFT calculations

The Raman and IR spectra of the biomolecule 5-aminoorotic acid in the solid state were simulated by a dimer and tetramer forms, with the special interest in the interactions that involve the NH and NH(2) groups. The unit cell expected in the crystal was simulated as a tetramer form by density functional calculations. They were performed to clarify wavenumber assignments of the experimentally observed bands in the spectra. Correlations with the molecule of uracil were made, and accurate scaling procedures deduced by us were employed in the calculated wavenumbers of 5-aminoorotic acid. Good reproduction of the experimental wavenumbers is obtained and the % error is very small in the majority of cases. This fact confirms our simplified solid state model. The scaling leads to a reassignment of the IR and Raman experimental bands. The NBO atomic charges and several thermodynamic parameters were also calculated.

Microsolvated transition state models for improved insight into chemical properties and reaction mechanisms

Over the years, several methods have been developed to effectively represent the chemical behavior of solutes in solvents. The environmental effects arising due to solvation can generally be achieved either through inclusion of discrete solvent molecules or by inscribing into a cavity in a homogeneous and continuum dielectric medium. In both these approaches of computational origin, the perturbations on the solute induced by the surrounding solvent are at the focus of the problem. While the rigor and method of inclusion of solvent effects vary, such solvation models have found widespread applications, as evident from modern chemical literature. A hybrid method, commonly referred to as cluster-continuum model (CCM), brings together the key advantages of discrete and continuum models. In this perspective, we intend to highlight the latent potential of CCM toward obtaining accurate estimates on a number of properties as well as reactions of contemporary significance. The objective has generally been achieved by choosing illustrative examples from the literature, besides expending efforts to bring out the complementary advantages of CCM as compared to continuum or discrete solvation models. The majority of examples emanate from the prevalent applications of CCM to organic reactions, although a handful of interesting organometallic reactions have also been discussed. In addition, increasingly accurate computations of properties like pK(a) and solvation of ions obtained using the CCM protocol are also presented.