A theoretical investigation of the mechanism for the reaction between a quebrachitol derivative and N3−

The reaction between a mesylated compound and sodium azide was previously studied experimentally at a temperature of 140 °C using dimethylformamide as a solvent. The product was assigned on the basis of the analysis of the NMR spectra. In this work semiempirical (AM1 and PM3), ab initio (Hartree–Fock and MP2) and density functional theory (BLYP functional) quantum mechanical calculations, using continuum models for describing the solvent effect, were carried out for this process to better understand the reaction mechanism. Three distinct mechanisms involving a carbocation and epoxide intermediates, and a transition-state structure for direct attack of the N₃ ⁻ species to the reactant were investigated. The theoretically calculated preferred reaction pathway passing through an epoxide intermediate agrees nicely with the experimental proposal, providing a good example of where theoretical calculations can be of great help to definitively elucidate the reaction mechanism. Received: 10 July 2001 / Accepted: 20 December 2001 / Published online: 8 April 2002     

A theoretical study of photoinduced electron transfer between tetracyanoethylene and arene

 Ab initio calculations have been performed to investigate the state transition in photoinduced electron transfer reactions between tetracyanoethylene and biphenyl as well as naphthalene. Face-to-face conformations of electron donor–acceptor (EDA) complexes were selected for this purpose. The geometries of the EDA complexes were determined by using the isolated optimized geometries of the donor and the acceptor to search for the maximum stabilization energy along the center-to-center distance. The correction of interaction energies for basis set superposition error was considered by using counterpoise methods. The ground and excited states of the EDA complexes were optimized with complete-active-space self-consistent-field calculations. The theoretical study of the ground state and excited states of the EDA complex in this work reveals that the S₁ and S₂ states of the EDA complexes are charge–transfer (CT) excited states, and CT absorption which corresponds to the S₀→S₁ and S₀→S₂ transitions arise from π−π^* excitation. On the basis of an Onsager model, CT absorption in dichloromethane was investigated by considering the solvent reorganization energy. Detailed discussions on the excited state and on the CT absorptions were made. Received: 30 April 2001 / Accepted: 18 October 2001 / Published online: 9 January 2002     

Computational insight into the effect of C(19) substituents on [1,7]-hydrogen shift in previtamin D

 B3LYP calculations in conjunction with natural bond orbital population analysis have been performed for a previtamin D model and corresponding transition structures for the [1,7]-hydrogen migration. In addition the 19,19-difluoro, 19-methoxy and 19-fluoro substituted analogs were investigated. The calculated activation barriers decrease in the following order: CHF₂>CH₃>CH₂OCH₃ (24.8, 23.5 and 20.1 kcal/mol). This is in qualitative agreement with experiments. It has been suggested that a decrease of the barrier by a 19-methoxy substituent and its increase by a 19,19-difluoro substituent are phenomena of different origin. In the case of 19-methoxy substitution, the effect is due to the charge redistribution in the triene system and the decrease of the C(19)–H bond energy. The effect of two fluorine substituents at C-19 on the activation barrier is suggested to originate from the combination and balance of several factors: electrostatic repulsion between the negative fluorine atom and the π-electron cloud over the conjugated system, an increase of the HOMO–LUMO gap, and geminal difluoro substitution affecting C–F and C–C bond energies. Received: 17 May 2002 / Accepted: 11 September 2002 / Published online: 14 February 2003     

