Physical analysis of the diatomic “chemical” energy components

 The recent “chemical energy component analysis” permits the total energy of a molecule to be presented approximately but to good accuracy as a sum of atomic and diatomic energy contributions. Here the diatomic energy components are further decomposed into terms of different physical origin: electrostatics (in point-charge approximation and the distributed charge corrections), exchange effects, diatomic overlap and atomic basis extension terms. This analysis may provide us with a deeper insight into the factors influencing both the chemical bonds and the nonbonded interatomic interactions. Received: 6 May 2002 / Accepted: 13 November 2002 / Published online: 19 March 2003 Acknowledgements. The authors are indebted to the Hungarian Scientific Research Fund for partial financial support (grant no. OTKA T29716). Correspondence to: I. Mayer e-mail: mayer@chemres.hu     

Semiempirical QM/MM potential with simple valence bond (SVB) for enzyme reactions. Application to the nucleophilic addition reaction in haloalkane dehalogenase

 We present a method for the correction of errors in combined QM/MM calculations using a semiempirical Hamiltonian for enzyme reactions. Since semiempirical models can provide a reasonable representation of the general shape of the potential energy surface for chemical reactions, we introduce a simple valence bond-like (SVB) term to correct the energies at critical points on the potential energy surface. The present SVB term is not a stand-alone potential energy function, but it is used purely for introducing small energy corrections to the semiempirical Hamiltonian to achieve the accuracy needed for modeling enzymatic reactions. We show that the present coupled QM-SVB/MM approach can be parameterized to reproduce experimental and ab initio results for model reactions, and have applied the PM3-SVB/MM potential to the nucleophilic addition reaction in haloalkane dehalogenase. In a preliminary energy minimization study, the PM3-SVB/MM results are reasonable, suggesting that it may be used in free energy simulations to assess enzymatic reaction mechanism. Received: 1 November 2001 / Accepted: 6 September 2002 / Published online: 19 February 2003 Contribution to the Proceedings of the Symposium on Combined QM/MM Methods at the 222nd National Meeting of the American Chemical Society, 2001 Correspondence to: Lakshmi S. Devi-Kesavan e-mail: kesavan@chem.umn.edu Acknowledgments. The work is partially supported by the NIH and the NSF.     

A thermodynamic and kinetic study of the formation of C<Subscript>20</Subscript> compounds encapsulating H,He and Ne atoms

 Binding energies of helium, neon and atomic hydrogen encapsulated inside a C₂₀ cage were calculated using an ab initio method at the B3LYP/6-31+G⋆ level of theory. The standard equilibrium constants for the reactions of noble-gas atoms going into the C₂₀ molecular cage have also been studied. The transition states for the reactions of C₂₀ with hydrogen and helium were further obtained with an ab initio method at the B3LYP/6-31+G⋆ level and the rate constants were estimated by using conventional transition-state theory. It was found that the hydrogen and helium atoms are extremely difficult to put into the C₂₀ cage. Once inside the cage, a helium atom can hardly get out, while a hydrogen atom can easily escape from the cage. The results are expected to enrich fullerene science and be helpful for fullerene applications such as storage. Received: 2 November 2002 / Accepted: 19 December 2002 / Published online: 30 April 2003 Correspondence to: R. Q. Zhang, e-mail: aprqz@cityu.edu.hk Acknowledgements. The work described in this paper was jointly supported by a grant from the City University of Hong Kong (project no. 7001222) and a grant from the Research Grants Council of the Hong Kong Special Administrative Region, China (project no. 9040633/CityU, 1011/01P].     

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     

