Natural bond orbital approach to the transmission of substituent effect through the fulvene and benzene ring systems

Electronic structure of 22 monosubstituted derivatives of benzene and exocyclically substituted fulvene with substituents: B(OH)₂, BH₂, CCH, CF₃, CH₃, CHCH₂, CHO, Cl, CMe₃, CN, COCH₃, CONH₂, COOH, F, NH₂, NMe₂, NO, NO₂, OCH₃, OH, SiH₃, SiMe₃ were studied theoretically by means of Natural Bond Orbital analysis. It is shown, that sum of π-electron population of carbon atoms of the fulvene and benzene rings, pEDA(F) and pEDA(B), respectively correlate well with Hammett substituent constants s_p⁺ \sigma_{\rm{p}}^{+} and aromaticity index NICS. The substituent effect acting on pi-electron occupation at carbon atoms of the fulvene ring is significantly stronger than in the case of benzene. Electron occupations of ring carbon atoms (except C1) in fulvene plotted against each other give linear regressions with high correlation coefficients. The same is true for ortho- and para-carbon atoms in benzene. Positive slopes of the regressions indicate similar for fulvene and benzene kind of substituent effect – mostly resonance in nature. Only the regressions of occupation at the carbon atom in meta- position of benzene against ortho- and para-positions gives negative slopes and low correlation coefficients.     

Modeling the effect of H-bonding interactions and molecular packing on the molecular structure of [Ag(ethylnicotinate)2]NO3 complex

The gas phase molecular structure of a single isolated molecule of [Ag(Etnic)₂NO₃];1 where Etnic = Ethylnicotinate was calculated using B3LYP method. The H-bonding interaction between 1 with one (complex 2) and two (complex 3) water molecules together with the dimeric formula [Ag(Etnic)₂NO₃]₂;4 and the tetrameric formula [Ag(Etnic)₂NO₃]₄;5 were calculated using the same level of theory to model the effect of intermolecular interactions and molecular packing on the molecular structure of the titled complex. The H-bond dissociation energies of complexes 2 and 3 were calculated to be in the range of 12.220–14.253 and 30.106–31.055 kcal?mol^?1, respectively, indicating the formation of relatively strong H-bonds between 1 and water molecules. The calculations predict bidentate nitrate ligand in the case of 1 and 2, leading to distorted tetrahedral geometry around the silver ion with longer Ag–O distances in case of 2 compared to 1, while 3 has a unidentate nitrate ligand leading to a distorted trigonal planar geometry. The packing of two [Ag(Etnic)₂NO₃] complex units; 4 does not affect the molecular geometry around Ag(I) ion compared to 1. In the case of 5, the two asymmetric units of the formula [Ag(Etnic)₂NO₃] differ in the bonding mode of the nitrate group, where the geometry around the silver ion is distorted tetrahedral in one unit and trigonal planar in the other. The calculations predicted almost no change in the charge densities at the different atomic sites except at the sites involved in the C–H?O interactions as well as at the coordinated nitrogen of the pyridine ring.

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.

A comparative theoretical investigation into the strength of the trigger-bond in the Na+, Mg2+ and HF complexes involving the nitro group of R–NO2 (R = –CH3, –NH2 and –OCH3) or the C = C bond of (E)-O2N–CH = CH–NO2

A comparative theoretical investigation into the change in strength of the trigger-bond upon formation of the Na⁺, Mg²⁺ and HF complexes involving the nitro group of RNO₂ (R?=? –CH₃, –NH₂, –OCH₃) or the C?=?C bond of (E)-O₂N–CH?=?CH–NO₂ was carried out using the B3LYP and MP2(full) methods with the 6-311++G**, 6-311++G(2df,2p) and aug-cc-pVTZ basis sets. Except for the Mg²⁺?π system with (E)-O₂N–CH?=?CH–NO₂ (i.e., C₂H₂N₂O₄?Mg²⁺), the strength of the trigger-bond X–NO₂ (X?=?C, N or O) was enhanced upon complex formation. Furthermore, the increment of bond dissociation energy of the X–NO₂ bond in the Na⁺ complex was far greater than that in the corresponding HF system. Thus, the explosive sensitivity in the former might be lower than that in the latter. For C₂H₂N₂O₄?Mg²⁺, the explosive sensitivity might also be reduced. Therefore, it is possible that introducing cations into the structure of explosives might be more efficacious at reducing explosive sensitivity than the formation of an intermolecular hydrogen-bonded complex. AIM, NBO and electron density shifts analyses showed that the electron density shifted toward the X–NO₂ bond upon complex formation, leading to a strengthened X–NO₂ bond and possibly reduced explosive sensitivity.

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.

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.

