G2, G3, and complete basis set calculations of the thermodynamic properties of cis- and trans-triazene

Following our recent study on triazane, we present a follow-up study on the thermodynamic properties of triazane’s unsaturated analog, triazene. We predict optimized structural parameters, vibrational frequencies, enthalpies of formation, enthalpies of combustion, specific enthalpies of combustion, and proton affinities. Our results indicate that the cis form of triazene has a specific enthalpy of combustion of −15.2 kJ g⁻¹ and the trans form has a specific enthalpy of combustion of −14.7 kJ g⁻¹. Figure Structures of cis- and trans-triazane, N₃H₃     

G2, G3, and complete basis set calculations on the thermodynamic properties of triazane

As a follow-up study to our study on tetrazane (N₄H₆), we present computed thermodynamic properties of triazane (N₃H₅). Calculated properties include optimized geometries, infrared vibrations, enthalpy of formation, enthalpy of combustion, and proton affinities. We have also mapped the potential energy surface as the molecule is rotated about the N-N bond. We have predicted a specific enthalpy of combustion for triazane of about -20 kJ g⁻¹. Figure Schematic diagram of the dielectric barrier discharge (left) and typical temporal profiles of voltage and current, as obtained from the simulations (right)     

Theoretical study of [ XN5 ]− (X=O, S, Se, Te) systems

A series of [XN₅]⁻ (X=O, S, Se, Te) compounds has been examined with ab initio and Density Functional Theory (DFT) methods. The five-membered nitrogen ring series of structures are global minima and may exist or be characterized due to their significant dissociation barriers (29.7–32.7 kcal mol⁻¹). Nucleus-independent chemical shifts (NICS) criteria and the presence of (4n+2) π-electrons confirmed that the five-membered nitrogen ring in their structures exhibits characteristics of aromaticity. Thus, the strong stability of the five-membered nitrogen ring structures may be attributed partially to their aromaticity.      

Theoretical 49Ti NMR chemical shifts

⁴⁹Ti chemical shifts for a total of 20 titanium complexes are reported, and several levels of theory are evaluated in order to identify a reliable approach for the calculation of titanium NMR data. The popular B3LYP/6–31G(d)//B3LYP/6–31G(d) proves to give very good agreement with experimental data over a range from 1,400 to −1,300 ppm. The MP2/6–31G(d)//MP2/6–31G(d) level computes even smaller average deviations but fails for TiI₄. This behavior together with its huge demand for computational resources requires careful handling of this theoretical level. In addition, NMR data for five titanium fulvene (or related) complexes are given. Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

Conformation and tautomerizm of the 2-methyl-4-pyridin-2′-yl-1,5-benzodiazepine molecule. An ab initio study

The protomeric tautomerizm and conformation of the 2-methyl-4-pyridin-2′-yl-1,5-benzodiazepine molecule were investigated, and its three neutral tautomers (B₁,B₂,B₃) and their rotamers (C₁,C₂,C₃) were considered. Full geometry optimizations were carried out at the HF/6-31G* and B3LYP/6-31G* levels in gas phase and in water. The tautomerization processes in water (ɛ = 78.54) were studied by using self-consistent reaction field theory. The calculation showed that the boat conformation is dominant for the seven-membered diazepine ring in all of the structures, even with different double bond positions. The calculated relative free energies (ΔG) showed that the tautomer C₁ was the most stable structure, and its conformer B₁ was the second most stable in the gas phase and in water. Figure 2-Methyl-4-pyridin-2′-yl-1,5-benzodiazepine     

Interactions in pseudorotoxanes based on crown ether-secondary ammonium motifs. A theoretical study

A theoretical analysis of the nature of the interactions in dibenzo[24]crown-8 (DB24C8)-n-dibutylammonium (DBM)—pseudorotaxane complex at the MP2 and DFT levels shows that the main contribution to the binding energy is the electrostatic interaction with moderate (20–25%) correlation stabilization. The total binding energy in the DB24C8-DBM complex represents a sum of the binding energies of two NH–O and one CH–O hydrogen bonds and the latter constitutes about 25% of the total interaction energy, giving the total binding energy of −41.2 kcal mol⁻¹ at the BHandHLYP/6-311++G** level. Deprotonation of the DB24C8-DBM complex reduces the binding energy by some 50 kcal mol⁻¹, giving metastable complexes DB24C8-DBA-1 or DB24C8-DBA-2, which will dissociate to give free crown ether and n-dibutylamine because of the strong exchange repulsion that prevails in neutral complexes. Figure Formation of DB24C8-DBM pseudorotoxane complex     

Adsorption and phase transitions in adsorbed systems: structural properties of CCl<Subscript>4</Subscript> layers adsorbed on a graphite surface

