Molecular recognition in different environments: β-cyclodextrin dimer formation in organic solvents.

Electrostatic and van der Waals interactions as well as entropy contribute to the energetics governing macromolecular complexation in biomolecules. Hydrogen bonds play a particularly important role in such interactions. Here we use molecular dynamics (MD) simulations to investigate the hydrogen bond (HB) orientations of free beta-cyclodextrin (β-CD) and head-to-head dimerization of β-CD monomers with and without guest molecules in different environments, namely, in 10 different solvents covering a wide range of polarity. Potentials of mean force for the dimer dissociation are derived from umbrella sampling simulations, allowing determination of the binding affinity between monomers. The HB orientations are in good agreement with available experimental data in water and dimethyl sulfoxide, yielding confidence in the force field used. HB exchanges at the secondary rim of β-CD are observed with a fast rate in water and with a low rate or even no exchange in other solvents. Orientational preferences of interglucopyranose HBs and their effects on the β-CD structure in these solvents are discussed in detail. Polar solvents with stronger HB accepting abilities can interrupt intermolecular HBs more easily, resulting in a less stable dimer. Guest molecules included in the channel-type cavity strengthen the binding affinity between two monomers to some extent, particularly in polar solvents. Formation of the head-to-head dimer is therefore solvent-dependent and guest-modulated. There is only limited correlation between the dimer binding energies and solvent properties like the dielectric constant. This implies that implicit solvent models will not be capable of predicting important properties like binding energy for other solvents than water without a complete reparameterization. This work provides a deeper comprehension on the properties of β-CD, and implications for the application of cyclodextrins in aqueous and nonaqueous media are discussed.

[1]  Igor V. Tetko,et al.  Prediction of n-Octanol/Water Partition Coefficients from PHYSPROP Database Using Artificial Neural Networks and E-State Indices , 2001, J. Chem. Inf. Comput. Sci..

[2]  P. Kollman,et al.  Atomic charges derived from semiempirical methods , 1990 .

[3]  Starr,et al.  Hydrogen-bond dynamics for the extended simple point-charge model of water , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

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

[5]  The role of intra- and intermolecular hydrogen bonds in the formation of β-cyclodextrin head-to-head and head-to-tail dimers. The results of ab initio and semiempirical quantum-chemical calculations , 2001 .

[6]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[7]  Wim E Hennink,et al.  Cyclodextrin-based polymeric materials: synthesis, properties, and pharmaceutical/biomedical applications. , 2009, Biomacromolecules.

[8]  Igor V. Tetko,et al.  Application of Associative Neural Networks for Prediction of Lipophilicity in ALOGPS 2.1 Program , 2002, J. Chem. Inf. Comput. Sci..

[9]  G. Cravotto,et al.  New asymmetrical per-substituted cyclodextrins (2-O-methyl-3-O-ethyl- and 2-O-ethyl-3-O-methyl-6-O-t-butyldimethylsilyl-beta-derivatives) as chiral selectors for enantioselective gas chromatography in the flavour and fragrance field. , 2010, Journal of chromatography. A.

[10]  M. Skiba,et al.  Synthesis of per-substituted hydrophilic and hydrophobic β-cyclodextrin derivatives , 2011 .

[11]  L. Morin-Allory,et al.  Structure and Thermodynamics of α-, β-, and γ-Cyclodextrin Dimers. Molecular Dynamics Studies of the Solvent Effect and Free Binding Energies† , 2002 .

[12]  L. Morin-Allory,et al.  α-, β-, and γ-Cyclodextrin Dimers. Molecular Modeling Studies by Molecular Mechanics and Molecular Dynamics Simulations , 2001 .

[13]  Serdar Kuyucak,et al.  Affinity and selectivity of ShK toxin for the Kv1 potassium channels from free energy simulations. , 2012, The journal of physical chemistry. B.

[14]  H. Ritter,et al.  Cyclodextrin-Modified Polyesters from Lactones and from Bacteria: An Approach to New Drug Carrier Systems , 2011 .

[15]  Cleber P. A. Anconi,et al.  Structure and stability of (alpha-CD)3 aggregate and OEG@(alpha-CD)3 pseudorotaxane in aqueous solution: a molecular dynamics study. , 2009, The journal of physical chemistry. B.

[16]  M. C. Feiters,et al.  SELECTIVE FUNCTIONALIZATION AND FLEXIBLE COUPLING OF CYCLODEXTRINS AT THE SECONDARY HYDROXYL FACE , 1995 .

[17]  Guangying Chen,et al.  Crystal Structures of β‐Cyclodextrin Inclusion Complexes with 7‐Hydroxycoumarin and 4‐Hydroxycoumarin and Substituent Effects on Inclusion Geometry , 2011 .

[18]  Y. Inoue,et al.  1H NMR Study of Intramolecular Hydrogen-Bonding Interaction in Cyclodextrins and Their Di-O-methylated Derivatives , 1988 .

