Theoretical ab initio study of the effects of methylation on structure and stability of G:C Watson-Crick base pair

Methylation of DNA occurs most readily at N(3), N(7), and O(6) of purine bases and N(3) and O(2) of pyrimidines. Methylated bases are continuously formed through endogenous and exogenous mechanisms. The results of a theoretical ab initio study on the methylation of G:C base pair components are reported. The geometries of the local minima were optimized without symmetry restrictions by the gradient procedure at DFT level of theory and were verified by energy second derivative calculations. The standard 6-31G(d) basis set was used. The single-point calculations have been performed at the MP2/6-31G(d,p), MP2/6-31++G(d,p), and MP2/6-311++G(2d,2p) levels of theory. The geometrical parameters, relative stability and counterpoise corrected interaction energies are reported. Also, using a variation-perturbation energy decomposition scheme we have found the vital contributions to the total interaction energy.     

Theoretical ab Initio Study of the Effects of Methylation on Structure and Stability of G:C Watson-Crick Base Pair

Abstract

Methylation of DNA occurs most readily at N(3), N(7), and O(6) of purine bases and N(3) and O(2) of pyrimidines. Methylated bases are continuously formed through endogenous and exogenous mechanisms. The results of a theoretical ab initio study on the methylation of G:C base pair components are reported. The geometries of the local minima were optimized without symmetry restrictions by the gradient procedure at DFT level of theory and were verified by energy second derivative calculations. The standard 6–31G(d) basis set was used. The single-point calculations have been performed at the MP2/6–31G(d,p), MP2/6–31₊₊G(d,p), and MP2/6–311₊₊G(2d,2p) levels of theory. The geometrical parameters, relative stability and counterpoise corrected interaction energies are reported. Also, using a variation-perturbation energy decomposition scheme we have found the vital contributions to the total interaction energy.     

Ab Initio Quantum Mechanics Analysis of Imidazole C-H…O Water Hydrogen Bonding and a Molecular Mechanics Forcefield Correction

Abstract

While it is well established that classical hydrogen bonds play an important role in enzyme structure, function and dynamics, the role of weaker, but ‘activated’ C-H donor hydrogen bonds is poorly understood. The most important such case involves histidine which often plays a direct role in enzyme catalysis and possesses the most acidic C-H donor group of the standard amino acids. In the present study, we obtained optimized geometries and hydrogen bond interaction energies for C-H…O hydrogen bonded complexes between methane, ethylene, benzene, acetylene, and imidazole with water at the MP2-FC/6-31++G(2d,2p) and MP2-FC/aug-cc-pVDZ//MP2-FC/6-31++G(2d,2p) levels of theory. A strong linear relationship is obtained between the stability of the various hydrogen bonded complexes and both separation distances for H…0 and C—O. In general, these calculations indicate that C-H…0 interactions can be classified as hydrogen bonding interactions, albeit significantly weaker than the classical hydrogen bonds, but significantly stronger than just van der Waals interactions. For instance, while the electronic energy of stabilization at the MP2-FC/aug-cc-pVDZ//MP2-FC/6-31++G(2d,2p) level of theory of a water C-H…O water hydrogen bond is 4.36 kcal/mol more stable than the methane C-H…O water interaction, the water-water hydrogen bond is only 2.06 kcal/mol more stable than the imidazole C_e?H…O water hydrogen bond. Neglecting this latter hydrogen bonding interaction is obviously unacceptable. We next compare the potential energy surfaces for the imidazole C_e?H…O water and imidazole N_d?H…O hydrogen bonded complexes computed at the MP2/6-31++G(2d,2p) level of theory with the potential energy surface computed using the AMBER molecular mechanics program and forcefields. While the Weiner et al and Cornell et al AMBER forcefields reasonably account for the imidazole N-H…O water interaction, these forcefields do not adequately account for the imidazole C_e?H…O water hydrogen bond. A forcefield modification is offered that results in excellent agreement between the ab initio and molecular mechanics geometry and energy for this C-H…O hydrogen bonded complex.     

