Experimental and theoretical study on the supramolecular complexes of 15-crown-5 with adrenaline

The electron transfer properties of supramolecular complexes of 15-crown-5 (15C5) with protonated adrenaline (PAd⁺) at different electrodes using cyclic voltammetry (CV) have been investigated in the article. The experimental results show that 15C5 will affect the electron transfer properties of adrenaline. The formed supramolecular complexes by ion-dipole and hydrogen bond interaction between PAd⁺ and 15C5 will slow down the diffusion ability of adrenaline and make it hard to donate electron and be oxidized.The interaction energies and NPA calculations for the supramolecular complexes of 15C5 with PAd⁺ at B3LYP/6-31+G(d) level have been performed. The calculational results confirm the experimental fact that 15C5 can form stable supramolecular complexes with PAd⁺.     

Hydrogen-bonding interactions in adrenaline–water complexes: DFT and QTAIM studies of structures, properties, and topologies

ωB97XD/6-311++G(d,p) calculations were carried out to investigate the hydrogen-bonding interactions between adrenaline (Ad) and water. Six Ad-H(2)O complexes possessing various types of hydrogen bonds (H-bonds) were characterized in terms of their geometries, energies, vibrational frequencies, and electron-density topology. Natural bond orbital (NBO) and quantum theory of atoms in molecules (QTAIM) analyses were performed to elucidate the nature of the hydrogen-bonding interactions in these complexes. The intramolecular H-bond between the amino and carboxyl oxygen atom of Ad was retained in most of the complexes, and cooperativity between the intra- and intermolecular H-bonds was present in some of the complexes. H-bonds in which hydroxyls of Ad/water acted as proton donors were stronger than other H-bonds. Both hydrogen-bonding interactions and structural deformation play important roles in the relative stabilities of the complexes. The intramolecular H-bond was broken during the formation of the most stable complex, which indicates that Ad tends to break the intramolecular H-bond and form two new intermolecular H-bonds with the first water molecule.     

Computational Study of the Sulfonylated Amino Acid Hydroxamates Binding to the Zinc Ion within the Active Site of Carbonic Anhydrase

Abstract

The B3LYP/6–311+G(d,p) method and three ONIOM extrapolation methods ONI-OM (B3LYP/6–311+G(d,p): AM1); ONIOM(B3LYP/6–311+G(d,p): MNDO); ONIOM (B3LYP/6–311+G(d,p): HF/3-21G(d)) were used to characterize the complexes of Zn²⁺ cation with anionic sulfonylated amino acid hydroxamates (RSO₂NH-AA-CON(-)OH), possessing an unsubstituted RSO₂NH—amino acyl moiety. According to the R moiety we distinguish between pentafluorophenyl and 4-methoxyphenyl derivates. The amino acid hydroxamates included in the study were the Gly, Ala, and Leu derivates. Of the inhibitors investigated, the weakest zinc affinity exhibits the pentafluorophenyl derivate with Gly amino acid and the strongest affinity the 4-methoxyphenyl derivate with Leu amino acid. The inhibitors form bidentate coordination bonds with the zinc cation by means of the sulfonyl oxygen and the ionized hydroxamate nitrogen atoms, respectively. The zinc affinities computed using the B3LYP/6–311 +G(d,p)//HF/6–31 +G(d,p) method are in very good agreement with the full density functional theory (DFT) B3LYP/6–311+G(d,p)//B3LYP/6- 311+G(d,p) method and this method can be adopted to model larger complexes of inhibitors with the active site of carbonic anhydrase.     

Molecular structure,vibrational spectra and HOMO,LUMO analysis of 4-piperidone by density functional theory and ab initio Hartree–Fock calculations

Vibrational frequencies and geometrical parameters of 4-piperidone (4-PID) in the ground state have been calculated by using the Hartree–Fock (HF) and density functional methods (B3LYP) with 6-311++G(d,p) and 6-311+G(3df,2p) basis sets. These methods are proposed as a tool to be applied in the structural characterisation of 4-PID (C₅H₉NO). The title molecule has C _s point group symmetry, thus providing useful support in the interpretation of experimental IR and Raman data. The DFT-B3LYP/6-311+G(3df,2p) calculations have been found more reliable than the ab initio HF/6-311++G(d,p) calculations for the vibrational study of 4-PID. The calculated highest occupied molecular orbital and lowest unoccupied molecular orbital energies show that charge transfer occurs within the molecule. The theoretical spectrograms for FT-IR and FT-Raman spectra of the title molecule have been constructed.     

