Low-energy conformers of pamidronate and their intramolecular hydrogen bonds: a DFT and QTAIM study

Extensive DFT and ab initio calculations were performed to characterize the conformational space of pamidronate, a typical pharmaceutical for bone diseases. Mono-, di- and tri-protic states of molecule, relevant for physiological pH range, were investigated for both canonical and zwitterionic tautomers. Semiempirical PM6 method were used for prescreening of the single bond rotamers followed by geometry optimizations at the B3LYP/6-31++G(d,p) and B3LYP/6-311++G(d,p) levels. For numerous identified low energy conformers the final electronic energies were determined at the MP2/6-311++G(2df,2p) level and corrected for thermal effects at B3LYP level. Solvation effects were also considered via the COSMO and C-PCM implicit models. Reasonable agreement was found between bond lengths and angle values in comparison with X-ray crystal structures. Relative equilibrium populations of different conformers were determined from molecular partition functions and the role of electronic, vibrational and rotational degrees of freedom on the stability of conformers were analyzed. For no level of theory is a zwitterionic structure stable in the gas-phase while solvation makes them available depending on the protonation state. Geometrically identified intramolecular hydrogen bonds were analyzed by QTAIM approach. All conformers exhibit strong inter-phosphonate hydrogen bonds and in most of them the alkyl-amine side chain is folded on the P-C-P backbone for further hydrogen bond formation.

DFT studies of the conversion of four mesylate esters during reaction with ammonia

The energetics of the Menshutkin-like reaction between four mesylate derivatives and ammonia have been computed using B3LYP functional with the 6-31+G** basis set. Additionally, MPW1K/6-31+G** level calculations were carried out to estimate activation barrier heights in the gas phase. Solvent effect corrections were computed using PCM/B3LYP/6-31+G** level. The conversion of the reactant complexes into ion pairs is accompanied by a strong energy decrease in the gas phase and in all solvents. The ion pairs are stabilized with two strong hydrogen bonds in the gas phase. The bifurcation at C2 causes a significant activation barrier increase. Also, bifurcation at C5 leads to noticeable barrier height differentiation. Both B3LYP/6-31+G** and MPW1K/6-31+G** activation barriers suggest the reaction 2 (2a?+?NH₃) to be the fastest in the gas phase. The reaction 4 is the slowest one in all environments.

Predicting the preferred conformations of luteolin-4′-O-β-D-glucoside in gas phase: a comparison of two computational approaches

A tree-step computational approach has been applied to determine the lowest-energy conformers of luteolin-4′-O-β-D-glucoside (L4′G). Fifty-seven starting structures of the L4′G have been built, and then by performing with density functional theory (DFT) optimizations and second-order Møller-Plesset (MP2) calculations, the preferred conformations of L4′G are predicted. In order to test the accuracy of the computational approach, a hybrid Monte-Carlo multiple minimum (MCMM)/quantum mechanical (QM) approach is applied to determine the favorable conformers of L4′G. The alternative classification is employed to put similar conformations into the same catalogue according to the dihedral angles among the luteolin rings, glycosidic dihedral angles, and the orientations of hydroxyl and hydroxymethyl groups. The low-energy conformations are located after the optimizations at the HF/6-31G(d) and B3LYP/6-311+G(d) levels. Compared with the hybrid MCMM/QM approach, the tree-step computational approach not only remains accurate but also saves a lot of computing resources.

A B3LYP and MP2(full) theoretical investigation into the strength of the C-NO2 bond upon the formation of the molecule-cation interaction between Na+ and the nitro group of nitrotriazole or its methyl derivatives

The changes of bond dissociation energy (BDE) in the C–NO₂ bond and nitro group charge upon the formation of the molecule-cation interaction between Na⁺ and the nitro group of 14 kinds of nitrotriazoles or methyl derivatives were investigated using the B3LYP and MP2(full) methods with the 6-311++G**, 6-311++G(2df,2p) and aug-cc-pVTZ basis sets. The strength of the C–NO₂ bond was enhanced in comparison with that in the isolated nitrotriazole molecule upon the formation of molecule-cation interaction. The increment of the C–NO₂ bond dissociation energy (ΔBDE) correlated well with the molecule-cation interaction energy. Electron density shifts analysis showed that the electron density shifted toward the C-NO₂ bond upon complex formation, leading to the strengthened C-NO₂ bond and the possibly reduced explosive sensitivity.

A B3LYP and MP2(full) theoretical investigation into the strength of the C–NO2 bond upon the formation of the intermolecular hydrogen-bonding interaction between HF and the nitro group of nitrotriazole or its methyl derivatives

