Theoretical study of the decomposition of ethyl and ethyl 3-phenyl glycidate

The mechanism of the decomposition of ethyl and ethyl 3-phenyl glycidate in gas phase was studied by density functional theory (DFT) and MP2 methods. A proposed mechanism for the reaction indicates that the ethyl side of the ester is eliminated as ethylene through a concerted six-membered cyclic transition state, and the unstable intermediate glycidic acid decarboxylates rapidly to give the corresponding aldehyde. Two possible pathways for glycidic acid decarboxylation were studied: one via a five-membered cyclic transition state, and the other via a four-membered cyclic transition state. The results of the calculations indicate that the decarboxylation reaction occurs via a mechanism with five-membered cyclic transition state.

Theoretical evaluation of flotation performance of carboxyl hydroxamic acids with different number of polar groups on the surfaces of diaspore (010) and kaolinite (001)

The adsorption behaviors of three carboxyl hydroxamic acids on diaspore (010) and kaolinite (001) have been studied by density functional theory (DFT) and molecular dynamics (MD) method. The results indicated that carboxyl hydroxamic acids could adsorb on diaspore surface by ionic bonds and hydrogen bonds, and adsorb on kaolinite surface by hydrogen bonds. The models of carboxyl hydroxamic acids adsorbed on diaspore and kaolinite surfaces are proposed.

Heavy periodane

The potential energy surface of the hypothetical NaMgAlSiPSCl system (heavy periodane) is exhaustively analyzed via the gradient embedded genetic algorithm (GEGA) in combination with density functional theory (DFT) computations. The electronegativity differences among the elements in both the second and third rows of the periodic table indicate that low-energy heavy periodane structures are obtained when highly electronegative and electropositive elements are bound together, but the global minimum of the heavy periodane system is completely different to its second-row analog (LiBeBCNOF).

Exploration of various electronic properties along the reaction coordinate for hydration of Pt(II) and Ru(II) complexes; the CCSD, MPx, and DFT computational study

In the study behavior of molecular electrostatic potential, averaged local ionization energy, and reaction electronic flux along the reaction coordinate of hydration process of three representative Ru(II) and Pt(II) complexes were explored using both post-HF and DFT quantum chemical approximations. Previously determined reaction mechanisms were explored by more detailed insight into changes of electronic properties using ωB97XD functional and MP2 method with 6–311++G(2df,2pd) basis set and CCSD/6–31(+)G(d,p) approach. The dependences of all examined properties on reaction coordinate give more detailed understanding of the hydration process.

Study on the role of SBA-15 in the oxidative dehydrogenation of n-butane over vanadia catalyst using density functional theory

The first step in the mechanism of n-butane oxidative dehydrogenation (ODH) on a V₄O₁₀ cluster and V₄O₁₀ supported SBA-15 is examined using DFT method. The activation and adsorption energies, oxidation state of V atoms are calculated. Over V₄O₁₀ the obtained results indicate that the activation of C-H bond of methylene group can occur at both the terminal and the bridging oxygen atoms with similar barrier (21.5–22.5 kcal mol^?1). The role of SBA-15 (with and without modification by Al) in n-butane adsorption step has been studied in detail. SBA-15 itself has mild effect on the reaction process, but the substitution of silicon atoms by aluminum atoms results in an active supporter for V₂O₅ in ODH reaction. In that, the ratio of Si/Al will decide the direction of initial interaction steps between n-butane and catalyst surface and it will result in the selectivity of the reaction products.

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).

Industrial applications of in silico ADMET

Quantitative structure activity relationship (QSAR) modeling has been in use for several decades now. One branch of it, in silico ADMET, became more and more important since the late 1990s as studies indicated that poor pharmacokinetics and toxicity were important causes of costly late-stage failures in drug development. In this paper we describe some of the available methods and best practice for the different stages of the in silico model building process. We also describe some more recent developments, like automated model building and the prediction probability. Finally we will discuss the use of in silico ADMET for “big data” and the importance and possible further development of interpretable models.

