Hydrogen bonds in galactopyranoside and glucopyranoside: a density functional theory study

Density functional theory calculations on two glycosides, namely, n-octyl-β-D-glucopyranoside (C₈O-β-Glc) and n-octyl-β-D-galactopyranoside (C₈O-β-Gal) were performed for geometry optimization at the B3LYP/6-31G level. Both molecules are stereoisomers (epimers) differing only in the orientation of the hydroxyl group at the C4 position. Thus it is interesting to investigate electronically the effect of the direction (axial/equatorial) of the hydroxyl group at the C4 position. The structure parameters of X-H???Y intramolecular hydrogen bonds were analyzed, while the nature of these bonds and the intramolecular interactions were considered using the atoms in molecules (AIM) approach. Natural bond orbital analysis (NBO) was used to determine bond orders, charge and lone pair electrons on each atom and effective non-bonding interactions. We have also reported electronic energy and dipole moment in gas and solution phases. Further, the electronic properties such as the highest occupied molecular orbital, lowest unoccupied molecular orbital, ionization energy, electron affinity, electronic chemical potential, chemical hardness, softness and electrophilicity index, are also presented here for both C₈O-β-Glc and C₈O-β-Gal. These results show that, while C₈O-β-Glc possess– only one hydrogen bond, C₈O-β-Gal has two intramolecular hydrogen bonds, which further confirms the anomalous stability of the latter in self-assembly phenomena.     

Acetylene chain reaction on hydrogenated boron nitride monolayers: a density functional theory study

Spin-polarized first-principles total-energy calculations have been performed to investigate the possible chain reaction of acetylene molecules mediated by hydrogen abstraction on hydrogenated hexagonal boron nitride monolayers. Calculations have been done within the periodic density functional theory (DFT), employing the PBE exchange correlation potential, with van der Waals corrections (vdW-DF). Reactions at two different sites have been considered: hydrogen vacancies on top of boron and on top of nitrogen atoms. As previously calculated, at the intermediate state of the reaction, when the acetylene molecule is attached to the surface, the adsorption energy is of the order of ?0.82 eV and ?0.20 eV (measured with respect to the energy of the non interacting molecule-substrate system) for adsorption on top of boron and nitrogen atoms, respectively. After the hydrogen abstraction takes place, the system gains additional energy, resulting in adsorption energies of ?1.52 eV and ?1.30 eV, respectively. These results suggest that the chain reaction is energetically favorable. The calculated minimum energy path (MEP) for hydrogen abstraction shows very small energy barriers of the order of 5 meV and 22 meV for the reaction on top of boron and nitrogen atoms, respectively. Finally, the density of states (DOS) evolution study helps to understand the chain reaction mechanism.

Open image in new window

Graphical abstract Acetylene chain reaction on hydrogenated boron nitride monolayers

     

Conformational preferences of 1,4,7-trithiacyclononane: a molecular mechanics and density functional theory study

Conformational preferences of 1,4,7-trithiacyclononane were studied using a highly efficient sampling technique based on local nonstochastic deformations and the MM2(91) force field. The results show that conformers that the molecule adopts in the crystal state were found to be low-energy conformers (LECs) within 5 kcal mol(-1) of the global minimum. A conformation with C1 symmetry was the global minimum and the C3 and C2 conformations were calculated to be 0.03 and 1.78 kcal mol(-1) higher in energy, respectively. The structures were further minimized using Density Functional Theory (DFT) calculations with two different functionals. The C2 and the C1 conformations were found to be LECs with the C3 conformation more than 4.0 kcal mol(-1) above the global minimum. The relative energies and structural ordering obtained using the BP86 functional are in agreement with the previously reported relative energies calculated using second-order Moller-Plesset (MP2) ab initio calculations. With the energy ordering being dependent on the molecular mechanics force field used, the approach of MM-->DFT (searching exhaustively the available conformational space at the MM level followed by generating the energy ordering through DFT calculations) appears to be appropriate for thiacrown ethers.     

The reaction mechanism of phenylethanolamine N-methyltransferase: A density functional theory study

