Evaluation of the chemical reactivity in lignin precursors using the Fukui function

The hydroxycinnamyl alcohols: p-coumarol, coniferol and sinapol are considered the basic units and precursors of lignins models. In this work, the specific reactivity of these molecules was studied. We investigate their intrinsic chemical reactivity in terms of the Fukui function, applying the principle of hard and soft acids and bases (HSAB) in the framework of the density functional theory (DFT). Comparisons of their nucleophilic, electrophilic and free radical reactivity show their most probably sites to form linkages among them. It is found that the most reactive sites, for reactions involving free radicals, are the carbons at the β-position in the p-coumarol and sinapol molecules, whilst the regions around the carbon-oxygen bond of the phenoxyl group are the most reactive in coniferol. Figure Isocontour plots for the free radical form of the Fukui function f ⁰ (r), showing the reactive sites toward electron-rich/poor reactants     

Theoretical studies on interaction of anticancer drugs (dacarbazine,procarbazine and triethylenemelamine) with normal (AT and GC) and mismatch (GG,CC, AA and TT) base pairs

This study is an attempt to gain a better understanding of the physicochemical interaction between novel anticancer drugs and DNA bases. We have employed quantum chemical tools to explore the interaction of a few anticancer drugs [namely procarbazine (PR), dacarbazine (DC) and triethylenemelamine (TR)] with isolated normal (GC and AT) and mismatch (AA, CC, GG and TT) base pairs. The molecular geometries, electronic structural stability, vibrational energies, chemical reactivity and other electronic properties were studied using MP2/6-311+G^**, B3LYP/6-311+G^** and M05-2X/6-311+G^** methods. The optimised geometries of the usual and mismatch base pairs are almost planar whereas the geometries of drug-interacting complexes deviate from planarity. The presence of steric hindrance and π-bond overlaps between C–C bonds in the complexes has distorted the planarity of the four- and five-member rings in the base pairs. Among the three drugs chosen, DC and PR bond well with normal and mismatch base pairs with large interaction energy. The electron density (ED) difference maps of the most stable GG–DC, GG–PR and GG–TR drug-interacting complexes show the information about sharing of ED and gain or loss of ED within the interacting molecules. The stabilisation energy of the charge transfer interaction between the relevant donor–acceptor orbital of GG–DC and GC–DC complexes has been found to be around 16 kcal/mol and GG–PR and GC–PR complexes has been found to be around 12 kcal/mol. But, for the GG–TR and GC–TR complexes, the stabilisation energy is found to be less than 6 kcal/mol. Moreover, the topological analysis of hydrogen bond network of DC and PR drug-interacting complexes have high electron and Laplacian density with structural stability at the bond critical points (BCPs), while compared TR drug-interacting complexes by atoms in molecules and natural bond orbital analysis. Finally, we may conclude that the drugs DC and PR are highly efficient drugs to target normal and mismatch base pair for control and inhibition of DNA replication.     

Interactions of anticancer drugs with usual and mismatch base pairs - density functional theory studies

The antitumor activity of a drug is associated with its molecular properties as well as its interactions with target molecules. The molecular structures of usual, mismatch base pairs and their drug (Hydroxyurea and 5-Fluorouracil) interacting complexes were studied using density functional theory methods. The two and three-body interaction energies have been used to analyze the influence of a drug on the stability of base pairs. The sharing of electron density between the interacting molecules is shown through electron density difference maps. The Atoms in Molecules theory and Natural Bond Orbital analysis have been performed to study the hydrogen bonds in the drug interacting complexes.     

Theoretical investigations on the hydrogen bonding of nitrile isomers with H2O,HF, NH3 and H2S

The hydrogen bonds formed by the interaction of nitriles with water, hydrogen fluoride, ammonia and hydrogen sulphide have been studied using B3LYP and second-order Møller–Plesset perturbation (MP2) methods and 6-311+ + G(d,p) basis set. The energies and structures of 80 hydrogen-bonded complexes between nitriles and small molecules were examined systematically using B3LYP and MP2 procedure. Categorisation of the hydrogen bonds involved in the various complexes led to an ordering of hydrogen bond donor and acceptor abilities for some functional groups. The interaction energies have been corrected for the basis set superposition error using Boy's counterpoise correction method. The Morokuma energy decomposition analysis reveals that the strong interactions are due to the attractive contributions from the electrostatic (ES), polarisation (PL) and charge transfer (CT) components. The topological parameters, electron density and Laplacian of electron density show excellent correlation with the hydrogen bond length. Natural bond orbital (NBO) analysis has also been performed to study the CT from proton acceptor to the antibonding orbital of the H–Y bond in the proton donor part of complexes. The frequency analysis of C–H…Y bond in the complexes indicates the blue-shifting nature largely in case of sp² hybridised carbon atom.     