Action of HCl on 3-hydroxypyrazolo(iso)quinolines to give 1-chloropyrazoles: evidence for an addition–elimination mechanism by ab initio calculations in gas phase and water

 Addition–elimination reactions involving a nucleophile and a remote leaving group [S^H _N(AE)^tele] are well-known under basic conditions, especially amongst electron-poor six-membered heterocycles, but are less commonly encountered for five-membered heterocycles and are rare under acidic conditions. Concentrated HCl converts 1-hydroxy-1H-pyrazolo[3,4-c] isoquinoline and 1-hydroxy-1H-pyrazolo[3,4-c]quinoline into 3-chloro-1H-pyrazolo[3,4-c]isoquinoline and 3-chloro-1H-pyrazolo[3,4-c]quinoline, respectively. However, apparently neither the isomeric 1-hydroxy-1H-pyrazolo[4,3-c](iso)-quinolines nor the parent 1-hydroxypyrazole undergo this reaction. Additionally, all these systems are refractory under basic conditions. We present a plausible mechanism for the reaction, involving the 3-addition of Cl^- to the diprotonated heterocycle, followed by the elimination of water. Calculations of the initial transition states and intermediates, using optimisation at B3LYP/6-311+G(d,p), including thermochemistry [HF/6-31+G(d)], and single-point Poisson–Boltzmann self-consistent reaction field determination of the free energy of solvation (Jaguar Poisson–Boltzmann self-consistent reaction field), support this mechanism and reproduce the observed order of reactivity, the addition step being 2–4 kcal less favourable for the isomeric 1-hydroxy-1H-pyrazolo[4,3-c](iso)quinolines and provide a rationalisation for the role of strong acid. Received: 27 June 2002 / Accepted: 6 September 2002 / Published online: 14 February 2003     

Ab initio and density functional theory studies on the mechanism of nucleophilic vinylic substitution of 4<Emphasis Type="Italic">H</Emphasis>-pyran-4-one and 2-methyl-4<Emphasis Type="Italic">H</Emphasis>-pyran-4-one with ammonia

 Nucleophilic vinylic substitutions of 4H-pyran-4-one and 2-methyl-4H-pyran-4-one with ammonia were calculated by the B3LYP method using the 6-31G(d,p) basis set. Bulk solvent effects of aqueous solution were estimated by the polarized continuum and Poisson–Boltzmann self-consistent reaction field models using the 6-311+G(d,p) basis set. In the gas phase different mechanisms were found for the two reaction systems calculated. The reaction of 4H-pyran-4-one proceeds through enol, whereas a feasible path for the less reactive 2-methyl-4H-pyran-4-one is the mechanism through a keto intermediate. Addition of ammonia in concert with proton transfer is the rate-determining step ofthe reaction. The mechanism proceeding either by a bimolecular nucleophilic substitution (S_N2) or by one involving a tetrahedral zwitterionic intermediate is shown to be unlikely in the gas phase or nonpolar solution. The effects of bulk solvent not only consist in a reduction of the various activation barriers by about 25–40 kJ mol⁻¹ but also in a change in the reaction mechanism. Received 26 May 2002 / Accepted 26 July 2002 / Published online: 14 February 2003     

Partial Hessian vibrational analysis: the localization of the molecular vibrational energy and entropy

The partial Hessian vibrational analysis (PHVA), in which only a subblock of the Hesssian matrix is diagonalized to yield vibrational frequencies for partially optimized systems, is extended to the calculation of vibrational enthalpy and entropy changes for chemical reactions. The utility of this method is demonstrated for various deprotonation reactions by reproducing full HVA values to within 0.1–0.4 kcal/mol, depending on the number atoms included in the PHVA. When combined with the hybrid effective fragment potential method [Gordon MS, et al. (2001) J Phys Chem A 105:293–307], the PHVA method can provide (harmonic) free-energy changes for localized chemical reactions in very large systems. Received: 21 September 2001 / Accepted: 30 October 2001 / Published online: 22 March 2002     

Multireference perturbation configuration interaction V. Third-order energy contributions in the Møller–Plesset and Epstein–Nesbet partitions

 A simple implementation of third-order perturbation theory applied to a multireference zero-order wavefunction is presented. Two different partitions of the Hamiltonian (M?ller–Plesset baricentric and Epstein–Nesbet) are considered. Two test cases, CH₂ and N₂, are examined. The third-order results are shown to be in good agreement in either partition and are generally an improvement with respect to the second-order results. The phenomenon of intruder states, absent in Epstein–Nesbet, appears to be magnified in the M?ller–Plesset partition. Received: 22 November 2001 / Accepted: 22 February 2002 / Published online: 3 May 2002     

