Linear free energy relationship for the anomeric effect: MP2, DFT and ab initio study of 2-substituted-1,4-dioxanes

The anomeric effect of 2-substituted 1,4-dioxane derivatives was calculated and compared with the values for substituted cyclohexane. The bond lengths, bond angles, torsion angles, and relative energies of axial and equatorial conformers of 2-substituted 1,4-dioxanes were calculated by the second-order Møller–Plesset (MP2), density functional theory (DFT/B3LYP), and Hartree–Fock (HF) methods using 6-31G^∗ basis set. The energy differences between the axial and equatorial conformers, endo and exo-anomeric effects, repulsive non-bond and H-bonding interactions were investigated. A linear free energy relationship (LFER) between calculated (MP2/6-31G^∗) anomeric effect and inductive substituent constants (σ_I) was obtained for 2-substituted-1,4-dioxanes (slope = 6.19 and r² = 0.967). The calculated energy differences indicate lower equatorial orientation for 2-substituted-1,4-dioxanes compared to the 2-substituted-tetrahydropyrans. The contribution of resonance, hyperconjugation, inductive, steric, hydrogen bonding, electrostatic interaction, and level of theory influences the anomeric effect.     

Complete basis set,hybrid-DFT study and NBO interpretation of conformational analysis of 2-methoxytetrahydropyran and its thiopyran and selenopyran analogues in relation to the anomeric effect

2-Methoxytetrahydropyran (1), -thiopyran (2) and -selenopyran (3) have been chosen as model compounds to investigate the origin of the anomeric effect (AE). The impacts of the hyperconjugation, electrostatic and steric interactions on the conformational preferences of compounds 1–3 have been analysed by means of complete basis set-4, hybrid-density functional theory (B3LYP/6-311+G**) based methods and natural bond orbital (NBO) interpretation. Both levels of theory showed that the axial conformations of compounds 1–3 are more stable than their equatorial conformations. The Gibbs free energy difference (G _eq–G _ax) values (i.e. ΔG _eq–ax) between the axial and equatorial conformations increase from compound 1 to compound 2 but decrease from compound 2 to compound 3. Based on the NBO results obtained, the AE associated with the electron delocalisation [i.e. Σ(endo-AE_eq + exo-AE_eq) ? Σ(endo-AE_ax + exo-AE_ax)] increase slightly from compound 1 to compound 2 but decrease from compound 2 to compound 3. Similar trend is also observed for the differences between the calculated total steric exchange energy values [i.e. Δ(TSEE)_eq–ax]. On the other hand, the calculated differences between the dipole moment values of the axial and equatorial conformations [i.e. Δ(μ_eq–μ_ax)] decrease from compound 1 to compound 3. These findings led to the proposal that the AE associated with the electron delocalisation (the hyperconjugation effect) is more significant for the explanation of the conformational preferences of compounds 1–3 than the electrostatic model. The correlations between the AE associated with the electron delocalisation, bond orders, TSEE, ΔG _eq–ax, dipole–dipole interactions, structural parameters and conformational behaviours of compounds 1–3 have been investigated.     

Effect of microsolvation on zwitterionic glycine: an ab initio and density functional theory study

The effect of microsolvation on zwitterionic glycine, considering both (-NH3(+)) as proton donor and (-COO(-)) as proton acceptor at correlated ab initio (MP2) level and density functional methods (B3LYP, PW91, MPW1PW91 and PBE) using 6-311++G** basis set has been reported. DFT methods have been employed so as to compare the performance/quality of different gradient-corrected correlation functionals (PW91, PBE), hybrid functionals (B3LYP, MPW1PW91) and to predict the near quantitative structural and vibrational properties, at reduced computational cost. B3LYP method outperforms among the different DFT methods for the computed hydrogen bond distances and found closer to the value obtained by correlated MP2 level, whereas MPW1PW91 and PBE methods shows very similar values but approximately 0.03 A less, compared to B3LYP method. MP2 calculation and single point CCSD(T)//MP2 calculation have been considered to decompose the interaction energy, including corrections for basis set superposition error (BSSE). Moreover, charge distribution analysis has also been carried out to understand the long raised questions, how and why the two body energies have significant contribution to the total binding energy.     

Valence isomerization of 2-phospha-4-silabicyclo[1.1.0]butane: a high-level ab initio study

