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.     

Quantum-chemical investigation of tautomerization ways of Watson-Crick DNA base pair guanine-cytosine

A novel physico-chemical mechanism of the Watson-Crick DNA base pair Gua.Cyt tautomerization Gua.Cyt*<---->Gua.Cyt<---->Gua*.Cyt (mutagenic tautomers of bases are marked by asterisks) have been revealed and realized in a pathway of single proton transfer through two mutual isoenergetic transition states with Gibbs free energy of activation 30.4 and 30.6 kcal/mol and they are ion pairs stabilized by three (N2H...N3, N1H...N4- and O6+H...N4-) and five (N2H...O2, N1H...O2, N1H...N3, O6+H...N4- and 06+H...N4-) H-bonds accordingly. Stable base pairs Gua-Cyt* and Gua*.Cyt which dissociate comparably easy into monomers have acceptable relative Gibbs energies--12.9 and 14.3 kcal/mol--for the explanation of the nature of the spontaneous transitions of DNA replication. Results are obtained at the MP2/6-311++G(2df,pd)//B3LYP/6-31 1++G(d,p) level of theory in vacuum approach.     

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.     

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.     

QM/MM study of the active site of free papain and of the NMA-papain complex.

Hybrid quantum mechanical/molecular mechanical (QM/MM) calculations using restricted and unrestricted Hartree-Fock and B3LYP ab initio (QM) and Amber force field (MM), respectively, have been applied to study the catalytic site of papain in both free and substrate bonded forms. Ab initio geometry optimizations have been performed for the active site of papain and the N-methyl-acetamide (NMA)-papain complex within the molecular mechanical treatment of the protein environment. A covalent tetrahedral intermediate structure could be obtained only when the amide N atom of the substrate molecule was protonated through a proton transfer from the His-159 in the catalytic site. Our results support the previous assumption that a proton transfer from His-159 to the amide N atom of the substrate occurs prior to or concerted with the nucleophilic attack of the Cys-25 sulfur atom to the carbonyl group of the substrate. The electron correlation effect will reduce the proton transfer barrier. Therefore, this proton transfer can be easily observed in the B3LYP/6-31G* calculations. The HF/6-31G* method overestimates the reaction barrier against this proton transfer. The sulfur atom of Cys-25 and the imidazole ring of His-159 are found to be coplanar in the free form of the enzyme. However, the rotation of the imidazole ring of His-159 was observed during the formation of the tetrahedral intermediate. Without the papain environment, the coplanar thiolate-imidazolium ion pair RS-...ImH+ is much less stable than the neutral form of RSH....Im. Within the protein environment, however, the thiolate-imidazolium ion pair becomes more stable than its neutral form by 4.1 and 0.4 kcal/mol in HF/6-31G* and B3LYP/6-31G* calculations, respectively. The barrier of proton transfer from S-H group of Cys-25 to the imidazole ring of His-159 was reduced from 22.0 kcal/mol to 15.2 kcal/mol by the protein environment in HF/6-31G* calculations. This barrier is found to be much smaller (2.5 kcal/mol) in B3LYP/6-31G* calculations.     

Molecular mechanisms of transitions induced by cytosine analogue: comparative quantum-chemical study

Using the simplest molecular models at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of the theory it has been shown for the first time that in addition to traditional incorporational errors caused by facilitated (compared with the canonical DNA bases cytosine (Cyt)) tautomerization of 6-(2-deoxy-beta-D-ribofuranosyl)-3,4-dihydro-6H,8H-pyrimido[4,5-c][1,2]oxazin-7-one (DCyt), this mutagen causes the replication errors, increasing one million times the population of mispair Gua.DCyt* (asterisk marked mutagenic tautomer) as compared with mispair Gua.Cyt*. It is also proved that DCyt in addition to traditional incorporational errors also induces similar errors by an additional mechanism - due to a facilitated tautomerization of the wobble base pair Ade.DCyt (compared to the same pair Ade.Cyt) to a mispair Ade.DCyt* which is quasirisomorphic Watson-Crick base pair. Moreover, the obtained results allowed interpreting non-inconsistently the existing experimental NMR data.     

