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     

Docking studies suggest ligand-specific δ-opioid receptor conformations

An automated docking procedure was used to study binding of a series of δ-selective ligands to three models of the δ-opioid receptor. These models are thought to represent the three ligand-specific receptor conformations. Docking results are in agreement with point mutation studies and suggest that different ligands—agonists and antagonists—may bind to the same binding site under different receptor conformations. Docking to different receptor models (conformations) also suggests that by changing to a receptor-specific conformation, the receptor may open or close different binding sites to other ligands. Figure  Ligands 5 (green) and 6 (orange) in bindingpocket BP1 of the R1 δ-opioid receptor model     

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     

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₃     

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.      

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     

Mechanistic study of palladium-catalyzed telomerization of 1,3-butadiene with methanol

The Pd-catalyzed telomerization in the presence of phosphine and carbene ligands has been computed. It is shown that the C–C coupling of the less stable complex A with one trans- and one cis-butadiene in syn orientation forms the most stable intermediate B and is favorable both kinetically and thermodynamically. Protonation of B leads to equilibrium of the two most stable isomers of intermediate C. The overall regioselectivity is favored thermodynamically.      

Probing the nature of hydrogen bonds in DNA base pairs

Energy decomposition analyses based on the block-localized wave-function (BLW-ED) method are conducted to explore the nature of the hydrogen bonds in DNA base pairs in terms of deformation, Heitler–London, polarization, electron-transfer and dispersion-energy terms, where the Heitler–London energy term is composed of electrostatic and Pauli-exchange interactions. A modest electron-transfer effect is found in the Watson–Crick adenine–thymine (AT), guanine–cytosine (GC) and Hoogsteen adenine-thymine (H-AT) pairs, confirming the weak covalence in the hydrogen bonds. The electrostatic attraction and polarization effects account for most of the binding energies, particularly in the GC pair. Both theoretical and experimental data show that the GC pair has a binding energy (−25.4 kcal mol⁻¹ at the MP2/6-31G** level) twice that of the AT (−12.4 kcal mol⁻¹) and H-AT (−12.8 kcal mol⁻¹) pairs, compared with three conventional N-H···O(N) hydrogen bonds in the GC pair and two in the AT or H-AT pair. Although the remarkably strong binding between the guanine and cytosine bases benefits from the opposite orientations of the dipole moments in these two bases assisted by the π-electron delocalization from the amine groups to the carbonyl groups, model calculations demonstrate that π-resonance has very limited influence on the covalence of the hydrogen bonds. Thus, the often adopted terminology “resonance-assisted hydrogen bonding (RHAB)” may be replaced with “resonance-assisted binding” which highlights the electrostatic rather than electron-transfer nature of the enhanced stabilization, as hydrogen bonds are usually regarded as weak covalent bonds. Figure Electron density difference (EDD) maps for the GC pair: a shows the polarization effect (isodensity 1.2×10⁻³ a.u.); b shows the charge transfer effect (isodensity 2×10⁻⁴ a.u.) Dedicated to Professor Paul von Ragué Schleyer on the occasion of his 75th birthday     

Theoretical study on the thermal decomposition of model compounds for Poly (dialkyl fumarate)

The thermal decomposition of model compounds for poly (dialkyl fumarate) was studied by using ab initio and density functional theory (DFT) calculations. To determine the most favorable reaction pathway of thermal decomposition, geometries, structures, and energies were evaluated for reactants, products, and transition states of the proposed pathways at the HF/6-31G(d) and B3LYP/6-31G(d) levels. Three possible paths (I, II and III) and subsequent reaction paths (IV and V) for the model compounds of poly (dialkyl fumarate) decomposition had been postulated. It has been found that the path (I) has the lowest activation energy 193.8 kJ mol⁻¹ at B3LYP/6-31G(d) level and the path (I) is considered as the main path for the thermal decomposition of model compounds for poly (dialkyl fumarate).      

A B3LYP and MP2 theoretical investigation into host-guest interaction between calix[4]arene and Li+ or Na+

