Extension of the PDDG/PM3 and PDDG/MNDO semiempirical molecular orbital methods to the halogens

The new semiempirical methods, PDDG/PM3 and PDDG/MNDO, have been parameterized for halogens. For comparison, the original MNDO and PM3 were also reoptimized for the halogens using the same training set; these modified methods are referred to as MNDO′ and PM3′. For 442 halogen‐containing molecules, the smallest mean absolute error (MAE) in heats of formation is obtained with PDDG/PM3 (5.6 kcal/mol), followed by PM3′ (6.1 kcal/mol), PDDG/MNDO (6.6 kcal/mol), PM3 (8.1 kcal/mol), MNDO′ (8.5 kcal/mol), AM1 (11.1 kcal/mol), and MNDO (14.0 kcal/mol). For normal‐valent halogen‐containing molecules, the PDDG methods also provide improved heats of formation over MNDO/d. Hypervalent compounds were not included in the training set and improvements over the standard NDDO methods with sp basis sets were not obtained. For small haloalkanes, the PDDG methods yield more accurate heats of formation than are obtained from density functional theory (DFT) with the B3LYP and B3PW91 functionals using large basis sets. PDDG/PM3 and PM3′ also give improved binding energies over the standard NDDO methods for complexes involving halide anions, and they are competitive with B3LYP/6‐311++G(d,p) results including thermal corrections. Among the semiempirical methods studied, PDDG/PM3 also generates the best agreement with high‐level ab initio G2 and CCSD(T) intrinsic activation energies for SN2 reactions involving methyl halides and halide anions. Finally, the MAEs in ionization potentials, dipole moments, and molecular geometries show that the parameter sets for the PDDG and reoptimized NDDO methods reduce the MAEs in heats of formation without compromising the other important QM observables. © 2003 Wiley Periodicals, Inc. J Comput Chem 1: 138–150, 2004

[1]  K M Merz,et al.  New developments in applying quantum mechanics to proteins. , 2001, Current opinion in structural biology.

[2]  J. Stewart.  Optimization of parameters for semiempirical methods II. Applications , 1989 .

[3]  Henry S. Rzepa,et al.  Ground states of molecules. 53. MNDO calculations for molecules containing chlorine , 1983 .

[4]  James J. P. Stewart,et al.  MOPAC: A semiempirical molecular orbital program , 1990, J. Comput. Aided Mol. Des..

[5]  Henry S. Rzepa,et al.  Ground states of molecules. 40. MNDO results for molecules containing fluorine , 2002 .

[6]  William L. Jorgensen,et al.  A Quantum Mechanical and Molecular Mechanical Method Based on CM1A Charges: Applications to Solvent Effects on Organic Equilibria and Reactions , 1998 .

[7]  J. Stewart.  Optimization of parameters for semiempirical methods I. Method , 1989 .

[8]  M. Dewar,et al.  Ground States of Molecules. 38. The MNDO Method. Approximations and Parameters , 1977 .

[9]  Nathalie Reuter,et al.  Frontier Bonds in QM/MM Methods: A Comparison of Different Approaches , 2000 .

[10]  W. Richards,et al.  Insights into Chorismate Mutase Catalysis from a Combined QM/MM Simulation of the Enzyme Reaction , 1995 .

[11]  J. W. Larson,et al.  Strong hydrogen bonding in gas-phase anions. An ion cyclotron resonance determination of fluoride binding energetics to Broensted acids from gas-phase fluoride exchange equilibrium measurements , 1983 .

[12]  Manuel F. Ruiz-López,et al.  Calibration of the Quantum/Classical Hamiltonian in Semiempirical QM/MM AM1 and PM3 Methods , 2000 .

[13]  James J. P. Stewart,et al.  Ground states of molecules. 67. MNDO Calculations for compounds containing iodine , 1983 .

[14]  Cristiano Ruch Werneck Guimarães,et al.  Contributions of conformational compression and preferential transition state stabilization to the rate enhancement by chorismate mutase. , 2003, Journal of the American Chemical Society.

[15]  M. Schwartz,et al.  Systematic errors in DFT calculations of haloalkane heats of formation , 2003, J. Comput. Chem..

[16]  Michael J. S. Dewar,et al.  Extension of AM1 to the halogens , 1988 .

[17]  Jill E. Gready,et al.  Hybrid Quantum and Molecular Mechanical (QM/MM) Studies on the Pyruvate to l-Lactate Interconversion in l-Lactate Dehydrogenase , 1997 .

[18]  ab Mikhail N. Glukhovtsev,et al.  Gas-Phase Non-Identity SN2 Reactions of Halide Anions with Methyl Halides: A High-Level Computational Study , 1995 .

[19]  Martin J. Field,et al.  A Hybrid-Potential Free-Energy Study of the Isomerization Step of the Acetohydroxy Acid Isomeroreductase Reaction , 2000 .

[20]  Manuel F. Ruiz-López,et al.  MODELING OF PEPTIDE HYDROLYSIS BY THERMOLYSIN. A SEMIEMPIRICAL AND QM/MM STUDY , 1998 .

[21]  William L. Jorgensen,et al.  Improved semiempirical heats of formation through the use of bond and group equivalents , 2002, J. Comput. Chem..

[22]  Charles W. Bauschlicher,et al.  A comparison of the accuracy of different functionals , 1995 .

[23]  Cristiano Ruch Werneck Guimarães,et al.  Investigation of solvent effects for the Claisen rearrangement of chorismate to prephenate: mechanistic interpretation via near attack conformations. , 2003, Journal of the American Chemical Society.

[24]  Shawn T. Brown,et al.  Assessment of Density Functional Theory for Model SN2 Reactions: CH3X + F- (X = F, Cl, CN, OH, SH, NH2, PH2) , 2001 .

[25]  J. W. Larson,et al.  Hydrogen bonding in gas-phase anions. An experimental investigation of the interaction between chloride ion and Broensted acids from ion cyclotron resonance chloride exchange equilibriums , 1984 .

[26]  V. A. Tikhomirov,et al.  A study of the carbon-halogen bond breaking in tert-butyl halides by the PM3 quantum chemical method , 1995 .

[27]  Walter Thiel,et al.  Extension of MNDO to d Orbitals: Parameters and Results for the Halogens , 1992 .

[28]  Eamonn F. Healy,et al.  Development and use of quantum mechanical molecular models. 76. AM1: a new general purpose quantum mechanical molecular model , 1985 .

[29]  M. Field,et al.  A hybrid potential reaction path and free energy study of the chorismate mutase reaction. , 2001, Journal of the American Chemical Society.

[30]  William L. Jorgensen,et al.  PDDG/PM3 and PDDG/MNDO: Improved semiempirical methods , 2002, J. Comput. Chem..

[31]  William L. Jorgensen,et al.  QM/MM Simulations for Diels−Alder Reactions in Water: Contribution of Enhanced Hydrogen Bonding at the Transition State to the Solvent Effect† , 2002 .

[32]  Peter Botschwina,et al.  Coupled cluster calculations for the SN2 reaction Cl− + CH3Br → ClCH3 + Br− , 2000 .

[33]  Kamaldeep K. Chohan,et al.  QM/MM studies show substantial tunneling for the hydrogen-transfer reaction in methylamine dehydrogenase. , 2001, Journal of the American Chemical Society.