Comparative theoretical study of the UV/Vis absorption spectra of styrylpyridine compounds using TD-DFT calculations

This study examined absorption properties of 2-styrylpyridine, trans-2-(m-cyanostyryl)pyridine, trans-2-[3-methyl-(m-cyanostyryl)]pyridine, and trans-4-(m-cyanostyryl)pyridine compounds based on theoretical UV/Vis spectra, with comparisons between time-dependent density functional theory (TD-DFT) using B3LYP, PBE0, and LC-ωPBE functionals. Basis sets 6–31G(d), 6–31G(d,p), 6–31+G(d,p), and 6–311+G(d,p) were tested to compare molecular orbital energy values, gap energies, and maxima absorption wavelengths. UV/Vis spectra were calculated from fully optimized geometry in B3LYP/6–311+G(d,p) in gas phase and using the IEFPCM model. B3LYP/6–311+G(d,p) provided the most stable form, a planar structure with parameters close to 2-styrylpyridine X-ray data. Isomeric structures were evaluated by full geometry optimization using the same theory level. Similar energetic values were found: ～4.5 kJ?mol^?1 for 2-styrylpyridine and ～1 kJ?mol^?1 for derivative compound isomers. The 2-styrylpyridine isomeric structure differed at the pyridine group N-atom position; structures considered for the other compounds had the cyano group attached to the phenyl ring m-position equivalent. The energy difference was almost negligible between m-cyano-substituted molecules, but high energy barriers existed for cyano-substituted phenyl ring torsion. TD-DFT appeared to be robust and accurate approach. The B3LYP functional with the 6–31G(d) basis set produced the most reliable λ_max values, with mean errors of 0.5 and 12 nm respect to experimental values, in gas and solution, respectively. The present data describes effects on the λ_max changes in the UV/Vis absorption spectra of the electron acceptor cyano substituent on the phenyl ring, the electron donor methyl substituent, and the N-atom position on the electron acceptor pyridine ring, causing slight changes respect to the 2-styrylpyridine title compound.     

Influence of N-H...O and O-H...O hydrogen bonds on the (17)O, (15)N and (13)C chemical shielding tensors in crystalline acetaminophen: a density functional theory study

A computational investigation was carried out to characterize the (17)O, (15)N and (13)C chemical shielding tensors in crystalline acetaminophen. We found that N-H...O and O-H...O hydrogen bonds around the acetaminophen molecule in the crystal lattice have different influences on the calculated (17)O, (15)N and (13)C chemical shielding eigenvalues and their orientations in the molecular frame of axes. The calculations were performed with the B3LYP method and 6-311++G(d, p) and 6-311+G(d) standard basis sets using the Gaussian 98 suite of programs. Calculated chemical shielding tensors were used to evaluate the (17)O, (15)N, and (13)C NMR chemical shift tensors in crystalline acetaminophen, which are in reasonable agreement with available experimental data. The difference between the calculated NMR parameters of the monomer and molecular clusters shows how much hydrogen-bonding interactions affect the chemical shielding tensors of each nucleus. The computed (17)O chemical shielding tensor on O(1), which is involved in two intermolecular hydrogen bonds, shows remarkable sensitivity toward the choice of the cluster model, whereas the (17)O chemical shielding tensor on O(2) involved in one N-H...O hydrogen bond, shows smaller improvement toward the hydrogen-bonding interactions. Also, a reasonably good agreement between the experimentally obtained solid-state (15)N and (13)C NMR chemical shifts and B3LYP/6-311++G(d, p) calculations is achievable only in molecular cluster model where a complete hydrogen-bonding network is considered. Moreover, at the B3LYP/6-311++G(d, p) level of theory, the calculated (17)O, (15)N and (13)C chemical shielding tensor orientations are able to reproduce the experimental values to a reasonably good degree of accuracy.     

Gas phase acidities of N-substituted amine-boranes