Topological analysis of tetraphosphorus oxides (P4O6+n (n = 0–4))

Quantum chemical calculations were used to analyze the chemical bonding and the reactivity of phosphorus oxides (P₄O_6+n (n?=?0–4)). The chemical bonding was studied using topological analysis such as atoms in molecules (AIM), electron localization function (ELF), and the reactivity using the Fukui function. A classification of the P-O bonds formed in all structures was done according to the coordination number in each P and O atoms. It was found that there are five P-O bond types and these are distributed among the five phosphorus oxides structures. Results showed that there is good agreement among the evaluated properties (length, bond order, density at the critical point, and disynaptic population) and each P-O bond type. It was found that regardless of the structure in which a P-O bond type is present the topological and geometric properties do not have a significant variation. The topological parameters electron density and Laplacian of electron density show excellent linear correlation with the average length of P-O bond in each bond type for each structure. From the Fukui function analysis it was possible to predict that from P₄O₆ until P₄O₈ the most reactive regions are basins over the P.     

A B3LYP and MP2(full) theoretical investigation into explosive sensitivity upon the formation of the molecule-cation interaction between the nitro group of 3,4-dinitropyrazole and H+, Li+, Na+, Be2+ or Mg2+

The explosive sensitivity upon the formation of molecule-cation interaction between the nitro group of 3,4-dinitropyrazole (DNP) and H⁺, Li⁺, Na⁺, Be²⁺ or Mg²⁺ has been investigated using the B3LYP and MP2(full) methods with the 6-311++G** and 6-311++G(2df,2p) basis sets. The bond dissociation energy (BDE) of the C3–N7 trigger bond has also been discussed for the DNP monomer and the corresponding complex. The interaction between the oxygen atom of nitro group and H⁺ in DNP…H⁺ is partly covalent in nature. The molecule-cation interaction and bond dissociation energy of the C3–N7 trigger bond follow the order of DNP…Be²⁺ > DNP…Mg²⁺ > DNP…Li⁺ > DNP…Na⁺. Except for DNP…H⁺, the increment of the trigger bond dissociation energy in comparison with the DNP monomer correlates well with the molecule-cation interaction energy, natural charge of the nitro group, electron density ρ _BCP(C3–N7), delocalization energy E ⁽²⁾ and NBO charge transfer. The analyses of atoms in molecules (AIM), natural bond orbital (NBO) and electron density shifts have shown that the electron density of the nitro group shifts toward the C3–N7 trigger bond upon the formation of the molecule-cation interaction. Thus, the trigger bond is strengthened and the sensitivity of DNP is reduced.     

Combined theoretical studies on solvation and hydrogen bond interactions in glycine tripeptide

The theoretical conformational analysis of glycine tripeptide (GT) has been carried out by molecular dynamics (MD) method in order to find minimum energy conformations. The MD studies on GT with water have been carried out for over 10 ns with a time step of 2 fs using fixed charge force field (AMBER ff03). By adding the solvation effect using water as a solvent, the GT conformers identified in this study exhibit α-helical conformation. Compared with the earlier reports, this MD study is able to identify the energetically favourable GT conformations. The obtained geometry of the five most stable GT conformations was optimised using the density functional theory method at B3LYP/6-311G** level of theory. Subsequently, the effects of solvation on the conformational characteristics of five most stable GT conformers with four water molecules (the number of water molecules in the first solvation shell of GT obtained from MD study) were investigated using the same method and the same level of theory. The effect of microsolvation on the fifth GT conformer has been studied with a cluster of 11 water molecules as the first hydration shell which generates folded structure. The interaction energies of all the complexes are calculated by correcting the basis set superposition error. The strong hydrogen bond mainly contributes to the interaction energies. The atoms in molecules theory and natural bond orbital analysis were used to study the origin of H-bonds. A good correlation between the structural parameters and the properties of charge density is found. NMR calculations show that the C = O carbons of the amine groups of the first and middle glycine fragments have maximum chemical shifts.     

