A three-layer ONIOM model for the outside binding of cationic porphyrins and nucleotide pair DNA

In this work we investigated the outside binding mode between a cationic porphyrin and a nucleotide pair of DNA, adenine-thymine and guanine-cytosine, in a supramolecular assembly. We used two structural models (semi-extended, extended) that differ in the size of porphyrin, two kinds of theoretical methods: a three layer ONIOM (B3LYP/6-31G(d)/PM3/UFF), and DFT B3LYP/6-31G(d,p), and three cationic porphyrins. ONIOM method was first tested on the semi-extended model that was calculated in four environments: gas phase, solution phase using an explicit solvent model (H₂O), in the presence of a sodium cation (Na⁺) and in both (H₂O + Na⁺). From interaction energy results, we found that the affinity of the cationic substituent by the adenine nucleotide is favored upon the thymine nucleotide. The extended model that considers the whole porphyrin was applied in the gas phase to the four nucleotides. All the cationic porphyrins showed affinity by the nucleotides in the order adenine > guanine > thymine > cytosine. The interaction energy values for outside binding showed a strong porphyrin-nucleotide interaction (≈-90 kcal?mol^-1), that slightly varies between the nucleotides suggesting that this kind of cationic porphyrin has a little selectivity for some of them. We also found that the effect of the nature of the cationic substituent (chain length) in the porphyrin on the outside binding is small (≈2–13 kcal?mol^-1). Coherence between the results showed that ONIOM is a useful tool to get a reasonable molecular geometry to be used as a starting point in calculations of density functional theory.

Key role of hydrazine to the interaction between oxaloacetic against phosphoenolpyruvic carboxykinase (PEPCK): ONIOM calculations

The interactions between oxaloacetic (OAA) and phosphoenolpyruvic carboxykinase (PEPCK) binding pocket in the presence and absence of hydrazine were carried out using quantum chemical calculations, based on the two-layered ONIOM (ONIOM2) approach. The complexes were partially optimized by ONIOM2 (B3LYP/6-31G(d):PM6) method while the interaction energies between OAA and individual residues surrounding the pocket were performed at the MP2/6-31G(d,p) level of theory. The calculated interaction energies (INT) indicated that Arg87, Gly237, Ser286, and Arg405 are key residues for binding to OAA with the INT values of ?1.93, ?2.06, ?2.47, and ?3.16 kcal mol^?1, respectively. The interactions are mainly due to the formation of hydrogen bonding interactions with OAA. Moreover, using ONIOM2 (B3LYP/6-31G(d):PM6) applied on the PEPCKHS complex, two proton transfers were observed; first, the proton was transferred from the carboxylic group of OAA to hydrazine while the second one was from Asp311 to Lys244. Such reactions cause the generation of binding strength of OAA to the pocket via electrostatic interaction. The orientations of Lys243, Lys244, His264, Asp311, Phe333, and Arg405 were greatly deviated after hydrazine incorporation. These indicate that hydrazine plays an important role in terms of not only changing the conformation of the binding pocket, but is also tightly bound to OAA resulting in its conformation change in the pocket. The understanding of such interaction can be useful for the design of hydrazine-based inhibitor for antichachexia agents.

DFT investigations of phosphotriesters hydrolysis in aqueous solution: a model for DNA single strand scission induced by N-nitrosoureas

DNA phosphotriester adducts are common alkylation products of DNA phosphodiester moiety induced by N-nitrosoureas. The 2-hydroxyethyl phosphotriester was reported to hydrolyze more rapidly than other alkyl phosphotriesters both in neutral and in alkaline conditions, which can cause DNA single strand scission. In this work, DFT calculations have been employed to map out the four lowest activation free-energy profiles for neutral and alkaline hydrolysis of triethyl phosphate (TEP) and diethyl 2-hydroxyethyl phosphate (DEHEP). All the hydrolysis pathways were illuminated to be stepwise involving an acyclic or cyclic phosphorane intermediate for TEP or DEHEP, respectively. The rate-limiting step for all the hydrolysis reactions was found to be the formation of phosphorane intermediate, with the exception of DEHEP hydrolysis in alkaline conditions that the decomposition process turned out to be the rate-limiting step, owing to the extraordinary low formation barrier of cyclic phosphorane intermediate catalyzed by hydroxide. The rate-limiting barriers obtained for the four reactions are all consistent with the available experimental information concerning the corresponding hydrolysis reactions of phosphotriesters. Our calculations performed on the phosphate triesters hydrolysis predict that the lower formation barriers of cyclic phosphorane intermediates compared to its acyclic counter-part should be the dominant factor governing the hydrolysis rate enhancement of DEHEP relative to TEP both in neutral and in alkaline conditions.

Degradation of polyvinyl alcohol under mechanothermal stretching

