Ionic interactions for substituted MCH1R inhibitors studied by pKa values

MCH1R inhibitors with the quinoline moiety having the aromatic amine and aliphatic amine chain were selected, and then the effect of substituents of the quinoline ring on the ionic interaction were studied by calculating pK_a values for these amines at the B3LYP/6-311++G(d,p)//B3LYP/6-31+G(d) level in the gas phase and in water. For substituent with C, N, O, and S atoms next to the quinoline ring, respectively, the pK_a values of aromatic amines are estimated to be 8.98, 12.19, 4.64, and 4.33 and those of the aliphatic amines are 12.65, 10.82, 9.94, and 11.55, respectively.     

An assessment to evaluate the validity of different methods for the description of some corrosion inhibitors

New research and development efforts using computational chemistry in studying an assessment of the validity of different quantum chemical methods to describe the molecular and electronic structures of some corrosion inhibitors were introduced. The standard and the highly accurate CCSD method with 6-311++G(d,p), ab initio calculations using the HF/6-31G++(d,p) and MP2 with 6-311G(d,p), 6-31++G(d,p), and 6-311++G(2df,p) methods as well as DFT method at the B3LYP, BP86, B3LYP*, M06L, and M062x/6-31G++(d,p) basis set level were performed on some triazole derivatives and sulfur containing compounds used as corrosion inhibitors. Quantum chemical parameters, such as the energy of the highest occupied molecular orbital energy (E_HOMO), the energy of the lowest unoccupied molecular orbital energy (E_LUMO), energy gap (ΔE), dipole moment (μ), sum of total negative charges (TNC), chemical potential (Pi), electronegativity (χ), hardness (η), softness (σ), local softness (s), Fukui functions (f ⁺,f ^?), electrophilicity (ω), the total energy change (?E_T) and the solvation energy (S.E), were calculated. Furthermore, the accuracy and the applicability of these methods were estimated relative to the highest accuracy and standard CCSD with 6-311++G(d,p) method. Good correlations between the quantum chemical parameters and the corresponding inhibition efficiency (IE%) were found.     

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.     

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.     

A density functional theory approach toward substituent effect in Meerwein-Eschenmoser-Claisen rearrangement

B3LYP/6-31 G(d) level of theory has been used for the examination of substituent effect in the concerted step of the Meerwein-Eschenmoser-Claisen rearrangement. In this regard, the effect of NO₂ and NH₂ groups in different positions has been investigated. The obtained results show that substituent effect is very sensitive to its position and configuration. Electron withdrawing and electron donating groups in different positions and various configurations show different and sometimes opposite results.     

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.     

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.     

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

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

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.

Synthesis,electronic, and spectral properties of novel geranylated chalcone derivatives: a theoretical and experimental study

Novel chalcone derivatives with different substituents attached to A and B-rings: hydroxyl, methoxyl, geranyl, and prenyl groups were synthesized. The obtained compounds were characterized by NMR, HRMS, UV-Vis, IR, and MS. The theoretical analysis was carried out in all the compounds using density functional theory (DFT) with the B3LYP, PBE0, and M06-2X functionals in combination with the 6-311G(d,p) Pople-type basis set. The excited state properties were calculated by time dependent density functional theory (TD-DFT) using the same methodology applied for the ground state properties. The calculated vertical absorption wavelengths (λ_max) in gas phase and in ethanol as a solvent are consistent with the experimental ones, being the TD-DFT:B3LYP/6-311G(d,p) and PCM-TD-DFT:PBE0/6-311G(d,p) the best methodologies for these calculations with good approximation to the experimental values. The calculated reorganization energies indicated that, the four chalcone derivatives present an electron transfer character due to the smaller registered values. From these parameters it is proposed that these show an n-type semiconductor character. The localization of the frontier orbitals (HOMO and LUMO) shows that only the compound containing a hydroxyl group on the A-ring displays a marked delocalization favoring the charge-transfer process in this system. The HOMO-LUMO gap energies indicate that the inclusion of different donor groups in the rings does not improve the obtained values for this property.

Open image in new window

Graphical Abstract Relationship between spectroscopic and geometrical properties of chalcones were carried out using quantum-chemical calculations and compared with experimental values.

     

Quantum chemical study of the structure and properties of citrinin

Detailed structures and electronic properties of three tautomeric forms of the toxin citrinin were investigated using several quantum calculation methods. Energetic preference of the predominant p- and o-quinone methide tautomeric forms is dependent on the method of calculation. A previously unstudied carboxylic acid enol tautomer was calculated to be surprisingly stable in vacuo, being within 2.5 kcal mol^? 1 at the B3LYP/6-311++G(2d,2p) level of theory. Despite differences in bond nature and connectivity of tautomers, the natural bond orbital analysis revealed that tautomeric forms share similar natural charges and natural electron configurations. Calculated bond lengths corresponded with experimentally observed values and assignments for the calculated infrared vibrational frequencies are reported.     