The rearrangements for 2-phospha-4-silabicyclo[1.1.0]butane, analogous to the valence isomerization of the hydrocarbons bicyclobutane, 1,3-butadiene, and cyclobutene, were studied at the (U)QCISD(T)/6-311+G**//(U)QCISD/6-31G* level of theory. The monocyclic 1,2-dihydro-1,2-phosphasiletes are shown to be the thermodynamically preferred product, in contrast to the isomerization of the hydrocarbons, which favors the 1,3-butadiene structure. Furthermore, an unprecedented direct isomerization pathway to the 1,2-dihydro-1,2-phosphasiletes was identified. This pathway is competitive with the isomerization via the open-chain butadienes and becomes favorable when electron-donating substituents are present on silicon. Figure 2-Phospha-4-silabicyclo[1.1.0]butane can isomerize directly into the more stable P,Si-cyclobutene via an unprecedented [sigma2s+sigma2a] process, which becomes favorable over the isomerization via the P,Si-butadiene when electron-donating substituents are present on silicon.     

Structural and electronic properties of diazonium functionalized (4, 4) single walled carbon nanotube: an ab initio study

In this work, ab initio density functional theory (DFT) calculations are performed to study the structural and electronic properties of diazonium reagent functionalized (4, 4) single-walled carbon nanotube (SWCNT). We find the aryl group covalently bonds with SWCNT and prefers to be perpendicular to the side wall of nanotube. It has a rotational barrier of 0.35 eV around the formed aryl-tube bond axis and should be thermodynamically stable at room temperature. Additionally, new peaks appeared around the Fermi energy in the density of state (DOS) due to the weak band dispersion. Increasing of the coverage of the functional group will result in significant upshift of the Fermi level.     

Metal ions in non-complementary DNA base pairs: an ab initio study of Cu(I), Ag(I), and Au(I) complexes with the cytosine-adenine base pair

 Ab initio calculations have been carried out to characterize the structure and energetics of a silver(I) complex with the cytosine-adenine DNA base pair and an aqua ligand in the coordination sphere of Ag. In addition, we have also studied analogous complexes with Cu(I) and Au(I), and structures in which adenine has been replaced by purine in order to investigate the structural role of the adenine amino group. The calculations revealed that all metal-modified structures are dominated by the metal-base interactions, while the water-metal ion interaction and many-body interligand repulsion are less important contributions. Nevertheless, the structural role of the water molecule in the complex is quite apparent and in agreement with an earlier crystallographic study. The metal-modified base pairs exhibit large conformational flexibility toward out-of-plane motions (propeller twist and buckle), comparable or, in some cases, even larger than that observed in the base pairs without metal ions. All structures have been optimized within the Hartree-Fock approximation, while interaction energies were evaluated with the inclusion of electron correlation. Received: 25 March 1999 / Accepted: 10 June 1999     

Direct ab initio dynamics studies of the hydrogen abstraction reactions of hydrogen atom with n-propyl radical and isopropyl radical

The kinetics of the hydrogen abstraction reactions of hydrogen atom with n-propyl radical and isopropyl radical were studied using the direct ab initio dynamics approach. BHandHLYP/cc-pVDZ method was employed to optimize the geometries of stationary points as well as the points on the minimum energy path (MEP). The energies of all the points for the two reactions were further refined at the QCISD(T)/cc-pVTZ level of theory. No barrier was found at the QCISD(T)/cc-pVTZ//BHandHLYP/cc-pVDZ level of theory for both reactions. The forward and reverse rate constants were evaluated with both canonical variational transition state theory (CVT) and microcanonical variational transition state theory ( VT) in the temperature range of 300–2,500 K. The fitted three-parameter Arrhenius expression of the calculated CVT rate constants at the QCISD(T)/cc-pVTZ//BHandHLYP/cc-pVDZ level of theory are k^CVT (n – C₃H₇)=1.68×10^–14 T^0.84 e^(319.5/T) cm³ molecule^–1 s^–1 and k^CVT (iso-C₃H₇)=4.99×10^–14 T^0.90 e^(159.5/T) cm³ molecule^–1 s^–1 for reactions of n-C₃H₇ + H and iso-C₃H₇ + H, respectively, which are in good agreement with available literature data. The variational effects were analysed.Figure Comparison of the calculated forward rate constants at the QCISD(T)/cc-pVTZ//BHandHLYP/cc-pVDZ level of theory and the available experimental and theoretical data of the reaction vs 1,000/T for the two reactions.     

The anomeric effect revisited. A possible role of the CH/n hydrogen bond