An assessment to evaluate the validity of different methods for the description of some corrosion inhibitors

New research and development efforts using computational chemistry in studying an assessment of the validity of different quantum chemical methods to describe the molecular and electronic structures of some corrosion inhibitors were introduced. The standard and the highly accurate CCSD method with 6-311++G(d,p), ab initio calculations using the HF/6-31G++(d,p) and MP2 with 6-311G(d,p), 6-31++G(d,p), and 6-311++G(2df,p) methods as well as DFT method at the B3LYP, BP86, B3LYP*, M06L, and M062x/6-31G++(d,p) basis set level were performed on some triazole derivatives and sulfur containing compounds used as corrosion inhibitors. Quantum chemical parameters, such as the energy of the highest occupied molecular orbital energy (E_HOMO), the energy of the lowest unoccupied molecular orbital energy (E_LUMO), energy gap (ΔE), dipole moment (μ), sum of total negative charges (TNC), chemical potential (Pi), electronegativity (χ), hardness (η), softness (σ), local softness (s), Fukui functions (f ⁺,f ^?), electrophilicity (ω), the total energy change (?E_T) and the solvation energy (S.E), were calculated. Furthermore, the accuracy and the applicability of these methods were estimated relative to the highest accuracy and standard CCSD with 6-311++G(d,p) method. Good correlations between the quantum chemical parameters and the corresponding inhibition efficiency (IE%) were found.     

An improved generalized AMBER force field (GAFF) for urea

We describe an improved force field parameter set for the generalized AMBER force field (GAFF) for urea. Quantum chemical computations were used to obtain geometrical and energetic parameters of urea dimers and larger oligomers using AM1 semiempirical MO theory, density functional theory at the B3LYP/6-31G(d,p) level, MP2 and CCSD ab initio calculations with the 6-311++G(d,p), aug-cc-pVDZ, aug-cc-pVTZ, and aug-cc-pVQZ basis sets, and with the CBS-QB3 and CBS-APNO complete basis set methods. Seven different urea dimer structures were optimized at the MP2/aug-cc-pVDZ level to obtain accurate interaction energies. Atomic partial charges were calculated at the MP2/aug-cc-pVDZ level with the restrained electrostatic potential (RESP) fitting approach. The interaction energies computed with these new RESP charges in the force field are consistent with those obtained from CCSD and MP2 calculations. The linear dimer structure calculated using the force field with modified geometrical parameters and the new RESP charge set agrees well with available experimental data.     

DFT studies on the intrinsic conformational properties of non-ionic pyrrolysine in gas phase

B3LYP/6-31G(d,p) level of theory is used to carry out a detailed gas phase conformational analysis of non-ionized (neutral) pyrrolysine molecule about its nine internal back-bone torsional angles. A total of 13 minima are detected from potential energy surface exploration corresponding to the nine internal back-bone torsional angles. These minima are then subjected to full geometry optimization and vibrational frequency calculations at B3LYP/6-31++G(d,p) level. Characteristic intramolecular hydrogen bonds present in each conformer, their relative energies, theoretically predicted vibrational spectra, rotational constants and dipole moments are systematically reported. Single point calculations are carried out at B3LYP/6-311++G(d,p) and MP2/6-31++G(d,p) levels. Six types of intramolecular H-bonds, viz. O…H–O, N…H-O, O…H–N, N…H–N, O…H–C and N…H–C, are found to exist in the pyrrolysine conformers; all of which contribute to the stability of the conformers. The vibrational frequencies are found to shift invariably toward the lower side of frequency scale corresponding to the presence of intramolecular H-bond interactions in the conformers.     

Concerning the solvent effect in the tautomerism of uracil, 5-fluorouracil, and thymine by density-functional theory and ab initio calculations

