Comprehensive conformational analysis of the nucleoside analogue 2'-beta-deoxy-6-azacytidine by DFT and MP2 calculations

A comprehensive conformational analysis of isolated 2'-beta-deoxy-6-azacytidine (d6AC), an analogue of therapeutically active 6-azacytidine (6AC), has been performed by means of ab initio calculations at the MP2/6-311++G(2df,pd)//DFT B3LYP/6-31G(d,p) level of theory. Among the 81 conformers located within a 7.83 kcal/mol Gibbs energy range at T = 298.15 K, 38 contain syn-oriented bases with respect to 2'-deoxyribose; the other conformers include anti-oriented bases. Energetic analysis of these conformers shows that conformational equilibrium of isolated d6AC at T = 298.15 K is shifted to syn conformation with a syn/anti ratio estimated as 61.4%:38.6%. As far as the sugar conformation is concerned, 40 conformers contain north (N) (with 0.3 degrees < or = P < or = 40.1 degrees), and the rest possess south (S) (with 157.1 degrees < or = P < or = 207.0 degrees) puckers, where P is the pseudorotational angle of the furanose ring. The S/N occupancy ratio is estimated as 80.2%:19.8% (T = 298.15 K). The two most stable conformers are energetically quasidegenerate and correspond to both C2'-endo/syn conformers differing only by orientation of the O3'H hydroxyl group. They are both stabilized by means of similar intramolecular H-bonds, i.e., O5'H...O2, C2'H2...O2, and C2'H2...O5'. As examined by AIM criteria, from 1 to 3 H-bonds per conformer were identified among 13 possible interactions: O5'H...O2, O5'H...N6, O3'H...O5', O5'H...O3', C1'H...O2, C2'H2...O2, C2'H2...O5', C3'H...O2, C3'H...N6, C5'H1...O2, C5'H2...O2, C5'H1...N6, and C5'H2...N6. The biological effect of d6AC is conceived as an inhibition of replicative DNA polymerase caused by an unusual orientation of the sugar residue against the base in the only A form DNA-like conformer.     

HF, MP2 and DFT calculations and spectroscopic study of the vibrational and conformational properties of N-diethylendiamine

The conformational stability and vibrational modes of the N-diethylendiamine organic cation (N-DD(2+)) were studied by experimental (Raman) spectroscopy combined with theoretical calculations. Various ab initio theories were used: Hartree-Fock (HF) theory, M?ller-Plesset second-order perturbation (MP2) theory and density functional theory (DFT). Three stable conformers of N-DD(2+), trans-trans, gauche-gauche and gauche-trans were calculated. A comparison between the computed structural parameters of the conformers at both levels of theory and the X-ray data was made. It is demonstrated that the N-DD(2+) cation adopts more probably the gauche-gauche conformation at room temperature. In order to make a more detailed interpretation of the low temperature phase transition of N-DDHP, the Raman spectra of N-DDHP were recorded at room and low temperature in the 200-3400 cm(-1) region. The vibrational frequencies of the different conformers of N-DD(2+) were also calculated using the DFT/B3LYP (6-31G(d)) level of theory. By comparison between the experimental and theoretical results, the conformational dynamic of the N-DD(2+) organic cation was confirmed. It is shown that the N-DD(2+) cation configuration changes from gauche-gauche conformer to gauche-trans conformer when decreasing the temperature.     

DFT,AIM, and NBO analyses of 1-methyl-2-thioxoimidazolidin-4-one tautomers and their complexes with iodine

Thioimidazoline derivatives can be used to treat hyperthyroidism due to their ability to make complexes with iodine. In this research designed to find new structures with the same ability, 1-methyl-2-thioxoimidazolidin-4-one (MTIO) and the structures of MTIO tautomers (5 tautomers), their isomers (total 9 isomers) and their complexes with iodine are optimized using the B3LYP method with two different basis sets to obtain their molecular parameters, relative energies, and vibrational frequencies. The relative energies show that in all tautomers and complexes, ketone and thione forms are more stable than enol and thienol forms, and also Z isomers are more stable than E isomers. Moreover, the NBO calculation is carried out for tautomers and complexes to obtain atomic charges and acceptor-donor interactions. These results confirm the ability of MTIO tautomers to form complexes and show that the planar complexes have more effective interaction than the perpendicular complexes. The essence and important complexation properties are also calculated and confirmed using the AIM analysis.     