The DFT-B3LYP and MP2 methods with 6-311G** and 6-311++G** basis sets have been applied to study the complexation energies of the host-guest complexes between the cone calix[4]arene and Li⁺ or Na⁺ on the B3LYP optimized geometries. A comparison of the complexation energies obtained from the MP2(full) with those from MP2(fc) method is also carried out. The result shows that it is essential to introduce the diffuse basis set into the geometry optimizations and complexation energy calculations of the alkali-metal cation-π interaction complexes of calix[4]arene, and the D _e values show a maximum of 21.13 kJ mol⁻¹ (14.45% of relative error) between the MP2(full)/6-311++G** and MP2(fc)/6-311++G** method. For Li⁺ cation, the complexation is mainly energetically stabilized by the lower rim/cation (namely O–Li⁺) interaction. However, binding energies and NBO analyses confirm that Na⁺ cation prefers to enter the calix[4]arene cavity and the cation-π interaction is predominant, which contradicts the previous low-level theoretical studies. Furthermore, the complexation with Li⁺ is preferred over that with Na⁺ by at least 12.70 kJ mol⁻¹ at MP2(full)/6-311++G**//B3LYP/6-311++G** level.      

Molecular model for host–guest interaction of tetraamino-<Emphasis Type="Italic">tert</Emphasis>-butylthiacalix[4]arene and tetraamino-<Emphasis Type="Italic">tert</Emphasis>-butylcalix[4]arene receptors with carboxylate and dicarboxylate guests: an ONIOM study

Geometry optimizations of tetraamino-tert-butylthiacalix[4]arene (tatbtc4a) and tetraamino-tert-butylcalix[4]arene (tatbc4a) complexes with acetate, oxalate, malonate, succinate, glutarate, adipate, and pimelate were carried out using the integrated MO:MO method. Thermodynamic quantities, preorganization energies and complexation energies of these complexes were obtained at the ONIOM(B3LYP/6-31G(d):AM1) level of theory. The relative stabilities of the tatbtc4a and tatbc4a complexes with carboxylate guests are reported. The complexes tatbtc4a/malonate and tatbc4a/oxalate were found to be the most stable species. The selectivity of the tatbtc4a receptor toward to malonate with respect to oxalate, in terms of selectivity coefficient, is 9.90×10². Figure Atom labeling of tatbtc4a/oxalate complex as a representative of host-guest system.     

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.     

Synchronous fluorescence determination and molecular modeling of 5-Aminosalicylic acid (5-ASA) interacted with human serum albumin

In this paper, we proposed a new method for the determination of either human serum albumin (HSA) or 5-Aminosalicylic acid (5-ASA) by synchronous fluorescence spectra and examined the interaction between them using the molecular modeling method under simulative physiological conditions. The optimum conditions of synchronous fluorometric determination of HSA were investigated and the method was successfully applied to the determination of 5-ASA added to serum, urine, and saliva samples. The linear range of the determination of HSA and 5-ASA were 1.60 – 414 μg mL⁻¹ and 0.76 –22.95 μg mL⁻¹, the detection limits were 0.552 μg mL⁻¹ and 0.38 μg mL⁻¹, respectively. In addition, the effect of various common ions on the determination of HSA with 5-ASA was also discussed at room temperature. Figure The salicylic acid moiety is located within the binding pocket. The ring of 5-ASA was inserted in the hydrophobic cavity of site I, and it is important to note that the residue ARG-218 and the trptophan residue of HSA (Trp214) are in close proximity to the ring of 5-ASA suggesting the existence of hydrophobic interaction between them.     

Aromatic C20F20 cage and its endohedral complexes X@C20F20 (X = H-, F-, Cl-, Br-, H, He)

The structure and stability of endohedral X@C₂₀F₂₀ complexes (X = H⁻, F⁻, Cl⁻, Br⁻, H, He) have been computed at the B3LYP level of theory. All complexes in I _h symmetry were found to be energy minimum structures. H⁻@C₂₀F₂₀ and F⁻@C₂₀F₂₀ complexes have negative inclusion energies, while other complexes have positive inclusion energies. Similarity between C₂₀F₂₀ and C₂₀H₂₀ has been found for X = H and He. On the basis of the computed nucleus independent chemical shift values at the cage center, both C₂₀F₂₀ and C₂₀F₂₀ are aromatic. Figure Endohedral X@C₂₀F₂₀ complexes     