Maxwell–Cartesian spherical harmonics in multipole potentials and atomic orbitals

 The nature of the Maxwell–Cartesian spherical harmonics S ⁽ⁿ⁾ _K and their relation to tesseral harmonics Y _nm is examined with the help of “tricorn arrays” that display the components of a totally symmetric Cartesian tensor of any rank in a systematic way. The arrays show the symmetries of the Maxwell–Cartesian harmonic tensors with respect to permutation of axes, the traceless properties of the tensors, the linearly independent subsets, the nonorthogonal subsets, and the subsets whose linear combinations produce the tesseral harmonics. The two families of harmonics are related by their connection with the gradients of 1/r, and explicit formulas for the transformation coefficients are derived. The rotational transformation of S ⁽ⁿ⁾ _K functions is described by a relatively simple Cartesian tensor method. The utility of the Maxwell–Cartesian harmonics in the theory of multipole potentials, where these functions originated in the work of Maxwell, is illustrated with some newer applications which employ a detracer exchange theorem and make use of the partial linear independence of the functions. The properties of atomic orbitals whose angular part is described by Maxwell–Cartesian harmonics are explored, including their angular momenta, adherence to an Uns?ld-type spherical symmetry relation, and potential for eliminating an angular momentum “contamination” problem in Cartesian Gaussian basis sets. Received: 9 July 2001 / Accepted: 7 September 2001 / Published online: 19 December 2001     

Calculations of the geometries, electronic structures and nonlinear second-order optical susceptibilities of spiroannelated quinone-type methanofullerenes

 The equilibrium geometries, electronic structures and UV–vis spectra of a series of spiroannelated quinone-type methanofullerenes have been determined by using Zerner's intermediate neglect of differential overlap method. The results show that between fullerene and the addend there exists a special interaction, “periconjugation”, which results in through-space orbital interactions. The calculated UV–vis spectra are in good agreement with experiments. On the basis of the electronic spectra, the β values are calculated. The results show that spiroannelated quinone-type methanofullerenes have quite large β values. We attribute the large β values to both the charge transfer from C₆₀ to benzoquinone and on the C₆₀ three-dimensional conjugated sphere. Received: 17 December 2000 / Accepted: 16 March 2001 / Published online: 13 June 2001     

Recent applications of density functional theory calculations to biomolecules

The purpose of this overview is to highlight the broad scope and utility of current applications of density functional theory (DFT) methods for the study of the properties and reactions of biomolecules. This is illustrated using examples selected from research carried out within our research group and in collaboration with others. The examples include the hyperfine coupling constants of amino acid radicals, the use of an amino acid as a chiral catalyst for the formation of carbon–carbon bonds in the aldol reaction, hydrogen-bond mediated catalysis of an aminolysis reaction, radiation-induced protein–DNA cross-links, and the mechanism by which an antitumor drug cleaves DNA. We demonstrate that DFT-based methods can be applied successfully to a broad range of problems that remain beyond the scope of conventional electron-correlation methods. Furthermore, we show that contemporary computational quantum chemistry complements experiment in the study of biological systems. Received: 19 December 2001 / Accepted: 8 April 2002 / Published online: 4 July 2002     

Molecular Schrödinger–Riccati calculations. Test for the hydrogen-ion molecule

 The hydrogen-ion molecule has been used as a test system for the application of the Schr?dinger–Riccati formulation to molecular calculations. Some of the points discussed are the characteristics (quasiconstancy of the local energies, size, number of points) of the sampling region to be chosen, the dependence on the starting function, the precision of the calculations, and the excellent behaviour of the predicted function (by comparison to an accurate function). Received: 1 May 2001 / Accepted: 25 July 2001 / Published online: 11 October 2001     