Complexation energies and acidities of 19 primary, secondary and tertiary amine-boranes were investigated using MP2/6-311+G(d,p) and B3LYP/6-311+G(d,p) methods. Gas phase acidities for free amines were also calculated. Acidity values for studied complexes range from 327.3 to 349.1 kcal mol^?1 and the most acidic are the ones with direct connection between deprotonation center and a π-system. Results obtained by both computational methods are in good agreement with each other and with known experimental data. Addition of BH₃ increases the acidity of amines by 30 to 50 kcal mol^?1. This enhancement effect was compared to the respective effect witnessed in phosphine-boranes and traced back to changes of charge delocalization on nitrogen. A question about the structural stability of several deprotonated amine-borane anions in the gas phase was also raised.     

Coordination and thermodynamics of stable Zn(II) complexes in the gas phase

Ab initio molecular orbital methods in combination with DFT calculations were used to study the structural and thermodynamic properties of 17 complexes containing zinc cation and four first-shell ligands as models of active site of metalloenzymes (e.g. angiotensin converting enzyme, thermolysin). The geometry of the complexes was relaxed by complete optimization by ab initio molecular orbital methods at Hertree-Fock level with 3-21G* basis set. Following single point calculation with tight SCF criteria at the B3LYP level with 6-311+G(2d,p) basis set was used to calculate accurate interaction enthalpies. The structure and thermodynamics of optimized complexes are discussed from the point of view of their biological importance.     

Quantum-chemical investigation of the structure and the antioxidant properties of α-lipoic acid and its metabolites

Quantum-chemical computations were used to investigate the structure-antioxidant parameter relationships of α-lipoic acid and its natural metabolites bisnorlipoic acid and tetranorlipoic acid in their oxidized and reduced forms. The enantiomers of lipoic and dihydrolipoic acid were optimized using the B3LYP/6-311+G(3df,2p), B3LYP/aug-cc-pVDZ and MP2(full)/6-31+G(d,p) levels of theory as isolated molecules and in the presence of water. The geometries of the metabolites and the values of their antioxidant parameters (proton affinity, bond dissociation enthalpy, adiabatic ionization potential, spin density, and the highest occupied molecular orbital energy) were calculated at the B3LYP/6-311+G(3df,2p) level of theory. The results obtained reveal similarities between these structures: a pentatomic, nonaromatic ring is present in the oxidized forms, while an unbranched aliphatic chain (as found in saturated fatty acids) is present in both the oxidized and the reduced forms. Analysis of the spin density and the highest occupied molecular orbital energy revealed that the SH groups exhibited the greatest electron-donating activities. The values obtained for the proton affinity, bond dissociation enthalpy and adiabatic ionization potential indicate that the preferred antioxidant mechanisms for α-lipoic acid and its metabolites are sequential proton loss electron transfer in polar media and hydrogen atom transfer in vacuum.     

Density functional theory study of N-H...O, O-H...O and C-H...O hydrogen-bonding effects on the 14N and 2H nuclear quadrupole coupling tensors of N-acetyl-valine

Hydrogen-bonding effects in the crystalline structure of N-acetyl-valine, NAV, were studied using the (14)N and (2)H quadrupole coupling tensors via density functional theory. The calculations were carried out at the B3LYP level with the 6-311++G(d,p) and 6-311+G(d) basis sets. The theoretical quadrupole coupling components and their relative orientation in the molecular frame axes at the nitrogen site are compared to experimental values. This nucleus is involved in a rather strong intermolecular O=CNH...O=CNH hydrogen bond, r(N-H...O(1))=2.04 A and angleN-H...O(1)=171.53 degrees. A reasonably good agreement between the experimentally obtained (2)H quadrupole coupling tensors and the B3LYP/6-311++G(d,p) calculations is achievable only in molecular model where a complete hydrogen-bonding network is considered.     

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.     

Structure and electronic properties of azadirachtin

We performed a combined DFT and Monte Carlo ¹³C NMR chemical-shift study of azadirachtin A, a triterpenoid that acts as a natural insect antifeedant. A conformational search using a Monte Carlo technique based on the RM1 semiempirical method was carried out in order to establish its preferred structure. The B3LYP/6-311++G(d,p), wB97XD/6-311++G(d,p), M06/6-311++G(d,p), M06-2X/6-311++G(d,p), and CAM-B3LYP/6-311++G(d,p) levels of theory were used to predict NMR chemical shifts. A Monte Carlo population-weighted average spectrum was produced based on the predicted Boltzmann contributions. In general, good agreement between experimental and theoretical data was obtained using both methods, and the ¹³C NMR chemical shifts were predicted highly accurately. The geometry was optimized at the semiempirical level and used to calculate the NMR chemical shifts at the DFT level, and these shifts showed only minor deviations from those obtained following structural optimization at the DFT level, and incurred a much lower computational cost. The theoretical ultraviolet spectrum showed a maximum absorption peak that was mainly contributed by the tiglate group.     