Absolute configuration of (−)-myrtenal by vibrational circular dichroism

The VCD spectrum of the monoterpene (−)-myrtenal (1) was compared with theoretical spectra using ab initio density functional theory (DFT) calculations at the B3LYP/6-31G(d,p), B3LYP/6-31G+(d,p), B3LYP/6-311G+(d,p), B3LYP/DGDZVP, and B3PW91/DGTZVP levels of theory. Conformational analysis of 1 indicated that the lowest energy conformer was s-trans-C2-C10, which contributes more than 98.5% to the total conformational population regardless of the employed level of theory. The use of a recently developed confidence level algorithm demonstrated that VCD spectra calculated for the main conformer, using the indicated hybrid functionals and basis set, gave no significant changes, from where it follows that B3LYP/DGDZVP calculations provide a superior balance between computer cost and VCD spectral accuracy. The DGDZVP basis set demanded around a quarter the time than the 6-311G+(d,p) basis set while providing similar results. The spectral comparison also provided evidence that the levorotatory enantiomer of myrtenal has the 1R absolute configuration.     

Density functional theory investigation of electrophilic addition reaction of bromine to tricyclo[4.2.2.22,5]dodeca-1,5-diene1

The geometry and the electronic structure of tricyclo[4.2.2.2^2,5]dodeca-1,5-diene (TCDD) molecule were investigated by DFT/B3LYP and /B3PW91 methods using the 6-311G(d,p) and 6-311++G(d,p) basis sets. The double bonds of TCDD molecule are syn-pyramidalized. The structure of π-orbitals and their mutual interactions for TCDD molecule were investigated. Potential energy surface (PES) of the TCDD-Br₂ system was studied by B3LYP/6-311++G(d,p) method and the configurations [molecular charge-transfer (CT) complex, transition states (TS1 and TS2), intermediate (INT) and product (P)] corresponding to the stationary points (minima or saddle points) were determined. Initially, a molecular CT-complex forms between Br₂ and TCDD. With a barrier of 22.336 kcal mol^-1 the CT-complex can be activated to an intermediate (INT) with energy 15.154 kcal mol^-1 higher than that of the CT-complex. The intermediate (INT) then transforms easily (barrier 5.442 kcal mol^-1) into the final, N-type product. The total bromination is slightly exothermic. Accompanying the breaking of Br-Br bond, C1-Br, C5-Br and C2-C6 bonds are formed, and C1 = C2 and C5 = C6 double bonds transform into single bonds. The direction of the reaction is determined by the direction of intramolecular skeletal rearrangement that is realized by the formation of C2-C6 bond.

Modulating weak intramolecular interactions through the formation of beryllium bonds: complexes between squaric acid and BeH2

The electronic structure of the two most stable isomers of squaric acid and their complexes with BeH₂ were investigated at the B3LYP/6-311?+?G(3df,2p)// B3LYP/6-31?+?G(d,p) level of theory. Squaric acid forms rather strong beryllium bonds with BeH₂, with binding energies of the order of 60 kJ?mol^?1. The preferential sites for BeH₂ attachment are the carbonyl oxygen atoms, but the global minima of the potential energy surfaces of both EZ and ZZ isomers are extra-stabilized through the formation of a BeH···HO dihydrogen bond. More importantly, analysis of the electron density of these complexes shows the existence of significant cooperative effects between the beryllium bond and the dihydrogen bond, with both becoming significantly reinforced. The charge transfer involved in the formation of the beryllium bond induces a significant electron density redistribution within the squaric acid subunit, affecting not only the carbonyl group interacting with the BeH₂ moiety but significantly increasing the electron delocalization within the four membered ring. Accordingly the intrinsic properties of squaric acid become perturbed, as reflected in its ability to self-associate.

Interference of H-bonding and substituent effects in nitro- and hydroxy-substituted salicylaldehydes

Two intramolecular interactions, i.e., (1) hydrogen bond and (2) substituent effect, were analyzed and compared. For this purpose, the geometry of 4- and 5-X-substituted salicylaldehyde derivatives (X = NO₂, H or OH) was optimized by means of B3LYP/6-311 + G(d,p) and MP2/aug-cc-pVDZ methods. The results obtained allowed us to show that substituents (NO₂ or OH) in the para or meta position with respect to either OH or CHO in H-bonded systems interact more strongly than in the case of di-substituted species: 4- and 3-nitrophenol or 4- and 3-hydroxybenzaldehyde by ∼31%. The substituent effect due to the intramolecular charge transfer from the para-counter substituent (NO₂) to the proton-donating group (OH) is ∼35% greater than for the interaction of para-OH with the proton-accepting group (CHO). The total energy of H-bonding for salicylaldehyde, and its derivatives, is composed of two contributions: ∼80% from the energy of H-bond formation and ∼20% from the energy associated with reorganization of the electron structure of the systems in question.     