Rydberg character of the higher excited states of free-base porphin

 The Rydberg character of the excited states of free-base porphin (FBP) has been investigated by the ab initio configuration interaction singles (CIS) method and the state-averaged complete-active-space self-consistent-field method. Double-zeta basis sets augmented with s, p, and d Rydberg functions and d polarization functions have been employed. Two types of molecular orbitals sets, the restricted Hartree–Fock molecular orbitals obtained for the ground state (¹A _g ) and for the cation state (²A _u ), have been used in the CIS calculations. All the calculations show that Rydberg-type excitations play important roles especially in the N bands. In this article we propose applying the model of a perturbed Rydberg series to interpret the excited states of FBP. By using this model, we have succeeded in analyzing the characteristics of the excited states as well as the experimental oscillator strengths, which have considerable magnitude even in the higher excited states. Received: 27 November 2000 / Accepted: 11 April 2001 / Published online: 27 June 2001     

Correlated electron-pair properties of the Li atom in momentum space

On the basis of multiconfiguration Hartree–Fock calculations, correlated electron-pair intracule (relative motion) and extracule (center-of-mass motion) properties are reported for the Li atom in momentum space. The present results are more accurate and consistent than those in the literature. Received: 10 September 2001 / Accepted: 11 December 2001 / Published online: 22 March 2002     

Nucleophilic Substitution at Heteroatoms—Identity Substitution Reactions at Phosphorus and Sulfur Centers: Do They Proceed in a Concerted (SN2) or Stepwise (A–E) Way?

The mechanisms of three selected identity substitution reactions at phosphorus and sulfur occurring with stereospecific inversion have been investigated using density functional theory (DFT). The first identity reaction between methoxyl anion and methyl ethylphenylphosphinate 1 reported in 1963 has been shown to proceed in a stepwise fashion according to the addition–elimination (A–E) mechanism involving formation of a pentacoordinate phosphorus intermediate (TBI-1). In contrast, the results of DFT studies of the identity chloride exchange reaction in (ethoxy)ethylphosphonochloridothionate 3 in acetone solution provided evidence that it proceeds synchronously according to the classical Ingold’s S_N2-P mechanism. DFT calculations of the methoxyl–methoxy exchange reaction at sulfur in methyl p-toluenesulfinate 4 catalyzed by trifluoroacetic acid in methanol revealed that it proceeds stepwise (A–E mechanism), involving the formation of the high-coordinate sulfurane intermediate. In both identity transesterification reactions, 1 and 4, the transiently formed trigonal bipyramidal intermediates with the two methoxyl groups occupying apical positions (TBI-1 and TBI-4) have higher free energy barriers for the Berry-type pseudorotation than those for direct decomposition to starting phosphinate and sulfinate ensuring stereospecific inversion of configuration at the phosphinyl and sulfinyl centers. Thus, the DFT method proved its usefulness in the distinction between both mechanisms that are often indistinguishable by kinetic measurements.     

Infrared spectra of CO in absorption and evaluation of radial functions for potential energy and electric dipolar moment

From quantum-chemical calculations of rotational g factor and new experimental measurements of strengths of lines in infrared spectra of vibration–rotational bands v′–0 in absorption, with 1≤v′≤4, of ¹²C¹⁶O, and from analysis of 16,947 frequencies and wave numbers assigned to pure rotational and vibration–rotational transitions within electronic ground state X ¹Σ⁺, including new measurements of band 4–0 of ¹²C¹⁶O, we evaluate radial functions for potential energy and electric dipolar moment, the latter both in polynomial form and as a rational function that has qualitatively correct behaviour under limiting conditions. Received: 8 November 2001 / Accepted: 5 February 2002 / Published online: 14 August 2002     

Probing O-dealkylation and deamination aging processes in tabun-conjugated AChE: a computational study

