Walking on the free energy hypersurface of the 18‐crown‐6 ion system using free energy derivatives

Molecular dynamics simulations have been used to calculate free energy derivatives with respect to nonbonded parameters for both ions and the host, 18‐crown‐6 in cation‐18‐crown‐6 complexes. The free energy derivatives have been used to examine the free energy hypersurface for such complexes and to ‘‘predict’’ the optimal cationic ligand for 18‐crown‐6 complexation in aqueous solution. This method has promise to enable the refinement of the intermolecular interaction parameters for both hosts and guests in aqueous solution. In particular, we show that the optimal ion size for 18‐crown‐6 is calculated to be between the sizes of Na+ and K+, consistent with the experimental data. More generally, it has potential to be a powerful tool in computer aided molecular design.

[1]  J. Dunitz,et al.  1,4,7,10,13,16‐Hexaoxacyclooctadecane: crystal structure at 100 K , 1980 .

[2]  H. Buschmann.  The macrocyclic and cryptate effects. 5. complexation of alkali ions by monocyclic and bicyclic ligands in methanol , 1985 .

[3]  J. Dunitz,et al.  Potassium thiocyanate complex of 1,4,7,10,13,16-hexaoxacyclooctadecane , 1974 .

[4]  P. Kollman,et al.  Molecular mechanics and dynamics studies of crown ether - cation interactions: free energy calculations on the cation selectivity of dibenzo-18-crown-6 and dibenzo-30-crown-10 , 1989 .

[5]  O. Alvarez,et al.  Ion-selective properties of a small ionophore in methanol studied by free energy perturbation simulations , 1992 .

[6]  Peter A. Kollman,et al.  A molecular-mechanics study of 18-crown-6 and its alkali complexes: an analysis of structural flexibility, ligand specificity, and the macrocyclic effect , 1982 .

[7]  T. P. Straatsma,et al.  Treatment of rotational isomers in free energy calculations. II. Molecular dynamics simulation study of 18‐crown‐6 in aqueous solution as an example of systems with large numbers of rotational isomeric states , 1989 .

[8]  J. Åqvist,et al.  Ion-water interaction potentials derived from free energy perturbation simulations , 1990 .

[9]  P. Kollman,et al.  Molecular dynamics studies of calixspherand complexes with alkali metal cations: calculation of the absolute and relative free energies of binding of cations to a calixspherand , 1992 .

[10]  W. C. Still,et al.  An 18-crown-6 derivative with only one conformation , 1993 .

[11]  W. L. Jorgensen,et al.  Modeling the complexation of substituted benzenes by a cyclophane host in water. , 1993, Proceedings of the National Academy of Sciences of the United States of America.

[12]  C. Riche,et al.  Enolates de l'acetylacetate d'ethyle. Structure et reactivite du sel de potassium en presence d'ether couronne et de cryptand , 1978 .

[13]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[14]  D. Beveridge,et al.  Free energy via molecular simulation: applications to chemical and biomolecular systems. , 1989, Annual review of biophysics and biophysical chemistry.

[15]  P. Kollman,et al.  Computational study of the structural, energetic, and acid-base properties of calix[4]arenes , 1990 .

[16]  P. Kollman,et al.  Free energy of association of the 18-crown-6:K+ complex in water: a molecular dynamics simulation , 1990 .

[17]  William L. Jorgensen,et al.  Enhanced view of structure and binding for cyclophane-arene complexes through joint theoretical and experimental study , 1992 .

[18]  J A McCammon,et al.  Theoretical calculation of relative binding affinity in host-guest systems. , 1986, Proceedings of the National Academy of Sciences of the United States of America.

[19]  M Karplus,et al.  Simulation analysis of the stability mutant R96H of T4 lysozyme. , 1991, Biochemistry.

[20]  H. K. Frensdorff.  Stability constants of cyclic polyether complexes with univalent cations , 1971 .

[21]  J. Lehn,et al.  A Polyfunctional Chiral Macrocyclic Polyether Derived from L‐(+)‐Tartaric Acid , 1975 .

[22]  I. Kuntz,et al.  Automated docking with grid‐based energy evaluation , 1992 .

[23]  David A. Pearlman,et al.  Free energy derivatives: A new method for probing the convergence problem in free energy calculations , 1994, J. Comput. Chem..

[24]  Chung F. Wong,et al.  Sensitivity analysis and principal component analysis in free energy calculations , 1991 .

[25]  I. Harada,et al.  Structures of 18-crown-6, 15-crown-5 and their metal complexes in methanol solution as studied by Raman spectroscopy , 1986 .

[26]  Chung F. Wong.  Systematic sensitivity analyses in free energy perturbation calculations , 1991 .

[27]  Wilfred F. van Gunsteren,et al.  An approximate but efficient method to calculate free energy trends by computer simulation: Application to dihydrofolate reductase-inhibitor complexes , 1993, J. Comput. Aided Mol. Des..

[28]  J. V. Eerden,et al.  Molecular dynamics of 18-crown-6 complexes with alkali-metal cations: calculation of relative free energies of complexation , 1988 .

[29]  P. Kollman,et al.  Crown ether-neutral molecule interactions studied by molecular mechanics, normal mode analysis, and free energy perturbation calculations. Near quantitative agreement between theory and experimental binding free energies , 1989 .

[30]  T. Lybrand,et al.  Molecular recognition in nonaqueous solvents: sodium ion, potassium ion, and 18-crown-6 in methanol , 1989 .

[31]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[32]  Ruifeng Liu,et al.  Comment on: Cation vibrational spectroscopy of trans and gauche n‐propylbenzene rotational isomers. Two‐color threshold photoelectron study and ab initio calculations , 1993 .

[33]  P. Kollman,et al.  An all atom force field for simulations of proteins and nucleic acids , 1986, Journal of computational chemistry.

[34]  Chung F. Wong,et al.  Sensitivity analysis of water thermodynamics , 1993 .

[35]  K. M. Doxsee,et al.  Unusual coordination number and geometry in a potassium 18-crown-6 complex , 1992 .

[36]  D. Moras,et al.  CHIRAL AND FUNCTIONALIZED FACE-DISCRIMINATED AND SIDE-DISCRIMINATED MACROCYCLIC POLYETHERS. SYNTHESES AND CRYSTAL STRUCTURES , 1981 .

[37]  Philip E. Gill,et al.  Practical optimization , 1981 .

[38]  R. E. Terry,et al.  Calorimetric titration study of the interaction of several uni- and bivalent cations with 15-crown-5, 18-crown-6, and two isomers of dicyclohexo-18-crown-6 in aqueous solution at 25.degree.C and .mu. = 0.1 , 1976 .

[39]  W. M. Haynes.  CRC Handbook of Chemistry and Physics , 1990 .

[40]  Peter A. Kollman,et al.  Many-body potential for molecular interactions , 1988 .

[41]  D. Cram.  Preorganization—From Solvents to Spherands , 1986 .

[42]  J. D. Lamb,et al.  A systematic study of the effect of macrocycle ring size and donor atom type on the log K, .DELTA.H, and T.DELTA.S of reactions at 25.degree.C in methanol of mono- and divalent cations with crown ethers , 1980 .

[43]  David A. Pearlman,et al.  Determining the contributions of constraints in free energy calculations: Development, characterization, and recommendations , 1993 .

[44]  Irwin D. Kuntz,et al.  A new technique to calculate low-energy conformations of cyclic molecules utilizing the ellipsoid algorithm and molecular dynamics: application to 18-crown-6 , 1988 .