Racemization of the gastrointestinal antisecretory chiral drug esomeprazole magnesium via the pyramidal inversion mechanism: A theoretical study

The pyramidal inversion mechanisms of the 6‐methoxy and the 5‐methoxy tautomers of (S)‐omeprazole were studied, employing ab initio and DFT methods. The conformational space of the model molecule (S)‐2‐[(3‐methyl‐2‐pyridinyl)methyl]sulfinyl‐1H‐benzimidazole was calculated, with respect to rotations around single bonds, at the B3LYP/6‐311G(d,p) level. All of the resulting conformations were used as starting points for full optimizations of (S)‐omeprazole, at B3LYP/6‐31G(d), B3LYP/6‐311G(d,p), B3LYP/6‐311++G(d,p), B3LYP/6‐311G(2df,2pd), MP2/6‐31G(d), and MP2/6‐311G(d,p) levels. Four distinct pathways were found for enantiomerization via the pyramidal inversion mechanism for each of the tautomers of (S)‐omeprazole. Each transition state, in which the sulfur, the oxygen and the two carbon atoms connected directly to the sulfur are in one plane, connects two diastereomeric minima. The enantiomerization is completed by free rotation around the sulfur–methylene bond, and around the methylene–pyridine ring bond. The effective Gibbs' free energy barrier for racemization ΔG of the two tautomers of (S)‐omeprazole are 39.8 kcal/mol (5‐methoxy tautomer) and 40.0 kcal/mol (6‐methoxy tautomer), indicating that the enantiomers of omeprazole are stable at room temperature (in the gas phase). The 5‐methoxy tautomer of (S)‐omeprazole was found to be slightly more stable than the 6‐methoxy tautomer, in the gas phase. The energy barrier (ΔG^?) for the(S,M) (S,P) diastereomerization of (S)‐omeprazole due to the rotation around the pyridine chiral axis was very low, 5.8 kcal/mole at B3LYP/6‐311G(d,p). Chirality 2010. © 2010 Wiley‐Liss, Inc.     

σ*-Aromaticity in Three-Membered Rings

The structures and energies of unsaturated three-membered rings of the general formula (CH)₂XY_n^m, with charge = m and n substituents Y at X (Y = fluorine, chlorine, bromine, hydrogen, X = phosphorus, silicon, nitrogen, carbon) are compared to their saturated analogs. The structures were optimized with B3LYP/6-311+G(2d,p) and at MP2/6-31+G(d), with single point energy calculations on the latter geometries at MP4SDTQ/6-31+G(d). The geometrical changes in bond lengths and angles, which correlate with substituent electronegativities, are discussed for the different ring systems. The relative stabilities of unsaturated and saturated rings are compared using isodesmic ring-opening reactions and homodesmic substituent-exchange reactions. C*-Aromaticity, a hyperconjugative effect found in the disubstituted rings, causes lowering of ring strain energies for the unsaturated rings and preference of unsaturated rings over saturated ones for the more electronegative substituents. For the mono-substituted ?-aromatic silacyclopropenes and cyclopropenes, a destabilization by more electronegative ligands is found. For the neutral rings monosubstituted at main group V atoms like the 1H-phosphirenes and also the isoelectronic negatively charged rings with main group IV atoms like the silacyclopropenium anions, no correlation of stabilization energies or geometrical changes with ligand electronegativity is found.     

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.     

Probing the antioxidant potential of phloretin and phlorizin through a computational investigation

The structures and energetics of two dihydrochalcones (phloretin and its glycoside phlorizin) were examined with density functional theory, using the B3LYP, M06-2X, and LC-ω PBE functionals with both the 6-311G(d,p) and 6-311 + G(d,p) basis sets. Properties connected to antioxidant activity, i.e., bond dissociation enthalpies (BDEs) for OH groups and ionization potentials (IPs), were computed in a variety of environments including the gas-phase, n-hexane, ethanol, methanol, and water. The smallest BDEs among the four OH groups for phloretin (three for phlorizin) were determined (using B3LYP/6-311 + G(d,p) in water) to be 79.36 kcal/mol for phloretin and 79.98 kcal/mol for phlorizin while the IPs (at the same level of theory) were obtained as 139.48 and 138.98 kcal/mol, respectively. By comparing with known antioxidants, these values for the BDEs indicate both phloretin and phlorizin show promise for antioxidant activity. In addition, the presence of the sugar moiety has a moderate (0-6 kcal/mol depending on functional) effect on the BDEs for all OH groups. Interestingly, the BDEs suggest that (depending on the functional chosen) the sugar moiety can lead to an increase, decrease, or no change in the antioxidant activity. Therefore, further experimental tests are encouraged to understand the substituent effect on the BDEs for phloretin and to help determine the most appropriate functional to probe BDEs for dihydrochalcones.