An investigation of the pyranose ring interconversion path of alpha-L-idose calculated using density functional theory

The interconversion pathways of the pyranose ring conformation of alpha-L-idose from a (4)C1 chair to other conformations were investigated using density functional calculations. From these calculations, four different ring interconversion paths and their transition state structures from the (4)C1 chair to other conformations, such as B(3,O), and (1)S3, were obtained. These four transition-state conformations cover four possible combinations of the network patterns of the hydroxyl group hydrogen bonds (clockwise and counterclockwise) and the conformations of the primary alcohol group (tg and gg). The optimized conformations, transition states, and their intrinsic reaction coordinates (IRC) were all calculated at the B3LYP/6-31G** level. The energy differences among the structures obtained were evaluated at the B3LYP/6-311++G** level. The optimized conformations indicate that the conformers of (4)C1, (2)S(O), and B(3,O) have similar energies, while (1)S3 has a higher energy than the others. The comparison of the four transition states and their ring interconversion paths, which were confirmed using the IRC calculation, suggests that the most plausible ring interconversion of the alpha-L-idopyranose ring occurs between (4)C1 and B(3,O) through the E3 envelope, which involves a 5.21 kcal/mol energy barrier.     

Catalytic mechanism of the inverting N-acetylglucosaminyltransferase I: DFT quantum mechanical model of the reaction pathway and determination of the transition state structure

The complex N-glycan structures on glycoproteins play important roles in cell adhesion and recognition events in metazoan organisms. A critical step in the biosynthetic pathway leading from high mannose to these complex structures includes the transfer of N-acetylglucosamine (GlcNAc) to a mannose residue by the inverting N-acetylglucosaminyltransferase I (GnT-I). The catalytic mechanism of this enzymatic reaction is explored herein using DFT quantum chemical methods. The computational model used to follow the reaction is based on the X-ray crystallographic structure of GnT-I and contains 127 atoms that represent fragments of residues critical for the substrate binding and catalysis. The mechanism of the catalytic reaction was monitored by means of a 2D potential energy map calculated as a function of predefined reaction coordinates at the B3LYP/6-31G** level. This potential energy surface revealed one transition state associated with a reaction pathway following a concerted mechanism. The reaction barrier was estimated, and the structure of the transition state was characterized at the B3LYP/6-311++G**// B3LYP/6-31G** level.     

Effect of carboxylate and Na ions on the tautomeric status of 9-methylguanine

Interactions of 9-methylguanine (m9Gua) with carboxylate ion of acetic acid (CH3COO-) and Na+ were studied by 1H NMR spectroscopy and ab initio quantum chemical calculations of the B3LYP/6-31++G(d,p) and B3LYP/6-311++G(d,p) levels of theory. Changes in the m9Gua 1H NMR spectrum in the presence of the equimolar amount of sodium acetate (NaAc), which in anhydrous DMSO dissociates to CH3COO- and Na+, were interpreted as a consequence of a complex formation of m9Gua in the amino-keto-N1H tautomeric form (m9GuaN1H) with carboxylate ions via two H-bonds involving amino and N1H-imino protons. Quantum chemical calculations of interactions of the m9GuaN1H ground-state tautomer and the m9GuaN3H high energy one with relative energy 20.01 kcal/mol show that the ground state tautomer forms the ground-state complex CH3COO-:mgGuaN1H, by 5.57 kcal/mol more stable than the CH3COO-:m9GuaN3H complex, and coordination of Na+ with the O6 and N7 atoms reduces this energy difference to 2.57 kcal/mol. Such a coordination of Na+ with tautomer m9GuaN3H therewith decreases its relative energy only to 13.31 kcal/mol. Non-additivity of the two ligands contributions to the 8-times reduction of the relative energy of the high energy tautomer in the CH3COO-:m9GuaN3H:Na+(O6,N7) triple complex was concluded, the role of CH3COO- being dominant. Besides, coordination with Na+ resulted in an iminoproton transfer from the base to CH3COO- in the triple complexes of both tautomers, according to calculations in vacuum. Biological significance of the results is noticed.     