Spectroscopic investigations and hydrogen bond interactions of 8-aza analogues of xanthine, theophylline and caffeine: a theoretical study

The structure, spectral properties and the hydrogen bond interactions of 8-aza analogues of xanthine, theophylline and caffeine have been studied by using quantum chemical methods. The time-dependent density functional theory (TD-DFT) and the singly excited configuration interaction (CIS) methods are employed to optimize the excited state geometries of isolated 8-azaxanthine, 8-azatheophylline tautomers and 8-azacaffeine in both the gas and solvent phases. The solvent phase calculations are performed using the polarizable continuum model (PCM). The absorption and emission spectra are calculated using the time-dependent density functional theory (TD-DFT) method. The results from the TD-DFT calculations reveal that the excitation spectra are red shifted relative to absorption in aqueous medium. These changes in the transition energies are qualitatively comparable to the experimental data. The examination of molecular orbital reveals that the molecules with a small H→L energy gap possess maximum absorption and emission wavelength. The relative stability and hydrogen bonded interactions of mono and heptahydrated 8-azaxanthine, 8-azatheophylline tautomers and 8-azacaffeine have been studied using the density functional theory (DFT) and Møller Plesset perturbation theory (MP2) implementing the 6-311++G(d,p) basis set. The formation of strong N-H…O bond has resulted in the highest interaction energy among the monohydrates. Hydration does not show any significant impact on the stability of heptahydrated complexes. The atoms in molecule (AIM) and natural bonding orbital (NBO) analyses have been performed to elucidate the nature of the hydrogen bond interactions in these complexes.

Iterative one-pot syntheses of chitotetroses

Rapid syntheses of chitotetrose derivatives were achieved in good yields using the newly developed reactivity independent iterative one-pot strategy. The protective groups on donors and acceptors were independently evaluated allowing matching of the two partners in glycosylation. No anomeric reactivity adjustments or intermediate purification were necessary thus significantly improving the overall synthetic efficiency. Only near stoichiometric amounts of building blocks were required for the assembly of target molecules further highlighting the potential of the iterative one-pot method in complex oligosaccharide synthesis.     

A new visualization scheme of chemical energy density and bonds in molecules

Covalent bond describes electron pairing in between a pair of atoms and molecules. The space is partitioned in mutually disjoint regions by using a new concept of the electronic drop region R_D, atmosphere region R_A, and the interface S (Tachibana in J Chem Phys 115:3497–3518, 2001). The covalent bond formation is then characterized by a new concept of the spindle structure. The spindle structure is a geometrical object of a region where principal electronic stress is positive along a line of principal axis of the electronic stress that connects a pair of the R_Ds of atoms and molecules. A new energy density partitioning scheme is obtained using the Rigged quantum electrodynamics (QED). The spindle structure of the stress tensor of chemical bond has been disclosed in the course of the covalent bond formation. The chemical energy density visualization scheme is applied to demonstrate the spindle structures of chemical bonds in H₂, C₂H₆, C₂H₄ and C₂H₂ systems.Figure Field theory of the energy density.     

Theoretical study of the interaction between 5-methylcytosine and acrylamide

The hydrogen-bonded complexes between 5-methylcytosine and acrylamide have been investigated using the density function theory (DFT) method. Five stable complexes have been found with no imaginary frequencies. Complex C3 is the most stable one with interaction energies of -69.01?kJ?mol(-1) corrected for basis set superposition error (BSSE). The charge change in the process of these complexes formation has also been examined. The atoms in molecules (AIM) theory and natural bond orbital (NBO) method have been performed to investigate the hydrogen bonds involved in all the complexes. The electron density and its corresponding Laplacian at the bond and ring critical points have been analyzed. In C3 complex, there is the largest stabilization energy (18.17?kJ?mol(-1)) between N11-H12 antibonding orbital and lone electron pair of O17. It can be seen that the hydrogen bonds play a crucial role in the stability of all the complexes between 5-methylcytosine and acrylamide. The theoretical results could provide helpful information for other researchers in further work.     