We present the results of simulations of a CCl₄ monolayer adsorbed on a graphite surface. The CCl₄ molecule was represented either by a shapeless superatom or by its atomic sites. The simulations were carried out over a large range of temperatures, from 20 K up to 340 K. We address the following problems: (1) the influence of molecular shape on the structure and stability of phases (particularly at low temperatures), and (2) the influence of the graphite corrugation on layer stability and mechanism of phase transitions. In particular, we discuss the possibility and conditions of the appearance of hexatic phase in the system. Figure Temperature dependence of Φ6 order parameter for CCl₄ monolayer adsorbed onsmooth and corrugated surfaces, in the spherical Lennard Jones (LJ) approximation.For comparison, the order parameter calculated for MacDonald’s five-site potential is also presented     

Theoretical study of hydrogenation of the doubly aromatic B 7 − cluster

We have studied the influence of hydrogenation on the relative stability of the low-lying isomers of the anionic B₇⁻ cluster, computationally. It is known that the pure-boron B₇⁻ cluster has a doubly (σ- and π-) aromatic C_6v (³A₁) quasi-planar wheel-type triplet global minimum (structure 1), a low-lying σ-aromatic and π-antiaromatic quasi-planar singlet C_2v (¹A₁) isomer 2 (0.7 kcal mol⁻¹ above the global minimum), and a planar doubly (σ- and π-) antiaromatic C_2v (¹A₁) isomer 3 (7.8 kcal mol⁻¹ above the global minimum). However, upon hydrogenation, an inversion in the stability of the species occurs. The planar B₇H₂⁻ (C_2v, ¹A₁) isomer 4, originated from the addition of two hydrogen atoms to the doubly antiaromatic B₇⁻ isomer 3, becomes the global minimum structure. The second most stable B₇H₂⁻ isomer 5, originated from the quasi-planar triplet wheel isomer 1 of B₇⁻, was found to be 27 kcal mol⁻¹ higher in energy. The inversion in stability occurs due to the loss of the doubly aromatic character in the wheel-type global minimum isomer (C_6v, ³A₁) of B₇⁻ upon H₂−addition. In contrast, the planar isomer of B₇⁻ (C_2v, ¹A₁) gains aromatic character upon addition of two hydrogen atoms, which makes it more stable. Figure The B₇H₂-global minimum structure and its σ-aromatic and π-antiaromatic MOs Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday.     

Ab initio analysis of the Cope rearrangement of germacrane sesquiterpenoids

The energetics of the Cope rearrangement of 17 germacrane sesquiterpenoids to their respective elemane forms have been calculated using both density functional theory (B3LYP/6-31G*) and post Hartee-Fock (MP2/6-31G**) ab initio methods. The calculations are in qualitative agreement with experimentally observed Cope rearrangements, but the two methods give slightly different results. MP2 calculations generally show more favorable elemene energies compared to the respective germacrenes (by around 3–4 kcal mol⁻¹) and smaller activation energies (by 2–3 kcal mol⁻¹). Additionally, neither method is accurate enough to consistently reproduce the germacrene/elemene equilibrium. Apparently, the generally small energy differences between the two forms in these sesquiterpenoids cannot be adequately reproduced at these levels of calculation. Figure The Cope rearrangement of the germacrane sesquiterpenoid bacchascandon to the elemane shyobunone     

The platinum-olefin binding energy in series of (PH<Subscript>3</Subscript>)<Subscript>2</Subscript>Pt(olefin) complexes - a theoretical study

Theoretical investigation of Pt(0)-olefin organometallic complexes containing tertiary phosphine ligands was focused on the strength of platinum-olefin electronic interaction. DFT theoretical study of electronic effects in a substantial number of ethylene derivatives was evaluated in terms of the Pt-olefin binding energy using MP2 correlation theory. Organometallics bearing coordinated olefins with general formula (R¹R²C = CR³R⁴)Pt(PH₃)₂ [R = various substituents] had been selected, including olefins containing both electron-donor substituents as well as electron-withdrawing groups. The stability of the corresponding complexes increases with a strengthening electron-withdrawal ability of the olefin substituents. Figure Representation of (CH₂ = CHR)Pt(PPh₃)₂ and the stability chart     

Theoretical study of spin-spin coupling across the hydrogen (O-H⋯N) bond in adenosine derivatives

The study of spin-spin coupling constants across hydrogen bond provides useful information about configuration of complexes. The interesting case of such interactions was observed as a coupling across an intramolecular hydrogen bond in 8-bromo-2′,3′-O-isopropylideneadenosine between the -CH₂OH (at 5″ proton) group and the nitrogen atom of adenine. In this paper we report theoretical investigations on the ^4h J _NH coupling across the H″-C-O-H···N hydrogen bond in adenosine derivatives in various solvent models. Figure Coupling constants in 8-bromo-2′,3′-O-isopropylideneadenosine Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A computational investigation of sulfur-containing heterocyclic components from the anal sac secretions of Mustela species

