A theoretical study of the cis-syn pyrimidine dimers in the gas phase and water cluster and a tautomer-bypass mechanism for the origin of UV-induced mutations

A quantum mechanical study of all cis-syn cyclobutane pyrimidine photodimers including the normal and rare tautomeric forms of bases has been performed using the ab initio method at HF/6-31G(d.p), MP2(fc)//HF/6-31G(d,p) and MP2(fc)/6-31G(d,p) levels. A puckering angle of the cyclobutyl ring and twist angle of pyrimidine rings with respect to each other is well described by these calculations. It is predicted that in the gas phase all photodimers containing the rare imino form of cytosine are more stable than those containing its normal form. The Monte Carlo simulations show that the dimer containing the imino form of cytosine is more stabilized by water cluster than that containing its amino forms. The possible biological significance stems from the fact that the cytosine in the dimer directs the incorporation of adenine in the complementary strand during replicative bypass. Data obtained point to the cytosine tautomerism as a possible mechanism for the origin of UV-induced mutation.     

Interactions Between Amino Groups in DNA. An Ab Initio Study and a Comparison with Empirical Potentials

Abstract

An ab initio quantum chemical analysis of the close amino group contacts, existing in many DNA crystal structures, is presented. The calculations are made at the Hartree-Fock (HF) level with medium 6–31G* and 6–31 G(NH₂*) basis sets as well as with inclusion of correlation energy using the second order Møller-Plesset theory (MP2) with the 6–31G* basis set. We demonstrate that the model system (methylamine dimer, cytosine dimer) amino groups are forced to adopt significantly non-planar geometry to stabilize their mutual interaction. Comparison is made with a representative set of empirical potentials including AMBER, CHARMM and GROMOS. The empirical potentials are not reliable enough to analyze the amino group contacts occurring in the DNA double helices. We propose that the mutual amino group interactions contribute to the conformational variability of the CpG and ApT B-DNA steps.     

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.     

Computational Study of the Sulfonylated Amino Acid Hydroxamates Binding to the Zinc Ion within the Active Site of Carbonic Anhydrase

Abstract

The B3LYP/6–311+G(d,p) method and three ONIOM extrapolation methods ONI-OM (B3LYP/6–311+G(d,p): AM1); ONIOM(B3LYP/6–311+G(d,p): MNDO); ONIOM (B3LYP/6–311+G(d,p): HF/3-21G(d)) were used to characterize the complexes of Zn²⁺ cation with anionic sulfonylated amino acid hydroxamates (RSO₂NH-AA-CON(-)OH), possessing an unsubstituted RSO₂NH—amino acyl moiety. According to the R moiety we distinguish between pentafluorophenyl and 4-methoxyphenyl derivates. The amino acid hydroxamates included in the study were the Gly, Ala, and Leu derivates. Of the inhibitors investigated, the weakest zinc affinity exhibits the pentafluorophenyl derivate with Gly amino acid and the strongest affinity the 4-methoxyphenyl derivate with Leu amino acid. The inhibitors form bidentate coordination bonds with the zinc cation by means of the sulfonyl oxygen and the ionized hydroxamate nitrogen atoms, respectively. The zinc affinities computed using the B3LYP/6–311 +G(d,p)//HF/6–31 +G(d,p) method are in very good agreement with the full density functional theory (DFT) B3LYP/6–311+G(d,p)//B3LYP/6- 311+G(d,p) method and this method can be adopted to model larger complexes of inhibitors with the active site of carbonic anhydrase.     

Post Hartree-Fock Studies of the Canonical Watson-Crick DNA Base Pairs: Molecular Structure and the Nature of Stability

Abstract

Gas-phase gradient optimization was carried out on the canonical Watson-Crick DNA base pairs using the second-order Møller-Plesset perturbation method at the 6–31G(d) and 6- 31G(d,p) basis sets. It is detected that full geometry optimization at the MP2 level leads to an intrinsically nonplanar propeller-twisted and buckled geometry of G-C and A-T base pairs; while HF and DFT methods predict perfect planar or almost planar geometry of the base pairs. Supposedly the nonplanarity of the pairs is caused by pyramidalization of the amino nitrogen atoms, which is underestimated by the HF and DFT methods. This justifies the importance of geometry optimization at the MP2 level for obtaining reliable prediction of the charge distribution, molecular dipole moments and geometrical structure of the base pairs. The Morokuma-Kitaura and the Reduced Variational Space methods of the decomposition for molecular HF interaction energies were used for investigation of the hydrogen bonding in the Watson-Crick base pairs. It is shown that the HF stability of the hydrogen-bonded DNA base pairs originates mainly from electrostatic interactions. At the same time, the calculated magnitude of the second order intramolecular correlation correction to the Coulomb energy showed that electron correlation reduces the contribution of the electrostatic term to the attractive interaction for the A-T and G-C base pairs. Polarization, charge transfer and dispersion interactions also make considerable contribution to the attraction energy of bases.     