Mechanical and thermal properties of polyvinyl alcohol (PVA) are characterized and analyzed using in situ X-ray photoelectron spectroscopy (XPS) and quantum chemistry calculations. It is found that the carbon peaks—commonly used as the reference for spectroscopic analysis—shift under mechanical and thermal stretching. Results also indicate that, at different temperatures and among the various functional groups present in PVA, the carbon in the C–O group is the most stable. Computational calculations showed that Hartree–Fock/10-31G (d) reproduces the binding energy of core carbon electrons with an accuracy of 95 %, which is enough to characterize bonds, allowing the results of the spectroscopic analysis to be corroborated.

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.

A new interaction mechanism of LiNH2 with MgH2: magnesium bond

Quantum chemical calculations were performed for LiNH₂–HMgX (X?=?H, F, Cl, Br, CH₃, OH, and NH₂) complexes to propose a new interaction mechanism between them. This theoretical survey showed that the complexes are stabilized through the combinative interaction of magnesium and lithium bonds. The binding energies are in the range of 63.2–66.5 kcal mol^?1, i.e., much larger than that of the lithium bond. Upon complexation, both Mg–H and Li–N bonds are lengthened. Substituents increase Mg-H bond elongation and at the same time decrease Li-N bond elongation. These cyclic complexes were characterized with the presence of a ring critical point and natural population analysis charges.

Theoretical study on the binding mechanism between N6-methyladenine and natural DNA bases

N6-methyladenine (m⁶A) is a rare base naturally occurring in DNA. It is different from the base adenine due to its N-CH₃. Therefore, the base not only pairs with thymine, but also with other DNA bases (cytosine, adenine and guanine). In this work, Møller-Plesset second-order (MP2) method has been used to investigate the binding mechanism between m⁶A and natural DNA bases in gas phase and in aqueous solution. The results show that N-CH₃ changed the way of N6-methyladenine binding to natural DNA bases. The binding style significantly influences the stability of base pairs. The trans-m⁶A:G and trans-m⁶A:C conformers are the most stable among all the base pairs. The existence of solvent can remarkably reduce the stability of the base pairs, and the DNA bases prefer pairing with trans-m⁶A to cis-m⁶A. Besides, the properties of these hydrogen bonds have been analyzed by atom in molecules (AIM) theory, natural bond orbital (NBO) analysis and Wiberg bond indexes (WBI). In addition, pairing with m⁶A decreases the binding energies compared to the normal Watson-Crick base pairs, it may explain the instability of the N6 site methylated DNA in theory.

Entropy versus aromaticity in the conformational dynamics of aromatic rings

Comparison of the results of Car-Parrinello molecular dynamics simulations of isolated benzene, pyrimidine and 1,2,4-triazine molecules reveals that the unusually low population of planar geometry of the benzene ring is caused by entropy effects despite its high aromaticity. The decrease in symmetry of the molecule results in smaller changes in entropy and Gibbs free energy due to out-of-plane deformations of the ring, leading to an increase in the population of planar geometry of the ring. This leads to differences in the topology of potential energy and Gibbs free energy surfaces.

Insight on the interaction of polychlorobiphenyl with nucleic acid–base

The interaction between one polychlorobiphenyl (3,3′,4,4′,-tetrachlorobiphenyl, coded PCB77) and the four DNA nucleic acid–base is studied by means of quantum mechanics calculations in stacked conformations. It is shown that even if the intermolecular dispersion energy is the largest component of the total interaction energy, some other contributions play a non negligible role. In particular the electrostatic dipole-dipole interaction and the charge transfer from the nucleobase to the PCB are responsible for the relative orientation of the monomers in the complexes. In addition, the charge transfer tends to flatten the PCB, which could therefore intercalate more easily between DNA base pairs. From these seminal results, we predict that PCB could intercalate completely between two base pairs, preferably between Guanine:Cytosine pairs.

Dependence of the optical absorption and Na+ binding energies of coumarin-crown ethers on the size and attachment position of ether ring: density functional investigation

The crowned coumarin complexes are well known compounds for their ion recognition abilities. They undergo photophysical changes upon cation binding. On the basis of density functional theory calculations, we examined the sodium cation (Na⁺) binding energies of coumarin-crown ethers based on 15-Crown-5 (15 C5) and 18-Crown-6 (18 C6) as well as the optical absorptions of coumarin-crown ethers based on 12-Crown-4 (12 C4), 15 C5 and 18 C6. We explored why the attachment of crown ether ring to coumarin affects the Na⁺ binding energies of coumarin-crown ethers and also why the optical absorption of coumarin is modified by the crown ethers. Our study reveals that the Na⁺ ion binding energies of coumarin-crown ethers depend strongly on the size of the crown ether ring and also on the attachment position of the ether ring on coumarin. These factors affect the intramolecular charge transfer and overall stability of the complexes. The absorptions of the coumarin and ether ring parts of coumarin-crown ether are red shifted from those of isolated coumarin and crown ether, respectively. The red-shift of the coumarin ester group absorption is much stronger depending on the attachment position of the ether ring to coumarin. The absorption intensity of the coumarin part in coumarin-crown ethers is reduced for the benzene group absorption, but is enhanced for the ester group absorption.