The changes of bond dissociation energy (BDE) in the C–NO₂ bond and nitro group charge upon the formation of the intermolecular hydrogen-bonding interaction between HF and the nitro group of 14 kinds of nitrotriazoles or methyl derivatives were investigated using the B3LYP and MP2(full) methods with the 6-311++G**, 6-311++G(2df,2p) and aug-cc-pVTZ basis sets. The strength of the C–NO₂ bond was enhanced and the charge of nitro group turned more negative in complex in comparison with those in isolated nitrotriazole molecule. The increment of the C–NO₂ bond dissociation energies correlated well with the intermolecular H-bonding interaction energies. Electron density shifts analyses showed that the electron density shifted toward the C–NO₂ bond upon complex formation, leading to the strengthened C–NO₂ bond and the possibly reduced explosive sensitivity.

Relationship between structure and entropy contributions in an anthraquinone mercapto derivative

The structural and thermodynamic properties of an anthraquinone derivative were studied by means of quantum-chemical calculations. Conformational analysis using ab initio and density functional theory methods revealed 14 low-energy conformers. In order to discuss similarities and differences in entropy of the conformers, the rotational and vibrational contributions to entropy were correlated with changes in conformer structure. The component of the moment of inertia perpendicular to the molecular plane gives significant input to ΔS _rot , whereas the largest contributions to the ΔS _vib have vibrations associated with the τ _S1C20 coordinate.

Electron density reactivity indexes of the tautomeric/ionization forms of thiamin diphosphate

The generation of the highly reactive ylide in thiamin diphosphate catalysis is analyzed in terms of the nucleophilicity of key atoms, by means of density functional calculations at X3LYP/6–31++G(d,p) level of theory. The Fukui functions of all tautomeric/ionization forms are calculated in order to assess their reactivity. The results allow to conclude that the highly conserved glutamic residue does not protonate the N1′ atom of the pyrimidyl ring, but it participates in a strong hydrogen bonding, stabilizing the eventual negative charge on the nitrogen, in all forms involved in the ylide generation. This condition provides the necessary reactivity on key atoms, N4′ and C2, to carry out the formation of the ylide required to initiate the catalytic cycle of ThDP- dependent enzymes. This study represents a new approach for the ylide formation in ThDP catalysis.

Zinc-, cadmium-, and mercury-containing one-dimensional tetraphenylporphyrin arrays: a DFT study

This work describes theoretical and experimental studies on glycerol esterification to obtain acetins focusing on the obtained isomers. The reaction of glycerol with acetic acid was carried out on Amberlyst 36 wet. Density functional theory calculations on the level of M06-2X functional and 6-311+G(d,p) basis set are carried out and the most stable structures of the reactants and products are located by considering a large number of conformers. The thermodynamics is discussed in terms of the calculated reaction Gibbs free energy. The AIM theory was used to characterize reactants and products. The glycerol esterification with acetic acid is found to be thermodynamically favored, with exothermal property. These agree well with experiments and allow us to explain the relative selectivity of products.

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.

Structural characterization of the (MeSH)4 potential energy surface

A random walk on the PES for (MeSH)₄ clusters produced 50 structural isomers held together by hydrogen-bonding networks according to calculations performed at the B3LYP/6–311++G** and MP2/6–311++G** levels. The geometric motifs observed are somewhat similar to those encountered for the methanol tetramer, but the interactions responsible for cluster stabilization are quite different in origin. Cluster stabilization is not related to the number of hydrogen bonds. Two distinct, well-defined types of hydrogen bonds scattered over a wide range of distances are predicted.

Assessing the quantum mechanical level of theory for prediction of linear and nonlinear optical properties of push-pull organic molecules

In this paper, we assessed the quantum mechanical level of theory for prediction of linear and nonlinear optical (NLO) properties of push-pull organic molecules. The electric dipole moment (μ), mean polarizability (〈α〉) and total static first hyperpolarizability (β_t) were calculated for a set of benzene, styrene, biphenyl and stilbene derivatives using HF, MP2 and DFT (31 different functionals) levels and over 71 distinct basis sets. In addition, we propose two new basis sets, NLO-V and aNLO-V, for NLO properties calculations. As the main outcomes it is shown that long-range corrected DFT functionals such as M062X, ωB97, cam-B3LYP, LC-BLYP and LC-ωPBE work satisfactorily for NLO properties when appropriate basis sets such as those proposed here (NLO-V or aNLO-V) are used. For most molecules with β ranging from 0 to 190 esu, the average absolute deviation was 13.2 esu for NLO-V basis sets, compared to 27.2 esu for the standard 6-31 G(2d) basis set. Therefore, we conclude that the new basis sets proposed here (NLO-V and aNLO-V), together with the cam-B3LYP functional, make an affordable calculation scheme to predict NLO properties of large organic molecules.