Theoretical study of geometrical and nonlinear optical properties of pyridinum N-phenolate betaine dyes

This paper presents an ab initio quantum chemical investigation of the geometrical structures and the non-linear optical properties (NLO) of three structural isomers of pyridinium N-phenolate betaine dye. The ground state geometrical parameters and the first-order hyperpolarizabilities were calculated using the Hartree-Fock (HF) as well as the second-order perturbation Møller-Pleset (MP2) method with the 6–31G, 6–31G(d), 6–31G(d,p), 6–31+G(d), 6–31++G(d,p), 6–311+G(d), aug-cc-PVDZ and the recently developed Z3PolX basis sets. Moreover, the first-order hyperpolarizability was calculated at the coupled cluster singles and doubles (CCSD/6–31+G(d)) level of theory. The analysis of the results of calculations for the investigated isomers indicates that there are important differences in their NLO activities. Additionally, it was shown that Z3PolX basis set works reasonable well for betaine dyes.

Ab initio Investigation of the Molecular Structure of Methyl Methoxymethyl Phosphonate,a Promising Nuclease-resistant Alternative of the Phosphodiester Linkage

Abstract

Conformational flexibility of the methyl methoxymethyl phosphonate anion (CH₃-O-PO₂- CH₂-O-CH₃)^?, a nuclease resistant alternative to the phosphodiester linkage in DNA, have been investigated by ab initio quantum mechanical calculations. The potential of backbone torsional degrees of freedom of methyl methoxymethyl phosphonate anion (MMP) was determined at the Hartree-Fock (HF) 3–21G* level using the adiabatic mapping technique. Energies, geometries, and effective atomic charges of different conformers were calculated at HF/6–31G* and MP2/6–31G* levels of theory. These were compared to the results obtained for dimethyl phosphate calculated at the same level. The impact on DNA structure from inserting a methylene group between phosphorus and oxygen of the nucleoside sugar moiety was examined via distance- and angle-constrained geometry optimizations. Due to its high flexibility, MMP has been shown to be compatible with both A and B forms of DNA.     

Studies on molecular structure and tautomerism of a vitamin B<Subscript>6</Subscript> analog with density functional theory

This work presents a computational study on the molecular structure and tautomeric equilibria of a novel Schiff base L derived from pyridoxal (PL) and o-phenylenediamine by using the density functional method B3LYP with basis sets 6-31 G(d,p), 6-31++G(d,p), 6-311 G(d,p) and 6-311++G(d,p). The optimized geometrical parameters obtained by B3LYP/6-31 G(d,p) method showed the best agreement with the experimental values. Tautomeric stability study of L inferred that the enolimine form is more stable than its ketoenamine form in both gas phase and solution. However, protonation of the pyridoxal nitrogen atom (LH) have accelerated the formation of ketoenamine form, and therefore, both ketoenamine and enolimine forms could be present in acidic media.     

Theoretical Investigation of the Molecular Structure of the πκ DNA Base Pair

Abstract

The structure of the nonclassical πκ base pair (7–methyl-oxoformycin … 2,4-diaminopyrimidine) was studied at the ab initio Hartree-Fock (HF) and MP2 levels using the 6–31G* and 6–31G** basis sets. The πκ base pair is bound by three parallel hydrogen bonds with the donor-acceptor-donor recognition pattern. Recently, these bases were proposed as an extension of the genetic alphabet from four to six letters (Piccirilli et al. Nature 343, 33(1990)). By the HF/6- 31G* method with full geometry optimization we calculated the 12 degree propeller twist for the minimum energy structure of this complex. The linearity of hydrogen bonds is preserved in the twisted structure by virtue of the pyramidal arrangement of the κ-base amino groups. The rings of both the π and κ molecules remain nearly planar. This nonplanar structure of the πκ base pair is only 0.1 kcal/mol more stable than the planar (Cs) conformation. The HF/6- 31G* level gas-phase interaction energy of πκ (—13.5 kcal/mol) calculated by us turned out to be nearly the same as the interaction energy obtained previously for the adenine-thymine base pair (—13.4 kcal/mol) at the same computational level. The inclusion of p-polarization functions on hydrogens, electron correlation effects (MP2/6–31G** level), and the correction for the basis set superposition error (BSSE) increase this energy to -14.0 kcal/mol.     

Structure,energetics and properties of some molecules with potent anti-HIV activity: a theoretical study

Density functional theory (DFT; B3LYP) and Hartree–Fock (HF; 3-21G, 6-31G(d) and 6-311G(d,p)) calculations with complete geometry optimisations are carried out in the ground state on five 6-aminoquinolone derivatives, which have been proved to be highly effective in inhibiting HIV replication, to study their structures, energetics and HOMO–LUMO correlation with physiological action. The gas-phase calculations and single-point polarisable continuum model water-phase calculations show that the molecules are highly effective in inhibiting HIV replication, which is in excellent agreement with the experiment. Structural features, energies, charge densities and HOMO–LUMO correlation have been found to substantiate the experimental findings. Compound 4 (pyrazine) shows some special features in DFT calculations which are not found in HF calculations. In the present series, HF results are more reliable as expected.     