[19]  Carl Caleman,et al.  GROMACS molecule & liquid database , 2012, Bioinform..

[20]  W. B. Almeida,et al.  Theoretical Study of the a-Cyclodextrin Dimer , 2005 .

[21]  Akira Harada,et al.  Cyclodextrin-based supramolecular polymers. , 2009, Chemical Society reviews.

[22]  C. Washington,et al.  Structural Behaviour of 2-Hydroxypropyl-β-Cyclodextrin in Water: Molecular Dynamics Simulation Studies , 2008, Pharmaceutical Research.

[23]  K. Yannakopoulou,et al.  β-Cyclodextrin trimers enclosing an unusual organization of guest: The inclusion complex β-cyclodextrin/4-pyridinealdazine , 2007 .

[24]  Christian Betzel,et al.  Topography of cyclodextrin inclusion complexes, part 20. Circular and flip-flop hydrogen bonding in .beta.-cyclodextrin undecahydrate: a neutron diffraction study , 1984 .

[25]  A. Coleman,et al.  Aggregation of cyclodextrins: An explanation of the abnormal solubility ofβ-cyclodextrin , 1992 .

[26]  D. van der Spoel,et al.  Large influence of cholesterol on solute partitioning into lipid membranes. , 2012, Journal of the American Chemical Society.

[27]  Y. Ishimaru,et al.  Synthesis of secondary face-to-face cyclodextrin dimers linked at each 2-position , 1997 .

[28]  Justin A. Lemkul,et al.  Assessing the stability of Alzheimer's amyloid protofibrils using molecular dynamics. , 2010, The journal of physical chemistry. B.

[29]  G. Torrie,et al.  Monte Carlo free energy estimates using non-Boltzmann sampling: Application to the sub-critical Lennard-Jones fluid , 1974 .

[30]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[31]  M. Costas,et al.  Similarities and differences between cyclodextrin-sodium dodecyl sulfate host-guest complexes of different stoichiometries: molecular dynamics simulations at several temperatures. , 2010, The journal of physical chemistry. B.

[32]  C. Cézard,et al.  Molecular dynamics studies of native and substituted cyclodextrins in different media: 1. Charge derivation and force field performances. , 2011, Physical chemistry chemical physics : PCCP.

[33]  Kenny B. Lipkowitz,et al.  Applications of Computational Chemistry to the Study of Cyclodextrins. , 1998, Chemical reviews.

[34]  A. Harada,et al.  STRUCTURES OF POLYROTAXANE MODELS , 1997 .

[35]  R. Couto,et al.  Mechanisms involved in the antinociceptive effects of 7-hydroxycoumarin. , 2011, Journal of natural products.

[36]  C. Mijangos,et al.  Controlling PVA Hydrogels with γ-Cyclodextrin , 2004 .

[37]  M. Pasquinelli,et al.  Molecular dynamics simulations of interactions between polyanilines in their inclusion complexes with β-cyclodextrins. , 2012, The journal of physical chemistry. B.

[38]  A. Stalcup,et al.  Cyclodextrins: a versatile tool in separation science. , 2000, Journal of chromatography. B, Biomedical sciences and applications.

[39]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

[40]  C. Chipot,et al.  How Do α-Cyclodextrins Self-Organize on a Polymer Chain? , 2012 .

[41]  C. Chipot,et al.  Can the anomalous aqueous solubility of beta-cyclodextrin be explained by its hydration free energy alone? , 2008, Physical chemistry chemical physics : PCCP.

[42]  G. N. Patey,et al.  The free energy of spheres with dipoles: Monte Carlo with multistage sampling , 1973 .

[43]  R J Woods,et al.  Quantum mechanical study of the nonbonded forces in water-methanol complexes. , 2001, The journal of physical chemistry. A.

[44]  J. F. Stoddart,et al.  Working mechanism for a redox switchable molecular machine based on cyclodextrin: a free energy profile approach. , 2010, The journal of physical chemistry. B.

[45]  Karl N. Kirschner,et al.  Solvent interactions determine carbohydrate conformation , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[46]  P. Wolschann,et al.  Molecular dynamics of β-CD in water/co-solvent mixtures , 2011 .

[47]  Bernhard Sellner,et al.  Molecular dynamics simulations of beta-cyclodextrin-aziadamantane complexes in water. , 2008, The journal of physical chemistry. B.

[48]  F. Hirayama,et al.  Cyclodextrin Drug Carrier Systems. , 1998, Chemical reviews.

[49]  S. Bandyopadhyay,et al.  Hydration properties of α-, β-, and γ-cyclodextrins from molecular dynamics simulations. , 2011, The journal of physical chemistry. B.

[50]  A. Maliniak,et al.  Glucose orientation and dynamics in alpha-, beta-, and gamma-cyclodextrins. , 2008, The journal of physical chemistry. B.

[51]  R. Bell,et al.  The electrostatic energy of dipole molecules in different media , 1931 .