Prodrugs for masking bitter taste of antibacterial drugs—a computational approach

DFT calculations for the acid-catalyzed hydrolysis of several maleamic acid amide derivatives revealed that the reaction rate-limiting step is determined on the nature of the amine leaving group. Further, it was established that when the amine leaving group was a secondary amine, acyclovir or cefuroxime moiety the tetrahedral intermediate formation was the rate-limiting step such as in the cases of acyclovir ProD 1- ProD 4 and cefuroxime ProD 1- ProD 4. In addition, the linear correlation between the calculated and experimental rates provided a credible basis for designing prodrugs for masking bitter taste of the corresponding parental drugs which have the potential to release the parent drug in a sustained release fashion. For example, based on the DFT calculated rates the predicted t_1/2 (a time needed for 50 % of the reactant to be hydrolyzed to products) for cefuroxime prodrugs, cefuroxime ProD 1- ProD 4, were 12 min, 18 min, 200 min and 123 min, respectively.

Nucleus-independent chemical shift criterion for aromaticity in π-extended tetraoxa[8]circulenes

Recently synthesized π-extended symmetrical tetraoxa[8]circulenes that exhibit electroluminescent properties were calculated at the density functional theory (DFT) level using the quantum theory of atoms in molecules (QTAIM) approach to electron density distribution analysis. Nucleus-independent chemical shift (NICS) indices were used to characterize the aromaticity of the studied molecules. The tetraoxa[8]circulene molecules were found to consist of two antiaromatic perimeters (according to the Hückel “4n” antiaromaticity rule) that include 8 and 24 π-electrons. Conversely, NICS calculations demonstrated the existence of a common π-extended system (distributed like a flat ribbon) in the studied tetraoxa[8]circulene molecules. Thus, these symmetrical tetraoxa[8]circulene molecules provide examples of diatropic systems characterized by the presence of induced diatropic ring currents.

A theoretical study on the hydrogen adducts of diamidocarbenes and diaminocarbenes

The hybrid-meta GGA DFT functional M06-2X was used to examine the potential of N,N′-diamidocarbenes for use as hydrogen storage materials. We previously discovered that borylene, which is isoelectronic with an Arduengo-type carbene, was a suitable candidate for a hydrogen storage material. We compared the capabilities of N,N′-diamidocarbenes and N-heterocyclic carbenes as hydrogen storage materials. The results indicate that diamidocarbenes are not suitable hydrogen storage materials because the removal of H₂ is more endothermic for diamidocarbenes than for diaminocarbenes.

Vibrational spectra of an RDX film over an aluminum substrate from molecular dynamics simulations and density functional theory

We report calculated vibrational spectra in the range of 0–3,500 cm^?1 of RDX (hexahydro-1,3,5-trinitro-1,3,5-triazine) molecules adsorbed on a model aluminum surface. A molecular film was modeled using two approaches: (1) density functional theory (DFT) was used to optimize a single RDX molecule interacting with its periodic images, and (2) a group of nine molecules extracted from the crystal structure was deposited on the surface and interacted with its periodic images via molecular dynamics (MD) simulations. In both cases, the molecule was initialized in the AAA conformer geometry having the three nitro groups in axial positions, and kept that conformation in the DFT examination, but some molecules were found to change to the AAE conformer (two nitro groups in axial and one in equatorial position) in the MD analysis. The vibrational spectra obtained from both methods are similar to each other, except in the regions where collective RDX intermolecular interactions (captured by MD simulations) are important, and compare fairly well with experimental findings.

Conformational entropy of a polymer chain grafted to rough surfaces

A polymer molecule (represented by a statistical chain) end-grafted to a topologically rough surface was studied by static MC simulations. A modified self-avoiding walk on a cubic lattice was used to model the polymer in an athermal solution. Different statistical models of surface roughness were applied. Conformational entropies of chains attached to uncorrelated Gaussian, Brownian, and fractional Brownian surfaces were calculated. Results were compared with the predictions of a simple analytical model of a macromolecule end-grafted to a fractal surface.