Insights into the strength and nature of carbene···halogen bond interactions: a theoretical perspective

Halogen-bonding, a noncovalent interaction between a halogen atom X in one molecule and a negative site in another, plays critical roles in fields as diverse as molecular biology, drug design and material engineering. In this work, we have examined the strength and origin of halogen bonds between carbene CH₂ and XCCY molecules, where X?=?Cl, Br, I, and Y?=?H, F, COF, COOH, CF₃, NO₂, CN, NH₂, CH₃, OH. These calculations have been carried out using M06-2X, MP2 and CCSD(T) methods, through analyses of surface electrostatic potentials V _S(r) and intermolecular interaction energies. Not surprisingly, the strength of the halogen bonds in the CH₂···XCCY complexes depend on the polarizability of the halogen X and the electron-withdrawing power of the Y group. It is revealed that for a given carbene···X interaction, the electrostatic term is slightly larger (i.e., more negative) than the dispersion term. Comparing the data for the chlorine, bromine and iodine substituted CH₂···XCCY systems, it can be seen that both the polarization and dispersion components of the interaction energy increase with increasing halogen size. One can see that increasing the size and positive nature of a halogen’s σ-hole markedly enhances the electrostatic contribution of the halogen-bonding interaction.

Penta- and heteropentadienyl ligands coordinated to beryllium

In this work we have performed a systematic study of new organometallic complexes containing penta- and heteropentadienyl (CH₂CHCHCHX, X?=?CH₂, O, NH, S) ligands coordinated to beryllium. Calculated complexes were studied using the density functional theory (PBE) in combination with the 6-311++G(3d,2p) basis set. The coordination number on the beryllium atom varies according to the type of ligand. Pentadienyl ligand shows hapticities η¹ and η⁵, while heteropentadienyl ligands display η¹ and η² hapticities. A Wiberg bond indices study was performed in order to get information about their bond orders.

Theoretical investigation into optical and electronic properties of 1,8-naphthalimide derivatives

A series of 1,8-naphthalimide derivatives has been designed to explore their optical, electronic, and charge transport properties as charge transport and/or luminescent materials for organic light-emitting diodes (OLEDs). The frontier molecular orbitals (FMOs) analysis have shown that the vertical electronic transitions of absorption and emission are characterized as intramolecular charge transfer (ICT) for electron-donating and aromatic groups substituted derivatives. However, the ICT character of the electron-withdrawing substituted derivatives is not significant. The calculated results show that their optical and electronic properties are affected by the substituent groups in 4-position of 1,8-naphthalimide. Our results suggest that 1,8-naphthalimide derivatives with electron-donating ?OCH₃ and ?N(CH₃)₂ (1 and 2), electron-withdrawing ?CN and?COCH₃ (3 and 4), 2-(thiophen-2-yl)thiophene (5), 2,3-dihydrothieno[3,4-b][1, 4]dioxine (6), 2-phenyl-1,3,4-oxadiazole (7), and benzo[c][1,2,5]thiadiazole (8) fragments are expected to be promising candidates for luminescent materials for OLEDs, particularly for 5 and 7. In addition, 3 and 7 can be used as promising hole transport materials for OLEDs. This study should be helpful in further theoretical investigations on such kind of systems and also to the experimental study for charge transport and/or luminescent materials for OLEDs.