A computational investigation of the sulfur-containing heterocyclic components (substituted thietanes and 1,2-dithiolanes) of Mustela anal sac secretions has been carried out. A cluster analysis of the chemical compositions of Mustela anal sac volatiles reveals little similarity with established phylogenetic relationships between members of the genus. Ab initio calculations [MP2/6–311++G(2df,2p)//B3LYP/6–311++G**] show the lowest-energy C₅H₁₀S isomeric thietane to be 2,2-dimethylthietane, which is also the most abundant of the Mustela thietanes. Similarly, 3,3-dimethyl-1,2-dithiolane is the lowest-energy C₅H₁₀S₂ compound. 2-n-Propylthietane is the highest-energy C₆H₁₂S compound, but the most abundant Mustela C₆H₁₂S compound produced, whereas cis-2-ethyl-4-methylthietane, the lowest-energy C₆H₁₂S thietane, has never been observed in Mustela anal sac secretions. A molecular docking analysis of the Mustela sulfur-containing heterocycles into both porcine and bovine odorant binding proteins reveals the interactions of the docked ligands with the proteins to be largely hydrophobic, and have binding energies generally lower than typical odorant molecules such as linalool or eugenol. Figure Mustela anal sac volatile components, 2,2-dimethylthietane and cis-3,4-dimethyl–1,2-dithiolane.     

The electronic structure of the C2H4O···2HF tri-molecular heterocyclic hydrogen-bonded complex: a theoretical study

The geometries of three isomers of the C₂H₄O···2HF tri-molecular heterocyclic hydrogen-bonded complex were examined through B3LYP/aug-cc-pVDZ calculations. Analysis of structural parameters, determination of CHELPG (charge electrostatic potential grid) intermolecular charge transfer, interpretation of infrared stretching modes, and Bader’s atoms in molecules (AIM) theory calculations was carried out in order to characterize the hydrogen bonds in each isomer of the C₂H₄O···2HF complex. The most stable structure was determined through the identification of hydrogen bonds between C₂H₄O and HF, (O···H), as well as in the hydrofluoric acid dimer, (HF^D–R···HF^D). However, the existence of a tertiary interaction (F^λ···H^α) between the fluoride of the second hydrofluoric acid and the axial hydrogen atoms of C₂H₄O was decisive in the identification of the preferred configuration of the C₂H₄O···2HF system. Figure Geometries of three isomers of the C₂H₄O···2HF tri-molecular heterocyclic hydrogen-bonded complex     

Molecular modeling and biophysical analysis of the c-MYC NHE-III<Subscript>1</Subscript> silencer element

G-Quadruplex and i-Motif-forming sequences in the promoter regions of several oncogenes show promise as targets for the regulation of oncogenes. In this study, molecular models were created for the c-MYC NHE-III₁ (nuclease hypersensitivity element III₁) from two 39-base complementary sequences. The NHE modeled here consists of single folded conformers of the polypurine intramolecular G-Quadruplex and the polypyrimidine intramolecular i-Motif structures, flanked by short duplex DNA sequences. The G-Quadruplex was based on published NMR structural data for the c-MYC 1:2:1 loop isomer. The i-Motif structure is theoretical (with five cytosine–cytosine pairs), where the central intercalated cytosine core interactions are based on NMR structural data obtained for a tetramolecular [d(A₂C₄)₄] model i-Motif. The loop structures are in silico predictions of the c-MYC i-motif loops. The porphyrin meso-tetra(N-methyl-4-pyridyl)porphine (TMPyP4), as well as the ortho and meta analogs TMPyP2 and TMPyP3, were docked to six different locations in the complete c-MYC NHE. Comparisons are made for drug binding to the NHE and the isolated G-Quadruplex and i-Motif structures. NHE models both with and without bound cationic porphyrin were simulated for 100 ps using molecular dynamics techniques, and the non-bonded interaction energies between the DNA and porphyrins calculated for all of the docking interactions. Figure Molecular models of the average structure of the final 20 ps of the molecular dynamics simulation of the c-MYC NHE-III₁ (nuclease hypersensitivity element III₁) “silencer” element. The G-Quadruplex structure is at the top-center, and the i-Motif is at the bottom-center of each picture. a “Rotation #1” of the G-Quadruplex, with the T15 loop at the top and rear and the G19/A20 loop at the top and front of the picture. b “Rotation #2” of the G-Quadruplex, with the T15 loop at the top and front of the image, and the G19/A20 loop at the front and adjacent to the G-Quadruplex/i-Motif interface     

Boron nitride cages from B12N12 to B36N36: square–hexagon alternants vs boron nitride tubes

The structures and stabilities of square–hexagon alternant boron nitrides (B _x N _x , x=12–36) vs their tube isomers containing octagons, decagons and dodecagons have been computed at the B3LYP density functional level of theory with the correlation-consistent cc-pVDZ basis set of Dunning. It is found that octagonal B₂₀N₂₀ and B₂₄N₂₄ tube structures are more stable than their square–hexagon alternants by 18.6 and 2.4 kcal mol⁻¹, respectively, while the square–hexagon alternants of other cages are more stable. Trends in stability as a function of cluster size are discussed.Figure The octagonal B₂₀N₂₀ and B₂₄N₂₄ tube structures are more stable than their square-hexagon alternant cagesDedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

Optimization of parameters for semiempirical methods V: Modification of NDDO approximations and application to 70 elements

Several modifications that have been made to the NDDO core-core interaction term and to the method of parameter optimization are described. These changes have resulted in a more complete parameter optimization, called PM6, which has, in turn, allowed 70 elements to be parameterized. The average unsigned error (AUE) between calculated and reference heats of formation for 4,492 species was 8.0 kcal mol⁻¹. For the subset of 1,373 compounds involving only the elements H, C, N, O, F, P, S, Cl, and Br, the PM6 AUE was 4.4 kcal mol⁻¹. The equivalent AUE for other methods were: RM1: 5.0, B3LYP 6–31G*: 5.2, PM5: 5.7, PM3: 6.3, HF 6–31G*: 7.4, and AM1: 10.0 kcal mol⁻¹. Several long-standing faults in AM1 and PM3 have been corrected and significant improvements have been made in the prediction of geometries. Figure Calculated structure of the complex ion [Ta₆Cl₁₂]²⁺ (footnote): Reference value in parenthesis Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Influence of different transition metals in phthalocyanines on their interaction energies with volatile organic compounds: an experimental and computational study

Experimental partition coefficients were determined for a series of volatile organic compounds (VOCs) (acetonitrile, n-butylamine, n-octane tetrachloroethene, and toluene) for the interaction with 2,3,9,10,16,17,23,24-octakis(octyloxy)-phthalocyaninato complexes, PcM(OR)₈, with varying central metal atoms [M=H₂ (metal-free), Ni, Pd, Cu, Zn]. Large partition coefficients for toluene were observed in the case of the nickel and palladium phthalocyanines, whereas for the corresponding zinc-containing compound, interaction with n-butylamine resulted in a high value for the partition coefficient. Interaction energies for model coordination complexes were obtained at the ab initio LMP2/ LACVP* level of theory. The interaction of various small volatiles with the various PcM(OR)₈ compounds was studied using the PM3 semiempirical Hamiltonian. Large values for interaction energies correspond to particularly strong partition coefficients, suggesting that coordination of the volatiles to the central metal dominates over the often discussed π-system stacking at the PcM(OR)₈’s. Figure: Chemical structure of the phthalocyanines investigated.     

Computational studies of the binding site of α<Subscript>1A</Subscript>-adrenoceptor antagonists

Aimed at achieving a good understanding of the 3-dimensional structures of human α_1A-adrenoceptor (α_1A-AR), we have successfully developed its homology model based on the crystal structure of β₂-AR. Subsequent structural refinements were performed to mimic the receptor’s natural membrane environment by using molecular mechanics (MM) and molecular dynamics (MD) simulations in the GBSW implicit membrane model. Through molecular docking and further simulations, possible binding modes of subtype-selective α_1A-AR antagonists, Silodosin, RWJ-69736 and (+)SNAP-7915, were examined. Results of the modeling and docking studies are qualitatively consistent with available experimental data from mutagenesis studies. The homology model built should be very useful for designing more potent subtype-selective α_1A-AR antagonists and for guiding further mutagenesis studies. Figure The superposition of β₂-AR crystal structure (gold ribbons) and α_1A-AR homology model (blue ribbons)     

ONIOM DFT/PM3 calculation on the interaction between STI-571 and abelson tyrosine kinase

The ONIOM2 (B3LYP/6–31G (d, p): PM3) and B3LYP/6–31G (d, p) methods were applied to investigate the interaction between STI-571 and abelson tyrosine kinase binding site. The complex of N-[4-methyl-3-(4-pyridin-3-yl-pyrimidin-2-ylamino)- phenyl]-benzamide (part of STI-571) and related 16 amino acid residues were found at B3LYP/6–31G (d, p) level to have hydrogen bonds and π....π stacking interaction, their binding energy via HAF optimization was −20.4 kcal mol⁻¹. The results derived from this study agreed well with the reported observation. Figure Optimized structure of STI-571 and Thr315 in abelson tyrosine kinase based on ONIOM2 method