DOCSLIB.ORG

Energy Analysis of Chemistry for Correct Insertion by DNA Polymerase ␤

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Energy Analysis of Chemistry for Correct Insertion by DNA Polymerase ␤ Energy analysis of chemistry for correct insertion by DNA polymerase ␤ Ping Lin†, Lars C. Pedersen‡, Vinod K. Batra‡, William A. Beard‡, Samuel H. Wilson‡, and Lee G. Pedersen†‡§ †Department of Chemistry, University of North Carolina, Chapel Hill, NC 27599; and ‡Laboratory of Structural Biology, National Institute of Environmental Health Sciences, National Institutes of Health, P.O. Box 12233, Research Triangle Park, NC 27709-2233 Communicated by Robert G. Parr, University of North Carolina, Chapel Hill, NC, July 19, 2006 (received for review March 8, 2006) X-ray crystallographic structures of human DNA polymerase ␤ with transition-state (TS) theory applied to an enzyme-catalyzed nonhydrolyzable analogs containing all atoms in the active site reaction (27) and an insertion rate constant of 10 sϪ1 measured required for catalysis provide a secure starting point for a theo- for gapped DNA (28), the apparent free energy of activation is retical analysis (quantum mechanics͞molecular mechanics) of the estimated to be 16.1 kcal (1 kcal ϭ 4.18 kJ)͞mol at room mechanism of chemistry without biasing of modeling assumptions temperature. Although there is some debate about the extent to as required in previous studies. These structures provide the basis which a conformational change or the chemical reaction is the for a detailed quantum mechanics͞molecular mechanics study of rate-limiting step (29), we may use this value for comparison with the path for the complete transfer of a monophosphate nucleoside our calculated reaction energy barrier. donor to the sugar acceptor in the active site. The reaction is largely The development of combined quantum mechanics (QM) and associative with the main energetic step preceded by proton molecular mechanics (MM) methods provides tools for investi- transfer from the terminal primer deoxyribose O3؅ to Asp-256. The gating the catalytic process for complex biological systems. The key residues that provide electrostatic stabilization of the transi- reactive center, where bonds break and form, is treated with an tion state are identified and compared with those identified by accurate QM approach, whereas the remainder of the enzyme mutational studies. environment is treated with empirical MM (30–32). There have been a number of successful applications of the QM͞MM quantum mechanics ͉ molecular mechanics ͉ reaction mechanism method in the studies of the mechanism of enzyme-catalyzed chemical reactions (refs. 33 and 34 and references therein). In ͞ NA polymerases play a central role in the replication and this study, the QM MM version of the ONIOM method devel- oped by the Morokuma group (35, 36) is used, with the electronic repair of DNA. High-fidelity proofreading-deficient poly- D coupling between the QM and MM regions implemented merases typically produce one error per million nucleotides through an electronic embedding approach. synthesized (1). DNA polymerase ␤ (pol ␤) is the simplest eukaryotic DNA polymerase and catalyzes the template- Results and Discussion directed nucleotidyl transfer reaction during the repair of ‘‘sim- We used a constrained dynamics simulation͞no-constraint min- ple’’ base lesions with moderate fidelity, producing approxi- imization procedure, which resulted in a full-atom model that mately one error per 3,000 nt synthesized during base excision well represents the initial structure before the reaction. The DNA repair (2). DNA pol ␤ has been characterized by x-ray rmsd of the backbone atoms between the minimized structure crystallography in various liganded states (3–8), and the dynam- and the crystal structure was 0.31 Å. The unconstrained mini- ics were characterized by NMR (9, 10) and time-resolved mized structure with the Amber ff99 force field was further fluorescence (11). A number of theoretical studies have been ͞ ␤ optimized at the ONIOM(B3LYP 3–21G*:Amber) and carried out with pol in an effort to understand details of ONIOM(B3LYP͞6–31G*:Amber) levels. A comparison of key catalysis (12–15). Additionally, computer modeling of confor- active-site distances between optimized structures and crystal mational changes deduced from x-ray crystallography before and structure is presented in Fig. 3, which is published as supporting after ligand binding and catalysis have provided insight into information on the PNAS web site. Overall, we found that the substrate discrimination (16–23). These structural and compu- geometry of the active site optimized by the QM͞MM hybrid tational studies can be interpreted in the wealth of kinetic and potential agrees well with the experimental values and signifi- ␤ site-directed mutagenesis studies available for pol (2). cantly better than if optimized with the classical force field. One A recent structural study of the precatalytic complex including of the most notable differences is in the coordination between ␤ Ј pol , DNA template, primer and dNTP analogue, 2 -deoxy- the Mg2ϩ (binding) and the O1A of the ␣-phosphate, which Ј ␣ ␤ uridine-5 -( , )-imido triphosphate, provides a ternary complex changes from a noncoordination distance at 2.79 Å (classical) to Ј structure including a hydroxyl group (O3 ) on the primer a distance at Ϸ2.1 Å (QM͞MM). The positions of some sur- terminus with two bound Mg ions (24). Both Mg ions are found rounding x-ray crystallographic waters and ions also changed to be octahedrally coordinated, with each ion bound firmly to a after simulation. ␣ nonbridging oxygen on the dNTP -phosphate and each with a We also investigated the interactions between the critical 3Ј distinct bound water molecule. The primer O3Ј atom, often hydroxyl group of the primer terminus and the surrounding missing in other x-ray crystal structures and necessary for the environment. The optimized structure with the Amber99 force nucleophilic attack, is securely bound in the inner coordination field suggests a H-bonding interaction involving the hydrogen sphere of the catalytic Mg ion. As all of the vital components for (H3T) of the terminal 3Ј hydroxyl group and the O5Ј on the the reaction are in place in this experimental structure, we have incoming nucleotide triphosphate, an interaction also suggested a secure initial reaction system that minimizes the chance of in an earlier study (15). However, such an initial configuration introducing modeling bias at the outset of a theoretical study. All families of DNA polymerases are believed to catalyze the nucleotidyl transfer reaction through a universal two-metal-ion Conflict of interest statement: No conflicts declared. mechanism (25). Therefore, the detailed reaction pathway we Abbreviations: pol ␤, polymerase ␤; QM, quantum mechanics; MM, molecular mechanics; may derive from this study should provide a general picture of TS, transition state. nucleotidyl transfer reactions. The nucleotide insertion reaction §To whom correspondence should be addressed. E-mail: [email protected]. catalyzed by pol ␤ occurs at a moderate rate (26). Based on © 2006 by The National Academy of Sciences of the USA 13294–13299 ͉ PNAS ͉ September 5, 2006 ͉ vol. 103 ͉ no. 36 www.pnas.org͞cgi͞doi͞10.1073͞pnas.0606006103 Downloaded by guest on October 1, 2021 of H3T does not provide a viable pathway for the critical deprotonation of the 3Ј hydroxyl group. Exploring the potential energy surface with the QM͞MM potential allows us to identify a second energy minimum in which H3T forms a H bond with OD2 of Asp-256. This minimum is Ϸ4 kcal͞mol lower in energy than the energy minimum of the structure with the H3T- O5Ј(dTTP) H bond. This configuration for H3T suggests a direct proton transfer pathway from the 3Ј hydroxyl primer terminus to Asp-256. Proton transfer to an ␣-phosphate oxygen suggested by a cluster model calculation (12) appears unlikely because such a proton transfer pathway will require large structural rearrange- ment of the active site, which is hindered by the surrounding residues in the initial structure. A water-mediated proton trans- fer pathway is also possible if water molecules are able to penetrate and interact with H3T. The proton will subsequently be transferred to the OD1 of Asp-256 (unpublished result) or migrate to the solution (37, 38). Reaction Mechanism. The deprotonation of the terminal hydroxyl group of the growing polynucleotide primer is believed to be an essential initial step for the success of the nucleotidyl transfer reaction. The proton transfer occurred automatically before the formation of O3Ј-P␣ bond when we drove the reaction forward by reducing the O3Ј-P␣ distance. A reverse exploration of the potential energy surface identifies a local minimum of the proton transfer product with an energy 3.5 kcal͞mol higher than the reactant state. The existence of the local minimum of the proton transfer product CHEMISTRY suggests the reaction proceeds by a two-step pathway, in which proton transfer reaction occurs before the nucleotidyl transfer reaction. Because the proton transfer reaction proceeds with a relatively lower energy barrier, the overall activation barrier is determined by the formation of O3Ј-P1A bond and the breaking of P1A-O3A bond. The rate-determining process in such a two-step reaction will have a lower activation barrier than the single-step process and therefore improves the catalytic proficiency. An illus- tration of the potential energy profile is presented in Fig. 4, which BIOCHEMISTRY is published as supporting information on the PNAS web site. The key active-site structures along the reaction path have been ob- tained in our calculations and are illustrated in Fig. 1. The transfer Fig. 1. Critical structures in the QM region determined along the reaction of the proton to Asp-256 significantly weakens the side-chain pathway. (a) Structures of ground state, local minimum state, TS, and product coordination with Mg2ϩ(catalytic) and enhances the ligand-metal state. Blue indicates the sugar of the primer terminus, green indicates Mg ions, interaction between O3Ј and Mg2ϩ(catalytic). The distance be- and yellow indicates phosphates. A nonmetal coordinated water molecule Ј 2ϩ included in the quantum atoms is deleted for clarity.
Load more