Hybrid density functional theory methods were used to investigate the reaction mechanism of human phenylethanolamine N-methyltransferase (hPNMT). This enzyme catalyzes the S-adenosyl-l-methionine-dependent conversion of norepinephrine to epinephrine, which constitutes the terminal step in the catecholamine biosynthesis. Several models of the active site were constructed based on the X-ray structure. Geometries of the stationary points along the reaction path were optimized and the reaction barrier and energy were calculated and compared to the experimental values. The calculations demonstrate that the reaction takes place via an S_N2 mechanism with methyl transfer being rate-limiting, a suggestion supported by mutagenesis studies. Optimal agreement with experimental data is reached using a model in which both active site glutamates are protonated. Overall, the mechanism of hPNMT is more similar to those of catechol O-methyltransferase and glycine N-methyltransferase than to that of guanidinoacetate N-methyltransferase in which methyl transfer is coupled to proton transfer.     

Structure-activity relationship study of flavone compounds with anti-HIV-1 integrase activity: a density functional theory study

Human immunodeficiency virus type-1 integrase (HIV-1 IN) is an essential enzyme for effective viral replication. Flavone compounds have been very much studied due to their activity during the inhibition process of HIV-1 IN. In this study, we employed density functional theory (DFT) using the B3LYP hybrid functional to calculate a set of molecular properties for 32 flavonoid compounds with anti-HIV-1 IN activity. The stepwise discriminant analysis (SDA), principal component analysis (PCA) and hierarchical cluster analysis (HCA) methods were employed to reduce dimensionality and investigate possible relationship between the calculated properties and the anti-HIV-1 IN activity. These analyses showed that the molecular hydrophobicity (ClogP), charge on atom 11 and electrophilic index (omega) are responsible for the separation between anti-HIV-1 IN active and inactive compounds.     

Adsorption of Mn atom on pristine and defected graphene: a density functional theory study

The functionalization of graphene with transition metals is of great interest due to its wide range of applications, such as hydrogen storage, spintronics, information storage, etc. Due to its magnetic property adsorption of Mn atom on graphene has a high consequence on the electronic properties of graphene. The increase in size of the graphene sheet with hydrogen termination has a high impact on the transformation of electronic properties of the graphene sheet. Hence in this work, we investigate the size as well as change in structural and electronic properties of pristine/defective graphene sheets on adsorption of Mn atom using density functional theory methods. From the results obtained a higher adsorption energy value of 3.04 eV is found for Mn adatom on the defected graphene sheet than the pristine, 1.85 eV. It is subject to the coverage effect which decreases on increasing number of carbon atoms. Moreover, a decrease in energy gap is observed in pristine and defected graphene sheets with a high number of carbon atoms. The density of states illustrates the significant effect for hydrogen termination in the conduction band of the Mn adsorbed graphene sheet with low carbon atoms.

Open image in new window

Graphical Abstract Mn adatom on graphene at different sites

     

Exchange interactions in a Fe5 complex: A theoretical study using density functional theory

Theoretical methods based on density functional theory have been employed to analyze the exchange interactions in an Fe₅ complex. The calculated exchange coupling constants are consistent with an S = 15/2 ground state and agree well with those reported previously for other Fe^III polynuclear complexes. The strongest antiferromagnetic interactions are those through a double oxo-benzotriazole bridging ligand, where the benzotriazole ligand seems to play a minor role and the exchange coupling constants present values typical of single oxo bridging systems.     

Electronic structure study using density functional theory in organic dendrimers

The electronic and structural properties of pyrrolic ring derivatives were studied using density functional theory (DFT) in terms of their application as organic semiconductor materials in photovoltaic devices. The B3LYP hybrid functional in combination with Pople type 6-31G(d) basis set with a polarization function was used in order to determine the optimized geometries and the electronic properties of the ground state, while transition energies and excited state properties were obtained from time-dependent (TD)-DFT with B3LYP/6-31G(d) calculation. The investigation of pyrrolic derivatives formed by the arrangement of several monomeric units revealed that three-dimensional (3D) conjugated architectures in which the combination of a triphenylamine (TPA) core with π-conjugated rings attached to the core, present the best geometric and electronic characteristics for use as an organic semiconductor material. The highest occupied molecular orbital (HOMO) − lowest unoccupied molecular orbital (LUMO) energy gap was decreased in 3D-structures that extend the absorption spectrum toward longer wavelengths, revealing a feasible intramolecular charge transfer process in these systems. All calculations in this work were performed using the Gaussian 03 W software package.     