First-principle DFT and MP2 modeling of infrared reflection-absorption spectra of oriented helical ethylene glycol oligomers

First-principle modeling is used to obtain a comprehensive understanding of infrared reflection absorption (RA) spectra of helical oligo(ethylene glycol) (OEG) containing self-assembled monolayers (SAMs). Highly ordered SAMs of methyl-terminated 1-thiaoligo(ethylene glycols) [HS(CH2CH2O)(n)CH3, n = 5, 6] on gold recently became accessible for systematic infrared analyses [Vanderah et al., Langmuir, 2003, 19, 3752]. We utilized the quoted experimental data to validate the first-principle modeling of infrared RA spectra of HS(CH2CH2O)(5,6)CH3 obtained by (i) DFT methods with gradient corrections (using different basis sets, including 6-311++G) and (ii) HF method followed by a M?ller-Plesset (MP2) correlation energy correction. In focus are fundamental modes in the fingerprint and CH-stretching regions. The frequencies and relative intensities in the calculated spectra for a single molecule are unambiguously identified with the bands observed in the experimental RA spectra of the corresponding SAMs. In addition to confirming our earlier assignment of the dominating peak in the CH-stretching region to CH2 asymmetric stretching vibrations, all other spectral features observed in that region have received an interpretation consistent (but not in all cases coinciding) with previous investigations. The obtained results provide an improved understanding of the orientation and conformation of the molecular building blocks within OEG-containing assemblies, which, in our opinion, is crucial for being able to predict the folding and phase characteristics and interaction of OEG-SAMs with water and proteins.     

MP2, DFT and ab initio calculations on thioxanthone

Density functional theory (DFT), HF and MP2 calculations have been carried out to investigate thioxanthone molecule using the standard 6-31+G(d,p) basis set. The results of MP2 calculations show a butterfly structure for thioxanthone. The calculated results show that the predicted geometry can well reproduce the structural parameters. The predicted vibrational frequencies were assigned and compared with experimental IR spectra. A good harmony between theory and experiment is found. The theoretical electronic absorption spectra have been calculated using CIS method. ¹³C and ¹H NMR of the title compound have been calculated by means of B3LYP density functional method with 6-31+G(d,p) basis set. The comparison of the experimental and the theoretical results indicate that density functional B3LYP method is able to provide satisfactory results for predicting NMR properties.     

Vibrational analysis,conformational stability,force constants,internal rotation barriers,MP2=full and DFT calculations of 1,3-dimethyluracil tautomers

The molecular structure of 1,3-dimethyluracil (C₆H₈N₂O₂; 1,3-DMU) is studied theoretically and experimentally using Gaussian 98 calculations and different spectroscopic techniques. The vibrational spectrum for 1,3-DMU in the solid phase is recorded in the IR range 4000-400 cm^–1. Initially, in order to get the most stable structure, twelve structures were proposed for the titled compound as a result of the internal rotation of CH₃ around C–N bonds and keto-enol tautomerism. The single point energy and frequency calculations are obtained by MP2 (Full) and DFT/B3LYP methods with the 6-31G(d) basis set using the Gaussian 98 computation package. After the complete relaxation of twelve isolated isomers, the (diketo) tautomer was the only favored structure owing to its low energy relative to the other isomers and the prediction of real frequencies. This interpretation is supported by the recorded infrared spectrum that shows the presence of only the diketo tautomer. Aided by the normal coordinate analysis and potential energy distributions, a confident vibrational assignment of the fundamental frequencies is calculated. The results are discussed herein and compared with similar molecules whenever possible.     

Structure and stability of thiourea with water, DFT and MP2 calculations

The relative stabilities of thiourea in water are investigated computationally by considering thiourea–water complexes containing up to 1–6 water molecules (CS(NH₂)₂(H₂O)_n=1–6) using density functional theory and MP2 ab initio molecular orbital theory. The results show that the thiourea complex is stable and has an unusually high affinity for incoming water molecules. The clusters are progressively stabilized by the addition of water molecules, as indicated by the increasing of the binding energy. The binding energy of the cluster to each H₂O molecule is about 33 kJ mol⁻¹ for n=1–5.The C–S bond, N–C bond distance, Mulliken populations and binding energy keep approximately constant as the clusters increase in size with an increasing number of H₂O molecules. As the solvation progresses, the C–S distance increases monotonically while the Mulliken populations on the C–S bond reduces monotonically with the addition of each H₂O molecule, indicating that the C–S bond of the thiourea unit in the clusters is de-stabilized with an increasing number of H₂O molecules. Charge transfers for the clusters are mainly found at N, S atoms of the thiourea.     