Interaction between NO 3 − and abscisic acid in the regulation of lateral root elongation in Zea mays L.

Lateral root (LR) elongation rate of 7–8-day maize seedlings depends on the availability of NO ₃ ^? , NO ₂ ^? , and abscisic acid (ABA) in an environment. Four-hour exposure to 0.01–1.5 mM NO ₂ ^? increases the relative LR elongation rate; in the case of NO ₂ ^? , the stimulation occurs only at an NO ₂ ^? concentration equal to 0.01 mM. Exogenous ABA (10^?6 M) inhibits the LR elongation process. In the case of a combined influence of NO ₃ ^? and ABA or NO ₂ ^? and ABA, the character of the response elongation reaction is different. The NO role in the regulation of LR elongation is discussed.     

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 investigations on the synthesis mechanism of cyanuric acid from NH3 and CO2

In the synthesis of cyanuric acid from NH₃ and CO₂, urea and isocyanic acid OCNH are two pivotal intermediates. Based on density functional theory (DFT) calculations, the synthesis mechanism of cyanuric acid from NH₃ + CO₂ was investigated systematically. Urea can be synthesized from NH₃ and CO₂, and cyanuric acid can be obtained from urea or NH₃ + CO₂. In the stepwise mechanism of cyanuric acid from urea or NH₃ + CO₂, the energy barriers are relatively high, and the condition of high pressure and temperature does not decrease the energy barriers. Our theoretical model shows that cyanuric acid is actually acquired from OCNH via a one-step cycloaddition reaction.

Modeling the scavenging activity of ellagic acid and its methyl derivatives towards hydroxyl, methoxy, and nitrogen dioxide radicals

The reaction mechanisms involved in the scavenging of hydroxyl (OH^·), methoxy (OCH₃ ^·), and nitrogen dioxide (NO₂ ^·) radicals by ellagic acid and its monomethyl and dimethyl derivatives were investigated using the transition state theory and density functional theory. The calculated Gibbs barrier energies associated with the abstraction of hydrogen from the hydroxyl groups of ellagic acid and its monomethyl and dimethyl derivatives by an OH· radical in aqueous media were all found to be negative. When NO₂ ^· was the radical involved in hydrogen abstraction, the Gibbs barrier energies were much larger than those calculated when the OH^· radical was involved. When OCH₃ ^· was the hydrogen-abstracting radical, the Gibbs barrier energies lay between those obtained with OH^· and NO₂ ^· radicals. Therefore, the scavenging efficiencies of ellagic acid and its monomethyl and dimethyl derivatives towards the three radicals decrease in the order OH^· >> OCH₃ ^· > NO₂ ^·. Our calculated rate constants are broadly in agreement with those obtained experimentally for hydrogen abstraction reactions of ellagic acid with OH· and NO₂· radicals.

Theoretical study of the triangular bonding complex formed by carbon tetrabromide, a halide, and a solvent molecule in the gas phase

MP2(full)/aug-cc-pVDZ(-PP) computations predict that new triangular bonding complexes (where X^? is a halide and H–C refers to a protic solvent molecule) consist of one halogen bond and two hydrogen bonds in the gas phase. Carbon tetrabromide acts as the donor in the halogen bond, while it acts as an acceptor in the hydrogen bond. The halide (which commonly acts as an acceptor) can interact with both carbon tetrabromide and solvent molecule (CH₃CN, CH₂Cl₂, CHCl₃) to form a halogen bond and a hydrogen bond, respectively. The strength of the halogen bond obeys the order CBr₄???Cl^? > CBr₄???Br^? > CBr₄???I^?. For the hydrogen bonds formed between various halides and the same solvent molecule, the strength of the hydrogen bond obeys the order C-H???Cl^? > C-H???Br^? > C-H???I^?. For the hydrogen bonds formed between the same halide and various solvent molecules, the interaction strength is proportional to the acidity of the hydrogen in the solvent molecule. The diminutive effect is present between the hydrogen bonds and the halogen bond in chlorine and bromine triangular bonding complexes. Complexes containing iodide ion show weak cooperative effects.

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.