A DFT study on the initial stage of thermal degradation of Poly(methyl methacrylate)/carbon nanotube system

DFT calculations, with VWN exchange correlation functional and double numeric basis set, were used to evaluate the energies required for the scission reactions taking place in the initial stage of the thermal degradation of Poly(methyl methacrylate) (PMMA) in the presence of a carbon nanotube (CNT). Side group and main chain scissions were investigated. The results averaged from five configurations of pure PMMA (DP?=?5) were used as references and compared to the results obtained for the five same configurations of PMMA grafted on three carbon nanotubes of similar diameter (1.49 nm). The bond dissociation energies (BDE) of main chain scission evaluated for grafted PMMA was 4 % less endothermic than for pure PMMA. These results seemed independent of the tested chirality (11,11); (12,10) and (16,5) of the carbon nanotubes. Comparisons with the BDE of the weakest bonds due to intrinsic defaults (head to head and unsaturated end chain) were performed.

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.

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.

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.

Nickel/zinc-catalyzed decarbonylative addition of anhydrides to alkynes: A DFT study

Density functional theory (DFT) was used to investigate the nickel- or nickel(0)/zinc- catalyzed decarbonylative addition of phthalic anhydrides to alkynes. All intermediates and transition states were optimized completely at the B3LYP/6-31+G(d,p) level. Calculated results indicated that the decarbonylative addition of phthalic anhydrides to alkynes was exergonic, and the total free energy released was ?87.6 kJ mol^?1. In the five-coordinated complexes M4a and M4b, the insertion reaction of alkynes into the Ni-C bond occurred prior to that into the Ni-O bond. The nickel(0)/zinc-catalyzed decarbonylative addition was much more dominant than the nickel-catalyzed one in whole catalytic decarbonylative addition. The reaction channel CA→M1'→T1'→M2'→T2'→M3a'→M4a'→T3a1'→M5a1' →T4a1'→M6a'→P was the most favorable among all reaction pathways of the nickel- or nickel(0)/zinc- catalyzed decarbonylative addition of phthalic anhydrides to alkynes. And the alkyne insertion reaction was the rate-determining step for this channel. The additive ZnCl₂ had a significant effect, and it might change greatly the electron and geometry structures of those intermediates and transition states. On the whole, the solvent effect decreased the free energy barriers.

Fluorescent sensors based on BODIPY derivatives for aluminium ion recognition: an experimental and theoretical study

Two BODIPY derivative sensors for metal ion recognition containing 10-(4-hydroxyphenyl) (L1) and 10-(3,4-dihydroxyphenyl) (L2) were synthesized in a one-pot reaction of benzaldehyde derivative and 2,4-dimethylpyrrole in the presence of trifluoroacetic acid as catalyst. The binding abilities between these sensors and 50 equivalents of Na⁺, K⁺, Ag⁺, Ca²⁺, Fe²⁺, Co²⁺, Ni²⁺, Cu²⁺, Zn²⁺, Cd²⁺, Pb²⁺, Al³⁺ and Cr³⁺ ions were studied using UV–vis and fluorescent spectroscopic methods. Of all the metal ions tested, Al³⁺ ion showed the greatest decrease in intensity in the spectra of the sensors, and therefore Al³⁺ ion forms the strongest complex. The binding abilities of BODIPY receptors with Na⁺, Ag⁺, Ca²⁺, Co²⁺, Ni²⁺, Cu²⁺, Zn²⁺ and Al³⁺ ions were also investigated using density functional theory (DFT) calculations at B3LYP/LanL2DZ theoretical level. The calculated results point to the same conclusion. DFT calculations also provided the HOMO–LUMO energy levels, which can explain the spectrum change upon complexation.