The mechanisms of the aging process of tabun-conjugated acetylcholinesterase were explored using density functional theory calculations. The free energy surfaces were calculated for O-dealkylation (C–O bond breaking) and deamination (P–N bond breaking) pathways for the aging process of tabun-conjugated acetylcholinesterase as suggested by mass and crystallographic studies. Initially, the calculations were performed using tabun-conjugated serine (SUN) molecule. O-dealkylation mechanism proceeds via one-step S_N2 type process, whereas the deamination process proceeds via two steps addition–elimination reaction at the phosphorus center of SUN molecule. The recent proposal of another deamination mechanism using human butyrylcholinesterase (hBChE) conjugated with N-mono methyl analogue of tabun (TA4) has also been explored (Nachon et al. in Chem Biol Interact 187:44–48, 2010). The rate-determining activation barrier calculated for this deamination mechanism (26.3 kcal/mol) was comparable with O-dealkylation process (26.9 kcal/mol) with B3LYP/6-31+G* level of theory. To examine the influence of catalytic residue His447, additional calculations were performed with imidazole group of His447 residue. The incorporation of imidazole group of catalytic residue His447 showed marked decrease in the free energies of activation for all the studied aging processes of tabun-inhibited serine. The aging mechanisms have been explored with TA4-inhibited serine, and calculated results showed that the deamination with the rearrangement process is markedly preferred in this case, which supports the Nachon et al. proposal based on the crystallographic studies.     

A correlation-consistent basis set for Fe

 A series of correlation-consistent basis sets are developed for Fe. Our best computed ⁵F–⁵D separation in the Fe atom is in excellent agreement with experiment. Our best estimate for the FeCO D ₀ value is in good agreement with experiment. The ⁵Σ⁻–³Σ⁻ separation in FeCO has an error of 3.6 kcal/mol; while the origin of this error is not clear, it is probably not due to the basis set. Received: 5 March 2001 / Accepted: 2 May 2001 / Published online: 9 August 2001     

An overlap criterion for selection of core orbitals

 An overlap criterion is defined that connects the identification of core orbitals in a molecular system, which can be problematic, to that in isolated atoms, which is well defined. This approach has been tested on a variety of troublesome systems that have been identified in the literature, including molecules containing third-row main-group elements, and is shown to remove errors of up to 100 kcal/mol arising from an inconsistent treatment of core orbitals at different locations on a potential-energy surface. For some systems and choices of core orbitals, errors as large as 19 kcal/mol can be introduced even when consistent sets of orbitals are frozen, and the new method is shown to identify these cases of substantial core–valence mixing. Finally, even when there is limited core–valence mixing, the frozen-core approximation can introduce errors of more than 5 kcal/mol, which is much larger than the presumed accuracy of models such as G2 and CBS-QB3. The source of these errors includes interatomic core–core and core–valence dispersion forces. Received: 31 August 2001 / Accepted: 11 October 2001 / Published online: 9 January 2002     

Crystal Structure and Conformation of a 10-Thiabilirubin

Summary.  A crystal structure determination of a bilirubin analog with a sulfur instead of a C(10)–CH₂ linking the two dipyrrinones is reported. Conformation-determining torsion angles and key hydrogen bond distances and angles are compared to those obtained from molecular dynamics calculations as well as to the corresponding data from X-ray determinations and molecular dynamics calculations of bilirubin. Like other bilirubins, the component dipyrrinones of the analog are present in the bis-lactam form with (Z)-configurated double bonds at C(4) and C(15). Despite the large differences in bond lengths and angles at –S–vs.–CH₂–, the crystal structure shows considerable similarity to bilirubin: both pigments adopt a folded, intramolecularly hydrogen-bonded ridge-tile conformation stabilized by six hydrogen bonds – although the interplanar angle of the ridge-tile conformation of the title compound is smaller (∼ 86°) than that of bilirubin (∼ 98°). The collective data indicate that even with long C–S bond lengths and a smaller C–S–C bond angle at the pivot point on the ridge-tile seam, intramolecular hydrogen bonding persists. Received August 16, 2001. Accepted September 12, 2001