Three-dimensional model of zeaxanthin binding PsbS protein associated with nonphotochemical quenching of excess quanta of light energy absorbed by the photosynthetic apparatus

A three-dimensional model of the PsbS protein was built with the help of homology-modeling methods. This protein is also known as CP22 and is associated with the protection of photosystem II of thylakoid from excess quanta of light energy absorbed by the photosynthetic apparatus. PsbS is reported to bind two molecules of zeaxanthin at low pH (<5.0) and is believed to be essential for rapid nonphotochemical quenching (qE) of chlorophyll a fluorescence in photosystem II. An attempt was made to explain the pH modulation of the conformation of protein through salt-bridges Glu⁻(122)-Lys⁺(113) and Glu⁻(226)-Lys⁺(217). Binding of two molecules of zeaxanthin in the three-dimensional model of PsbS is postulated. The molecular mechanism of photoprotection by PsbS is explained through the model. 1 Backbone structure of the PsbS protein with two molecules of all trans zeaxanthin (ZEX). Residues Glu 90, 122, 194, 226 and Lys 113, 217 are shown. The figure is drawn with RASMOL (Molecular Visualization Program, RasMol V2.6, Roger Sayle, Glaxo Wellcome Research and Development, Stevenage, Hertfordshire, UK) Electronic Supplementary Material Supplementary material is available for this article at     

Ab initio and DFT study on the electrophilic addition reaction of bromine to tetracyclo[5.3.0.02,6.03,10]deca–4,8–diene

The electronic and geometric structures of tetracyclo[5.3.0.0^2,6.0^3,10]deca-4,8-diene (hypostrophene) have been investigated by ab initio and DFT/B3LYP methods using the 6-31G* and 6-311G* basis sets. The double bonds of hypostrophene are endo-pyramidalized. The cationic intermediates and products formed in the addition reaction have been investigated using the HF/6-311G*, HF/6-311G**, and B3LYP/6-311G* methods. The bridged bromonium cation was more stable than the U-type cation. Considering that the bridged cation does not isomerize to the less stable U-type cation, it is not possible for the U-type product to be obtained in the reaction. The bridged bromonium cation transformed into the more stable N-type cation and the N-type product was obtained via this cation. The thermodynamic stability of the exo, exo and exo, endo isomers of the N-type dibromide molecule were almost identical. The N-type product was 16.6 kcal mol⁻¹ more stable than the U-type product. Figure General energy diagram of the hypostrophene–bromine (HS–Br₂) system (kcal mol⁻¹) (MP2/6-311G*//HF/6-311G*)     

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.     

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     

Molecular surface electrostatic potentials and anesthetic activity

General anesthetics apparently act through weak, noncovalent and reversible interactions with certain sites in appropriate brain proteins. As a means of gaining insight into the factors underlying anesthetic potency, we have analyzed the computed electrostatic potentials V _S(r) on the surfaces of 20 molecules with activities that vary between zero and high. Our results are fully consistent with, and help to interpret, what has been observed experimentally. We find that an intermediate level of internal charge separation is required; this is measured by Π, the average absolute deviation of V _S(r), and the approximate window is 7 < Π < 13 kcal mol⁻¹. This fits in well with the fact that anesthetics need to be lipid soluble, but also to have some degree of hydrophilicity. We further show that polyhalogenated alkanes and ethers, which include the most powerful known anesthetics, have strong positive potentials, V _S,max, associated with their hydrogens, chlorines and bromines (but not fluorines). These positive sites may impede the functioning of key brain proteins, for example by disrupting their normal hydrogen-bond patterns. It has indeed been recognized for some time that the most active polyhalogenated alkanes and ethers contain hydrogens usually in combination with chlorines and/or bromines. Figure The computed HF/6-31G* electrostatic potential, in kcal mol⁻¹, on the 0.001 electrons/bohr³ surface of halothane, CF₃CHBrCl. The color ranges are: red, more positive than 25; yellow, between 15 and 25; green between 0 and 15; blue, between −10 and 0. The strongly positive (red) potential is due to the hydrogen; the yellow and green positive regions at the right are on the bromine surface Proceedings of “Modeling Interactions in Biomolecules II”, Prague, September 5th–9th, 2005.