Ab initio MO calculations were carried out at the MP2/6-311++G(d,p) level to investigate the conformational energy of 2-substituted oxanes and 1,3-dioxanes. It has been found that the Gibbs free energies of the axial conformers are smaller than those of the corresponding equatorial conformers in every case when the 2-substituent Z is electron withdrawing (OCH(3), F, Cl, Br). The difference in Gibbs energy between the equatorial and axial conformers DeltaG(eq-ax) increases from Z=OCH(3) to F, Cl, and then to Br. In the axial conformers, the interatomic distance between Z and the axial C-H, separated by four covalent bonds, has been found to be appreciably shorter than the van der Waals distance, suggesting the importance of the five-membered CH/n (CH/O or CH/halogen) hydrogen bond in stabilizing these conformations. Natural bonding orbital (NBO) charges of the relevant atoms have been shown to be different between the two conformers: more positive for H and more negative for C in the axial conformers than in the corresponding equatorial conformers. In view of the above findings, we suggest that the CH/n hydrogen bond plays an important role in stabilizing the axial conformation in 2-substituted oxanes and 1,3-dioxanes, and by implication, in the anomeric effect in carbohydrate chemistry.     

Ab initio study of phenyl benzoate: structure,conformational analysis,dipole moment,IR and Raman vibrational spectra

The molecular structure (bond distances and angles), conformational properties, dipole moment and vibrational spectroscopic data (vibrational frequencies, IR and Raman intensities) of phenyl benzoate were calculated using Hartree–Fock (HF), density functional (DFT), and second order Møller–Plesset perturbation theory (MP2) with basis sets ranging from 6-31G* to 6-311++G**. The theoretical results are discussed mainly in terms of comparisons with available experimental data. For geometric data, good agreement between theory and experiment is obtained for the MP2, B3LYP and B3PW91 levels with basis sets including diffuse functions. The B3LYP/6-31+G* theory level estimates the shape of the experimental functions for phenyl torsion around the Ph–O and Ph–C bonds well, but reproduces the height of the rotational barriers poorly. The B3LYP/6-31+G* harmonic force constants were scaled by applying the scaled quantum mechanical force field (SQM) technique. The calculated vibrational spectra were interpreted and band assignments were reported. They are in excellent agreement with experimental IR and Raman spectra.Figure Calculated and experimental (GED) potential energy functions for torsional motion of phenyl benzoate relative to the minimum value. a The potential function for torsion about the O₃–C₄ bond. b The potential function for torsion about the C₂–C₁₀ bond.     

Synthesis, crystal structures and ab initio studies of some heptaaza manganese(II) macrocyclic Schiff-base complexes with two 2-aminoethyl pendant arms

Two new heptaaza Mn(II) macrocyclic bis(pendant donor) Schiff-base complexes, [MnL⁶]²⁺ and [MnL⁷]²⁺, have been prepared via Mn(II) templated [1 + 1] cyclocondensation of 2,6-diacetylpyridine with an appropriate branched hexaamine, and characterized by X-ray crystallography and spectroscopic methods. The ligands are 16- and 17-membered pentaaza macrocycles with two 2-aminoethyl pendant arms [L⁶ = 2,14-dimethyl-6,10-bis(2-aminoethyl)-3,6,10,13,19-pentaazabicyclo[13.3.1]nonadeca-1(19),2,13,15,17-pentaene and L⁷ = 2,15-dimethyl-6,11-bis(2-aminoethyl)-3,6,11,14,20-pentaazabicyclo[14.3.1]eicosa-1(20),2,14,16,18-pentaene]. The crystal structures of [MnL⁷]²⁺ and another related complex, [MnL⁵]²⁺ [L⁵ = 2,13-dimethyl-6,9-bis(2-aminoethyl)-3,6,9,12,18-pentaazabicyclo[12.3.1]octadeca-1(18),2,12,14,16-pentaene] were determined from X-ray diffraction data and the structure of each complex indicates that the geometry of the coordination sphere is that of a slightly distorted pentagonal bipyramid with the Mn(II) ion located within a pentaaza macrocycle and two pendant amines coordinating on opposite sides. The increased macrocycle ring size in [MnL⁷]²⁺ introduces an unusual seven-membered chelate ring and leads to long Mn-N bond lengths within the macrocycle, but to shorter pendant amine-Mn(II) bond lengths. Ab initio HF-MO calculations produce structures for the two complexes which are in close agreement with experimental findings and predict a similar structure for [MnL⁶]²⁺, but with the longest pendant amine-Mn(II) bond lengths of the three structures.     