Recommended publications
  • UNITS This Appendix Explains Some of the Abbreviations1•2 Used For
    APPENDIX: UNITS This appendix explains some of the abbreviations1•2 used for units in this book and gives conversion factors to SI units and atomic units: length 1 a0 = 1 bohr = 0.5291771 X 10-10 m 1 A= 1 Angstrom= lo-10 m = 1.889727 ao mass 1 me = 1 atomic unit of mass = mass of an electron 9.109534 X 10-31 kg = 5.485803 X 10-4 U 1 u 1 universal atomic mass unit = one twelfth the mass of a 12c atom 1.6605655 x lo-27 kg = 1822.887 me time 1 t Eh 1 = 1 atomic unit of time = 2.418884 x l0-17 s 1 s = 1 second = 4.134137 x 1016 t/Eh temperature 1 K = 1 Kelvin amount of substance 1 mol = 1 mole 6.022045 x 1023 atoms, molecules, or formula units energy 1 cm-1 = 1 wavenumber 1 kayser 1.986477 x lo-23 J 4.556335 x 10-6 Eh 857 858 APPENDIX: UNITS 1 kcal/mol = 1 kilocalorie per mole 4.184 kJ/mol = 1.593601 x 10-3 Eh 1 eV 1 electron volt = 1.602189 x lo-19 J 3.674902 X 10-2 Eh 1 Eh 1 hartree = 4.359814 x lo-18 J Since so many different energy units are used in the book, it is helpful to have a conversion table. Such a table was calculated from the recommended values of Cohen and Taylor3 for the physical censtants and is given in Table 1. REFERENCES 1. "Standard for Metric Practice", American Society for Testing and Materials, Philadelphia (1976).
    [Show full text]
  • A Mathematical and Computational Review of Hartree-Fock SCF
    February 1, 2008 5:24 Molecular Physics introSCF Molecular Physics, Vol. 00, No. 00, DD Month 200x, 1–64 A mathematical and computational review of Hartree-Fock SCF methods in Quantum Chemistry Pablo Echenique∗ and J. L. Alonso †‡ †‡ Departamento de F´ısica Te´orica, Universidad de Zaragoza, † Pedro Cerbuna 12, 50009, Zaragoza, Spain. Instituto de Biocomputaci´on y F´ısica de los Sistemas Complejos (BIFI), ‡ Edificio Cervantes, Corona de Arag´on 42, 50009, Zaragoza, Spain. (May 2007 ) We present here a review of the fundamental topics of Hartree-Fock theory in Quantum Chemistry. From the molecular Hamiltonian, using and discussing the Born-Oppenheimer ap- proximation, we arrive to the Hartree and Hartree-Fock equations for the electronic problem. Special emphasis is placed in the most relevant mathematical aspects of the theoretical deriva- tion of the final equations, as well as in the results regarding the existence and uniqueness of their solutions. All Hartree-Fock versions with different spin restrictions are systematically extracted from the general case, thus providing a unifying framework. Then, the discretization of the one-electron orbitals space is reviewed and the Roothaan-Hall formalism introduced. This leads to a exposition of the basic underlying concepts related to the construction and selection of Gaussian basis sets, focusing in algorithmic efficiency issues. Finally, we close the review with a section in which the most relevant modern developments (specially those related to the design of linear-scaling methods) are commented and linked to the issues discussed. The whole work is intentionally introductory and rather self-contained, so that it may be useful for non experts that aim to use quantum chemical methods in interdisciplinary applications.
    [Show full text]
  • Hydrogen Fluoride and the Thermochemistry of Fluorine
    Jmxutm ' *•* ft 9 fi HI I '! fi£ HH 4 • acL I V'.. ^HH 9 HI jamH ; <v IfaflonaJ Bureau of Standard', Library, E-01 Admin. Bldg. KITED STATES URTMENT OF MAR 2 8 1970 MMERCE LICATION NBS TECHNICAL NOTE 513 8E~ -, I <^sk? Hydrogen Fluoride and the Thermochemistry of Fluorine U.S. :partment OF COMMERCE National Bureau of Standards - NATIONAL BUREAU OF STANDARDS The National Bureau of Standards ' was established by an act of Congress March 3, 1901. Today, in addition to serving as the Nation's central measurement laboratory, the Bureau is a principal focal point in the Federal Government for assuring maximum application of the physical and engineering sciences to the advancement of technology in industry and commerce. To this end the Bureau conducts research and provides central national services in four broad program areas. These are: (1) basic measurements and standards, (2) materials measurements and standards, (3) technological measurements and standards, and (4) transfer of technology. The Bureau comprises the Institute for Basic Standards, the Institute for Materials Research, the Institute for Applied Technology, the Center for Radiation Research, the Center for Computer Sciences and Technology, and the Office for Information Programs. THE INSTITUTE FOR BASIC STANDARDS provides the central basis within the United States of a complete and consistent system of physical measurement; coordinates that system with measurement systems of other nations; and furnishes essential services leading to accurate and uniform physical measurements throughout the Nation's scientific community, industry, and com- merce. The Institute consists of an Office of Measurement Services and the following technical divisions: Applied Mathematics—Electricity—Metrology—Mechanics—Heat—Atomic and Molec- ular Physics—Radio Physics -—Radio Engineering -'—Time and Frequency -—Astro- physics -—Cryogenics.
    [Show full text]
  • One-Dimensional Transport with Equilibrium Chemistry (OTEQ) — a Reactive Transport Model for Streams and Rivers
    One-Dimensional Transport with Equilibrium Chemistry (OTEQ) — A Reactive Transport Model for Streams and Rivers Techniques and Methods Book 6, Chapter B6 U.S. Department of the Interior U.S. Geological Survey Cover. Upper workings of the Pennsylvania Mine in the headwaters of Peru Creek, Colorado. Photograph by Robert L. Runkel, U.S. Geological Survey, September 2009. One-Dimensional Transport with Equilibrium Chemistry (OTEQ): A Reactive Transport Model for Streams and Rivers By Robert L. Runkel Toxic Substances Hydrology Program Techniques and Methods Book 6, Chapter B6 U.S. Department of the Interior U.S. Geological Survey U.S. Department of the Interior KEN SALAZAR, Secretary U.S. Geological Survey Marcia K. McNutt, Director U.S. Geological Survey, Reston, Virginia: 2010 For product and ordering information: World Wide Web: http://www.usgs.gov/pubprod Telephone: 1–888–ASK–USGS For more information on the USGS—the Federal source for science about the Earth, its natural and living resources, natural hazards, and the environment: World Wide Web: http://www.usgs.gov Telephone: 1–888–ASK–USGS Any use of trade, product, or firm names is for descriptive purposes only and does not imply endorsement by the U.S. Government. Although this report is in the public domain, permission must be secured from the individual copyright owners to reproduce any copyrighted materials contained within this report. Suggested citation: Runkel, R.L., 2010, One-dimensional transport with equilibrium chemistry (OTEQ) — A reactive transport model for streams and rivers: U.S. Geological Survey Techniques and Methods Book 6, Chapter B6, 101 p. iii Contents Abstract .
    [Show full text]
  • Neutrons & Catalysis
    NEUTRONS & CATALYSIS by Juergen Eckert and Phillip J. Vergamini atalysis—the ability of some ture of chemicals, could not function catalysts, is still. as it always has been, substances to alter the rate of without catalysts. It has been estimated, largely empirical. chemical reactions without for example, that catalysts are involved In the last few years, however, so- being consumed-was rec- at some point in the production of 60 to phisticated new analytical and compu- c tational techniques have helped stimu- ognized more than 150 years ago and 70 percent of all industrial chemicals. has been applied on an industrial scale Yet the store of knowledge about how late a renaissance in catalysis research. since the beginning of this century. catalysts work is surprisingly small. The Powerful economic forces have mo- Modem industrial chemistry, especially search for a catalyst for a particular re- tivated the study of catalysis as well: petroleum processing and the manufac- action, or for ways to improve existing the need for new sources of energy and Neutrons and Catalysis chemicals, changes in the availability Unfortunately, these details of struc- of’ raw materials, potential restrictions ture and dynamics cannot easily be de- on the availability of noble-metal cat- termined in a “real-world” situation— alysts, and the desire for new products that is, during an actual catalytic reac- have pointed up the need for a clearer tion. Catalytic processes usually pro- understanding of catalytic processes. ceed under conditions that preclude the Much of the research into cataly- direct application of many powerful an- sis is directed toward metals.
    [Show full text]
  • 1 I. Units and Conversions 1 A. Atomic Units 1 B. Energy Conversions 2 II
    1 CONTENTS I. Units and Conversions 1 A. Atomic Units 1 B. Energy Conversions 2 II. Approximation Methods 2 A. Semiclassical quantization 2 B. Time-independent perturbation theory 4 C. Linear variational method 7 D. Orthogonalization 11 1. Gram-Schmidt orthogonalization 11 2. Diagonalization of the overlap matrix 12 E. MacDonald’s Theorem 14 F. DVR method for bound state energies 16 References 21 I. UNITS AND CONVERSIONS A. Atomic Units Throughout these notes we shall use the so-called Hartree atomic units in which mass is reckoned in units of the electron mass (m = 9.1093826 10−31 kg) and distance in terms e × of the Bohr radius (a =5.299175 10−2 nm). In this system of units the numerical values 0 × of the following four fundamental physical constants are unity by definition: Electron mass m • e Elementary charge e • reduced Planck’s constant ~ • Coulomb’s constant 1/4πǫ • 0 The replacement of these constants by unity will greatly simplify the notation. 2 It is easiest to work problems entirely in atomic units, and then convert at the end to SI units, using Length (Bohr radius) 1a =5.299175 10−2 nm = 0.5291775 A˚ • 0 × Energy (Hartree) 1E =4.35974417 10−18 J • h × B. Energy Conversions Atomic (Hartree) units of energy are commonly used by theoreticians to quantify elec- tronic energy levels in atoms and molecules. From an experimental viewpoint, energy levels are often given in terms of electron volts (eV), wavenumber units, or kilocalories/mole (kcal/mol). From Planck’s relation hc E = hν = λ The relation between the Joule and the kilocalorie is 1kcal = 4.184 kJ Thus, 1 kcal/mole is one kilocalorie per mole of atoms, or 4.184 103 J divided by Avogadro’s × number (6.022 1023) = 6.9479 10−21 J/molecule.
    [Show full text]
  • Exercise 1 Q.2 the Gas Is Cooled Such That It Loses 65 J of Heat
    Chemistry | 4.43 Solved Examples JEE Main/Boards WAB= P A( V B −= V A ) nR( TB − T A ) −−11 Q.1 Calculate the increase in internal energy of 1 kg of = (2.0 mol )( 8.3 J K mol) (200K) = 3320 J. water at 100° C when it is converted into steam at the The work done by the gas during the process CA is zero same temperature and at 1 atm (100 kPa). The density as the volume remains constant. From (i), of water and steam are 1000 kg m-3 and 0.6 kg m-3 +=−=− respectively. 3320 J WBC 1200 J 4520 J The latent heat of vaporization of water = 2.25 × 106 J kg−1 Q.3 A mole of a monoatomic ideal gas at 1 atm and 273 Sol: Mass and density of water and steam is provided K is allowed to expand adiabatically against a constant so from the given data calculate the increase in volume. pressure of 0.395 bar until equilibrium is reached. Pressure term is given so by using pressure volume (a) What is the final temperature? relationship calculate the work done by the system. Change in internal energy can be calculated using heat (b) What is the final volume? and work relation. (c) How much Work is done by the gas? The volume of 1 kg of water (d) What is the change in internal energy? 11 = m33 and of 1 kg of steam= m 1000 0.6 Sol: Let the initial and final volumes of the gas be 1V 3 and V2m respectively.
    [Show full text]
  • Free Energies of Formation and Aqueous
    Awrican Mineralogist Vol. 57, pp. 1163-1789 (1972) / FREE ENERGIESOF FORMATION AND AQUEOUS SOLUBILITIESOF ALUMINUM HYDROXIDES AND OXIDE HYDROXIDES AT 25"C GnoncnA. P.e,nrs,Departntant of ApptiedEarth Scienaes,StanJord Uniu ersity, Stanl ord, Cal'it ornia I d305 AssrBAc{ comparison of the standard free energies of formation at 25'c of corundum, diaspore, boehmite, gibbsite, and bayerite estimated from carefully selected uq,ruor6 solubility data with those of calorimetric origin .permits selection of a set which is more consistent with stabilities observed in field and laboratory than the calorimetric set alone. The selection process involves acceptance of -378'2 and -1160 kcal per mole for AG'1 of aALo" and Al"* (aq) respectively and assignmentof -311'0 kcal per mole for the AG'1 of AI(OE); (aq)' Selected values for the solids are: kcal/Mole diaspore -219.5+ 0.5 boehmite -218.7+0.2 gibbsite -275.3+ 0.2 bayerite -274.6+ O.L Solubility-pE curves based on these asignments agree with independent solubility data only if it is assumed that slight supersaturation rezults during many experiments and leads, upon aging toward equilibrium, to precipitation of gibbsite or bayerite, depending on pE. One set of data which suggeststhis inter- pretation is discussed. Isrnopucrrow There is convincing evidencein phase equilibrium investigations that corundum is unsiable with respectto gibbsite in the presenceof of water aL25" and onebar pressure(Kennedy, 1959). Unfortunately, tabulated free energiesof formation (Wagmanet aI., L968,and Robie and Waldbaum, 1968) predict the opposite.During an attempt to find free energiesof formation with which to calculatesolubilities of the aluminum hydroxides,discovery of this discrepancyrevealed the need for critical appraisalof availabledata.
    [Show full text]
  • High-Temperature Relative Enthalpies and Related Thermodynamic Properties of Cui
    PLEASE DO NOT REMOVE FRCM LIBRARY 9074 IRli I Bureau of Mines Report of Investigations/1987 High-Temperature Relative Enthalpies and Related Thermodynamic Properties of CuI By M. J. Ferrante, R. V. Mrazek, and R. R. Brown UNITED STATES DEPARTMENT OF THE INTERIOR Report of Investigations 9074 High-Ten1perature Relative Enthalpies and Related Thermodynamic Properties of CuI By M. J. Ferrante, R. V. Mrazek, and R. R. Brown UNITED STATES DEPARTMENT OF THE INTERIOR Donald Paul Hodel, Secretary BUREAU OF MINES Robert C. Horton, Director library of Congress Cataloging in Publication Data: Ferrante, M. J. (Michael John), 1930- High-temperature relative enthalpies and related thermodynamic properties of CuI. (Report of investigations; 9074) Bibliography: p. 9-10. Supt. of Docs. no.: I 28.23: 9074. 1. Cuprous iodide-Thermal properties. 2. High temperatures. 1. Mrazek, R. V. (Robert V.). II. Brown, R. R. (Robert R.). III. Title. IV. Series: Report of investigations (United States. Bureau of Mines) ; 9074. TN23.U43 [QD181.I11 622 s [546'.6522J 86-607924 CONTENTS Abstract ....•..........•..................•...................•...............•. 1 Introduction ..•..............................................•.•...•............ 2 Materials .............................•...•..................................... 2 Experimental work and resul t s ............••..•.•..........•.•....... iii • G ••••••••• 2 Discussion. 4 References. 9 ILLUSTRATION 1. Variation of mean heat capacity with temperature for CuI(c.t) .••.........••• 7 TABLES 1. Experimental relative
    [Show full text]
  • Ab Initio MO Calculation and Monte Carlo Simulation
    http://www.paper.edu.cn J. Phys. Chem. A 2003, 107, 7945-7951 7945 Effects of Substituents and Solvents on the Reactions of Iminophosphorane with Formaldehyde: Ab Initio MO Calculation and Monte Carlo Simulation Ying Xue and Chan Kyung Kim* Department of Chemistry, Inha UniVersity, Inchon 402-751, Korea ReceiVed: May 20, 2003; In Final Form: August 4, 2003 Ab initio molecular orbital method and Monte Carlo (MC) simulation with free energy perturbation (FEP) techniques have been used to study the aza-Wittig reaction of iminophosphoranes (H3PdNH) with formaldehyde (H2CO) in the gas phase and in three different solvents: water, methanol, and tetrahydrofuran (THF). The optimized structures and thermodynamic properties of stationary points for the title reaction system in the gas phase were calculated at the MP2/6-31G** level of theory. The effects of substituents on the reactivity of iminophosphorane were discussed. This aza-Wittig reaction is more favorable for XdH and CH3 than for XdCl in the gas phase. The potential energy profiles along the minimum energy path in the gas phase and in three solvents were obtained. The solvent effects on the H3PdNH + H2CO reaction increase in the order water ≈ methanol > THF, suggesting that the protic polar solvents are more suitable for the aza-Wittig reaction. 1. Introduction performed ab initio calculations at the MP2/6-31G** level of theory for the reactions of X PdNH with H CO (X ) HorCl) The Wittig reaction of methylenephosphorane with a carbonyl 3 2 to examine the identity of INT1 and INT2 and to gain some compound is one of the most important organic transformations information about these reactions in the gas phase.
    [Show full text]
  • A Mathematical and Computational Review of Hartree-Fock SCF Methods in Quantum Chemistry
    October 18, 2007 13:12 Molecular Physics introSCF View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Digital.CSIC Molecular Physics, Vol. 00, No. 00, DD Month 200x, 1–64 A mathematical and computational review of Hartree-Fock SCF methods in Quantum Chemistry Pablo Echenique∗ and J. L. Alonso †‡ †‡ Departamento de F´ısica Te´orica, Universidad de Zaragoza, † Pedro Cerbuna 12, 50009, Zaragoza, Spain. Instituto de Biocomputaci´on y F´ısica de los Sistemas Complejos (BIFI), ‡ Edificio Cervantes, Corona de Arag´on 42, 50009, Zaragoza, Spain. (May 2007 ) We present here a review of the fundamental topics of Hartree-Fock theory in Quantum Chemistry. From the molecular Hamiltonian, using and discussing the Born-Oppenheimer ap- proximation, we arrive to the Hartree and Hartree-Fock equations for the electronic problem. Special emphasis is placed in the most relevant mathematical aspects of the theoretical deriva- tion of the final equations, as well as in the results regarding the existence and uniqueness of their solutions. All Hartree-Fock versions with different spin restrictions are systematically extracted from the general case, thus providing a unifying framework. Then, the discretization of the one-electron orbitals space is reviewed and the Roothaan-Hall formalism introduced. This leads to a exposition of the basic underlying concepts related to the construction and selection of Gaussian basis sets, focusing in algorithmic efficiency issues. Finally, we close the review with a section in which the most relevant modern developments (specially those related to the design of linear-scaling methods) are commented and linked to the issues discussed.
    [Show full text]