DFT calculation of four new potential agents muscarinic of bispyridinium type: structure, synthesis, biological activity, hydration, and relations with the potents W84 and DUO-3O

Biomolecules of 2-Thiouracil, 4-Thiouracil and 2,4-Dithiouracil: A DFT Study of the Hydration, Molecular Docking and Effect in DNA:RNAMicrohelixes

The molecular structure of 2-thiouracil, 4-thiouracil and 2,4-dithiouracil was analyzed under the effect of the first and second hydration shell by using the B3LYP density functional (DFT) method, and the results were compared to those obtained for the uracil molecule. A slight difference in the water distribution appears in these molecules. On the hydration of these molecules several trends in bond lengths and atomic charges were established. The ring in uracil molecule appears easier to be deformed and adapted to different environments as compared to that when it is thio-substituted. Molecular docking calculations of 2-thiouracil against three different pathogens: Bacillus subtilis, Escherichia coli and Candida albicans were carried out. Docking calculations of 2,4-dithiouracil ligand with various targeted proteins were also performed. Different DNA: RNA hybrid microhelixes with uridine, 2-thiouridine, 4-thiouridine and 2,4-dithiouridine nucleosides were optimized in a simple model with three nucleotide base pairs. Two main types of microhelixes were analyzed in detail depending on the intramolecular H-bond of the 2'-OH group. The weaker Watson-Crick (WC) base pair formed with thio-substituted uracil than with unsubstituted ones slightly deforms the helical and backbone parameters, especially with 2,4-dithiouridine. However, the thio-substitution significantly increases the dipole moment of the A-type microhelixes, as well as the rise and propeller twist parameters.

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.

DFT computations on vibrational spectra: Scaling procedures to improve the wavenumbers

The performance of ab initio and density functional theory (DFT) methods in calculating the vibrational wavenumbers in the isolated state was analyzed. To correct the calculated values, several scaling procedures were described in detail. The two linear scaling equation (TLSE) procedure leads to the lowest error and it is recommended for scaling. A comprehensive compendium of the main scale factors and scaling equations available to date for a good accurate prediction of the wavenumbers was also shown. Examples of each case were presented, with special attention to the benzene and uracil molecules and to some of their derivatives. Several DFT methods and basis sets were used. After scaling, the X3LYP/DFT method leads to the lowest error in these molecules. The B3LYP method appears closely in accuracy, and it is also recommended to be used. The accuracy of the results in the solid state was shown and several additional corrections are presented.

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.

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.