[52]  J. Marcy,et al.  Cyclodextrin inclusion complex formation and solid-state characterization of the natural antioxidants alpha-tocopherol and quercetin. , 2009, Journal of agricultural and food chemistry.

[53]  Michael H. Abraham,et al.  Scales of solute hydrogen-bonding: their construction and application to physicochemical and biochemical processes , 2010 .

[54]  A. Harada,et al.  Cyclodextrin-based molecular machines. , 2001, Accounts of chemical research.

[55]  H. Dodziuk.  Rigidity versus flexibility. A review of experimental and theoretical studies pertaining to the cyclodextrin nonrigidity , 2002 .

[56]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[57]  Dong Wang,et al.  Molecular dynamics simulations of the supramolecular assembly between an azobenzene-containing surfactant and α-cyclodextrin: role of photoisomerization. , 2012, The journal of physical chemistry. B.

[58]  G. Gallo,et al.  Conformation of O-methylated amylose and cyclodextrins☆ , 1968 .

[59]  Carl Caleman,et al.  Force Field Benchmark of Organic Liquids: Density, Enthalpy of Vaporization, Heat Capacities, Surface Tension, Isothermal Compressibility, Volumetric Expansion Coefficient, and Dielectric Constant , 2011, Journal of chemical theory and computation.

[60]  Paul J van Maaren,et al.  Thermodynamics of hydrogen bonding in hydrophilic and hydrophobic media. , 2006, The journal of physical chemistry. B.

[61]  Bert L. de Groot,et al.  g_wham—A Free Weighted Histogram Analysis Implementation Including Robust Error and Autocorrelation Estimates , 2010 .

[62]  L. Onsager.  Electric Moments of Molecules in Liquids , 1936 .

[63]  T. Tan,et al.  Investigation of the inclusions of puerarin and daidzin with beta-cyclodextrin by molecular dynamics simulation. , 2010, The journal of physical chemistry. B.

[64]  K. Harata,et al.  Crystal structure of the cyclomaltohexaose (alpha-cyclodextrin) complex with isosorbide dinitrate. Guest-modulated channel-type structure. , 2002, Carbohydrate research.

[65]  A molecular dynamics study of the inclusion complexes of C60 with some cyclodextrins. , 2010, The journal of physical chemistry. B.

[66]  Carl Caleman,et al.  Atomistic simulation of ion solvation in water explains surface preference of halides , 2011, Proceedings of the National Academy of Sciences.

[67]  Carl Caleman,et al.  Organic molecules on the surface of water droplets--an energetic perspective. , 2012, Physical chemistry chemical physics : PCCP.

[68]  J. Szejtli.  Introduction and General Overview of Cyclodextrin Chemistry. , 1998, Chemical reviews.

[69]  T. Aree,et al.  Crystal structure of -cyclodextrinbenzoic acid inclusion complex , 2003 .

[70]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[71]  W. Mattice,et al.  Inclusion Complexes of Chain Molecules with Cycloamyloses. 2. Molecular Dynamics Simulations of Polyrotaxanes Formed by Poly(ethylene glycol) and α-Cyclodextrins , 1997 .

[72]  Y. Takashima,et al.  Double-threaded dimer and supramolecular oligomer formed by stilbene modified cyclodextrin: effect of acyl migration and photostimuli. , 2011, The Journal of organic chemistry.

[73]  Robert J. Woods,et al.  Solvated ensemble averaging in the calculation of partial atomic charges , 2001, J. Comput. Chem..

[74]  Cleber P. A. Anconi,et al.  Ab initio calculations on low-energy conformers of alpha-cyclodextrin. , 2007, Journal of Physical Chemistry A.

[75]  Michael H. Abraham,et al.  Linear solvation energy relationships. 23. A comprehensive collection of the solvatochromic parameters, .pi.*, .alpha., and .beta., and some methods for simplifying the generalized solvatochromic equation , 1983 .

[76]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[77]  T. Tan,et al.  Insight into the structural deformations of beta-cyclodextrin caused by alcohol cosolvents and guest molecules. , 2012, The journal of physical chemistry. B.

[78]  B. L. de Groot,et al.  Potentials of mean force and permeabilities for carbon dioxide, ammonia, and water flux across a Rhesus protein channel and lipid membranes. , 2010, Journal of the American Chemical Society.

[79]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[80]  K. Yannakopoulou,et al.  Similar modes of inclusion in complexes of β-cyclodextrin with sulfonylurea hypoglycemic drugs , 2010 .

[81]  Yong Chen,et al.  Cooperative binding and multiple recognition by bridged bis(beta-cyclodextrin)s with functional linkers. , 2006, Accounts of chemical research.

[82]  A. Maliniak,et al.  Molecular Properties Related to the Anomalous Solubility of β-Cyclodextrin , 2004 .

[83]  W. Saenger.  Crystal packing patterns of cyclodextrin inclusion complexes , 1984 .