Computation of the bond dissociation enthalpies and free energies of hydroxylic antioxidants using the ab initio Hartree–Fock method

Abstract

Introduction

A new method for calculating theoretical bond dissociation enthalpy (BDE) and bond dissociation free energy (BDFE) of hydroxylic antioxidants is forwarded. BDE and BDFE may be understood as activation energies accompanying the formation of transition states, which may undergo downhill homolytic dissociation. The new method does not involve the complete fission of O–H bonds.

Method

Theoretical gas phase BDE values were calculated with the ab initio unrestricted Hartree–Fock (UHF) method, as changes in enthalpy between ground singlet states (GS) and triplet dissociative states (DS). Similarly, gas phase BDFEs were estimated from the corresponding changes in Gibbs free energy. The results were then compared with reliable experimental reports.

Results

The proposed theoretical approach of BDE and BDFE determination was tested using 10 simple phenols, 5 flavonoids, and l-ascorbic acid derivatives. The agreement between our calculated gas phase results and the adopted experimental values were generally within 0.5 kcal mol^?1, with a very few exceptions.

Discussion

Generally, steric interactions as well as intramolecular hydrogen bonding involving the dissociating OH group should be minimized in the GS. The DS are both electronically and vibrationally exited transition states. They have one unpaired electron on the carbon atom, which bears the homolytically dissociating OH group and are second order saddle points with a fixed <C–O–H bond angel of 180°.

Conclusion

It was concluded that ab initio UHF was well suited for the estimation of gas phase BDE and BDFE. The method presented has a good potential for application across a range of hydroxylic antioxidants. Currently, work is underway to extend its application in other class of antioxidants.     

Relative roles of K region and bay region towards determining the carcinogenic potencies of polycyclic aromatic hydrocarbons

Using the quantum mechanical MINDO/3 method a bond reactivity index based on superdelocalization of electrons in a compound has been formulated in order to estimate the relative carcinogenic properties of polycyclic aromatic hydrocarbons. The values of the index calculated for the K region as well as the bay region suggest that a knowledge of the reactivities of both regions is necessary for a reasonable estimation of the relative carcinogenic potencies of PAHs. However, no direct relationship between the two reactivity indices was observed in the calculations. It was found that the results for the bay-region reactivity index correlated well with the relative carcinogenic potencies of the molecules.     

Applications of natural orbitals for chemical valence in a description of bonding in conjugated molecules

Natural orbitals for chemical valence (NOCV) were used to describe bonding in conjugated pi-electron molecules. The 'single' C-C bond in trans-1,3-butadiene, 1,3-butadiene-1,1,4,4-tetra-carboxilic acid, 1,3,5,7-octatetraene, and 11-cis-retinal was characterized. In the NOCV framework, the formation of the sigma-bond appears as the sum of two complementary charge transfer processes from each vinyl fragment to the bond region, and partially to the other fragment. The formation of the pi-component of the bond is described by two pairs of NOCV representing the transfer of charge density from the neighboring 'double' C-C bonds. The NOCV eigenvalues and the related fragment-fragment bond multiplicities were used as quantitative measures of the sigma- and pi- contributions. The sigma-component of the 'single' C-C bonds appears to be practically constant in the systems analyzed, whereas the pi-contributions increase from butadiene (ca. 7.5%) to retinal (ca. 14%).     

Theoretical studies on the interaction between the nitrile-based inhibitors and the catalytic triad of Cathepsin K

Computational studies on the interaction of novel inhibitor compounds with the Cathepsin K protease have been performed to study the inhibition properties of the inhibitor compounds. The quantum chemical calculations have been performed to analyze the molecular geometries, structural stability, reactivity, nature of interaction, and the charge transfer properties using B3LYP level of theory by implementing 6-311g(d,p) basis set. The calculated C–S and N–H…N bond lengths of the inhibitor-triad complexes are found to agree well with the previous literature results. The chemical reactivity of the inhibitors and catalytic triad are analyzed through frontier molecular orbital analysis and found that the inhibitors are subjected to nucleophilic attack by the catalytic triad. The nature of inhibition of the inhibitor compounds is examined using the quantum theory of Atoms in Molecules analysis and found to be partially covalent. The NBO stabilization energy for the Cys – inhibitor are found to be most stable than the other interactions. The molecular dynamic simulations were performed to study the influence of dynamic of the active site on the QM results. The many body decomposition interaction energy calculated for the final results of MD simulation reveals that the dynamic of the active site induces significant changes in the interaction energy and occupancy of H-bonds plays a major role in the stabilizing the active site inhibitor interactions. The present study reveals that the inhibitor compounds can inhibit the proteolytic activity of the proteases on binding with the catalytic active site.     