New insights into the stability of alkenes and alkynes, fluoro-substituted or not: a DFT, G4, QTAIM and GVB study

Many undergraduate organic chemistry books do not agree with the order of relative stability of alkenes towards hydrogenation reactions. Although they ascribe the stability of alkenes to the number and spatial position of the alkyl groups attached to the vinyl carbon atoms, results from the quantum theory of atoms in molecules indicate that the influence of an alkyl substituent on the stability of unsaturated hydrocarbons arises from the slight removal of electron density of the π bond, not from donation of their charge density to unsaturated carbon atoms as stated in many text books. There is an inverse relation between delocalization index—the number of shared electrons between two atoms, or Wiberg bond index of C=C bond—and the number of methyl groups attached to the vinyl carbon atoms. Electron withdrawing groups (EWGs) attached to unsaturated carbon atoms of alkenes and alkynes have two different behaviors: slight EWGs (alkyl groups) stabilize unsaturated carbon atoms while the strong EWG destabilizes the unsaturated carbon atoms. Generalized valence bond theory was also used to study the ambiguous behavior of fluorine substituents bonded to vinyl carbon atoms.

Enhancing and modulating the intrinsic acidity of imidazole and pyrazole through beryllium bonds

The structure and electronic properties of the complexes formed by the interaction of imidazole and pyrazole with different BeXH(BeX₂) (X = H, Me, F, Cl) derivatives have been investigated via B3LYP/6?311+G(3df,2p)//B3LYP/6?31+G(d,p) calculations. The formation of these azole:BeXH(BeX₂) complexes is accompanied by a dramatic enhancement of the intrinsic acidity of the azole, as the deprotonated azole is much more stable after the aforementioned interaction. Most importantly, the increase in acidity is so large that the azole:BeXH or azole:BeX₂ complexes behave as NH acids, which are stronger than typical oxyacids such as phosphoric acid and oxalic acid. Interestingly, the increase in acidity can be tuned through appropriate selection of the substituents attached to the Be atom, permitting us to modulate the electron-accepting ability of the BeXH or BeX₂ molecule.

Modeling the activity of glutathione as a hydroxyl radical scavenger considering its neutral non-zwitterionic form

Glutathione is an immensely important antioxidant, particularly in the central nervous system. The scavenging mechanism of glutathione towards the OH radical was studied theoretically, considering its neutral, non-zwitterionic form relevant to acidic media. Gibbs free barrier and released energies involved in hydrogen abstraction from the different sites of glutathione by an OH radical were studied at the B3LYP/6-31G(d,p), B3LYP/AUG-cc-pVDZ, M06/AUG-cc-pVDZ, M06-2X/AUG-cc-pVDZ levels of density functional theory. Solvation in bulk aqueous media was also studied at all these levels of theory employing the polarizable continuum model. Our study shows that a hydroxyl radical can abstract a hydrogen atom easily from glutathione. Thus, glutathione is shown to be an efficient scavenger of OH radicals, which is in agreement with the results of previous studies.

Computationally designed prodrugs of statins based on Kirby’s enzyme model

DFT calculations at B3LYP/6-31G(d,p) for intramolecular proton transfer in Kirby’s enzyme models 1–7 demonstrated that the reaction rate is dependent on the distance between the two reacting centers, r_GM, and the hydrogen bonding angle, α, and the rate of the reaction is linearly correlated with r_GM and α. Based on these calculation results three simvastatin prodrugs were designed with the potential to provide simvastatin with higher bioavailability. For example, based on the calculated log EM for the three proposed prodrugs, the interconversion of simvastatin prodrug ProD 3 to simvastatin is predicted to be about 10 times faster than that of either simvastatin prodrug ProD 1 or simvastatin ProD 2. Hence, the rate by which the prodrug releases the statin drug can be determined according to the structural features of the promoiety (Kirby’s enzyme model).

Theoretical design of donor-acceptor conjugated copolymers based on furo-, thieno-, and selenopheno[3,4-c] thiophene-4,6-dione and benzodithiophene units for organic solar cells