A computational investigation on the antioxidant potential of myricetin 3,4′-di-<Emphasis Type="Italic">O</Emphasis>-<Emphasis Type="Italic">α</Emphasis>-<Emphasis Type="Italic">L</Emphasis>-rhamnopyranoside

In this work, we present a computational study on the antioxidant potential of myricetin 3,4\(^{\prime }\)-di-O-α-L-rhamnopyranoside (Compound M). A density functional theory (DFT) approach with the B3LYP and LC-ωPBE functionals and with both the 6-311G(d,p) and 6-311+G(d,p) basis sets was used. The focus of the investigation was on the structural and energetic parameters including both bond dissociation enthalpies (BDEs) and ionization potentials (IPs), which provide information on the potential antioxidant activity. The properties computed were compared with BDEs and IPs available in the literature for myricetin, a compound well known for presenting antioxidant activity (and the parent molecule of the compound of interest in the present work). Myricetin 3,4\(^{\prime }\)-di-O-α-L-rhamnopyranoside presented the lowest BDE to be 79.13 kcal/mol (as determined using B3LYP/6-311G(d,p) in water) while myricetin has a quite similar value (within 3.4 kcal/mol). IPs computed in the gas phase [B3LYP/6-311G(d,p)] are 157.18 and 161.4 kcal/mol for myricetin 3,4\(^{\prime }\)-di-O-α-L-rhamnopyranoside and myricetin, respectively. As the values of BDEs are considerably lower than the ones probed for IPs (in the gas phase or in any given solvent environment), the hydrogen atom transfer mechanism is preferred over the single electron transfer mechanism. The BDEs obtained suggest that myricetin 3,4\(^{\prime }\)-di-O-α-L-rhamnopyranoside can present antioxidant potential as good as the parent molecule myricetin (a well-known antioxidant). Therefore, experimental tests on the antioxidant activity of Compound M are encouraged.     

Theoretical investigation on the structure and thermodynamic properties of the 2,4-dinitroimidazole complex with methanol

The structure and thermodynamic properties of the 2, 4-dinitroimidazole complex with methanol were investigated using the B3LYP and MP2(full) methods with the 6-31++G(2d,p) and 6-311++G(3df,2p) basis sets. Four types of hydrogen bonds [N–H?O, C–H?O, O–H?O (nitro oxygen) and O–H?π] were found. The hydrogen-bonded complex having the highest binding energy had a N–H?O hydrogen bond. Analyses of natural bond orbital (NBO) and atoms-in-molecules (AIM) revealed the nature of the intermolecular hydrogen-binding interaction. The changes in thermodynamic properties from monomers to complexes with temperatures ranging from 200.0 to 800.0 K were investigated using the statistical thermodynamic method. Hydrogen-bonded complexes of 2,4-dinitroimidazole with methanol are fostered by low temperatures.

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.     

Bis-dibenzo[a.i]fluorenylidene, does it exist as stable 1,2-diradical?

The geometries, energies, and electronic properties of the two possible configurations of bis-[dibenzo[a.i]fluorenylidene] were investigated theoretically by density functional theory DFT B3LYP at the UB3LYP/6-311?+?G(2d,p) // UB3LYP/6-31?+?G(d,p) level of theory. According to the performed calculations, it was found that the singlet is 3.4?kcal?mol^-1 lower in energy compared to triplet state at room temperature. This gap is compared with those of other alkenes like ethylene, (61.9?kcal?mol^-1) tetra-tert-butyethylene, (6.4?kcal?mol^-1) and bis-fluorenylidene (19.5?kcal?mol^-1). These results confirm the experimental findings of the paramagnetic properties determined by Franzen and Joschek. The low singlet-triplet gap in the case of bis-[dibenzo[a.i]fluorenylidene] is the result of a steric destabilization of the singlet due to strain and stabilization of the triplet electronic state by delocalization of each free electron within each aromatic moiety. This correlates with the special electronic structure of the triplet state of this compound, where facial interaction of two hydrogen atoms lying close to the lobes of each p-orbital occupied with a single electron at the distorted double bond in the triplet electronic state.

