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.

Microsolvation of Mg2+, Ca2+: strong influence of formal charges in hydrogen bond networks

A stochastic exploration of the quantum conformational spaces in the microsolvation of divalent cations with explicit consideration of up to six solvent molecules [Mg (H ₂ O) _n )]²⁺, (n?=?3, 4, 5, 6) at the B3LYP, MP2, CCSD(T) levels is presented. We find several cases in which the formal charge in Mg²⁺ causes dissociation of water molecules in the first solvation shell, leaving a hydroxide ion available to interact with the central cation, the released proton being transferred to outer solvation shells in a Grotthus type mechanism; this particular finding sheds light on the capacity of Mg²⁺ to promote formation of hydroxide anions, a process necessary to regulate proton transfer in enzymes with exonuclease activity. Two distinct types of hydrogen bonds, scattered over a wide range of distances (1.35–2.15 Å) were identified. We find that in inner solvation shells, where hydrogen bond networks are severely disturbed, most of the interaction energies come from electrostatic and polarization+charge transfer, while in outer solvation shells the situation approximates that of pure water clusters.

Calixarene building block bis(2-hydroxyphenyl)methane (2HDPM) and hydrogen-bonded 2HDPM-H2O complex in electronic excited state

Intramolecular and intermolecular hydrogen bonding in electronic excited states of calixarene building blocks bis(2-hydroxyphenyl)methane (2HDPM) monomer and hydrogen-bonded 2HDPM-H₂O complex were studied theoretically using the time-dependent density functional theory (TDDFT). Twenty-four stable conformations (12 pairs of enantiomers) of 2HDPM monomer have been found in the ground state. From the calculation results, the conformations 1a and 1b which both have an intramolecular hydrogen bond are the most stable ones. The infrared spectra of 2HDPM monomer and 2HDPM-H₂O complex in ground state and S₁ state were calculated. The stretching vibrational absorption band of O₂???H₃ group in the monomer and complex disappeared in the S₁ state. At the same time, a new strong absorption band appeared at the C=O stretching region. From the calculation of bond lengths, it indicates that the O₂???H₃ bond is significantly lengthened in the S₁ state. However, the C₁???O₂ bond is drastically shortened upon electronic excitation to the S₁ state and has the characteristics of C=O band. Furthermore, the intramolecular hydrogen bond O₂???H₃?·?·?·?O₄ of the 2HDPM monomer and the intermolecular hydrogen bonds O₂???H₃?·?·?·?O₇ and O₇???H₉?·?·?·?O₄ of 2HDPM-H₂O complex are all shortened and strengthened in the S₁ state.

Conformation-dependent conductance through a molecular break junction

Ab initio molecular dynamics simulations have been performed of a gold—1,4-benzenedithiol (BDT)—gold nanojunction under mechanical stress. For three different pulling rates between 10 and 40 m s^-1, it is found that the nanowire always ruptures between the second and third Au atom from the thiol sulfur. Larger rupture forces and longer extensions are required at higher pulling rates and vice versa. The electrical conductance was calculated along a pulling trajectory using the DFT-NEGF method to study the effect of thermal and stress-induced structural changes on the electrical transport properties. While the mechanically induced stretching of the junction is seen to lower the time-averaged conductance, thermal conformational changes are capable of altering the conductance by one order of magnitude. No single geometric quantity could be identified as the main contributor to the conductance fluctuations. Small modulations, however, can be explained in terms of C=C double bond vibrations in the BDT molecule. The dependence of the conductance on different geometric variables has further been investigated systematically by performing constrained geometry optimizations along a number of angle and dihedral coordinates. The largest changes in the conductance are observed when the Au-S-C angle and the Au-S-C-C dihedral are simultaneously constrained.

Deamination features of 5-hydroxymethylcytosine,a radical and enzymatic DNA oxidation product