Structural, electronic and magnetic effects of Al-doped niobium clusters: a density functional theory study

The application of the ab initio stochastic search procedure with Saunders "kick" method has been carried out for the elucidation of global minimum structures of a series of Al-doped clusters, Nb(n)Al (1 ≤ n ≤ 10). We have studied the structural characters, growth behaviors, electronic and magnetic properties of Nb(n)Al by the density functional theory calculations. Unlike the previous literature reported on Al-doped systems where ground state structures undergo a structural transition from the Al-capped frame to Al-encapsulated structure, we found that Al atom always occupies the surface of Nb(n)Al clusters and structural transition does not take place until n = 10. Note that the fragmentation proceeds preferably by the ejection of an aluminum atom other than niobium atom. According to the natural population analysis, charges always transfer from aluminum to niobium atoms. Furthermore, the magnetic moments of the Nb(n)Al clusters are mainly located on the 4d orbital of niobium atoms, and aluminum atom possesses very small magnetic moments.     

Unsaturation in binuclear cyclopentadienyliron carbonyl thiocarbonyls: Four-electron donor bridging thiocarbonyl groups versus metal-metal multiple bonds

Theoretical studies on the binuclear cyclopentadienyliron carbonyl thiocarbonyl derivatives Cp₂Fe₂(CO)₂(μ-CS)(μ-CO) and Cp₂Fe₂(CO)₂(μ-CS)₂ indicate that the trans and cis isomers are nearly degenerate in energy, consistent with experiment. Structures with bridging CS groups are of lower energy than corresponding structures with bridging CO groups. The corresponding unbridged Cp₂Fe₂(CS)(CO)₃ and Cp₂Fe₂(CS)₂(CO)₂ isomers are predicted to lie 11 and 16 kcal/mol, respectively, above their global minima, indicating increasing activation energies for the cis/trans interconversion as bridging CO groups are replaced by bridging CS groups. The unsaturated species Cp₂Fe₂(μ-CS)(μ-CO)₂ and Cp₂Fe₂(μ-CS)₂(μ-CO) are predicted to have triply bridged triplet spin state structures with FeFe double bonds of lengths 2.26 Å, analogous to the experimentally known triplet (Me₅C₅)₂Fe₂(μ-CO)₃. However, low-lying singlet Cp₂Fe₂(CS)(CO)₂ and Cp₂Fe₂(CS)₂(CO) structures with four-electron donor bridging η²-μ-CS groups and formal Fe-Fe single bonds are also found. The lowest lying Cp₂Fe₂(CS)(CO) and Cp₂Fe₂(CS)₂ structures have two bridging groups and very short FeFe distances of ∼2.14 Å, suggesting formal triple bonds. Several higher energy four-electron donor η²-μ-CS bridged structures are also found for Cp₂Fe₂(CS)(CO) and Cp₂Fe₂(CS)₂. In addition, singlet and triplet structures are found for Cp₂Fe₂(CS)₂ in which the two CS ligands have coupled to form a bridging SCCS group with a carbon-carbon bond. Only a η²-μ-CS bridged singlet structure is predicted for Cp₂Fe₂(CS), rather than the normal bridged structure with a FeFe quadruple bond such as that predicted for the carbonyl analog Cp₂Fe₂(CO).     

A density functional theory study of the decomposition mechanism of nitroglycerin