Introducing a novel method based on the imperialistic competitive algorithm to determine fluorine intermolecular potential from ab initio calculations and calculation of some properties via MD simulations

In this work, the possibility of obtaining an accurate site-site potential model suitable for use in molecular dynamics (MD) simulations of fluorine from ab initio calculations has been explored. The exploration was made on ab initio calculations. To reduce the ab initio pair potentials into a site-site potential, a higher significance was assigned to the configuration which is more stable. For this purpose, the imperialistic competitive algorithm (ICA) was implemented as a powerful optimisation tool. The calculated second virial coefficients were compared to the experimental values to test the quality of the presented intermolecular potential. The relative error for the calculated second virial coefficient ranged from 0.1 to 5.6%. MD simulations were used to evaluate the ability of the proposed intermolecular potential function. The relative error for the MD simulations ranged from 0.5 to 5.2%. The results are in good agreement with experimental data.     

New macrocyclic Schiff base complexes incorporating a phenanthroline unit: Part 1; Template synthesis of three cadmium(II) complexes and crystal structure, NMR and ab initio studies

Complexes of three Cd(II)-containing macrocyclic Schiff base complexes containing a phenanthroline ligand (L) of the type [CdLⁿ(Cl)]⁺ (n=2,3,4), have been prepared via [1+1] cyclocondensation of 2,9-dicarboxaldehyde-1,10-phenanthroline and a number of linear triamines via a metal-templated reaction and coordination features have been examined. The ligands, L, are 16-, 17-, and 18-membered pentaaza macrocycles and all the complexes incorporate a 1,10-phenanthroline unit as an integral part of their cyclic structure. The complexes have been characterized by a variety of methods including IR, ¹H, ¹³C, DEPT, COSY(H,H) and HMQC(H,C) NMR studies and MALDI mass spectrometry. The polymeric structure of was determined by X-ray crystallography, which showed that the complex cation consisted of a pentagonal bipyramidally coordinated Cd(II) ion. The seven-coordinated Cd(II) ion is ligated by the five nitrogen atoms of the macrocycle in the equatorial plane and has two bridging chloride ligands in the axial positions resulting in a ribbon of such complex ions. Supporting ab initio HF-MO calculations have been undertaken using the standard 3-21G^∗ and 6-31G^∗ basis sets.     

Quantum simulation for muoniated and deuterated methyl radicals in implicit water solvent: combined ab initio path integral molecular dynamics and the polarizable continuum model simulation study

We have combined ab initio path integral molecular dynamics (PIMD) simulation and the polarizable continuum model (PCM) method to efficiently incorporate solvent effects into nuclear quantum fluctuation of molecular systems. Our combined ab initio PIMD–PCM simulation was applied to muoniated and deuterated methyl radical immersed in implicit water solvent to gain information on solvent and isotope effects from one simulation run. We found that solvent effects lead to the bond elongation and a decrease in the magnitude of isotropic hyperfine coupling constants. These are consistent with the trends in conventional static calculations and experiments. In addition, the performance of cavity models (universal force field, united atom specified for Kohn–Sham and these hybrid models) and the conservation of the PIMD–PCM Hamiltonian were accessed. We confirmed that solvent effects on nuclear quantum fluctuation are efficiently computed using our combined simulation of quantum solute in implicit solvent.     

Solvent effects on the conformational preferences of model peptoids. MP2 study

The influence of aqueous environment on the main‐chain conformation (ω₀, ?, and ψ dihedral angles) of two model peptoids: N‐acetyl‐N‐methylglycine N’‐methylamide (Ac‐N(Me)‐Gly‐NHMe) ( 1 ) and N‐acetyl‐N‐methylglycine N’,N’‐dimethylamide (Ac‐N(Me)‐Gly‐NMe₂) ( 2 ) was investigated by MP2/6‐311++G(d,p) method. The Ramachandran maps of both studied molecules with cis and trans configuration of the N‐terminal amide bond in the gas phase and in water environment were obtained and all energy minima localized. The polarizable continuum model was applied to estimate the solvation effect on conformation. Energy minima of the Ac‐N(Me)‐Gly‐NHMe and Ac‐N(Me)‐Gly‐NMe₂ have been analyzed in terms of the possible hydrogen bonds and C = O dipole attraction. To validate the theoretical results obtained, conformations of the similar structures gathered in the Cambridge Crystallographic Data Centre were analyzed. Obtained results indicate that aqueous environment in model peptoids 1 and 2 favors the conformation F (? and ψ = ?70º, 180º), and additionally significantly increases the percentage of structures with cis configuration of N‐terminal amide bond in studied compounds. Copyright © 2014 European Peptide Society and John Wiley & Sons, Ltd.     