The tautomerism of uracil, 5-fluorouracil, and thymine has been investigated in the gas phase and in solution. Electron correlation effects were included in ab initio computations at the MP2 level, and DFT calculations were performed using the B3LYP level. Full geometry optimizations were conducted at the HF/6-31G**, HF/6-31+G**, and B3LYP/6-31+G** levels. Single-point MP2/6-31+G** calculations were performed on the HF/6-31+G** optimized geometries. The influence of the solvent was examined from self-consistent reaction field calculations performed with )=2.21 (1,4-dioxane) and )=78.54 (water). The calculated relative free energies ((G) indicate that substitution of uracil at the position group does not change the relative free energy order of the uracil tautomers in the gas phase and in 1,4-dioxane (except at the MP2 level) whereas this ordering changes in water. Attachment of a fluorine atom changes the relative free energy order of uracil tautomers in the gas phase and in solution.     

Evaluation of several economical computational methods for geometry optimization of phosphorus acid derivatives

Several economical methods for geometry optimisation, applicable to larger molecules, have been evaluated for phosphorus acid derivatives. MP2/cc-pVDZ and B3LYP/6-31+G(d) geometry optimizations are used as reference points, results from geometry optimizations for other methods and their subsequent single point energy calculations are compared to these references. The geometries from HF/MIDI! optimizations were close to those of the references and subsequent single point energies with B3LYP/6-31+G(d,p) or EDF1/6-31+G(d) gave a mean average deviation (MAD) of less than 0.5 kcal mol-1 from those obtained with the reference geometries.     

The reaction of the aminoboranylidene-iminoborane isomerization: a CASSCF direct dynamics study

The barrier and the potential-energy surface of the isomerization from aminoboranylidene (BNH₂) to iminoborane (HBNH) have been studied using complete active space self-consistent field (CASSCF) with the 6−31+G(d, p) basis set and higher-level energy methods. The rate constants of the isomerization reaction are reported by employing the direct ab initio dynamics method. The geometries of all the stationary points were optimized using the B3LYP and CCSD methods with the cc-pVTZ and cc-pVQZ basis sets. The information along the intrinsic reaction coordinate (IRC) was also calculated at the CASSCF/6−31+G (d,p) level of theory. The energies were refined at the G3, G3MP2, G3MP2B3, CBS-Q, CBS-QB3, and two high-level (HL) methods based on the geometries optimized using CASSCF/6-31+G(d,p). The rate constants were evaluated using conventional transition-state theory (TST), canonical variational transition-state theory (CVT), and canonical variational transition-state theory with small curvature tunneling correction (CVT/SCT) and conventional transition-state theory with Eckart tunneling correction (TST/Eckart). According to the calculated results, we conclude that the tunneling effect is very important to this isomerization reaction.     

Theoretical potential functions and vibrational analysis for halocarbonyl azides CXO-NNN (X=F,Cl and Br)

The structural stability of halocarbonyl azides CXO-NNN (X=F, Cl and Br) was investigated by DFT and MP2 calculations using the 6-311++G** basis set. From the calculations, the molecules were found to have an s-cis<--> s-trans conformational equilibrium with cis being the lower -energy form. Full energy optimizations were carried out for the transition states and the minima at the B3LYP/6 -311++G** and MP2/6 -311++G** levels, from which the rotational barriers were calculated to be of the order 8-10 kcal x mol(-1). The vibrational frequencies were computed at the DFT -B3LYP level and the vibrational assignments for the normal modes of the stable conformers were made on the basis of normal coordinate calculations.     

EVALUATION OF SEVERAL ECONOMICAL COMPUTATIONAL METHODS FOR GEOMETRY OPTIMISATION OF PHOSPHORUS ACID DERIVATIVES

Several economical methods for geometry optimisation, applicable to larger molecules, have been evaluated for phosphorus acid derivatives. MP2/cc-pVDZ and B3LYP/6-31+G(d) geometry optimisations are used as reference points, results from geometry optimisations for other methods and their subsequent single point energy calculations are compared to these references. The geometries from HF/MIDI! optimisations were close to those of the references and subsequent single point energies with B3LYP/6-31+G(d,p) or EDF1/6-31+G(d) gave a mean average deviation (MAD) of less than 0.5 kcal mol^?1 from those obtained with the reference geometries.     

A theoretical study of cooperative and anticooperative effects on hydrogen-bonded clusters of water and the cyanuric acid

Ab initio and density functional calculations are used to analyse the interaction between a molecule of the cyanuric acid and one, two and three molecules of water at B3LYP/6-311++ G(d,p) and MP2/6-311++ G(d,p) computational levels. Also, the cooperative effect (CE) in terms of the stabilisation energy of clusters is calculated and discussed. Depending on the geometry of clusters under study, the cooperative, non- or anti-CE was found with an increasing cluster size. Red shifts of N–H and C = O stretching frequencies illustrate a good dependence on the CE. The atoms in molecules theory is used to analyse the CE on topological parameters.     

MP2, density functional theory, and molecular mechanical calculations of C-H···π and hydrogen bond interactions in a cellulose-binding module-cellulose model system

Exploring non-covalent interactions, such as C-H···π stacking and classical hydrogen bonding (H-bonding), between carbohydrates and carbohydrate-binding modules (CBMs) is an important task in glycobiology. The present study focuses on intermolecular interactions, such as C-H?π (sugar-aromatic stacking) and H-bonds, between methyl β-d-glucopyranoside and l-tyrosine—a proxy model system for a cellulose-CBM complex. This work has made use of various types of quantum mechanics (QM) and molecular mechanics (MM) methods to determine which is the most accurate and computationally efficient. The calculated interaction potential energies ranged between −24 and −38 kJ/mol. The larger interaction energy is due to H-bonding between the phenyl hydroxyl of tyrosine and the O4 of the sugar. Density functional theory (DFT) methods, such as BHandHLYP and B3LYP, exaggerate the H-bond. Although one of the MM methods (viz. MM+) considered in this study does maintain the C-H?π stacking configuration, it underestimates the interaction energy due to the loss of the H-bond. When the O-H bond vector is in the vicinity of O4 (O-H?O4 ≈ 2 Å, e.g., in the case of MP2/6-31G(d)), the torsional energy drops to a minimum. For this configuration, natural bond orbital (NBO) analysis also supports the presence of this H-bond which arises due to orbital interaction between one lone pair of the sugar O4 and the σ∗(O-H) orbital of the phenyl group of tyrosine. The stabilization energy due to orbital delocalization of the H-bonded system is ∼13 kJ/mol. This H-bond interaction plays an important role in controlling the CH/π interaction geometry. Therefore, the C-H?π dispersive interaction is the secondary force, which supports the stabilization of the complex. The meta-hybrid DFT method, M05-2X, with the 6-311++G(d,p) basis set agrees well with the MP2 results and is less computationally expensive. However, the M05-2X method is strongly basis set dependent in describing this CH/π interaction. Computed IR spectra with the MP2/6-31G(d) method show blue shifts for C1-H, C3-H, and C5-H stretching frequencies due to the C-H?π interaction. However, the M05-2X/6-311++G(d,p) method shows a small red shift for the C1-H stretching region and blue shifts for the C2-H and C3-H stretches. For the aromatic tyrosine C_δ1-C_ε1 and C_δ2-C_ε2 bonds in the complex, the calculated IR spectra show red shifts of 12 cm⁻¹ (MP2/6-31G(d)) and 5 cm⁻¹ (M05-2X/6-311++G(d,p)). This study also reports the upfield shifts of computed ¹H NMR chemical shifts due to the C-H?π interaction.     

The influence of the peptide bond on the conformation of amino acids: a theoretical and FT-IR matrix-isolation study of N-acetylproline

A combined experimental matrix-isolation FT-IR and theoretical study has been performed to investigate the conformational behavior of N-acetylproline. The conformational landscape of N-acetylproline was explored using successively higher computational methods, i.e. HF, DFT(B3LYP) and finally MP2. The exploration resulted in 10 conformations with a relative energy difference smaller than 22 kJ.mol(-1) at the HF/3-21G level of theory. These conformations led to six different conformations after DFT(B3LYP) optimizations. Further optimization at the MP2/6-31++G** level of theory resulted in the same six conformations, all of them with an energy difference smaller than 11.5kJ.mol(-1). One conformation with an intramolecular H-bond was found which was energetically the most favorable conformation. The vibrational and thermodynamical features were calculated using the DFT and MP2 methodologies. In the experimental matrix-isolation FT-IR spectrum, the most stable conformation was dominant and at least two non-H-bonded conformations could be identified. An experimental rotamerization constant between the H-bonded and the other non-H-bonded conformations was estimated and appeared to agree reasonably well with the theoretical MP2 predictions. Some new spectral features of N-acetylproline compared to proline were discovered which might be used to discriminate between the acetylated and non-acetylated form.     

Density functional computations of enantioselective alkynylation of aldehyde catalyzed by chiral zinc(II)-complexes

The enantioselective alkynylation of aldehyde catalyzed by chiral zinc(II)-complexes was studied by means of the density functional theory (DFT). All the structures were optimized completely at the B3LYP/6-31G(d,p) level. To obtain more exact energies, single-point energy calculations at B3LYP/6-31+G(d,p) level were carried out on the B3LYP/6-31G(d,p) geometries. As shown, this enantioselective alkynylation was endothermic. The chirality-determining step for the alkynylation was the formation of the catalyst–ethanol complexes and the transition states for this step involved a six-membered ring. The dominant products predicted theoretically were of (R)-chirality, in good agreement with experiment.      

Low-energy conformers of pamidronate and their intramolecular hydrogen bonds: a DFT and QTAIM study

Extensive DFT and ab initio calculations were performed to characterize the conformational space of pamidronate, a typical pharmaceutical for bone diseases. Mono-, di- and tri-protic states of molecule, relevant for physiological pH range, were investigated for both canonical and zwitterionic tautomers. Semiempirical PM6 method were used for prescreening of the single bond rotamers followed by geometry optimizations at the B3LYP/6-31++G(d,p) and B3LYP/6-311++G(d,p) levels. For numerous identified low energy conformers the final electronic energies were determined at the MP2/6-311++G(2df,2p) level and corrected for thermal effects at B3LYP level. Solvation effects were also considered via the COSMO and C-PCM implicit models. Reasonable agreement was found between bond lengths and angle values in comparison with X-ray crystal structures. Relative equilibrium populations of different conformers were determined from molecular partition functions and the role of electronic, vibrational and rotational degrees of freedom on the stability of conformers were analyzed. For no level of theory is a zwitterionic structure stable in the gas-phase while solvation makes them available depending on the protonation state. Geometrically identified intramolecular hydrogen bonds were analyzed by QTAIM approach. All conformers exhibit strong inter-phosphonate hydrogen bonds and in most of them the alkyl-amine side chain is folded on the P-C-P backbone for further hydrogen bond formation.

Stability Order of Basic Peptide Conformations Reflected by Density Functional Theory

Ab initio density functional methods at the B3LYP/6-311+G(d, p) and 6-31G(d) levels were performed on several basic peptide conformations representing typical elements of secondary structure (β-sheets, β- and γ-turns). The results are compared with those from Hartree-Fock and MP2 correlation energy calculations. Whereas the geometries of the structures are well described at all approximation levels, there are considerable discrepancies of the stability orders. Contrary to the Hartree-Fock calculations, the correlation energy methods provide the more compact structures with intramolecular hydrogen bonds distinctly favoured over extended conformations when comparing the energy differences. However, due to considerable compensation of correlation energy and entropy contributions, the stability order at the Gibbs free energy level closely corresponds to that at the Hartree-Fock level.