The detailed decomposition mechanism of nitroglycerin (NG) in the gas phase was studied by examining reaction pathways using density functional theory (DFT) and canonical variational transition state theory combined with a small-curvature tunneling correction (CVT/SCT). The mechanism of NG autocatalytic decomposition was investigated at the B3LYP/6-31G(d,p) level of theory. Five possible decomposition pathways involving NG were identified and the rate constants for the pathways at temperatures ranging from 200 to 1000 K were calculated using CVT/SCT. There was found to be a lower energy barrier to the β-H abstraction reaction than to the α-H abstraction reaction during the initial step in the autocatalytic decomposition of NG. The decomposition pathways for CHOCOCHONO₂ (a product obtained following the abstraction of three H atoms from NG by NO₂) include O–NO₂ cleavage or isomer production, meaning that the autocatalytic decomposition of NG has two reaction pathways, both of which are exothermic. The rate constants for these two reaction pathways are greater than the rate constants for the three pathways corresponding to unimolecular NG decomposition. The overall process of NG decomposition can be divided into two stages based on the NO₂ concentration, which affects the decomposition products and reactions. In the first stage, the reaction pathway corresponding to O–NO₂ cleavage is the main pathway, but the rates of the two autocatalytic decomposition pathways increase with increasing NO₂ concentration. However, when a threshold NO₂ concentration is reached, the NG decomposition process enters its second stage, with the two pathways for NG autocatalytic decomposition becoming the main and secondary reaction pathways.     

Asparagine deamidation: pH-dependent mechanism from density functional theory

Asparagine deamidation is a decisive event in chemotherapy-induced apoptosis and a major obstacle in the formulation of monoclonal antibodies. Despite the importance of deamidation, little is known about the elementary reactions involved. B3LYP/6-31+G(d,p)/COSMO-RS calculations were used to obtain stable structures and transition states for a network of reactions. Calculated rate constants were incorporated into a kinetic model of the pH dependence and compared to a pseudo-steady-state model. At low pH, the calculations show that deamidation occurs by direct acid-catalyzed hydrolysis to aspartate. At neutral to basic pH, deamidation proceeds by the initial formation of a tetrahedral intermediate. The intermediate can be converted to succinimide by two pathways and three rate-determining steps that shift in relative importance with pH. The calculated pH-dependent rate constant qualitatively agrees with the experimental pH dependence. The rate-determining transition state structures may help to understand chemotherapy-induced apoptosis and improve protein formulations.     

A density functional theory study on oxygen reduction reaction on nitrogen-doped graphene

Nitrogen (N)-doped carbons reportedly exhibit good electrocatalytic activity for the oxygen reduction reaction (ORR) of fuel cells. This work provides theoretical insights into the ORR mechanism of N-doped graphene by using density functional theory calculations. All possible reaction pathways were investigated, and the transition state of each elementary step was identified. The results showed that OOH reduction was easier than O–OH breaking. OOH reduction followed a direct Eley–Rideal mechanism (the OOH species was in gas phase, but H was chemisorbed on the surface) with a significantly low reaction barrier of 0.09 eV. Pathways for both four-electron and two-electron reductions were possible. The rate-determining step of the two-electron pathway was the reduction of O₂ (formation of OOH), whereas that of the four-electron pathway was the reduction of OH into H₂O. After comparing the barriers of the rate-determining steps of the two pathways, we found that the two-electron pathway was more energetically favored than the four-electron pathway.     

Antioxidant activity of phenolic and related compounds: a density functional theory study on the O-H bond dissociation enthalpy