Influence of diffuse and polarization functions on the second-order Møller–Plesset optimized dihedral angle of biphenyl

 Using 6-31G and 6-311G basis sets to which diffuse and polarization functions were added in a stepwise fashion (a total of 16 basis sets), Hartree–Fock (HF), MP2 and B3LYP geometry optimizations were performed on biphenyl. With the MP2 method, diffuse functions raise the dihedral angle φ, for example, from 46.3° for 6-31G to 54.1° for 6-311++G, while polarization functions lower it, for example, from 54.1° for 6-311++G to 42.1° for 6-311++G(2d,2p). For a single set of polarization functions, φ(MP2) lies close to or above φ(HF) (44–47°), but for a double set it is below φ(HF) and is close to B3LYP values (38–42°) which show little basis set dependence. The most reliable value for φ, 42.1° [MP2/6-311++G(2d,2p)], is expected to increase slightly by adding more diffuse functions. The corresponding best calculated energy barrier at 0° (coplanar conformation) is 2.83 kcal/mol, much higher than the experimental estimate (1.4 ± 0.5 kcal/mol). The barrier at 90° is 1.82 kcal/mol, in line with the experimental estimate (1.6 ± 0.5 kcal/mol) and with previous theoretical results. Received: 9 September 2002 / Accepted: 15 November 2002 / Published online: 1 April 2003 Correspondence to: Friedrich Grein e-mail: fritz@unb.ca Acknowledgement. The author would like to thank NSERC (Canada) for financial support.     

Transition states from empirical force fields

 This is an overview of the use of empirical force fields in the study of reaction mechanisms. Empirical-valence-bond-type methods (including reactive force field and multiconfigurational molecular mechanics) produce full reaction surfaces by mixing, in the simplest case, known force fields describing reactants and products. The SEAM method instead locates approximate transition structures by energy minimization along the intersection of the component force fields. The transition-state force-field approach (including Q2MM) designs a new force field mimicking the transition structure as an energy minimum. The scope and applicability of the various methods are compared. Received: 17 April 2002 / Accepted: 26 July 2002 / Published online: 4 November 2002 Correspondence to: P.-O. Norrby e-mail: pon@kemi.dtu.dk     

A DFT study of pyramidalized alkenes: 7-oxasesquinorbornenes and 7,7′-dioxasesquinorbornenes

  DFT calculations of 7′–oxasesquinorbornenes and 7,7′-dioxasesquinorbornenes using the B3LYP/6–31G* method are reported. All the investigated structures (syn- and anti- derivatives) showed significant non-planarity of the central double bond, with the exception of those anti-derivatives possessing symmetrical structures. The influence of the replacement of the methylene groups at position 7- of the norbornene fragment with oxygen and the introduction of second and third (peripheral) double bonds and benzene rings on the molecular and electronic structures of these molecules have also been investigated. Received: 11 November 2002 / Accepted: 6 June 2002 / Published online: 29 April 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     

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     

Role of twisting and sliding on the solvation of a stacked cytosine dimer: an ab initio study

Ab initio calculations with inclusion of correlation effects at the MP2/6-31G* level have been used to predict the interaction energy of stacked cytosine dimer (C/C) as a function of twisting and sliding in the gas phase. Systematic calculations have also been carried out on the solvation free energies of various rotated and translated C/C dimers using a polarized continuum model approach at the HF/6-31G* level with a view to probe the role of various degrees of freedom on the free energy of solvation of the C/C dimer. The interaction energy of the C/C dimer decreases upon changing from a parallel to an antiparallel conformation in the gas phase. The 180°-rotated conformation has been found to be the most stable arrangement when compared to other rotated positions. The rotated and translated dimers exhibit lower solvation free energy than the parallel conformation. The decrease in the dipole moment upon rotation from the parallel to the antiparallel conformation indicates the cancellation of charge distribution upon rotation in the z direction of one cytosine base with respect to the other. The calculation reveals that the present approach could not yield association energy, ΔΔG _Asso, in a solvent medium. This may be due to the fact that in the case of floppy molecules the contribution from translational, rotational and vibrational free energies plays a significant role in the calculation of ΔΔG _Asso. Received: 13 December 2001 / Accepted: 25 March 2002 / Published online: 13 June 2002     

Gradient of the ZAPT2 energy

 An explicit expression for the analytical first derivative of the Z-averaged perturbation theory taken to second order energy, due to Lee and Jayatilaka, is presented for application to high-spin systems described by a restricted open-shell Hartree–Fock wavefunction. The use of frozen core orbitals is incorporated into the derivation. Received: 23 April 2001 / Accepted: 31 August 2001 / Published online: 9 January 2002     