Geometric consequences of electron delocalization were studied for all possible adenine tautomers in aqueous solution by means of ab initio methods {PCM(water)//DFT(B3LYP)/6-311+G(d,p)} and compared to those in the gas phase {DFT(B3LYP)/6-311+G(d,p)}. To measure the consequences of any type of resonance conjugation (π-π, n-π, and σ-π), the geometry-based harmonic oscillator model of electron delocalization (HOMED) index, recently extended to the isolated (DFT) and hydrated (PCM//DFT) molecules, was applied to the molecular fragments (imidazole, pyrimidine, 4-aminopyrimidine, and purine) and also to the whole tautomeric system. For individual tautomers, the resonance conjugations and consequently the bond lengths strongly depend on the position of the labile protons. The HOMED indices are larger for tautomers (or their fragments) possessing the labile proton(s) at the N rather than C atom. Solvent interactions with adenine tautomers slightly increase the resonance conjugations. Consequently, they slightly shorten the single bonds and lengthen the double bonds. When going from the gas phase to water solution, the HOMED indices increase (by less than 0.15 units). There is a good relation between the HOMED indices estimated in water solution and those in the gas phase for the neutral and ionized forms of adenine. Subtle effects, being a consequence of intramolecular interactions between the neighboring groups, are so strongly reduced by solvent that the relation between the HOMED indices and the relative energies for the neutral adenine tautomers seems to be better in water solution than in the gas phase.

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.

Theoretical studies on the tautomerism of tetrazole selenone

The tautomerism of all possible forms of tetrazole selenone (A–G), induced by proton transfer, was studied, theoretically, in different environments including gas phase, continuum solvent and microsolvated environment with one or two explicit water or ammonia molecules. The calculations were performed using two different levels of theory including mPW2PLYP and DFT-B3LYP. The 6-311++G(d,p) basis set was used for C, H, O and N and the standard relativistic effective core pseudo potential LANL2DZ basis set was used for Se atom. It was found that the tetrazole selenone, in the form of A, is the most stable isomer in all of the environments considered in this work. The kinetics of proton transfer reaction was studied in both gas and solvent environments and it was concluded that the activation energy of the reaction increases with going from the gas phase to polar solvents. Moreover, the proton transfer reaction assisted by one or two water or ammonia molecules was investigated and it was found that the activation energy significantly reduces.     

Cooperativity between fluorine-centered halogen bonds: investigation of substituent effects

This article analyzes the substitution effects on cooperativity between fluorin-centered halogen bonds in NCF?·?·?·?NCF?·?·?·?NCX and CNF?·?·?·?CNF?·?·?·?CNX complexes, where X?=?H, F, Cl, CN, OH, and NH₂. These effects are investigated theoretically in terms of geometric and energetic features of the complexes, which are computed by ab initio methods. The topological analysis, based on the quantum theory of atoms in molecules (QTAIM), is used to characterize the interactions and analyze their enhancement with varying electron density at bond critical points. It is found that the complexes with electron-donating groups exhibit a strong cooperativity, while a much weaker cooperativity occurs in the NCF?·?·?·?NCF?·?·?·?NCCN and CNF?·?·?·?CNF?·?·?·?CNCN trimers. An excellent correlation is found between the cooperative energy in the ternary complexes and the calculated three-body interaction energies. The energy decomposition analysis (EDA) indicates that the electrostatic and dispersion effects play a main role in the cooperativity of fluorine-centered halogen bonding.

Sequence selectivity of azinomycin B in DNA alkylation and cross-linking: a QM/MM study

Azinomycin B—a well-known antitumor drug—forms cross-links with DNA through alkylation of purine bases and blocks tumor cell growth. This reaction has been modeled using the ONIOM (B3LYP/6-31?+?g(d):UFF) method to understand the mechanism and sequence selectivity. ONIOM results have been checked for reliability by comparing them with full quantum mechanics calculations for selected paths. Calculations reveal that, among the purine bases, guanine is more reactive and is alkylated by aziridine ring through the C10 position, followed by alkylation of the epoxide ring through the C21 position of Azinomycin B. While the mono alkylation is controlled kinetically, bis-alkylation is controlled thermodynamically. Solvent effects were included using polarized-continuum-model calculations and no significant change from gas phase results was observed.

Effects of trimethylaluminium and tetrakis(ethylmethylamino) hafnium in the early stages of the atomic-layer-deposition of aluminum oxide and hafnium oxide on hydroxylated GaN nanoclusters

We calculate the interactions of two atomic layer deposition (ALD) reactants, trimethylaluminium (TMA) and tetrakis(ethylmethylamino) hafnium (TEMAH) with the hydroxylated Ga-face of GaN clusters when aluminum oxide and hafnium oxide, respectively, are being deposited. The GaN clusters are suitable as testbeds for the actual Ga-face on practical GaN nanocrystals of importance not only in electronics but for several other applications in nanotechnology. We find that TMA spontaneously interacts with hydroxylated GaN; however it does not follow the atomic layer deposition reaction path unless there is an excess in potential energy introduced in the clusters at the beginning of the optimization, for instance, using larger bond lengths of various bonds in the initial structures. TEMAH also does not interact with hydroxylated GaN, unless there is an excess in potential energy. The formation of a Ga—N(CH₃)(CH₂CH₃) bond during the ALD of HfO₂ using TEMAH as the reactant without breaking the Hf—N bond could be the key part of the mechanism behind the formation of an interface layer at the HfO₂/GaN interface.

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.

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.

Revealing substitution effects on the strength and nature of halogen-hydride interactions: a theoretical study

A quantum chemistry study was carried out to investigate the strength and nature of halogen bond interactions in HXeH···XCCY complexes, where X = Cl, Br and Y = H, F, Cl, Br, CN, NC, C₂H, CH₃, OH, SH, NH₂. Examination of the electrostatic potentials V(r) of the XCCY molecules reveals that the addition of substituents has a significant effect upon the most positive electrostatic potential on the surface of the interacting halogen atom. We found that the magnitude of atomic charges and multipole moments depends upon the halogen atom X and is rather sensitive to the electron-withdrawing/donating power of the remainder of the molecule. An excellent correlation was found between the most positive electrostatic potentials on the halogen atom and the interaction energies. For either HXeH···ClCCY or HXeH···BrCCY complexes, an approximate linear correlation between the interaction energies and halogens multipole moments are established, indicating that the electrostatic and polarization interactions are responsible for the stability of the complexes. According to energy decomposition analysis, it is revealed that the electrostatic interactions are the major source of the attraction in the HXeH···XCCY complexes. Furthermore, the changes in the electrostatic term are mainly responsible for the dependence of interaction energy on the halogen atom.

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.

Competition between halogen, dihalogen and hydrogen bonds in bromo- and iodomethanol dimers

O-H…X and O-H…O H-bonds as well as C-X…X dihalogen and C-X…O halogen bonds have been investigated in halomethanol dimers (bromomethanol dimer, iodomethanol dimer, difluorobromomethanol…bromomethanol complex and difluoroiodomethanol…iodomethanol complex). Structures of all complexes were optimized at the counterpoise-corrected MP2/cc-pVTZ level and single-point energies were calculated at the CCSD(T)/aug-cc-pVTZ level. Energy decomposition for the bromomethanol dimer complex was performed using the DFT-SAPT method based on the aug-cc-pVTZ basis set. OH…O and OH…X H-bonds are systematically the strongest in all complexes investigated, with the former being the strongest bond. Halogen and dihalogen bonds, being of comparable strength, are weaker than both H-bonds but are still significant. The strongest bonds were found in the difluoroiodomethanol…iodomethanol complex, where the O-H…O H-bond exceeds 7 kcal mol^-1, and the halogen and dihalogen bonds exceed 2.5 and 2.3 kcal mol^-1, respectively. Electrostatic energy is dominant for H-bonded structures, in halogen bonded structures electrostatic and dispersion energies are comparable, and, finally, for dihalogen structures the dispersion energy is clearly dominant.

Quantitative characterization of protein tertiary motifs

Dimerization of 2-naphthalenecarbonitrile (2-NpCN) mediated by cucurbit[8]uril (CB[8]) has been investigated employing the density functional theory. Different structures of 2-NpCN dimers were generated by combining monomers in anti-head-to-head (A), anti-head-to-tail (B) and syn-head-to-tail (C) fashion. All these dimeric structures possess rigid cube-like architecture. On confinement within the CB[8] dimer A turns out to be the lowest energy structure. Calculated ¹H NMR spectra revealed that the 2-NpCN dimer exhibits large shielding for aromatic protons consistent with the experiment. The protons attached to cubane moiety on the other hand, led to down-field signals. Dimerization mediated with CB[8] cavitand is further accompanied by the frequency up-shift (blue shift) of methylene stretching vibration in its infrared spectra.

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.

N(4)-Tolyl-2-acetylpyridine thiosemicarbazones and their platinum(II,IV) and gold(III) complexes: cytotoxicity against human glioma cells and studies on the mode of action

Complexes [Au(2Ac4oT)Cl][AuCl₂] (1), [Au(H_py2Ac4mT)Cl₂]Cl·H₂O (2), [Au(H_py2Ac4pT)Cl₂]Cl (3), [Pt(H2Ac4oT)Cl]Cl (4), [Pt(2Ac4mT)Cl]·H₂O (5), [Pt(2Ac4pT)Cl] (6) and [Pt(L)Cl₂OH], L = 2Ac4mT (7), 2Ac4oT (8), 2Ac4pT (9) were prepared with N(4)-ortho- (H2Ac4oT), N(4)-meta- (H2Ac4mT) and N(4)-para- (H2Ac4pT) tolyl-2-acetylpyridine thiosemicarbazone. The cytotoxic activities of all compounds were assayed against U-87 and T-98 human malignant glioma cell lines. Upon coordination cytotoxicity improved in 2, 5 and 8. In general, the gold(III) complexes were more cytotoxic than those with platinum(II,IV). Several of these compounds proved to be more active than cisplatin and auranofin used as controls. The gold(III) complexes probably act by inhibiting the activity of thioredoxin reductase enzyme whereas the mode of action of the platinum(II,IV) complexes involves binding to DNA. Cells treated with the studied compounds presented morphological changes such as cell shrinkage and blebs formation, which indicate cell death by apoptosis induction.     

Computational investigation on microsolvation of the osmolyte glycine betaine [GB (H2O)1-7]

The preferential interactions of glycine betaine (GB) with solvent components and the effect of solvent on its stability have been examined. In particular, the microsolvation of organic osmolyte and widely important osmoprotectant in nature as glycine betaine has been reported by using M06 method. A number of configurations (b_X (a-z)) of the clusters for one to seven water molecules (×?=?1-7) have been considered for the microsolvation. Structures of stable conformers are obtained and denoted as b1a, b2a, b3a, b4a, b5a, b6a and b7a. It is observed from the interaction energy difference (?E) that only seven water molecules can be accommodated in the first solvation shell to stabilize GB. It is also observed that the calculated relative energy using M06 is in close agreement with calculations at the MP2 level of theory.