The (α-1,6) glycosidic bond of isomaltose: a tricky system for theoretical conformational studies

Stable conformations of β-isomaltose (α-d-glucopyranosyl-(1→6)-β-d-glucose) in gas-phase and aqueous solution are investigated in this study using quantum mechanical calculations. Conformational maps are calculated at HF/6-31G(d,p) level and lower energy structures are sampled in the most stable regions. Entropic and thermal corrections are considered and the Boltzmann population is obtained for conformers that are representative of the 18 most stable regions found on the potential energy surface. B3LYP/6-31+G(d,p) and B3LYP/6-311+G(2d,2p) calculations are used in conformational samplings. Solvation effects are considered through the polarizable continuum model approach. Hydroxymethyl group orientations are investigated for the most stable conformers. The influence of electronic correlation and solvation on the glycosidic linkage preference (TG, GT, and GG) and hydroxymethyl group orientation (tg, gt, and gg) are discussed. Heteronuclear spin coupling constants (³J_C,H) along the glycosidic linkage are calculated and comparison with other theoretical results and experiments is used to validate the obtained structures.     

Tautomeric energetics of xanthine oxidase substrates: xanthine, 2-oxo-6-methylpurine, and lumazine

The stability of the tautomers of each of the three important substrates of xanthine oxidase, xanthine, 2-oxo-6-methylpurine, and lumazine, was examined by quantum mechanical calculations. The geometries of these tautomers were optimized at the AM1, Hartree-Fock (HF/6-31G), and hybrid Hartree-Fock/density functional theory (B3LYP/6-31G(d)) levels of theory. The single point energies of some of the more stable tautomers for each of the substrates were calculated at the B3LYP/6-311 +G(2d,p) level of theory. The Conductor Polarized Continuum Model (CPCM) was used to evaluate the solvent effects on the relative stabilities of these tautomers. The calculations clearly identify the lowest energy tautomeric form for xanthine and lumazine. On the other hand, there appear to be three tautomers for 2-oxo-6-methylpurine, with only minor energetic differences in vacuo. In water, however, only one of them predominates. The lowest energy tautomers presumably represent the predominant tautomeric forms at the molybdenum center of xanthine oxidase during catalysis. Implications of these computational results are discussed in the context of enzyme catalysis.     

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.     

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.     

Molecular conformational analyses of dehydroalanine analogues

The conformational preferences of dehydroalanine (ΔAla) were examined through ab initio calculations. The geometries of the minimum energy conformations for N-formyldehydro alanilamide and N-acetyl-N′-methylamide of dehydroalanine were determined by gradient optimization at the HF/6-31G* level, and correlation corrections were examined with MP2 single-point energy calculations. Furthermore, HF/3-21G ab initio geometry optimizations were performed on nine conformations of the model tripeptide N-acetyl-N′-methylamide of didehydroalanine. The results indicate that the C₅ is the lowest energy conformation at all levels of theory. However, the relative energy of the helix conformation decreases when the number of ΔAla residues in the peptide chain increases. On the other hand, significant variations of the geometry upon conformational change were observed for the three compounds investigated. These results permit to extract important conformationally dependent geometry trends. The results of this study were compared to x-ray diffraction data on single crystals of dehydroalanine-containing peptides. © 1995 John Wiley & Sons, Inc.     

Comparative theoretical studies of energetic pyrazole-pyridine derivatives

The pyrazole-pyridine derivatives were optimized to obtain their molecular geometries and electronic structures at the DFT-B3LYP/6-31G(d,p) and DFT-B3P86/6-31G(d,p) levels. Molecular mechanics (MM) calculations were performed for the title compounds. Heats of formation (HOFs) were predicted through designed isodesmic reactions. Detonation performance was evaluated by using the Kamlet-Jacobs equations based on the calculated densities and heats of formation. The thermal stability of the title compounds was investigated via the bond dissociation energies (BDEs). The simulation results reveal that the compound with one pyrazole ring that is fully nitro-substituted performs similarly to the famous explosive HMX, and the compound with two pyrazole rings that are fully nitro-substituted outperforms HMX. According to the quantitative standard of energetics and stability as high energy density materials (HEDMs), the compound with two pyrazole rings that are fully nitro-substituted essentially satisfies this requirement.     

Zwitterionic conformers of pyrrolysine and their interactions with metal ions—a theoretical study