Molecular Structure and Properties of Protonated and Methylated Derivatives of Cytosine

Abstract

Ab initio techniques were used to determine the effects of protonation and methylation on cytosine's molecular geometry, molecular affinities, relative stability, and structural rigidity. The geometries of local minima were fully optimized by the gradient procedure at DFT and MP2 levels of theory with the medium size 6–31G(d,p) basis set. The results of energetic analysis indicate that N(3)-methyl-cytosine and C(5)-methyl-cytosine are the most stable derivatives for monocationic and neutral species, respectively. The structural rigidity of each species was assessed by an analysis of normal out-of-plane frequencies, the amplitudes, and by the contribution of internal coordinates to the potential energy distributions. The obtained evidence suggests that methylation increases the overall structural flexibility of cytosine and that all molecules in this study populate a non-planar conformation 50% of the time.     

Structural and vibrational spectroscopic elucidation of sulpiride in solid state

The study on the conformational and vibrational behaviors of sulpiride molecule which is known as a neuroleptic or antipsychotic drug that is widely used clinically in the treatment of schizophrenic or depressive disorders is an important scientific and practical task. In here, a careful enough study of monomer and dimeric forms of sulpiridine {5-(aminosulfonyl)-N-[(1-ethyl-2-pyrrolidinyl) ethyl]-2-methoxy-benzamide (C₁₅H₂₃N₃O₄S)} is undertaken by density functional theory (DFTB3LYP) method with the B3LYP/6-31G(d,p) basis set. The conformations of free molecule were searched by means of torsion potential energy surfaces scan studies through dihedral angles D₁ (8?N, 18C, 20C, 23?N), D₂ (18C, 20C, 23?N, 25C) and D₃ (28C, 30C, 41S, 44?N) in electronically ground state, employing 6-31G basic set. The final geometrical parameters for the obtained stable conformers were determined by means of geometry optimization, carried out at DFT/B3LYP/6-31G(d,p) theory level. Afterwards, the possible dimer forms of the molecule were formed and their energetically preferred conformations were investigated. Moreover, the effect of basis set superposition error on the structure and energy of the three energetically favourable sulpiride dimers has been determined. The optimized structural parameters of the most stable monomer and three low energy dimer forms were used in the vibrational wavenumber calculations. Raman and IR (4000–400?cm^?1) spectra of sulpiride have been recorded in the solid state. The assignment of the bands was performed based on the potential energy distribution data. The natural bond orbital analysis has been performed on both monomer and dimer geometries in order to elucidate delocalization of electron density within the molecule. The predicted frontier molecular orbital energies at DFT/B3LYP/6-31G(d,p) theory level show that charge transfer occurs within the molecule. The first-order hyperpolarizability (β₀) and related properties (μ and α) of the title molecule were also calculated.     

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.     

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.     

Post Hartree-Fock studies of the canonical Watson-Crick DNA base pairs: molecular structure and the nature of stability

Gas-phase gradient optimization was carried out on the canonical Watson-Crick DNA base pairs using the second-order M?ller-Plesset perturbation method at the 6-31G(d) and 6-31G(d,p) basis sets. It is detected that full geometry optimization at the MP2 level leads to an intrinsically nonplanar propeller-twisted and buckled geometry of G-C and A-T base pairs; while HF and DFT methods predict perfect planar or almost planar geometry of the base pairs. Supposedly the nonplanarity of the pairs is caused by pyramidalization of the amino nitrogen atoms, which is underestimated by the HF and DFT methods. This justifies the importance of geometry optimization at the MP2 level for obtaining reliable prediction of the charge distribution, molecular dipole moments and geometrical structure of the base pairs. The Morokuma-Kitaura and the Reduced Variational Space methods of the decomposition for molecular HF interaction energies were used for investigation of the hydrogen bonding in the Watson-Crick base pairs. It is shown that the HF stability of the hydrogen-bonded DNA base pairs originates mainly from electrostatic interactions. At the same time, the calculated magnitude of the second order intramolecular correlation correction to the Coulomb energy showed that electron correlation reduces the contribution of the electrostatic term to the attractive interaction for the A-T and G-C base pairs. Polarization, charge transfer and dispersion interactions also make considerable contribution to the attraction energy of bases.     

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.     

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.     

Theoretical Ab initio study of the effects of methylation on the nature of hydrogen bonding in A:T base pair

We report the results of a theoretical ab initio study of methylation in Watson-Crick A:T base pairs. Equilibrium geometries were obtained without symmetry restrictions by the gradient procedure at DFT level of theory with the standard 6-31G(d) basis set. Each local minima was verified by energy second derivative calculations. Single-point calculations for the DFT geometries 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 stabilities and counterpoise corrected interaction energies are reported. In addition, using a variation-perturbation energy decomposition scheme, we have found the important contributions to the total interaction energy.