Structure and the energy of base pairing in non-natural bases of nucleic acids: the azaguanine-cytosine and azaadenine-thymine base pairs

Watson-Crick optimized geometries and the energies of base pairing for the natural pairs of nucleic bases: adenine-thymine (AT) and guanine-cytosine (GC) have been recalculated by ab initio methods in order to compare results to those found for the non-natural azaadenine-thymine (AAT) and azaguanine-cytosine (AGC) pairs. Geometry optimizations carried out at the HF/6-31G** level and energies obtained at MP2/6-31G**, show that AAT and AGC have hydrogen bonding patterns similar to the natural AT and GC and that the interaction energies (DeltaH0int) for the former are ca. 7 kcal/mol more stable than the latter. Accordingly, the pairs based on azapurines would be favored with respect to the natural pairs. Some possible explanations why nature does not use extensively the azabases in base pairing are given.     

Enantiomerization of Environmentally Significant Overcrowded Polychlorinated Biphenyls (PCBs)

The rotational barriers of overcrowded PCBs are predicted by ab initio methods including electron correlation, thus settling the controversy between theory and experiment. For 2,2′,3,3′,4,6′-hexachlorobiphenyl ( PCB 132 ), an enantiomerization barrier of 185 kJ/mol is calculated by B3LYP/6-31G*, in excellent agreement with the experimental data. Chirality 9:350–353, 1997. © 1997 Wiley-Liss, Inc.     

Structural, IR, and EPR studies of the bis(methoxyacetato)diaquo-copper(II) complex

The structure, stability, and the IR, and EPR spectroscopic properties of bis(methoxyacetato)diaquo-copper(II) were studied both experimentally using FT-IR and theoretically using B3LYP/6-31G**, B3LYP/6-311G, BWP91/6-31G** methods. The same approaches were used to calculate the harmonic frequencies and to compare them to the experimental solid state values. The g-tensors are calculated using the NMR/GIAO computational method.     

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.     

Estimation of strength in different extra Watson-Crick hydrogen bonds in DNA double helices through quantum chemical studies

It was shown earlier, from database analysis, model building studies, and molecular dynamics simulations that formation of cross-strand bifurcated or Extra Watson-Crick hydrogen (EWC) bonds between successive base pairs may lead to extra rigidity to DNA double helices of certain sequences. The strengths of these hydrogen bonds are debatable, however, as they do not have standard linear geometry criterion. We have therefore carried out detailed ab initio quantum chemical studies using RHF/6-31G(2d,2p) and B3LYP/6-31G(2p,2d) basis sets to determine strengths of several bent hydrogen bonds with different donor and acceptors. Interaction energy calculations, corrected for the basis set superposition errors, suggest that N-H...O type bent EWC hydrogen bonds are possible along same strands or across the strands between successive base pairs, leading to significant stability (ca. 4-9 kcal/mol). The N-H...N and C-H...O type interactions, however, are not so stabilizing. Hence, consideration of EWC N-H...O H-bonds can lead to a better understanding of DNA sequence directed structural features.     

Ab initio and DFT study of the inner mechanism and dynamic stereochemistry of electrophilic addition reaction of bromine to bisbenzotetracyclo[6.2.2.23,6.02,7]tetradeca-4,9,11,13-tetraene

