Ab initio calculations on the thermodynamic properties of azaborospiropentanes

Following our recent studies of the thermodynamic properties of azaspiropentane and borospiropentane, in consideration of their usefulness as new potential high energy materials, we follow up with ab initio calculations on the thermodynamic properties of azaborospiropentanes. Properties reported in this study include optimized structural parameters, vibrational frequencies, enthalpies of formation, specific enthalpies of combustion, proton affinities, and hydride affinities. Our results indicate that azatriborospiropentane gives off most energy when combusted, as evidenced by its specific enthalpy of combustion of about −52 kJ per gram. Figure Optimized geometry for R-azatriborospiropentane (10) Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The design of organic catalysis for epoxidation by hydrogen peroxide

The potential of various organic species to catalyze epoxidation of ethene by hydrogen peroxide is explored with B3LYP/6-31G* DFT calculations. Electronic Supplementary Materials Supplementary material is available for this article at http://dx.doi.org/10.1007/s00894-005-0044-4.Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

Density functional computations of enantioselective alkynylation of aldehyde catalyzed by chiral zinc(II)-complexes

The enantioselective alkynylation of aldehyde catalyzed by chiral zinc(II)-complexes was studied by means of the density functional theory (DFT). All the structures were optimized completely at the B3LYP/6-31G(d,p) level. To obtain more exact energies, single-point energy calculations at B3LYP/6-31+G(d,p) level were carried out on the B3LYP/6-31G(d,p) geometries. As shown, this enantioselective alkynylation was endothermic. The chirality-determining step for the alkynylation was the formation of the catalyst–ethanol complexes and the transition states for this step involved a six-membered ring. The dominant products predicted theoretically were of (R)-chirality, in good agreement with experiment.      

A computational approach to the synthesis of 1,3,5-thiadiazinane-2-thiones in aqueous medium: theoretical evidence for water-promoted heterocyclization

Based on experimental evidence and DFT studies, a probable cyclization route to 1,3,5-thiadiazinanes-2-thiones in aqueous medium is proposed. Experimental facts suggest the formation of a {[hydroxymethyl (substituted) carbamothioyl] sulfanyl}methanol intermediate via reaction of dithiocarbamate (DTC) and formaldehyde. Nucleophilic addition of glycine to this intermediate generates an adduct that undergoes intramolecular heterocyclization via an S_N2 reaction. Computational calculations predict an active role of water in the reaction mechanism that promotes intramolecular cyclization. Figure Energy profile of the proposed reaction mechanism for the synthesis of thiadiazinane-2-thione ring 11 in aqueous medium from a (hydroxymethylcarbamothioyl)sulfanylmethanol intermediate, 9     

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₃     

Structure prediction of a multi-domain EF-hand Ca2+ binding protein by PROPAINOR

PROPAINOR is a new algorithm developed for ab initio prediction of the 3D structures of proteins using knowledge-based nonparametric multivariate statistical methods. This algorithm is found to be most efficient in terms of computational simplicity and prediction accuracy for single-domain proteins as compared to other ab initio methods. In this paper, we have used the algorithm for the atomic structure prediction of a multi-domain (two-domain) calcium-binding protein, whose solution structure has been deposited in the PDB recently (PDB ID: 1JFK). We have studied the sensitivity of the predicted structure to NMR distance restraints with their incorporation as an additional input. Further, we have compared the predicted structures in both these cases with the NMR derived solution structure reported earlier. We have also validated the refined structure for proper stereochemistry and favorable packing environment with good results and elucidated the role of the central linker. Figure The predicted 3D Structure of EhCaBP with bound Ca²⁺ ions (CaBP-0). In the structure, α-helices are shown in pink and the β-strands in yellow. Ca²⁺ ions are depicted as fluorescent green balls. Some of the residues in the calcium-binding loops are depicted in space-fill representation.      

Synthesis, DFT calculations, linear and nonlinear optical properties of binuclear phthalocyanine gallium chloride

The axially substituted binuclear GaCl/GaCl phthalocyanine 1 with an unsymmetrical pattern of substitution has been prepared and its nonlinear optical (NLO) properties determined. The resulting binuclear complex retains approximately the same transition energies of monomeric (RO)₈PcGaCl as far as the linear optical spectrum is concerned, although 1 has a double concentration of central atoms per molecule and an enlarged conjugated ligand. The lack of significant spectral shifts in passing from mononuclear to binuclear complexes has been rationalized theoretically by means of density functional theory calculations. The purpose of the present study is to determine whether binuclearity affects the optical limiting behavior of 1 with respect to monomeric (RO)₈PcGaCl in the NLO regime determined by nanosecond laser pulses. Figure Bis axially substituted binuclear phthalocyanine: synthesis, DFT calculations and NLO properties Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday.     

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)     