A DFT/ab initio study of hydrogen bonding and conformational preference in model cellobiose analogs using B3LYP/6-311++G**

A series of beta-cellobiose analogs were studied at the B3LYP/6-311++G** level of theory to isolate and understand how the various electronic components of the beta-(1-->4)-linked disaccharide, cellobiose, contribute to the energetic stability of the molecule in vacuo. Previous studies on beta-cellobiose (see accompanying paper) showed that the most energetically stable conformation was that in which the dihedral angle phi (phi(H)) was 'flipped' by approximately 180 degrees relative to the 'normal' form. From our examination of eight sets of structures in which different combinations of functional hydroxyl and hydroxymethyl groups were removed, it was determined that only beta-cellobiose and one other analog (analog 7, beta-xylobioside), an analog in which both hydroxymethyl groups were removed but the exocyclic hydroxyl groups retained, can form a 'cooperative' hydrogen-bonding network. Only in these two molecules did we find continuous synergistic 'communication' through hydrogen bonding from one sugar moiety to the other. This 'cooperative' hydrogen bonding energetically stabilizes the 'flipped' conformation of beta-cellobiose and beta-xylobioside, while the other analogs studied were unable to form a 'cooperative' grouping of hydrogen bonds and thus were more stable in their 'normal' conformational state.     

Ab initio computational insight into the ion-pair SN2 reaction of lithium isothiocyanate and methyl fluoride in the gas phase and in acetone solution

The ion-pair S_N2 reaction LiNCS + CH₃F with two mechanisms, inversion and retention, was investigated at the MP2(full)/6-311+G**//HF/6-311+G** level in the gas phase and in acetone solution. All HF-optimized structures were confirmed by vibrational frequency analysis. Based on IRC analyses, eight possible reaction pathways in the title reaction are proposed. The inversion mechanism through a six-membered-ring transition-state structure is the most favorable. Methyl thiocyanate should form preferentially in the gas phase and the more stable methyl isothiocyanate will be the main product in CH₃COCH₃. The retardation of the reaction in CH₃COCH₃ solution was attributed to the differences in the solvation free energies in the separated reactants and transition structures. All of the theoretical results are consistent with the experiment.     

Chiral bias of amyloid fibrils revealed by the twisted conformation of Thioflavin T: an induced circular dichroism/DFT study

Since it was implicated in a number of neurodegenerative conditions, such as Alzheimer disease, formation of beta-sheet-rich protein fibrils (amyloids) has been drawing a lot of attention. One of elusive aspects of amyloidogenesis concerns the mechanisms of specific binding of molecules such as Congo red, or Thioflavin T by amyloid fibrils. A comprehensive understanding of these docking interactions is needed, however, for the sake of furthering biochemical studies and developing molecular, pharmacological strategies preventing proliferation of amyloids in vivo. Through the application of circular dichroism, here we show that upon binding to insulin fibrils, a twisted conformation is enforced in molecules of Thioflavin T, manifested in a strong negative Cotton effect around 450 nm, which is supported by density functional theory-based calculations. This finding may lead to circular dichroism of Thioflavin T becoming a new diagnostic technique for protein fibrils, complementary to fluorescence spectroscopy.     

DFT study of a substrate and inhibitors of 1-deoxy-2-xylulose-5-phosphate reductoisomerase – the potential novel target molecule for anti-malaria drug development

1-Deoxy-2-xylulose-5-phosphate (DOXP) reductoisomerase is a novel target for developing anti-malaria drugs. The determination of structural and electronic properties of the inhibitor molecules is of crucial importance for analyzing the interactions between DOXP-reductoisomerase and its inhibitors. Geometry-optimizations and single point calculations at the B3LYP/3-21G*//B3LYP/3-21G** and B3LYP/3-21G*//MP2/3-21G** levels were performed to determine the structures and charge distributions of an enzyme substrate (1-deoxy-D-xylulose 5-phosphate) and the two inhibitors (fosmidomycin and FR-900098). The theoretically derived bond lengths are in excellent agreement with the corresponding experimental values reported for similar structures. Partial charges and dipole moments are assigned using the Mulliken and natural population analyses. The calculated structures and partial charge distributions can readily be used for the further development of biologically active inhibitors of DOXP-reductoisomerase as well as parameters for docking experiments.Electronic Supplementary Material available.     

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.