The inner mechanism and dynamic stereochemistry of electrophilic addition of bromine to bisbenzotetracyclo[6.2.2.2(3,6).0(2,7)]tetradeca-4,9,11,13-tetraene(BBTT) molecule have been investigated by the methods of quantum chemistry. The structure of the BBTT molecule has been studied by ab initio and DFT/B3LYP methods using the 6-31G(d) and 6-311G(d) basis sets. The double bonds of BBTT molecule are endo-pyramidalized. The structure and stability of the cationic intermediates and products of the addition reaction have been investigated by HF/6-311G(d), HF/6-311G(d,p), B3LYP/6-311G(d) and B3LYP/6-311++G(2d,p)//B3LYP/6-311G(d) methods. The bridged bromonium cation isomerized into the more stable nonclassical delocalized N- and U-type cations and the difference between the stability of these cations is small. For the determination of the direction of addition reaction and the stereochemistry of the products, the stability of nonclassical delocalized N- and U-type ions and the structure of their cationic centres play a vital role. Since the cationic centre of the N-type ion is in interaction with the benzene ring from the exo face, the nucleofilic attackof the bromide anion to this centre occurs from the endo face and the exo,endo-isomer of the N-type product is obtained. The attack of bromide anion, towards the cationic centre of U-type ion from the endo face is sterically hindered by the hydrogen atom therefore the attack occurs from the exo face, which interacts with the benzene ring and the more stable exo,exo-isomer of U-type product is formed. Although, the U-type cation was 2.232 kcal mol(-1) more stable than the N-type cation, the U-type product was 0.587 kcal mol(-1) less stable than the N-type product.     

Quantum mechanical and NMR spectroscopy studies on the conformations of the hydroxymethyl and methoxymethyl groups in aldohexosides

The potential energy surfaces of the hydroxymethyl and methoxymethyl groups in methyl hexopyranosides have been extensively studied, employing quantum mechanical calculations and high resolution NMR data. The structure and energy of the C-5-C-6 rotamers were calculated at the B3LYP level of the density functional theory (DFT). For all, geometry optimizations were carried out for 264 conformers of 16 methyl D-gluco- and methyl D-galactopyranoside derivatives 1-16 at the B3LYP/6-31G** level. For all calculated minima, single-point calculations were performed at the B3LYP/6-311++G** level. Solvent effects were considered using a self-consistent reaction field method. Values of the vicinal coupling constants 3J(H-5-H-6R), 3J(H-5-H-6S), 3J(C-4-H-6R), and 3J(C-4-H-6S) for methyl D-glucopyranosides, methyl D-galactopyranosides and their 6-O-methyl derivatives 9-16 were measured in two solvents, methanol and water. The calculated gg, gt, and tg rotamer populations of the hydroxymethyl and methoxymethyl groups in 9-16 agreed well with experimental data. The results clearly showed that the population of gg, gt, and tg rotamers is sensitive to solvent effects. It was concluded that the preference of rotamers in 1-16 is due to the hydrogen bonding and solvent effects.     

Configurational stability of chiral lithiated cyclopropylnitriles: a density functional study

Chiral, configurationally stable lithiated nitriles would be valuable intermediates for asymmetric carbon-carbon bond-forming reactions. To gain insight into the design of such species, Walborsky's attempted enantioselective deprotonation/trapping reactions of a chiral cyclopropylnitrile were studied computationally up to the MP2(fc)/6-31+G* and B3LYP/6-31+G* levels. Investigation of cyclopropylnitrile/LiNH(2) deprotonation transition structures demonstrated a significant (20-23 kcal/mol) kinetic preference for N-lithiation, and a facile (4-6 kcal/mol barrier) "conducted tour" racemization pathway for the N-lithiated nitrile product. Addition of a model directing group (formyl) to the beta-carbon of the cyclopropyl ring is predicted to significantly favor C-lithiation over N-lithiation, both kinetically and thermodynamically. Thus, chiral beta-Lewis base substituted cyclopropylnitriles may serve as precursors to chiral, configurationally stable organolithium reagents.     

H-bonded complexes between acetylacetone and two molecules of methanol: HF and DFT level study

Five stable H-bonded complexes (supersystems) between acetylacetone and two methanol molecules were investigated at the B3LYP and HF levels of theory using the 6-311G** and 6-11++G** basis sets. The most stable complex was found as the one with the highest relative bonding and interaction energies. All vibrational frequencies resulting from calculations with the 6-311++G** basis set were compared with the recorded IR spectrum of acetylacetone/methanol mixture in a molar ratio 1:2.