B3LYP/6-311++G** geometry-optimization study of pentahydrates of alpha- and beta-D-glucopyranose

Five water molecules were placed in 37 different configurations around alpha- and beta-D-glucopyranose in the gt, gg, and tg conformational states, and the glucose-water complexes were geometry optimized using density functionals at the B3LYP/6-311++G** level of theory. The five water molecules were organized in space and energy minimized using an empirical potential, AMB02C, and then further geometry optimized using DFT algorithms to minimum energy positions. Electronic energy, zero point vibrational energy, enthalpy, entropy, stress energy on glucose and the water cluster, hydrogen-bond energy, and relative free energy were obtained for each configuration using thermodynamic procedures and an analytical Hessian program. The lowest energy complex was that of a clustering of water molecules around the 1- and 6-hydroxyl positions of the beta-gt anomer. Configurations in which the water molecules created a favorable network completely around and under glucose were found to have low energy for both alpha and beta anomers. Calculation of the alpha/beta anomeric ratio using the zero point corrected energy gave, approximately 32/68%, highly favoring the beta anomer in agreement with the experimental approximately 36/64% value. This ratio is better than the approximately 50/50% ratio found in our previous monohydrate study. An approximate hydroxymethyl population was obtained by noting average relative energies among the three conformational states, gt, gg, and tg. In the beta anomer complexes the gt conformation was favored over the gg state, while in the alpha anomer complexes the gg state was favored over the gt conformation, with the tg conformations all being of higher energy making little or no contribution to the rotamer population. Some geometry variances, found between glucose in vacuo and glucose after interaction with water molecules, are described and account for some observed C-5-C-6 bond length anomalies reported by us previously for the vacuum glucose structures.     

Morphology and reactivity characteristics of char biomass particles

In this work, 10 different biomasses were selected which included directly grown energy crops, industrial waste material and different wood types. Each biomass was sieved into six different size fractions and pyrolysed in a fixed bed furnace preheated to 1000 °C to produce a char residue. Intrinsic reactivity during burnout was measured using a non-isothermal thermogravimetric method. Scanning electron microscopy and oil immersion microscopy were used to characterise the morphology of the products. Char morphology was summarised in terms of degree of deformation, internal particle structure and wall thickness. Intrinsic reactivity corresponded directly with these morphology groupings showing a significant correlation between char morphotypes, char reactivity and the initial biomass material.     

Role of Molecular Electronic Properties of Some Novel Antimalarial Cyclic Peroxy Ketals in Relation to Potency and Potential Toxicity

Calculated molecular electronicproperties of 20 recently reported cyclic peroxy ketals have been rationalized with their in vitro antimalarial activity with the overall goal to guide the design for safer and effective peroxide-containing antimalarial agents. Stereoelectronicproperties were calculated on the optimized geometry of each compound using the ab initio 3-21G* quantum chemical basis set. Potency appears to be related to electronic properties rather than structural properties such as bond lengths and bond angles though an aliphatic cyclic ring seems to be a structural requirement for potent activity. Electronic properties such as differences in molecular polarity, in electrostatic potential profiles about the peroxide bond and the aromatic ring, in peroxide bond strength, and in the LUMO orbital energy of the molecules are all associated with potency. The three-dimensionalisopotential profile beyondthe van der Waals surface at –10 kcal/mol for the more potent analogs has a distinct large negative potential region by the aromatic ring extending to the methoxy moiety, suggesting a site for initial recognition interaction with the receptor away from the peroxide bond. The HOMO and LUMO isodensity surfaces for all molecules are located on the aromatic ring. The peroxide bond strength of the compounds is around 100 kcal/mol greater than the peroxide-containing clinically used antimalarial artemisinin compounds. In addition, density of the peroxy ketals also appears related to potent activity. The above features of the cyclic peroxy ketals are consistent with these compounds being less potent than the artemisinin compounds, but at the same time are less likely to be as neurotoxic as the artemisinins.