Testing electronic structure methods for describing intermolecular H · · · H interactions in supramolecular chemistry

In this article a wide variety of computational approaches (molecular mechanics force fields, semiempirical formalisms, and hybrid methods, namely ONIOM calculations) have been used to calculate the energy and geometry of the supramolecular system 2‐(2′‐hydroxyphenyl)‐4‐methyloxazole (HPMO) encapsulated in β‐cyclodextrin (β‐CD). The main objective of the present study has been to examine the performance of these computational methods when describing the short range H · · · H intermolecular interactions between guest (HPMO) and host (β‐CD) molecules. The analyzed molecular mechanics methods do not provide unphysical short H · · · H contacts, but it is obvious that their applicability to the study of supramolecular systems is rather limited. For the semiempirical methods, MNDO is found to generate more reliable geometries than AM1, PM3 and the two recently developed schemes PDDG/MNDO and PDDG/PM3. MNDO results only give one slightly short H · · · H distance, whereas the NDDO formalisms with modifications of the Core Repulsion Function (CRF) via Gaussians exhibit a large number of short to very short and unphysical H · · · H intermolecular distances. In contrast, the PM5 method, which is the successor to PM3, gives very promising results. Our ONIOM calculations indicate that the unphysical optimized geometries from PM3 are retained when this semiempirical method is used as the low level layer in a QM:QM formulation. On the other hand, ab initio methods involving good enough basis sets, at least for the high level layer in a hybrid ONIOM calculation, behave well, but they may be too expensive in practice for most supramolecular chemistry applications. Finally, the performance of the evaluated computational methods has also been tested by evaluating the energetic difference between the two most stable conformations of the host(β‐CD)‐guest(HPMO) system. © 2003 Wiley Periodicals, Inc. J Comput Chem 25: 99–105, 2004

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

[2]  A. Bondi.  van der Waals Volumes and Radii , 1964 .

[3]  Alistair P. Rendell,et al.  Comparison of semiempirical and ab initio QM decomposition analyses for the interaction energy between molecules , 2002 .

[4]  J. Mayoral,et al.  SOLVENT EFFECTS ON DIELS-ALDER REACTIONS: A SEMI-EMPIRICAL STUDY , 1995 .

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

[6]  Norman L. Allinger,et al.  Conformational analysis. 130. MM2. A hydrocarbon force field utilizing V1 and V2 torsional terms , 1977 .

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

[8]  T. Halgren.  MMFF VII. Characterization of MMFF94, MMFF94s, and other widely available force fields for conformational energies and for intermolecular‐interaction energies and geometries , 1999, Journal of computational chemistry.

[9]  G. Chang,et al.  Macromodel—an integrated software system for modeling organic and bioorganic molecules using molecular mechanics , 1990 .

[10]  Claude Millot,et al.  Improving description of hydrogen bonds at the semiempirical level: water–water interactions as test case , 2000 .

[11]  J. Baker.  An algorithm for the location of transition states , 1986 .

[12]  János G. Ángyán,et al.  The origin of the problems with the PM3 core repulsion function , 1997 .

[13]  J. Pople,et al.  Self‐Consistent Molecular‐Orbital Methods. I. Use of Gaussian Expansions of Slater‐Type Atomic Orbitals , 1969 .

[14]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[15]  J. M. Lluch,et al.  The photoinduced intramolecular proton transfer in 2-(2′-hydroxyphenyl)-4-methyloxazole embedded in β-cyclodextrin , 2002 .

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

[17]  W. Goddard,et al.  UFF, a full periodic table force field for molecular mechanics and molecular dynamics simulations , 1992 .

[18]  D. Truhlar,et al.  Direct dynamics calculations with NDDO (neglect of diatomic differential overlap) molecular orbital theory with specific reaction parameters , 1991 .

[19]  M. Frisch,et al.  Using redundant internal coordinates to optimize equilibrium geometries and transition states , 1996, J. Comput. Chem..

[20]  T. Halgren.  Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94 , 1996, J. Comput. Chem..

[21]  A. Sabljić,et al.  Dual-Level Direct Dynamics of the Hydroxyl Radical Reaction with Ethane and Haloethanes: Toward a General Reaction Parameter Method , 1999 .

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

[23]  A. Douhal,et al.  Probing Hydrophobic Nanocavities in Chemical and Biological Systems with a Fluorescent Proton‐Transfer Dye , 1999 .

[24]  Gábor I. Csonka,et al.  Analysis of the core‐repulsion functions used in AM1 and PM3 semiempirical calculations: Conformational analysis of ring systems , 1993, J. Comput. Chem..

[25]  Manuel F. Ruiz-López,et al.  Basic ideas for the correction of semiempirical methods describing H-bonded systems , 2000 .

[26]  J. Szejtli.  Cyclodextrins and their inclusion complexes , 1982 .

[27]  I. Csizmadia,et al.  Ab initio and density functional study of the conformational space of 1C4 α‐L‐fucose , 1997 .

[28]  David F. Shanno,et al.  An example of numerical nonconvergence of a variable-metric method , 1985 .