In this work, a series of donor-acceptor (D-A) copolymers (PBDTFPD(Pa1), PBDTTPD (Pa2) and PBDTSePD(Pa3)) were selected and theoretically investigated using O3LYP/6-31G(d), PBE0/6-31G(d), TD-O3LYP/6-31G(d)//O3LYP/6-31G(d) and periodic boundary conditions methods. The calculated results go well with the available experimental data of highest occupied and lowest unoccupied molecular orbital (HOMO/LUMO) energy levels and band gaps. A series of conjugated polymers (Pb1?~?Pb3) comprised of electron-deficient benzodithiophene and electron-rich furo-, thieno-, and selenopheno[3,4-c]thiophene-4,6-dione were further designed and studied. Compared with Pa1-Pa3, the designed polymers of Pb1?~?Pb3 show better performances with smaller band gaps, lower HOMO energy levels, red shift of absorption spectra, and larger open circuit voltage (Voc). For investigated polymers (Pa1, Pa2, Pa3, Pb1, Pb2, Pb3), the power conversion efficiencies (PCEs) of ~6.1 %, ~7.2 %, ~7.9 %, ~8.0 %, ~9.5 % and ~9.0 % are predicted by Scharber diagrams when they are used in combination with PC60BM as an acceptor. The results illustrate that these designed polymers which turn the electron-withdrawing capability in D-A conjugated polymers are expected to turn into highly efficient donor materials for organic solar cells.

Prodrugs of fumarate esters for the treatment of psoriasis and multiple sclerosis—a computational approach

Density functional theory (DFT) calculations at B3LYP/6-31 G (d,p) and B3LYP/6-311?+?G(d,p) levels for the substituted pyridine-catalyzed isomerization of monomethyl maleate revealed that isomerization proceeds via four steps, with the rate-limiting step being proton transfer from the substituted pyridinium ion to the C=C double bond in INT1. In addition, it was found that the isomerization rate (maleate to fumarate) is solvent dependent. Polar solvents, such as water, tend to accelerate the isomerization rate, whereas apolar solvents, such as chloroform, act to slow down the reaction. A linear correlation was obtained between the isomerization activation energy and the dielectric constant of the solvent. Furthermore, linearity was achieved when the activation energy was plotted against the pK _a value of the catalyst. Substituted-pyridine derivatives with high pK _a values were able to catalyze isomerization more efficiently than those with low pK _a values. The calculated relative rates for prodrugs 1–6 were: 1 (406.7), 2 (7.6?×?10⁶), 3 (1.0), 4 (20.7), 5 (13.5) and 6 (2.2?×?10³). This result indicates that isomerizations of prodrugs 1 and 3–5 are expected to be slow and that of prodrugs 2 and 6 are expected to be relatively fast. Hence, prodrugs 2 and 3–5 have the potential to be utilized as prodrugs for the slow release of monomethylfumarate in the treatment of psoriasis and multiple sclerosis.

Role of gold in a complex cascade reaction involving two electrocyclization steps

Quantum chemical computations (B3LYP/LACVP**) were applied to assess the impact of Au(I) complexation on activation barriers for sequential electrocyclization reactions (one a 1,2-dihydroazete ring-opening and another a pentadienyl cation ring-closure) proposed to occur during a complex reaction cascade that converts alkynes and imines to cyclopentenimines.

Regioselectivity in Sonogashira synthesis of 6-(4-nitrobenzyl)-2-phenylthiazolo[3,2-b]1,2,4-triazole: a quantum chemistry study

In the present research, the experimentally observed regioselectivity in Sonogashira synthesis of 6-(4-nitrobenzyl)-2-phenylthiazolo[3,2-b]1,2,4triazole has been modeled by means of density functional theory (DFT) employed to investigate the structural and thermochemical aspects of this synthesis in the gas and solution phases. Comparison of our calculated structural parameters of the title compound with the available X-ray crystallographical data demonstrate a reliable agreement. Then, the effect of two different solvents, DMF and ethanol, are examined via polarized continuum model calculations, showing a significant decrease in the computed values of the reaction enthalpy and free energy changes compared with the gas phase results. We have also considered two tautomeric structures of the intermediate species that it seems the mode of its intermolecular cyclization has an important role in regioselectivity of the final products. Moreover, all obtained results in the gas and solution phases also confirm that the synthesis of the title compound is thermodynamically more favorable than the other regioisomeric product. We also discuss the thermodynamical feasibility of this reaction at higher temperatures. Finally, we concentrate on the survey of substituent effect by choosing electron-withdrawing and electron-donating groups on the aryl iodide. Our calculated thermochemical data in the gas and solution phases indicate that the use of electron-withdrawing moieties is more favorable thermodynamically than electron-donating ones which has been previously concluded via the experimental elucidations.