Molecular models of protein kinase 6 from Plasmodium falciparum

Cyclin-dependent kinases (CDKs) have been identified as potential targets for development of drugs, mainly against cancer. These studies generated a vast library of chemical inhibitors of CDKs, and some of these molecules can also inhibit kinases identified in the Plasmodium falciparum genome. Here we describe structural models for Protein Kinase 6 from P. falciparum (PfPK6) complexed with Roscovitine and Olomoucine. These models show clear structural evidence for differences observed in the inhibition, and may help designing inhibitors for PfPK6 generating new potential drugs against malaria. Figure Ribbon diagram of PfPK6 complexed with a roscovitine and b olomoucine     

An ab initio quantum mechanical drug designing procedure: application to the design of balanced dual ACE/NEP inhibitors

This article describes in a sequential fashion how ab initio quantum mechanical methods can be applied to study the pharmacophoric features of drugs. It also describes how accurate drug–receptor interaction calculations can guide the careful design of balanced dual inhibitors, which form an important class of second generation drugs. As an example, the authors have chosen balanced inhibitors of angiotensin converting enzyme/neutral endopeptidase (ACE/NEP) as modern antihypertensive drugs. A unified, accurate, in silico design approach is presented, encompassing all steps from pharmacophore derivation to complete understanding of mechanistic aspects leading to drug design.      

Alpha particle chemistry. On the formation of stable complexes between He2+ and other simple species: implications for atmospheric and interstellar chemistry

The possibility that stable complexes may be formed between alpha particles (He²⁺) and small molecules is investigated using QCISD quantum mechanical calculations. Implications for their presence in the terrestrial atmosphere and/or in interstellar space are discussed. Figure Optimized structure of a stable H₂OHe²⁺ complex     

Ab initio and DFT study of the aromaticity of some Fulvalenes derived from Methylidenecyclopropabenzene

Methylidencyclopropabenzene (MCPB) 1 and Fulvalenes 2–4 are molecules of special interest due to the relation between structure and aromaticity. The aim of this work was to analyze this relation and to quantify the aromaticity in 1–4 using different methods. Magnetic properties are directly related with aromaticity; here we studied the magnetic susceptibility and the anisotropy of the magnetic susceptibility. Nucleus indepedent chemical shift (NICS) and the anisotropy of the induced current density (ACID) were also employed. Tools of very different nature, geometric indexes HOMA and Bird, were determinated too for 1–4. All of these measures were found to be in agreement. Figure Both spatial NICS and ACID plot allow to show the aromaticity/antiaromaticity of a ring     

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.      

Study of a ligand complexed with Cdk2/Cdk4 by computer simulation

Cyclin-dependent kinases (Cdks) play important roles in the regulation of the cell cycle. Their inhibitors have entered clinical trials to treat cancer. Very recently, Davis et al. (Nat Struct Biol 9:745–749, 2002) have found a ligand NU6102, which has a high affinity with cyclin-dependent kinase 2 (K _i =6 nM) but a low affinity with cyclin-dependent kinase 4 (K _i =1,600 nM). To understand the selectivity, we use homology modeling, molecular docking, molecular dynamics and free-energy calculations to analyze the interactions. A rational 3D model of the Cdk4–NU6102 complex is built. Asp86 is a key residue that recognizes NU6102 more effectively with Cdk2 rather than Cdk4. Good binding free energies are obtained. Energetic analysis reveals that van der Waals interaction and nonpolar contributions to solvent are favorable in the formation of complexes and the sulfonamide group of the ligand plays a crucial role for binding selectivity between Cdk2 and Cdk4. Figure Two-dimensional representative for the interacting model of NU6102 complexed with the Cdk4 from a predicted structure by LIGPLOT.      

Role of hydrophobic interactions and salt-bridges in β-hairpin folding

β-Hairpins are the simplest form of β-sheets which, due to the presence of long-range interactions, can be considered as tertiary structures. Molecular dynamics simulation is a powerful tool that can unravel whole pathways of protein folding/unfolding at atomic resolution. We have performed several molecular dynamics simulations, to a total of over 250 ns, of a β-hairpin peptide in water using GROMACS. We show that hydrophobic interactions are necessary for initiating the folding of the peptide. Once formed, the peptide is stabilized by hydrogen bonds and disruption of hydrophobic interactions in the folded peptide does not denature the structure. In the absence of hydrophobic interactions, the peptide fails to fold. However, the introduction of a salt-bridge compensates for the loss of hydrophobic interactions to a certain extent. Figure Model of b-hairpin folding: Folding is initiated by hydrophobic interactions (Brown circles). The folded structure, once formed, is stabilized by hydrogen bonds (red lines) and is unaffected by loss of hydrophobic contacts     