Second-order Møller–Plesset perturbation theory with terms linear in the interelectronic coordinates and exact evaluation of three-electron integrals

 The second-order correlation energy of M?ller–Plesset perturbation theory is computed for the neon atom using a wave function that depends explicitly on the interelectronic coordinates (MP2-R12). The resolution-of-identity (RI) approximation, which is invoked in the standard formulation of MP2-R12 theory, is largely avoided by rigorously computing the necessary three-electron integrals. The basis-set limit for the second-order correlation energy is reached to within 0.1 mE _h. A comparison with the conventional RI-based MP2-R12 method shows that only three-electron integrals over s and p orbitals need to be computed exactly, indicating that the RI approximation can be safely used for integrals involving orbitals of higher angular momentum. Received: 9 May 2001 / Accepted: 31 October 2001 / Published online: 9 January 2002     

Dipole–reaction field interaction model for the solvent reorganization energy and its application to the benzoquinone–benzoquinone anion radical system

 Based on the spherical cavity approximation and the Onsager model, a dipole–reaction field interaction model has been proposed to elucidate the solvent reorganization energy of electron transfer (ET). This treatment only needs the cavity radius and the solute dipole moment in the evaluation of the solvent reorganization energy, and fits spherelike systems well. As an application, the ET reaction between p-benzoquinone and its anion radical has been investigated. The inner reorganization energy has been calculated at the level of MP2/6–31+G, and the solvent reorganization energies of different conformations have been evaluated by using the self-consistent reaction field approach at the HF/6–31+G level. Discussions have been made on the cavity radii and the values are found to be reasonable when compared with the experimental ones of some analogous intramolecular ET reactions. The ET matrix element has been determined on the basis of the two-state model. The fact that the value of the ET matrix element is about 10 times larger than RT indicates that this ET reaction can be treated as an adiabatic one. By invoking the classical Marcus ET model, a value of 4.9 × 10⁷M⁻¹s⁻¹ was obtained for the second-order rate constant, and it agrees quite well with the experimental one. Received: 19 October 2001 / Accepted: 17 January 2002 / Published online: 3 May 2002     

Signatures of reactive resonance: three case studies

 Recent advances in the study of short-lived reactive resonances are reviewed. Special emphasis is given to addressing the question of how reactive resonances might be observed in molecular beam scattering experiments. Three case studies are presented for simple triatomic systems that are believed to exhibit resonance phenomena: F + HD → D + HF, F + H₂ → H + HF, and H + HD → D + H₂. It is seen that reactive resonances do strongly influence collision observables, but in a different way for each case. At this stage, there does not appear to be a unique resonance signature that can be applied to all reactions. Received: 2 March 2002 / Accepted: 2 June 2002 / Published online: 17 October 2002 Acknowledgements. For their valuable input into the research discussed here, we are grateful to our collaborators, K. Liu, S.H. nk;Lee, F. Dong, D.E. Manolopoulos, D, Skouteris, X. Yang, S. Harich, D. Dai, and C.C. Wang. This work was supported by the National Science Foundation. Correspondence to: R. T. Skodje e-mail: skokje@spot.colorado.edu     

NMR chemical shifts in the low–pH form of a-chymotrypsin. A QM/MM and ONIOM–NMR study

 The relationship between hydrogen bonding and NMR chemical shifts in the catalytic triad of low-pH α-chymotrypsin is investigated by combined use of the effective fragment potential [(2001) J Phys Chem A 105:293] and ONIOM–NMR [(2000) Chem Phys Lett 317:589] methods. Our study shows that while the His57 N^δ1−H bond is stretched by a relatively modest amount (to about 1.060 ?) this lengthening, combined with the polarization due to the molecular environment, is sufficient to explain the experimentally observed chemical shifts of 18.2 ppm. Furthermore, the unusual down-field shift of H^ɛ1 (9.2 ppm) observed experimentally is reproduced and shown to be induced by interactions with the C=O group of Ser214 as previously postulated. The free-energy cost of moving H^δ1 from His57 to Asp102 is predicted to be 5.5 kcal/mol. Received: 26 September 2001 / Accepted: 6 September 2002 / Published online: 21 January 2003 Contribution to the Proceedings of the Symposium on Combined QM/MM Methods at the 222nd National Meeting of the American Chemical Society, 2001 Correspondence to: J. H. Jensen e-mail: jan-jensen@uiowa.edu Acknowledgements. This work was supported by a Research Innovation Award from the Research Corporation and a type G starter grant from the Petroleum Research Fund. The calculations were performed on IBM RS/6000 workstations obtained through a CRIF grant from the NSF (CHE-9974502) and on supercomputers at the National Center for Supercomputer Applications at Urbana-Champaign. The authors are indebted to Visvaldas Kairys for help with the CHARMM program, and to Daniel Quinn for many helpful discussions.     