Revealing the physical nature and the strength of charge‐inverted hydrogen bonds by SAPT(DFT), MP2, SCS‐MP2, MP2C,and CCSD(T) methods

The physical nature of charge‐inverted hydrogen bonds in H₃XH YH₃ (X = Si, Ge; Y = Al, Ga) dimer systems is studied by means of the SAPT(DFT)‐based decomposition of interaction energies and supermolecular interaction energies based on MP2, SCS‐MP2, MP2C, and CCSD(T) methods utilizing dimer‐centered aug‐cc‐pCVnZ (n = D, T, Q) basis sets as well as an extrapolation to the complete basis set limit. It is revealed that charge‐inverted hydrogen bonds are inductive in nature, although dispersion is also important. Computed interaction energies form the following relation: . It is confirmed that the aug‐cc‐pCVDZ basis set performs poorly and that very accurate values of interaction and dispersion energies require basis sets of at least quadrupole‐ζ quality. Considerably large binding energies suggest potential usefulness of charge‐inverted hydrogen bonds as an important structural motif in molecular binding. Terminology applying to σ‐ and π‐hole interactions as well as to triel and tetrel bonds is discussed. According to this new terminology the charge‐inverted hydrogen bond would become the first described case of a hydride‐triel bond. © 2017 Wiley Periodicals, Inc.     

X-Ray, MP2 and DFT studies of the structure and vibrational spectra of homarinium chloride

The effects of hydrogen bonding, inter- and intramolecular electrostatic interactions on the structure of homarinium chloride, HOMH·Cl, in the crystal and its isolated molecule have been studied by X-ray diffraction, FT-IR, Raman, ¹H and ¹³C NMR spectroscopies, and by the MP2 and DFT theoretical methods. In the crystal, the Cl⁻ anion is connected with protonated homarine via the O–HCl⁻ hydrogen bond of the length of 2.937(4) Å, and two N⁺Cl⁻ intermolecular electrostatic interactions. In the isolated molecule, according to the MP2 and B3LYP calculations, the Cl⁻ anion is engaged in a shorter hydrogen bond (O–HCl⁻ of 2.811–2.861 Å) and in one type of intramolecular electrostatic interactions. The calculated bond lengths and bond angles at the MP2 and B3LYP levels of theory are in good agreement with the X-ray data, except the conformation of the COOH group, which is cis (syn) in the crystal and trans (anti) in the isolated molecule. The tentative assignments for the experimental solid state vibrational spectra of HOMH·Cl and HOMD·Cl have been made on the basis of the B3LYP/6-31G(d,p) calculated frequencies and intensities. The effect of quaternization of picolinic acid on the chemical shifts of the ring protons and carbons is analyzed.     

MP2 and DFT theoretical studies of the geometry, vibrational and electronic absorption spectra of 2-aminopyrimidine

The conformational stability, relative IR intensities and harmonic vibrational wave numbers of 2-aminopyrimidine in the electronically ground state have been investigated by a comparison of the experimental and calculated spectroscopy. Two ab initio theories, Möller–Plesset second-order perturbation theory and density functional theory DFT-B3LYP methods, were used at different extended basis sets. The vibrational wave numbers calculated with DFT method were scaled using two different methods. The solvent effects on the structural parameters were also extensively studied. A comparison between the computed structural parameters of the investigated compound and the experimental data was shown. A good agreement between the theoretical results and the experimental data was found.     

Parallel computation of the MP2 energy on distributed memory computers

A parallel distributed implementation of the second-order Møller-Plesset perturbation theory method, widely used in quantum chemistry, is presented. Parallelization strategy and performance for the HONDO quantum chemistry program running on a network of Unix computers are also discussed. Superlinear speedups are obtained through a combined use of the CPU and memory of the different processors. Performance for standard and direct algorithms are presented and discussed. A superdirect algorithm that eliminates the communication bottleneck during the integral transformation step is also proposed. © 1995 by John Wiley & Sons, Inc.     