A total of 16 pyrrolysine conformers in their zwitterionic forms are studied in gas and simulated aqueous phase using a polarizable continuum model (PCM). These conformers are selected on the basis of our study on the intrinsic conformational properties of non-ionic pyrrolysine molecule in gas phase [Das and Mandal (2013) J Mol Model 19:1695?1704]. In aqueous phase, the stable zwitterionic pyrrolysine conformers are characterized by full geometry optimization and vibrational frequency calculations using B3LYP/6-311++G(d,p) level of theory. Single point calculations are also carried out at MP2/6-311++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. The calculated relative energy range of the conformers at B3LYP/6-311++G(d,p) level is 5.19 kcal mol^?1 whereas the same obtained by single point calculations at MP2/6-311++G(d,p) level is 4.58 kcal mol^?1. A thorough analysis reveals that four types of intramolecular H-bonds are present in the conformers; all of which play key roles in determining the energetics and in imparting the observed conformations to the conformers. The vibrational frequencies are found to shift invariably toward the lower side of frequency scale corresponding to the presence of the H-bonds. This study also points out that conformers with diverse structural motifs may differ in their thermodynamical stability by a narrow range of relative energy. The effects of metal coordination on the relative stability order and structural features of the conformers are examined by complexing five zwitterionic conformers of pyrrolysine with Cu⁺² through their carboxylate groups. The interaction enthalpies and Gibbs energies, rotational constants, vibrational frequencies and dipole moments of the metal complexes calculated at B3LYP level are also reported. The zwitterionic conformers of pyrrolysine are not stable in gas phase; after geometry optimization they are converted to the non-ionic forms.     

Study of the conformational profile of the cyclohexane analogs of L-phenylalanine.

The conformational profile of the conformationally constrained cyclohexane analogs of phenylalanine (1-amino-2-phenylcyclohexanecarboxylic acids, c6Phe) was assessed using computational methods. For this purpose, the conformational space of the N-acetyl methylamide derivatives of the stereoisomers (2S,3R)c6Phe and (2S,3S)c6Phe was explored by computing their respective Ramachandran maps, and low-energy minima were characterized at molecular mechanics level by means of the AMBER program, using the parm94 force field set of parameters. In order to assess the performance of the molecular mechanics calculations, each of the low-energy conformations was also investigated further at the ab initio level. Accordingly, the molecular mechanics geometries were used as starting conformations to perform full geometry optimizations at the Hartree-Fock level, using a 6-31G(d) basis set. Analysis of the results revealed that the cyclohexane structure directly induces some restrictions on the backbone, and constrains the orientation of the aromatic side-chain to two narrow regions for each stereoisomer. The conformational profile of these amino acids is then explained on the grounds of the interaction between the rigidly held phenyl ring and the main chain NH and CO groups. The results obtained are in good accordance with the experimental observations.     

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.     

Predictive modelling of topology and loop variations in dimeric DNA quadruplex structures

We have used a combination of simulated annealing (SA), molecular dynamics (MD) and locally enhanced sampling (LES) methods in order to predict the favourable topologies and loop conformations of dimeric DNA quadruplexes with T2 or T3 loops. This follows on from our previous MD simulation studies on the influence of loop lengths on the topology of intramolecular quadruplex structures [P. Hazel et al. (2004) J. Am. Chem. Soc., 126, 16405–16415], which provided results consistent with biophysical data. The recent crystal structures of d(G4T3G4)2 and d(G4BrUT2G4) (P. Hazel et al. (2006) J. Am. Chem. Soc., in press) and the NMR-determined topology of d(TG4T2G4T)2 [A.T. Phan et al. (2004) J. Mol. Biol., 338, 93–102] have been used in the present study for comparison with simulation results. These together with MM-PBSA free-energy calculations indicate that lateral T3 loops are favoured over diagonal loops, in accordance with the experimental structures; however, distinct loop conformations have been predicted to be favoured compared to those found experimentally. Several lateral and diagonal loop conformations have been found to be similar in energy. The simulations suggest an explanation for the distinct patterns of observed dimer topology for sequences with T3 and T2 loops, which depend on the loop lengths, rather than only on G-quartet stability.     

Theoretical study of the NMR chemical shift of Xe in supercritical condition

In this work we investigate the level of theory necessary for reproducing the non-linear variation of the ¹²⁹Xe nuclear magnetic resonance (NMR) chemical shift with the density of Xe in supercritical conditions. In detail we study how the ¹²⁹Xe chemical shift depends under supercritical conditions on electron correlation, relativistic and many-body effects. The latter are included using a sequential-QM/MM methodology, in which a classical MD simulation is performed first and the chemical shift is then obtained as an average of quantum calculations of 250 MD snapshots conformations carried out for Xe _n clusters (n =?2 ? 8 depending on the density). The analysis of the relativistic effects is made at the level of 4-component Hartree-Fock calculations (4c-HF) and electron correlation effects are considered using second order Møller-Plesset perturbation theory (MP2). To simplify the calculations of the relativistic and electron correlation effects we adopted an additive scheme, where the calculations on the Xe _n clusters are carried out at the non-relativistic Hartree-Fock (HF) level, while electron correlation and relativistic corrections are added for all the pairs of Xe atoms in the clusters. Using this approach we obtain very good agreement with the experimental data, showing that the chemical shift of ¹²⁹Xe in supercritical conditions is very well described by cluster calculations at the HF level, with small contributions from relativistic and electron correlation effects.     

Study of the conformational profile of the norbornane analogues of phenylalanine.