Comments on relativistic basis sets

 Adding the tight and diffuse Gaussian-type functions (GTFs), Faegri's variationally determined double-zeta-quality basis sets for molecular relativistic calculations are examined. An example atom is Cm. When the tight s-type GTF is added the total energy increases, whereas when diffuse GTFs are added the total energy decreases. The reasons for these findings are clarified. It is also pointed out that not only the Faegri's sets but also other variationally determined basis sets would show similar behavior so far as the expansion terms are not sufficient. Received: 22 July 2002 / Accepted: 21 October 2002 / Published online: 31 January 2003 Correspondence to: H. Tatewaki e-mail: htatewaki@nsc.nagoya.cu.ac.jp     

Long-range and short-range Coulomb correlation effects as simulated by Hartree–Fock,local density approximation,and generalized gradient approximation exchange functionals

 Exchange functionals used in density functional theory (DFT) are generally considered to simulate long-range electron correlation effects. It is shown that these effects can be traced back to the self-interaction error (SIE) of approximate exchange functionals. An analysis of the SIE with the help of the exchange hole reveals that both short-range (dynamic) and long-range (nondynamic) electron correlation effects are simulated by DFT exchange where the local density approximation (LDA) accounts for stronger effects than the generalized gradient expansion (GGA). This is a result of the fact that the GGA exchange hole describes the exact exchange hole close to the reference electron more accurately than the LDA hole does. The LDA hole is more diffuse, thus leading to an underestimation of exchange and stronger SIE effects, where the magnitude of the SIE energy is primarily due to the contribution of the core orbitals. The GGA exchange hole is more compact, which leads to an exaggeration of exchange in the bond and the nonbonding region and negative SIE contributions. Partitioning of the SIE into intra-/interelectronic and individual orbital contributions makes it possible to test the performance of a given exchange functional in different regions of the molecule. It is shown that Hartree–Fock exchange always covers some long-range effects via interelectronic exchange while self-interaction-corrected DFT is lacking these effects. Received: 25 May 2002 / Accepted: 7 October 2002 / Published online: 21 January 2003 Correspondence to: E. Kraka e-mail: kraka@theoc.gn.se Acknowledgements. This work was supported financially by the Swedish Natural Science Research Council (NFR). Calculations were done on the supercomputers of Nationellt Superdatorcentrum (NSC), Link?ping, Sweden. The authors thank the NSC for a generous allotment of computer time.     

Aconityl-derived polymers for biomedical applications. Modeling study of cis–trans isomerisation

 Soluble polymers have been prepared that are designed to undergo enhanced rates of hydrolysis at pH values less than that observed in blood circulation. The degradable element in the polymer mainchain is derived from cis-aconityl acid and is defined by a carboxylic acid pendent functionality (C-4) that is cis across a double bond to an amide at C-1 in the polymer mainchain. While degradation studies in vitro have confirmed these polymers do undergo enhanced rates of degradation at acidic pH values, there is also increasing evidence that during the degradation process the double bond isomerises to the trans configuration and thus prevents the full degradation of a polymer. From a molecular modelling perspective we are seeking to understand the propensity for this cis–trans isomerisation and the mechanism of this cis–trans isomerisation is discussed. Received: 29 April 2002 / Accepted: 6 September 2002 / Published online: 14 February 2003