Berberine alkaloid: Quantum chemical study of different forms by the DFT and MP2 methods

The stable structures and electronic properties for the berberine cation as well as possible ammonium, carbinol and amino-aldehyde forms of protoberberine salts in the presence of hydroxyl ions were investigated by the B3LYP/6-31G(d,p) and MP2/6-31++G(d,p) methods. The geometry optimizations by both methods lead to the nonplanar propeller-twisted and buckled structure for the all forms. The obtained bond lengths and bond angles agree with the experimental values. The comparison of total energies elucidates that the amino-aldehyde form is the most preferable tautomer in gas phase, while the carbinol form is less stable. The least stable tautomer is the ammonium form.     

DFT and MP2 study of simple and water-assisted tautomerism in amidrazones

DFT and MP2 methods were used to calculate structural parameters, vibrational modes, solvent effect, and energetic properties of amidrazones. Amidrazones can be presented by three tautomeric forms and six isomers. All tautomers and transition states were optimized at the B3LYP/6-311++g** and MP2/6-311++G** levels of theory. The relative stabilities of amidrazone isomers in the gas phase were found to be as 1Z > 1E > 2E > 2Z > 3E > 3Z > TS(1–2) > TS(1–3). The calculated energy differences between E and Z isomers are very low and between different tautomers are nearly low, but the energy barriers for tautomerism interconversions at the gas phase are high. The kinetic and thermodynamic data in solvents (chloroform, tetrahydrofuran, acetone, and water) are nearly similar to those in the gas phase but their rate constants are slightly less than those in the gas phase. Moreover, equilibrium and rate constants of intermolecular tautomerism in presence of 1–3 molecules of water were calculated. Computed energy barriers show that the barrier energy of water-assisted tautomerism is very lower than that in simple tautomerism and also calculated binding energies show that water can stabilize transition states more than tautomers. Therefore, this water-assisted tautomerism can be performed fast, especially with the assistance of two molecules of water.     

Molecular structure and vibrational raman and infrared spectra of bromochlorofluoromethane and its silicon and germanium analogs: quantum-mechanical DFT and MP2 calculations

DFT(B3LYP) and MP2 calculations with the 6-311G(2d, 2p)-type basis set have been carried out for the prediction of molecular parameters (bond distances, bond angles, rotational constants, and dipole moments) and vibrational Raman and infrared spectra (harmonic wavenumbers, absolute intensities, Raman scattering activities, and depolarization ratios) of bromochlorofluoromethane (HCBrCIF) and its silicon and germanium analogs (HSiBrClF and HGeBrCIF). The predicted geometry and vibrational Raman and infrared spectra of HCBrClF agree well with the available experimental data for this molecule and their deuterated derivatives. This agreement allows one to believe that the predicted molecular parameters and vibrational spectra of HSiBrClF, HGeBrClF, and their deuterated derivatives will guide their future experimental studies.     

AlCl_3催化歧化制备二甲基二氯硅烷的DFT和MP2研究(II)

研究一甲基三氯硅烷与三甲基氯硅烷分别与三氯化铝作用并得到中间体后,中间体分别与一甲基三氯硅烷和三甲基氯硅烷作用的机理,反应的IRC(Intrinsicreactioncoordinate)曲线,各关键原子沿IRC的变化曲线及产物的理论产率.结果发现,中间体趋于与一甲基三氯硅烷作用生成二甲基二氯硅烷,且三甲基氯硅烷的转化率高于一甲基三氯硅烷.实验结果与计算结果吻合.     

胸腺嘧啶-BH3构型与性质的DFT和MP2研究

Geometries and combination energies are predicated at B3LYP / 6-31G（d）and MP2 / 6-31G（d）level for thymine-BH3 complexes and 5 geometries have been obtained. Then single point energy calculations using larger basis sets（6-311 + G（2df）and aug-cc-pVDZ）and vibrational analysis and natural bond orbital analysis are carried out on the 5 optimized conformers. The outcome indicates that the conformers with the boron atom combined with O directly are relatively stable ones,（a）and（b）,with the combination energies of 90. 4 and 88. 0 kJ / mol （B3LYP / 6-31G（d）,BSSE corrected）. The fact is that the nitrogen atom offers electron to the empty atomic orbital of boron which produces the conformers（c）and（d）. Only one conformer is found which is formed because two carbon atoms offer π electron to the empty orbital of boron. The charge transference exists in all the conformers. The combination energies have a good line relation with their charge transference. The calculated results show that when the complex forms their IR spectrum moved to the red side and the frequency shifts are relative to the stabilities of the complexes.     