TCNE-aniline charge transfer complex: ab initio and TDDFT investigations in gas phase

The geometric and electronic structure of tetracyanoethylene (TCNE)-aniline (donor-acceptor type) complex has been investigated in gas phase using ab initio and time dependent density functional theory calculations. Both the above calculations predict a composed structure for the complex, in which the interacting site is a C≡N and C=C bond center in the TCNE and, –NH₂ and π-electrons of aniline. The N atom of aniline is oriented toward the TCNE molecule. The charge transfer transition energy, estimated by calculating the ground-to-excited state transition electric dipole moments of the complex, agree well with the reported experimental value in chloroform medium. TCNE-aniline at ground state. TCNE-aniline at excited state     

Structures of inclusion complexes of halogenbenzoic acids and α-cyclodextrin based on AM1 calculations

Semi-empirical quantum mechanics calculations using AM1 (Austin Method 1) were carried out for various host-guest combinations of α-cyclodextrin and mono-halogen benzoic acids. The energetically favorable inclusion structures were identified. The AM1 results show that α-cyclodextrin complexes with mono-halogen benzoic acid acids (where the halogen is chlorine, bromide, iodine) as guest compounds are more stable in the “head first” position than in the “tail-first” position for meta and para isomers while ortho mono-halogen benzoic acids complexes with α-cyclodextrin are more stable in “tail-first” position. The calculated structures were found to be in good agreement with those obtained from crystalographic databases.      

Elongation cutoff technique: low-order scaling SCF method

The elongation cutoff technique at restricted Hartree-Fock (HF) level of theory in conventional type of calculations, i.e., with two electron integrals (TEI) stored on a disc, is presented for two model systems. It is demonstrated that the number of TEI in the elongation cutoff calculations increases linearly with the system size thus, allowing to extend the conventional type of calculations to bigger systems. The step CPU (central processing unit) time in the elongation cutoff calculations is much lower than in the HF reference calculations. Such behavior reduces significantly the prefactor in the quadratic scaling relation. The total CPU time in the elongation calculation is about 40% lower than in the conventional HF calculations or comparable to direct type of calculations with the quantum fast multipoles method employed. It is shown that by introducing the interaction radius one can obtain linear scaling in the SCF calculations. Figure: The structure of density matrix and total CPU timings for polyglycine clusters in the elongation cutoff calculations. The structure of density matrix and total CPU timings for polyglycine clusters in the elongation cutoff calculations     

Aromaticity and electronic properties of Heterosuperbenzene (Heterohexabenzocoronene)

A global electrophilicity parameter and the aromaticity of some heterocyclic polyaromatic hydrocarbons were evaluated on the basis of DFT calculations. The substitution of carbon atoms by nitrogen atoms dramatically changes the global electrophilicity of the molecules, with the fully substituted molecule being the most electrophilic with a reactivity very close to that of fullerene. Figure Fully substituted heterohexabenzocoronene Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Theoretical study of the novel sandwich compound [Au<Subscript>3</Subscript>Cl<Subscript>3</Subscript>Tr<Subscript>2</Subscript>]<Superscript>2+</Superscript>

A theoretical study of a sandwich compound with a metal monolayer sheet between two aromatic ligands is presented. A full geometry optimization of the [Au₃Cl₃Tr₂]²⁺ (1) compound, which is a triangular gold(I) monolayer sheet capped by chlorines and bounded to two cycloheptatrienyl (Tr) ligands was carried out using perturbation theory at the MP2 computational level and DFT. Compound (1) is in agreement with the 18–electron rule, the bonding nature in the complex may be interpreted from the donation interaction coming from the Tr rings to the Au array, and from the back-donation from the latter to the former. NICS calculations show a strong aromatic character in the gold monolayer sheet and Tr ligands; calculations done with HOMA, also report the same aromatic behavior on the cycloheptatrienyl fragments giving us an insight on the stability of (1). The Au –Au bond lengths indicate that an intramolecular aurophilic interaction among the Au(I) cations plays an important role in the bonding of the central metal sheet. Figure (a) Ground state geometry of complex 1; (b) Top view of compound 1 and Wiberg bond orders computed with the MP2/B1 computational method; (c) Lateral view of compound 1 and NICS values calculated with the MP2/B1 method; the values in parenthesis were obtained at the VWN/TZP level