Probing mechanism of metal catalyzed hydrolysis of Thymidylyl (3′-O, 5′-S) thymidine phosphodiester derivatives

Hydrolysis of nucleic acids is of fundamental importance in biological sciences. Kinetic and theoretical studies on different substrates wherein the phosphodiester bond combined with alkyl or aryl groups and sugar moiety have been the focus of attention in recent literature. The present work focuses on understanding the mechanism and energetics of alkali metal (Li, Na, and K) catalyzed hydrolysis of phosphodiester bond in modeled substrates including Thymidylyl (3′-O, 5′-S) thymidine phosphodiester (Tp-ST) (1), 3′-Thymidylyl (1-trifluoroethyl) phosphodiester (Tp-OCH₂CF₃) (2), 3′-Thymidylyl (o-cholorophenyl) phosphodiester (Tp-OPh(o-Cl)) (3) and 3′-Thymidylyl(p-nitrophenyl) phosphodiester (Tp-OPh(p-NO₂)) (4) employing density functional theory. Theoretical calculations reveal that the reaction follows a single-step (A_ND_N) mechanism where nucleophile attack and leaving group departure take place simultaneously. Activation barrier for potassium catalyzed Tp-ST hydrolysis (12.0 kcal mol^?1) has been nearly twice as large compared to that for hydrolysis incorporating lithium or sodium. Effect of solvent (water) on activation energies has further been analyzed by adding a water molecule to each metal ion of the substrate. It has been shown that activation barrier of phosphodiester hydrolysis correlates well with basicity of leaving group.

Small-conductance chloride channels in human peripheral T lymphocytes

During whole-cell patch-clamp recording from normal (nontransformed) human T lymphocytes a chloride current spontaneously activated in >98% of cells (n > 200) in the absence of applied osmotic or pressure gradients. However, some volume sensitivity was observed, as negative pressure pulses reduced the current. With iso-osmotic bath and pipette solutions the peak amplitude built up (time constant ≈23 sec at room temperature), a variable-duration plateau phase followed, then the current ran down spontaneously (time constant ≈280 sec). The anion permeability sequence, calculated from reversal potentials was I^?, Br^? > NO ₃ ^? , Cl^? > CH₃SO ₃ ^? , HCO ₃ ^? > CH₃COO^? > F^? > aspartate, gluconate, SO ₄ ^2? and there was no measurable monovalent cation permeability. The Cl^? current was independent of time during long voltage steps and there was no evidence of voltage-dependent gating; however, the current showed intrinsic outward rectification in symmetrical Cl^? solutions. The conductance of the channels underlying the whole-cell current was calculated from fluctuation analysis, using power-spectral density and variance-vs.-mean analysis. Both methods yielded a single channel conductance of about 0.6 pS at ?70 mV (close to the normal resting potential of T lymphocytes). The power spectral density function was best fit by the sum of two Lorentzian functions, with corner frequencies of 30 and 295 Hz, corresponding to mean open times of 0.54 and 5.13 msec. The pharmacological profile included rapid block by external application of flufenamic acid (50 μm), 5-nitro-2-(3-phenylpropylamino)-benzoic acid (NPPB, 100 μm, [6,7-dichloro-2-cyclopentyl-2,3-dihydro-2-methyl-1-oxo-1H-inden-5-y1) oxy] acetic acid (IAA-94, 250 μm) or 100 μm 1,9-dideoxyforskolin. The stilbene derivatives DIDS (4,4′-diisothiocyano-2,2′ di-sulphonic acid stilbene, 500 μm) and SITS (4-acetamido-4′-isothiocyano-2, 2′-disulphonic acid stilbene, 500 μm) prevented buildup of Cl^? current after a 30-min preincubation at 500 μm. When tested in a mitogenic assay, DIDS, flufenamic acid, NPPB and IAA-94 all inhibited T-cell proliferation, suggesting a physiological function in addition to the observed volume sensitivity.