A comparison of Hartree—Fock,MP2, and DFT results for the HCN dimer and crystal

A number of hydrogen-bond related quantities—geometries, interaction energies, dipole moments, dipole moment derivatives, and harmonic vibrational frequencies—were calculated at the Hartree—Fock, MP2, and different DFT levels for the HCN dimer and the periodic HCN crystal. The crystal calculations were performed with the Hartree—Fock program CRYSTAL92, which routinely allows an a posteriori electron-correlation correction of the Hartree—Fock obtained lattice energy using different correlation-only functionals. Here, we have gone beyond this procedure by also calculating the electron-correlation energy correction during the structure optimization, i.e., after each CRYSTAL92 Hartree—Fock energy evaluation, the a posteriori density functional scheme was applied. In a similar manner, we optimized the crystal structure at the MP2 level, i.e., for each Hartree—Fock CRYSTAL92 energy evaluation, an MP2 correction was performed by summing the MP2 pair contributions from all HCN molecules within a specified cutoff distance. The crystal cell parameters are best reproduced at the Hartree—Fock and the nongradient-corrected HF + LDA and HF + VWN levels. The BSSE-corrected MP2 method and the HF + P91, HF + LDA, and HF + VWN methods give lattice energies in close agreement with the ZPE-corrected experimental lattice energy. The (HCN)₂ dimer properties are best reproduced at the MP2 level, at the gradient-corrected DFT levels, and with the B3LYP and BHHLYP methods. © 1996 John Wiley & Sons, Inc.     

Electronic and infrared spectroscopy of jet-cooled (+/-)-cis-1-amino-indan-2-ol hydrates

The role of conformational isomerism in molecular interaction has been studied using the example of jet-cooled complexes of (+/-)-cis-1-amino-indan-2-ol with water. The two formerly evidenced conformers of (+/-)-cis-1-amino-indan-2-ol easily form hydrates and dihydrates, which have been studied by means of laser-induced fluorescence and IR/UV double resonance spectroscopy, as well as ab initio calculations. All the 1 : 1 and 1 : 2 complexes with water evidenced in this work involve "ring" structures, in which the water monomer or dimer acts as an acceptor from the NH(2) and a donor to the OH groups of (+/-)-cis-1-amino-indan-2-ol. However, the water lies externally to the indan frame in the hydrates of conformer I of (+/-)-cis-1-amino-indan-2-ol, which possesses axial NH(2) and equatorial OH groups, and above it for the hydrates with the less stable conformer II, with equatorial NH(2) and axial OH groups. Consequently, the different steric constraints which exist in the two conformers result in different hydrogen bond topologies, with an additional OH[dot dot dot]pi interaction for the hydrates of conformer II.     

Peptide models XXIII. Conformational model for polar side‐chain containing amino acid residues: A comprehensive analysis of RHF,DFT, and MP2 properties of HCO‐L‐SER‐NH2

Geometric and energetic properties of a diamide of serine, HCO‐NH‐L ‐CH(CH₂OH)CO‐NH₂, are investigated by standard methods of computational quantum chemistry. Similarly to other amino acid residues, conformational properties of HCO‐L ‐Ser‐NH₂ can be derived from the analysis of its E=E(ϕ,ψ;χ₁,χ₂) hypersurface. Reoptimization of 44 RHF/3‐21G conformers at the RHF/6‐311++G** level resulted in 36 minima. For all conformers, geometrical properties, including variation of H‐bond parameters and structural shifts in the torsional space, are thoroughly investigated. Results from further single‐point energy calculations at the RHF, DFT, and MP2 levels, performed on the entire conformational data set, form a database of 224 energy values, perhaps the largest set calculated so far for any single amino acid diamide. A comprehensive analysis of this database reveals significant correlation among energies obtained at six levels of ab initio theory. Regression parameters provide an opportunity for extrapolation in order to predict the energy of a conformer at a high level by doing explicit ab initio computations only for a few selected conformers. The computed conformational and relative energy data are compared with structural and occurrence results derived from a nonhomologous protein database incorporating 1135 proteins. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 626–655, 2000