Electric field effect on the zigzag (6,0) single-wall BC2N nanotube for use in nano-electronic circuits

We have analyzed the effect of external electric field on the zigzag (6,0) single-wall BC₂N nanotube using density functional theory calculations. Analysis of the structural parameters indicates that the nanotube is resistant against the external electric field strengths. Analysis of the electronic structure of the nanotube indicates that the applied parallel electric field strengths have a much stronger interaction with the nanotube with respect to the transverse electric field strengths and the nanotube is easier to modulate by the applied parallel electric field. Our results show that the properties of the nanotube can be controlled by the proper external electric field for use in nano-electronic circuits.

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.

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.

Thermal, electronic and ductile properties of lead-chalcogenides under pressure

Fully relativistic pseudo-potential ab-initio calculations have been performed to investigate the high pressure phase transition, elastic and electronic properties of lead-chalcogenides including the less known lead polonium. The calculated ground state parameters, for the rock-salt structure show good agreement with the experimental data. PbS, PbSe, PbTe and PbPo undergo a first-order phase transition from rock-salt to CsCl structure at 19.4, 15.5, 11.5 and 7.3 GPa, respectively. The elastic properties have also been calculated. The calculations successfully predicted the location of the band gap at L-point of Brillouin zone and the band gap for each material at ambient pressure. It is observed that unlike other lead-chalcogenides, PbPo is semi-metal at ambient pressure. The pressure variation of the energy gap indicates that these materials metalize under pressure. The electronic structures of these materials have been computed in parent as well as in high pressure B2 phase.

Cytotoxic effect and molecular docking of 4-ethoxycarbonylmethyl-1-(piperidin-4-ylcarbonyl)-thiosemicarbazide—a novel topoisomerase II inhibitor

The preliminary cytotoxic effect of 4-ethoxycarbonylmethyl-1-(piperidin-4-ylcarbonyl)-thiosemicarbazide hydrochloride (1)—a potent topoisomerase II inhibitor—was measured using a MTT assay. It was found that the compound decreased the number of viable cells in both estrogen receptor-positive MCF-7 and estrogen receptor-negative MDA-MB-231breast cancer cells, with IC₅₀ values of 146?±?2 and 132?±?2 μM, respectively. To clarify the molecular basis of the inhibitory action of 1, molecular docking studies were carried out. The results suggest that 1 targets the ATP binding pocket.

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.

Non-covalent interactions – QTAIM and NBO analysis

MP2(full)/6-311++G(3df,3pd) calculations were carried out on complexes linked through various non-covalent Lewis acid – Lewis base interactions. These are: hydrogen bond, dihydrogen bond, hydride bond and halogen bond. The quantum theory of ´atoms in molecules´ (QTAIM) as well as the natural bond orbitals (NBO) method were applied to analyze properties of these interactions. It was found that for the A-H…B hydrogen bond as well as for the A-X…B halogen bond (X designates halogen) the complex formation leads to the increase of s-character in the A-atom hybrid orbital aimed toward the H or X atom. In opposite, for the A…H-B hydride bond, where the H-atom possesses negative charge, the decrease of s-character in the B-atom orbital is observed. All these changes connected with the redistribution of the electron charge being the effect of the complex formation are in line with Bent´s rule. The numerous correlations between energetic, geometrical, NBO and QTAIM parameters were also found.

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.

Influence of quantum dots on the aromaticity of thiosalicylic acid

When ligands are coordinated to quantum dots (QDs), the ring current of the ligand strongly influences the applications of the QDs, for example in solar cell technology. The Raman spectrum of the ligand can be used to probe and identify ions or measure ion concentrations. Here, we investigated, using a theoretical method, the aromaticities and Raman spectra of CdTe, CdSe, and CdS QDs coordinated with thiosalicylic acid ligands. We found that the aromaticity of the benzene ring in free thiosalicylic acid increased when it was used as a QD ligand. The ring currents of the benzene rings in the CdTe–ligand, CdSe–ligand, and CdS–ligand systems were stronger than the ring current of the benzene ring in free thiosalicylic acid; in other words, the QDs influence the ring current—they enhance the electron transfer rate of the benzene ring. We also discovered that the CdTe–ligand and CdSe–ligand systems have stronger ring currents than the CdS–ligand system. The high electronegativity and vacant d orbital of the sulfur atom influence the ring current of the ligand in the CdS–ligand system. Further, the Raman spectrum of free thiosalicylic acid was different from the spectra of the ligands in the QD–ligand systems: the Raman spectra of COO^? in each QD–ligand system was enhanced compared with that of the COO^? in free thiosalicylic acid.