The conformational profile of the eight stereoisomeric 2-amino-3-phenylnorbornane-2-carboxylic acids (2-amino-3-phenylbicyclo[2.2.1]heptane-2-carboxylic acids) has been assessed by computational methods. These molecules constitute a series of four enantiomeric pairs that can be considered as rigid analogues of either L- or D-phenylalanine. The conformational space of their N-acetyl methylamide derivatives has been explored within the molecular mechanics framework, using the parm94 set of parameters of the AMBER force field. Local minimum energy conformations have been further investigated at the ab initio level by means of the Hartree-Fock and second order Moller-Plesset perturbation energy calculations using a 6-31G(d) basis set. The results of the present work suggest that the bulky norbornane structure induces two kinds of conformational constraints on the residues. On one hand, those of a steric nature directly imposed by the bicycle on the peptide backbone and, on the other hand, those that limit the orientations attainable by the phenyl ring which, in turn, reduces further the flexibility of the peptide backbone. A comparative analysis of the conformational profile of the phenylnorbornane amino acids with that of the norbornane amino acids devoid of the beta-phenyl substituent suggests that the norbornane system hampers the residue to adopt extended conformations in favour of C7-like structures. However, the bicycle itself does not impart a clear preference for any of the two possible C7 minima. It is the aromatic side chain, which is forced to adopt an almost eclipsed orientation, that breaks this symmetry introducing a marked preference for a single region of the (phi, psi) conformational space in each of the phenylalanine norbornane analogues investigated.     

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.

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.     

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.     

Conformational behaviour of the architectural units of peptides and proteins. Assessment of current understanding by ab initio quantum mechanical methods and refinement of the dipeptide energy surface.

The conformational energy surfaces of analogues of the dipeptide unit of polypeptides and proteins are calculated by ab initio methods using extended basis sets.The calculations are not particularly sensitive to the choice of (extended) basis set.The calculations are shown to support a particular empirical method parameterized with respect to crystal data. Non-hydrogen bonded conformations agree to within 3 kcal mol^?1, even for conformations in which quite considerable degrees of atomic overlap occur.Hydrogen bonded conformations, are, however, in less satisfactory agreement and it is the ab initio calculations which appear to be at fault.A simple correction is applied to the ab initio energy for hydrogen bonded conformations, and with the use of the empirical energy surface a full quantum mechanical conformational energy map is interpolated for the alanyl dipeptide.The effect of flexibility in the peptide backbone is taken into account, and supports recent empirical findings that distortions in valence angles must be considered in calculations of the conformational behaviour of peptides.     

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.     

Ab initio computational study of beta-cellobiose conformers using B3LYP/6-311++G**

The molecular structure of 27 conformers of beta-cellobiose were studied in vacuo through gradient geometry optimization using B3LYP density functionals and the 6-311++G** basis set. The conformationally dependent geometry changes and energies were explored as well as the hydrogen-bonding network. The lowest electronic energy structures found were not those suggested from available crystallographic and NMR solution data, where the glycosidic dihedral angles fall in the region (phi, psi) approximately (40 degrees, -20 degrees ). Rather, 'flipped' conformations in which the dihedral angles are in the range (phi, psi) approximately (180 degrees, 0 degrees ) are energetically more stable by approximately 2.5 kcal/mol over the 'experimentally accepted' structure. Further, when the vibrational free energy, deltaG, obtained from the calculated frequencies, is compared throughout the series, structures with (phi, psi) in the experimentally observed range still have higher free energy ( approximately 2.0 kcal/mol) than 'flipped' forms. The range of bridging dihedral angles of the 'normal' conformers, resulting from the variance in the phi dihedral is larger than that found in the 'flipped' forms. Due to this large flat energy surface for the normal conformations, we surmise that the summation of populations of these conformations will favor the 'normal' conformations, although evidence suggests that polar solvent effects may play the dominant role in providing stability for the 'normal' forms. Even though some empirical studies previously found the 'flipped' conformations to be lowest in energy, these studies have been generally discredited because they were in disagreement with experimental results. Most of the DFT/ab initio conformations reported here have not been reported previously in the ab initio literature, in part because the use of less rigorous theoretical methods, i.e. smaller basis sets, have given results in general agreement with experimental data, that is, they energetically favored the 'normal' forms. These are the first DFT/ab initio calculations at this level of theory, apparently because of the length and difficulty of carrying out optimizations at these high levels.     

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.     

Vibrational assignments and derived potential energy distributions for tri- and difluoromethyl ketene by density functional calculations

The structures of 3,3,3-trifluoromethyl ketene and 3,3-difluoromethyl ketene were studied by utilizing ab initio calculations with the 6-311++G** basis set at the (B3LYP) Density Functional level. Full optimization was performed for both molecules in their ground and transition states. Energy optimization of the systems under investigation shows that trifluoromethyl ketene exists only in the cis conformation (fluorine atom eclipses the ketene group). Difluoromethyl ketene was predicted to have two stable conformations: the cis (hydrogen atom eclipses the ketene group) and the gauche (fluorine atom eclipses the ketene group) form. The conformational stability of the molecules was found to be governed mainly by electrostatic and molecular orbital interactions. The vibrational frequencies were computed and complete assignments were provided on the basis of normal coordinate calculations and comparison with similar molecules. The potential energy distributions (PED) among symmetry coordinates were derived for the stable conformations of the two molecules.