Computational study on C−H…π interactions of acetylene with benzene, 1,3,5-trifluorobenzene and coronene

Meta-hybrid density functional theory calculations using M06-2X/6-31+G(d,p) and M06-2X/6-311+G(d,p) levels of theory have been performed to understand the strength of C?H^…π interactions of two possible types for benzene-acetylene, 1,3,5-trifluorobenzene-acetylene and coronene-acetylene complexes. Our study reveals that the C?H^...π interaction complex where acetylene located above to the center of benzene ring (classical T-shaped) is the lowest energy structure. This structure is twice more stable than the configuration characterized by H atom of benzene interacting with the π-cloud of acetylene. The binding energy of 2.91 kcal/mol calculated at the M06-2X/6-311+G(d,p) level for the lowest energy configuration (1A) is in very good agreement with the experimental binding energy of 2.7?±?0.2 kcal/mol for benzene-acetylene complex. Interestingly, the C?H^...π interaction of acetylene above to the center of the aromatic ring is not the lowest energy configuration for 1,3,5-trifluorobenzene-acetylene and coronene-acetylene complexes. The lowest energy configuration (2A) for the former complex possesses both C?H^...π interaction and C?H^...F hydrogen bond, while the lowest energy structure for the coronene-acetylene complex involves both π-π and C?H^...π interactions. C?H stretching vibrational frequencies and the frequency shifts are reported and analyzed for all of the configurations. We observed red-shift of the vibrational frequency for the stretching mode of the C-H bond that interacts with the π-cloud. Acetylene in the lowest-energy structures of the complexes exhibits significant red-shift of the C?H stretching frequency and change in intensity of the corresponding vibrational frequency, compared to bare acetylene. We have examined the molecular electrostatic potential on the surfaces of benzene, 1,3,5-trifluorobenzene, coronene and acetylene to explain the binding strengths of various complexes studied here.     

Studies on molecular structure and tautomerism of a vitamin B<Subscript>6</Subscript> analog with density functional theory

This work presents a computational study on the molecular structure and tautomeric equilibria of a novel Schiff base L derived from pyridoxal (PL) and o-phenylenediamine by using the density functional method B3LYP with basis sets 6-31 G(d,p), 6-31++G(d,p), 6-311 G(d,p) and 6-311++G(d,p). The optimized geometrical parameters obtained by B3LYP/6-31 G(d,p) method showed the best agreement with the experimental values. Tautomeric stability study of L inferred that the enolimine form is more stable than its ketoenamine form in both gas phase and solution. However, protonation of the pyridoxal nitrogen atom (LH) have accelerated the formation of ketoenamine form, and therefore, both ketoenamine and enolimine forms could be present in acidic media.     

Theoretical study of the pre- and post-translational effects of adenine and thymine tautomers and methyl derivatives

The study of pre-translational effects (ionization, tautomerization) and post-translational effects (methylation) of adenine and thymine has only recently been the focus of some studies. These effects can potentially help regulate gene expression as well as potentially disrupt normal gene function. Because of this wide array of roles, greater insight into these effects in deoxyribonucleic acids (DNA) are paramount. There has been considerable research of each phenomenon (tautomerization, methylation and ionization) individually. In this work, we attempt to shed light upon the pre-translational effects and post translational effects of adenine and thymine by investigating the electron affinities (EAs) and ionization potentials (IPs) of the major and minor tautomers and their methyl derivatives. We performed all calculations using the density functional theory (DFT) B3LYP functional accompanied with 6-311G(d,p), 6-311+G(d,p) and 6-311++G(df,pd) basis sets. Our results reveal that the thymine tautomer has a higher EA and IP than the adenine tautomers. The higher EA suggests that an electron that attaches to the AT base pair would predominately attach to the thymine instead of adenine. The higher IP would suggest that an electron that is removed from the AT base pair would be predominately removed from the adenine within the base pair. Understanding how tautomerization, ionization and methylation differences change effects, discourages, or promotes one another is lacking. In this work, we begin the steps of integrating these effects with one another, to gain a greater understanding of molecular changes in DNA bases.