We report here on calculations at the hybrid DFT/HF (B3-LYP/6-31G(d, p)) level of the O-H bond dissociation enthalpy (O-H BDE) of phenylpropenoic acids (caffeic, ferulic, p-coumaric and cinnamic) and phenolic acids and related compounds (gallic, methylgallate, vanillic and gentisic) in order to gain insight into the understanding of structure-antioxidant activity relationships. The results were correlated and discussed mainly on the basis of experimental data in a companion work (Galato D, Giacomelli C, Ckless K, Susin MF, Vale RMR, Spinelli A. Antioxidant capacity of phenolic and related compounds: correlation among electrochemical, visible spectroscopy methods and structure-antioxidant activity. Redox Report 2001; 6: 243-250). The O-H BDE values showed remarkable dependence on the hydroxyl position in the benzene ring and the existence of additional interaction due to hydrogen bonding. For parent molecules, the experimental antioxidant activity (AA) order was properly obeyed only when intramolecular hydrogen bonding was present in the radicalized structures of o-dihydroxyl moieties. In structurally related compounds, excellent correlation with experimental data was in general observed (0.64 < rho < 0.99). However, it is shown that excellent correlation can also be obtained for this series of compounds considering p-radicalized structures which were not stabilized by intramolecular hydrogen bonding, but this had no physical meaning. These findings suggested that the antioxidant activity evaluation of phenolic and related compounds must take into consideration the characteristics of each particular compound.     

Inter- versus intra-molecular cyclization of tripeptides containing tetrahydrofuran amino acids: a density functional theory study on kinetic control

Density functional B3LYP method was used to investigate the preference of intra- and inter-molecular cyclizations of linear tripeptides containing tetrahydrofuran amino acids. Two distinct model pathways were conceived for the cyclization reaction, and all possible transition states and intermediates were located. Analysis of the energetics indicate intermolecular cyclization being favored by both thermodynamic and kinetic control. Geometric and NBO analyses were performed to explain the trends obtained along both the reaction pathways. Conceptual density functional theory-based reactive indices also show that reaction pathways leading to intermolecular cyclization of the tripeptides are relatively more facile compared to intramolecular cyclization.     

Structures and formation mechanisms of aquo/hydroxo oligomeric beryllium in aqueous solution: a density functional theory study

The structures and formation mechanisms of a wide variety of aquo/hydroxo oligomeric beryllium clusters were investigated using density functional theory. The structural parameters of beryllium clusters were found to vary regularly with the stepwise substitution of bound water molecules in the inner coordination sphere by hydroxyl groups. According to the Gibbs free energies deduced from SMD solvation model computations, unhydrolyzed oligomeric beryllium species are the most favorable products of polymerization, independent of the degrees of hydrolysis of the reactants. Simulation of the formation processes of oligomeric beryllium showed that polymerization, in essence, involves the nucleophilic attack of a terminal hydroxyl group in one BeO₄ tetrahedron on the beryllium center in another BeO₄ tetrahedron, leading to the bridging of two BeO₄ tetrahedrons by a hydroxyl group.     

Potassium permeation through the KcsA channel: a density functional study

We present a theoretical study on structural and electronic aspects of K+ permeation through the binding sites of the KcsA channel's selectivity filter. Density functional calculations are carried out on models taken from selected snapshots of a molecular dynamics simulation recently reported [FEBS Lett. 477 (2000) 37]. During the translocation process from one binding site to the other, the coordination number of the permeating K+ ion turns out to decrease and K+ ion polarizes significantly its ligands, backbone carbonyl groups and a water molecule. K+-induced polarization increases significantly at the transition state (TS) between the two binding sites. These findings suggest that polarization effects play a significant role in the microscopic mechanisms regulating potassium permeation.     

A density functional theory for Yukawa chain fluids in a nanoslit

A weighted density functional theory is developed for Yukawa chain fluids confined in a nanoslit. The excess free-energy functional is separated into repulsive and attractive contributions. A simple Heaviside function is used as the weighting function to calculate the weighted density in both contributions. The excess free-energy functional of repulsive interaction is calculated by the equation of state developed by Liu et al., while the contribution to excess free-energy functional by attractive interaction is calculated using the statistical associating fluids theory for chain molecules with attractive potentials of variable range. For pure fluids, the predicted density profiles near the nanoslit wall are in good agreement with simulations. The effect of cut-off introduced in the weighting function for the attractive part is examined; in addition, the surface excess and partition coefficient are calculated. The density profiles are also predicted for mixtures of two Yukawa chain fluids with different chain lengths, hard-core diameters, fluid–fluid and wall–fluid interactions. This work reveals that it is important to decompose the excess free-energy functional into repulsive and attractive contributions, and a simple weighting function can be used for both contributions.