Revisiting the internal conformational dynamics and solvation properties of cyclodextrins

Abstracts

Molecular dynamics simulations were used to investigate the internal conformational dynamics and solvation properties of the three natural cyclodextrins, <FONT FACE=Symbol>a-, b</FONT>- and gamma-cyclodextrin in aqueous solution at room temperature. These glucose-derived oligosacharides present a molecular structure that confers them the ability to complex guest molecules and change their physicochemical properties. The structural behavior of cyclodextrins in solution is crucial for their complexation abilities. Analyses of the obtained trajectories show that inter-glucose hydrogen bonds between the secondary hydroxyl groups are present in solution, but show a very dynamical character where alternative hydrogen bonds to water molecules can be formed. Despite the lower hydrophilicity of the cyclodextrins inner-cavities, they were found to be solvated and the number of water molecules inside the cavity roughly doubles per glucose unit added to the ring. The residence times for water molecules inside the cavities are inversely proportional to the cavity size.

cyclodextrins; molecular dynamics; solvation; flexibility

Simulações de dinâmica molecular são utilizadas para investigar a dinâmica conformacional interna e propriedades de solvatação das três ciclodextrinas naturais, alfa-, beta- e gama-ciclodextrina em água e temperatura ambiente. Estes oligossacarídeos derivados da glicose apresentam uma estrutura molecular peculiar que lhes confere a habilidade de formar complexos de inclusão e modificar as propriedades físico-químicas das moléculas complexadas. As características estruturais das ciclodextrinas em solução são determinantes para a complexação. As análises das trajetórias obtidas mostram que ligações de hidrogênio secundárias entre resíduos de glicose adjacentes estão presentes em solução e apresentam um caráter altamente dinâmico, onde ligações de hidrogênio alternativas podem ser formadas com moléculas de água. Apesar do baixo caráter hidrofílico, as cavidades apresentam-se solvatadas e o número de moléculas de água dentro da cavidade é aproximadamente duplicada por unidade de glicose adicionada ao macrociclo. Os tempos de residência para as moléculas de água dentro das cavidades são inversamente proporcionais ao seu respectivo tamanho.

ARTICLE

Revisiting the internal conformational dynamics and solvation properties of cyclodextrins

Cristina S. Pereira^I; André F. de Moura^II; Luiz Carlos G. Freitas^II; Roberto D. Lins^III,*

^ILaboratory of Physical Chemistry, Swiss Federal Institute of Technology Zürich (ETHZ) CH-8093 Zürich, Switzerland

^IIDepartamento de Química, Universidade Federal de São Carlos, 13565-905 São Carlos-SP, Brazil

^IIIComputational Biology and Bioinformatics, Pacific Northwest National Laboratory, Richland, 99352, WA, USA

ABSTRACT

Molecular dynamics simulations were used to investigate the internal conformational dynamics and solvation properties of the three natural cyclodextrins, a-, b- and g-cyclodextrin in aqueous solution at room temperature. These glucose-derived oligosacharides present a molecular structure that confers them the ability to complex guest molecules and change their physicochemical properties. The structural behavior of cyclodextrins in solution is crucial for their complexation abilities. Analyses of the obtained trajectories show that inter-glucose hydrogen bonds between the secondary hydroxyl groups are present in solution, but show a very dynamical character where alternative hydrogen bonds to water molecules can be formed. Despite the lower hydrophilicity of the cyclodextrins inner-cavities, they were found to be solvated and the number of water molecules inside the cavity roughly doubles per glucose unit added to the ring. The residence times for water molecules inside the cavities are inversely proportional to the cavity size.

Keywords: cyclodextrins, molecular dynamics, solvation, flexibility

RESUMO

Simulações de dinâmica molecular são utilizadas para investigar a dinâmica conformacional interna e propriedades de solvatação das três ciclodextrinas naturais, a-, b- e g-ciclodextrina em água e temperatura ambiente. Estes oligossacarídeos derivados da glicose apresentam uma estrutura molecular peculiar que lhes confere a habilidade de formar complexos de inclusão e modificar as propriedades físico-químicas das moléculas complexadas. As características estruturais das ciclodextrinas em solução são determinantes para a complexação. As análises das trajetórias obtidas mostram que ligações de hidrogênio secundárias entre resíduos de glicose adjacentes estão presentes em solução e apresentam um caráter altamente dinâmico, onde ligações de hidrogênio alternativas podem ser formadas com moléculas de água. Apesar do baixo caráter hidrofílico, as cavidades apresentam-se solvatadas e o número de moléculas de água dentro da cavidade é aproximadamente duplicada por unidade de glicose adicionada ao macrociclo. Os tempos de residência para as moléculas de água dentro das cavidades são inversamente proporcionais ao seu respectivo tamanho.

Introduction

Cyclodextrins (CDs) have been extensively studied since their discovery in the late nineteenth century.¹ Their peculiar properties and behavior have interested scientists in many different fields.² These cyclic oligosaccharides, obtained from starch by glucanotransferase degradation, are formed by a certain number of a-D-glucopyranose units linked by a(1 ® 4) glycosidic bonds. This arrangement produces a truncated cone shape with a central cavity relatively hydrophobic and an external surface relatively hydrophilic.^1-3 Naturally occurring CDs, known as a-, b- and g-CD, present six, seven and eight glucose units, respectively.^1-3 The structural features of these saccharides allow them to form inclusion complexes with a great variety of compounds.³ Due to their ability to form host-guest complexes, resulting in significant changes in the physicochemical properties of the guest molecules,^2,3 CDs have been extensively applied in the pharmaceutical field,^4,5 industry,^6,7 separation methodologies^8,9 and as enzyme models for the elucidation of the molecular recognition phenomenon.¹⁰ However, a complete understanding of the driving forces and mechanisms involved in the complexation phenomenon remains elusive.³

Computer simulation methods can provide a detailed picture of the structural aspects controlling molecular behavior, therefore contributing to understand the chemistry of the CDs and the details of the complexation mechanisms.¹¹ In fact, computer based methodologies have been previously applied in order to characterize CDs and the complexation phenomenona.¹¹ Among these methods, geometry optimization protocols have been extensively applied, as a result of its popularization in the last decades.^11-28 Some of these studies have aided the interpretation of nuclear magnetic resonance (NMR) and other experimental techniques in the context of CD complexes.^23-28 However, conformational analyses of energetically stable conformations of the CDs^12-14 and their complexes^15-22 have been the main targeted subjects.

Molecular dynamics (MD) simulations are more time-consuming but provide a better connection between atomistic constitution and macroscopic properties when compared to pure geometry optimizations methods.²⁹ In addition to structural^{30-32,35,36,39,41} and solvation^33,34,37-43 aspects, dynamic information can be monitored. The latter is essential to the understanding of the properties and mechanisms involved in host-guest complexation.^44-53 Even with the considerable number of studies using MD to investigate CDs, long simulations of fully-hydrated CDs presenting details about their intrinsic structural features are rarely reported.^41,43 In this contribution, a comparative detailed description of the CDs internal dynamics and solvation behavior has been addressed by means of explicit-solvent MD simulations.

Methods

Molecular dynamics simulations

All the simulations were performed using the Amber96 force field⁵⁴ in conjunction with a modified (optimized charge set, as described below) version of the GLYCAM98⁵⁵ force field, implemented in the TINKER program.⁵⁶ Newton's equations of motion were integrated based on the Verlet algorithm⁵⁷ using a 1-fs time step. The RATTLE⁵⁸ method was applied to constrain all bond lengths involving hydrogen atoms. The simulations were carried out in the isothermal-isobaric (NPT) ensemble. The temperature was kept at 298 K by coupling the system to a heat bath⁵⁹ with a relaxation time of 0.1 ps. The pressure was kept at 1 atm by coupling the system to a pressure bath via isotropic coordinate scaling⁵⁹ with a relaxation time of 0.5 ps. A cutoff of 1.0 nm was applied to evaluate non-bonded interactions. The Particle Mesh Ewald method (PME)⁶⁰ was employed to account for long range electrostatic interactions. At least 2 ns MD simulations were performed for each CD at 298 K and 1 atm after an initial 100 ps period of heating and equilibration at the same temperature and constant volume. The atomic coordinates were saved every 0.1 ps and the analyses were carried out over the entire production runs, unless otherwise specified.

Molecular systems

A central problem dealing with carbohydrate simulations consists on the parameterization of a force field to correctly describe conformational preferences and intermolecular interactions simultaneously. Several efforts have been made to set up force fields to reproduce the specific features of carbohydrates.^55,61-68

Modeling of CDs, just like other oligosaccharides, involves a particular problem regarding the electrostatic representation. Charge sets in carbohydrate force fields have been often derived based either on small fragments^62,65 or on the glucose molecule.⁵⁵ However, CDs are oligomers of glucose linked by a(1 ® 4) glycosidic bonds. The groups involved in the glycosidic linkage present a different chemical environment that is expected to alter its electrostatic properties. Aiming to remedy this problem, atomic partial charges were re-optimized to account for the specific molecular characteristics of the CDs. Hartree-Fock calculations at the 6-31G* level on an a(1 ® 4) glucose trimer and a restrained electrostatic potential (RESP)⁶⁹ fit based on the optimized geometry were performed using the NWChem4.1⁷⁰ program to compute the partial charges. Values of 0.005 au for the harmonic and 0.001 au for the hyperbolic restraints, respectively, were used for all RESP fittings. The final charge set derived for the constituent atoms of the CDs molecules is compared to the corresponding charges of the GLYCAM98⁵⁵ force field in Table 1. The numbering used to represent the atoms are displayed in Figure 1, which shows a schematic drawing of the a-CD molecule.

The two charge sets are relatively similar. As expected, the most significant differences are found in the atoms near the glycosidic bond. Charge derivation by means of RESP⁶⁹ fitting tends to maintain the dipole moment of a molecule close to its quantum-mechanical value and should lead to a realistic description of intermolecular interactions during the simulation. The charge-derivation procedure was performed according to the Amber96/GLYCAM98^54,55 force fields philosophy in order to preserve compatibility.

The optimized geometry of a modeled a-CD using the re-optimized force field presented excellent structural agreement with an ab initio optimized structure via HF/3-21G* basis set (root mean square atomic position deviation (RMSD) of 0.047 nm between the two structures). The empirical energy minimization was performed using the steepest-descent method implemented in the Tinker package⁵⁶ and the ab initio HF/3-21G* calculation was performed using the Gaussian98 package.⁷¹ In both calculations an energy gradient of 10^-4 kcal mol^-1 Å^-1 was used as the convergence criterion.

The descriptive models for a-, b- and g-CD have been created using the derived point charge set along with the Amber96/GLYCAM98^54,55 force field parameters. The systems consisted of one CD molecule in a cubic box containing previously equilibrated TIP3P⁷² water molecules. The initial coordinates of the CDs were obtained from the Cambridge Crystallographic Database (http://www.ccdc.cam.ac.uk/). In the simulations of a- and b-CD the initial dimensions of the box containing 1450 water molecules were 3.53 × 3.53 × 3.53 nm³. The g-CD simulation involved 1600 water molecules resulting in an initial box of 3.70 × 3.70 × 3.70 nm³. The energy of the initial configuration was further minimized in order to relax unfavorable energetic interactions. All atoms of the systems were represented explicitly, CD molecules were considered fully flexible, with constraints only in the bonds involving hydrogen atoms.

Analysis

Time-dependent RMSD was evaluated for the three CD molecules during the simulations. Two different atom sets were considered: i) all solute atoms and ii) ring backbone, excluding hydrogen atoms and lateral hydroxymethyl and hydroxyl groups. The coordinates from the geometry-optimized structures were used as reference structures for the fitting. The corresponding distributions for all atoms were also calculated. The time series of the two glycosidic torsion angles C1_n-O4_n+1-C4_n+1-C3 _n+1 [y], and O5_n-C1_n-O4_n+1-C4 _n+1 [f] (Figure 1) were evaluated (n refers to any given glucose residue in the macrocycle) for all glucose residues of each CD molecule. The internal hydrogen bonds (H-bonds) between secondary hydroxyl groups (OH2, OH3) in the CD molecules (referred throughout the text as secondary H-bonds) were evaluated considering the time series of the distances between the atoms O2_n...O3_n+1, O2_n...HO3_n+1 and O3_n+1...HO2_n of neighbor glucose residues during the simulation (n refers to any given glucose residue in the macrocycle). A 0.2 and 0.3 nm cutoff distances were used as H-bond definition criteria for hydrogen-oxygen and oxygen-oxygen pairs, respectively. Radial distribution functions, g(r), were evaluated for water around all the oxygen atoms of the CD molecules. The solvation properties of the cavities have been characterized by a number of analyses. The center of the cavity was defined evaluating the geometric center of all the atoms comprising the CD molecule. Radial distribution functions, g(r), and the number of water molecules (coordination number) were calculated as a function of the distance from the center of the cavity. Additionally, the probability of finding water molecules as a function of the distance from the center of the CDs was calculated. From these analyses an estimation of the average radius of the cavity for each CD molecule has been drawn. Residence time or correlation time for the water molecules inside the cavity was calculated during the simulation as a function of an arbitrary radial cutoff (in the range 0.05 to 1.0 nm from the center of the cavities). The distance of some selected water molecules from the cavity center as well as the number of water molecules inside the region defined by the cavity were calculated as a function of the simulation time.

Results and Discussion

Structure of the CDs

After the initial equilibration period of 100 ps no significant drifts were observed in the energy (potential, kinetic and total energy) as well as in the density of the systems (data not shown). The RMSD analyses, performed for a-, b- and g-CD during the simulations, indicated no significant difference among the structural fluctuations of the three molecules in solution, as can be seen from the RMSD time evolution (Figure 2, left panel). The Gaussian-like distributions of the corresponding RMSDs for all atoms (Figure 2, right panel) show that all three simulations have reached structural equilibrium. Both all-atoms and backbone RMSDs present the same profile with a small shift in amplitude indicating that the structural changes in CDs come mainly from the fluctuations in the ring backbone. Positional variations of the hydroxyl and hydroxymethyl groups, normally highly-mobile groups in free glucose, here seem to play a less central role for the overall structural fluctuations of these molecules.

The phenomenon of conformational flexibility in CDs is still nowadays matter of debate.⁷³ For a long while CD molecules have been considered conformationally quite rigid, but this view has undergone revisions.^73,74 Supporting the latter is the ability of CDs to accommodate themselves to the spatial and electronic character of a wide variety of guests. It is hard to conceive such structural adaptation in naturally highly-rigid structures. The highly symmetrical representation of these cyclic carbohydrates is in contrast to their physicochemical plasticity.^1-6,73 Most of the discussions about CDs flexibility have been based on crystallographic information obtained from solid state studies or gas-phase analyses, as discussed in the review by Dodziuk.⁷³ The presence of H-bonds between the neighbor secondary hydroxyl groups OH2 and OH3 is generally considered a factor that enhance the rigidity of these cyclic molecules.^10,75 In studies concerning CDs conformational flexibility, the inter-glucose H-bonds are normally considered to be the dominant factor for the macrocycle rigidity, having influence even on the conformational space defined by y and f.⁷⁵ In order to probe these assumptions, we have performed analyses of both glycosidic torsion angles (y and f) (Table 2) and inter-glucose H-bonds as a function of time (Figure 3).

The glycosidic torsion angles (Table 2) are clearly more flexible than the ring torsion angles, whose values do not deviate too much from the mean values characterizing the stable chair conformation (⁴C₁) of the glucose residues (data not shown for conciseness). The RMSD analyses have shown that structural fluctuations in the CD molecules were defined basically by the fluctuations in the ring backbone (Figure 2). Since deviations from the average values of the internal torsion angles in the glucose residues are very small, it can be concluded that the structural flexibility in CD molecules comes mainly from fluctuations of the glycosidic torsion angles. The average values for y and f of the three CD molecules are shown in Table 2 and are in agreement with the calculated mean values obtained from a statistical analysis of these angles in CD crystal structures. ⁷⁶ Interestingly, average standard deviations for b-CD arguably display a smaller internal motion amplitude than its counter parts. While such observation relies on rather small differences over the analyzed final 2 ns, a recent study combining pulse-field-gradient-spin NMR with MD simulations has reached the same conclusion.⁴³

Distances between the atoms O2_n...HO3_n+1 and O3_n+1...HO2_n in each glucose unit of a-, b- and g-CD, as a function time, are illustrated in Figure 3 (a, b and c, respectively). Inter-glucose H-bonds between the secondary hydroxyl groups (OH2 and OH3) are present most of the time in the different glucose residues of the three CD molecules. However, it is important to notice the dynamical behavior of these H-bonds, with apparently no preferential direction (H-bonds are formed between O2...HO3 as frequently as between O3...HO2). This dynamical behavior, so called flip-flop phenomenon, has been already observed experimentally in the crystalline state of b-CD^77,78 and reported in previous MD simulations.³² Our findings are in agreement with references 32 and 77-78 and suggest that the flip-flop phenomenon may be an intrinsic property of these cyclic molecules. It is worth to point that the dynamics of these interactions are probably influenced by the nature of the solvent. In fact, solvation effects on the conformation of carbohydrates are a well-known phenomenon (widely discussed in references 67 and 68). Very likely, when an internal H-bond is broken almost simultaneously there is the formation of an H-bond with water, hence restoring the energetic equilibrium of the system. Therefore, H-bonds between water molecules and polar groups of the CDs would act contributing to increase its conformational diversity. Furthering the assumption, solvents with little or no possibility for intermolecular H-bonds, might hinder such degrees of freedom to certain preferred configurations.

Hydration and interaction with water

Radial distribution functions, g(r), for water around the oxygen atoms of a-, b- and g-CD molecules are shown in Figure 4 (a, b and c, respectively). g(r) for atoms O2, O3, O5 and O6 present structural features characteristic of h-bonding to water molecules (correlation peaks at 0.18 nm and 0.28 nm for hydrogen and oxygen atoms of water, respectively). These findings confirm our previous statement that energetically stable arrangements in these systems are also possible when H-bonds are formed between water and CD. Nonetheless, the correlation peaks for the water hydrogen around the oxygen atoms of the CDs sometimes present intensities smaller than one, showing that the backbone hinder water molecules from freely approaching the hydrophilic sites. This hindrance becomes so effective for the O4 atoms that no significant correlation was observed between these sites and the surrounding water molecules. The effect is also evident for O5 (located inside the ring structure) and O2 and O3, reflecting the presence of H-bonds between OH2 and OH3 groups.

The g(r) functions calculated for water around the cavity center of the three CDs (Figure 5a) show two regions of correlation, one closer to the center, indicating the existence of water molecules inside the cavity, and another from 0.85 to 1.1 nm, corresponding to the hydration shell molecules outside the cavity. To fully describe the correlation between water and CDs, coordination number and radial probabilities have also been evaluated. The radial probability of finding water around the geometric center of the CD molecules increases until this distance approaches the inner surface of the cavity (Figure 5b, inset graph), since the CD backbone may be pictured as an excluded volume to the water molecules. Accordingly, the probability must increase again after the distance becomes larger than the outer radius of the CD (Figure 5b). Therefore, the first minimum in the radial probability of finding a water molecule around the geometric center of the CD molecules may be considered as a parameter to evaluate the size of their cavities. The radii found in this way were 0.25, 0.29 and 0.37 nm for a-, b- and g-CD, respectively (Figure 5b, inset graph), in good agreement with experimental data obtained from hydrated crystals (0.24-0.27 nm, 0.30-0.33 nm, 0.38-0.42 nm for a-, b- and g-CD, respectively⁷⁵). Note that the predicted radii are towards the lower end value measured in the crystal. Partial solvation and intrinsic flexibility restraints of the crystal are expected to favor intra-molecular interactions in the CDs, by fully-solvating these molecules an increase in the entropy of the hydroxyl groups would probably lead to an average slightly smaller cavity radius. Coordination numbers show that the larger the cavity radius, the more water molecules are found inside it (Figure 5c).

The number of water molecules inside the cavity as a function of the time in the last 1 ns of each simulation (Figure 6, left panels) shows the correlation between the cavity size and the number of guest water molecules. In a-CD, most of the time there are 2 water molecules inside the cavity while in b-CD this number increases to an average of 4. In g-CD a further increase to an average of 8 water molecules inside the cavity is observed.

Given the geometric center of two molecules, it is possible to evaluate the autocorrelation function that determines the probability of a pair of two such molecules to still exist after a time lag t, given that it existed at time t₀.⁷⁹ The autocorrelation function integral is a characteristic time, which may be regarded as a residence time of the molecular dimer. Residence times for water as a function of an arbitrarily-defined radial cutoff from the center of the cavity of the CD molecules are shown in Figure 7. As can be seen, a-CD presents a different behavior when compared to b- and g-CD. Water molecules present a residence time of approximately 98 ps in the spherical region defining the cavity of a-CD (ca. 0.25 nm). This residence time is effectively longer than the residence time observed for water in the cavity of b-CD (27 ps) and g-CD (25 ps) defined approximately by the spherical radii 0.29 nm and 0.37 nm, respectively. These results indicate that the water molecules entering in the cavity of a-CD spend, on average, more time than they would do in the cavities of b- and g-CD. This finding may be related with the size of the cavities, i.e., bigger cavities permit a higher flux of water molecules, as shown in Figures 5c and 6. a-CD cavity holds a smaller number of water molecules that can eventually remain trapped inside the cavity, giving rise to the elevated residence time. It is worth to note that the trapped molecules continue to exhibit relatively high rotational and translational motion inside of the cavity. Probably due to such mobility, we were unable to identify any specific interactions responsible for this observation. The small likelihood of g-CD to trap water molecules within the cavity is also evident and agree with a recently reported experimental study.⁸⁰ The dynamics of the water molecules through the CDs cavities is better appreciated in Figure 6, right panels. The distance between the four water molecules presenting the longest residence times inside the cavity and the geometric center of the CDs is plotted as a function of time during the final 1 ns of each simulation (only the four most stable dimers are shown, i.e., the number of water molecules in the cavities at any given time is actually larger than the one shown in the plot). a-CD presents stable dimmers most of the time, while b- and g-CD have some periods in which more transient dimers were found (Figure 6, right panels).

Conclusions

MD simulations were used to investigate the internal conformational dynamics and solvation of the three natural CDs, a-, b- and g-CD, in water at room temperature. A modified version of Amber96/GLYCAM98 force field was developed by the re-optimization of the partial charges to account for some of the specific molecular characteristics of the CDs. The results obtained have shown that these molecules present high conformational freedom in solution. Overall, flexibility seems to be defined by fluctuations in the glycosidic torsion angles (f and y) and interaction with water. The ring torsion angles do not scatter too much from the average values defining the stable ⁴C₁ chair conformation for the glucose residues. Secondary H-bonds are formed in solution and present a fast dynamic flip-flop character. This observation agrees with other studies (experimental and theoretical) on the structure of CDs either in solution or in the solid state. The molecule of g-CD presents in solution a more symmetric structure when compared to a-, b-CD, as shown by the persistent chain of secondary H-bonds present in this molecule. We speculate that due to solvation effects, the secondary H-bonds do not confer as much rigidity to the CD molecules as in the solid state. If that would be the case, interactions between the CDs and high-polar solvents, such as water, may increase molecular flexibility allowing dynamical rearrangements of the structure with minimum energetic costs. Additional simulations of the CDs in low-dielectric solvents would be required to asses such assumption. The CDs cavities are known to be relatively hydrophobic, however the present simulation results show that the cavities are able to host water molecules, whose number and residence times depend on the size of the cavity. In summary, conformational variability in the molecular structure of CDs is assisted by interactions with water molecules suggesting that flexibility is strongly dependent on the chemical environment. The formation of H-bonds with neighbor glucose residues can significantly affect their conformational changes, and certainly the intrinsic flexibility of the CD molecules has a strong influence on the complexation mechanisms. Additionally, the findings presented here offer ground basis information for the novel design of carbohydrate-derived cyclic compounds and insights into their complexation mechanisms.

Acknowledgments

This work has received partial financial support from FAPESP (grants 99/12359-4, 01/08360-9 and 01/05891-3). CSP is grateful to CAPES for the award of an exchange-student fellowship and to the Swiss Science National Foundation for financial support (grant 21-105397). The NWChem 4.0 computational chemistry package for massively parallel computers used in this study was developed by the High Performance Computational Chemistry group, Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory. Battelle Memorial Institute operates Pacific Northwest National Laboratory for the U.S. Department of Energy.

References

1. Szejtli, J.; Pure Appl. Chem. 2004, 76, 1825.

2. Szejtli, J.; Cyclodextrins and their Inclusion Complexes, Akademiai Klado: Budapest, 1982.

3. Liu, L.; Guo, Q. X.; J. Inclusion Phenom. Macrocyclic Chem. 2002, 42, 1.

4. Challa, R.; Ahuja, A.; Ali, J.; Khar, R. K.; AAPS PharmSciTech. 2005, 6, E329.

5. Singh, M.; Sharma, R.; Banerjee, U. C.; Biotech. Adv. 2002, 20, 341.

6. Duchêne, D.; Cyclodextrins and their Industrial Uses, Editions de Santé: Paris, 1987.

7. Hedges, A. R.; Chem. Rev. 1998, 98, 2035.

8. Schneiderman, E.; Stalcup, A. M.; J. Chromatogr., B 2000, 745, 83.

9. Fanali, S.; J. Chromatogr., A 2000, 875, 89.

10. Szejtli, J. In Comprehensive Supramolecular Chemistry, Cyclodextrins; Szejtli, J.; Osa, T., eds.; Pergamon Elsevier: Oxford, 1996; vol. 3, pp. 189.

11. Lipkowitz, K. B.; Chem. Rev. 1998, 98, 1829.

12. Lipkowitz, K. B.; J. Org. Chem. 1991, 56, 6357.

13. Myles, A. M. C.; Barlow, D. J.; France, G.; Lawrence, M. J.; Biochim. Biophys. Acta 1994, 1199, 27.

14. Reinhardt, R.; Richter, M.; Mager, P. P.; Carbohydr. Res. 1996, 291, 1.

15. Bonnet, P.; Jaime, C.; Morin-Allory, L.; J. Org. Chem. 2001, 66, 689.

16. Sánchez-Ruiz, X.; Ramos, M.; Jaime, C.; J. Mol. Struct. 1998, 442, 93.

17. Beier, T.; Höltje, H. -D.; J. Chromatogr. B. 1998, 708, 1.

18. Hassonville, S. H.; Dive, G.; Evrard, B.; Barillaro, V.; Bertholet, P.; Delattre, L.; Piel, G.; J. Drug Deliv. Sci. Tech. 2004, 14, 357.

19. Dodziuk, H.; Lukin, O.; Nowinski, K. S.; THEOCHEM 2000, 503, 221.

20. Koca, J.; Ludin, M.; Perez, S.; Imberty, A.; J. Mol. Graphics Modell. 2000, 182, 108.

21. Lipkowitz, K. B.; J. Chromatogr. A. 2001, 906, 417.

22. Lino, A. S. C.; Takahata, Y.; Jaime, C.; THEOCHEM 2002, 594, 207.

23. Salvatierra, D.; Jaime, C.; Virgili, A.; Sánchez-Ferrando, F.; J. Org. Chem. 1996, 61, 9578.

24. Hartell, M. G.; Hicks, R.; Bhattacharjee, A. K.; Koser, B. W.; Carvalho, K.; Van Hamont, J. E.; J. Pharm. Sci. 2004, 93, 2076.

25. El-Barghouthi, M. I.; Masoud, N. A.; Al-Kafawein, J. K.; Zughul, M. B.; Badwan, A. A.; J. Inclusion Phenom. Macrocyclic Chem. 2005, 53, 15.

26. Zheng, Y.; Haworth, I. S.; Zuo, Z.; Chow, M. S. S.; Chow, A. H. L.; J. Pharm. Sci. 2005, 94, 1079.

27. Wen, X. H.; Liu, Z. Y.; Zhu, T. Q.; Chem. Phys. Lett. 2005, 405, 114.

28. Bratu, I.; Gavira-Vallejo, J. M.; Hernanz, A.; Biopolymers 2005, 77, 361.

29. Rapaport, D. C.; The Art of Molecular Dynamics Simulation, Cambridge University Press: Cambridge, 1995.

30. Koehler, J. E. H.; Saenger, W.; van Gunsteren, W. F.; Eur. Biophys. J. 1987, 15, 197.

31. Prabhakaran, M.; Harvey, S. C.; Biopolymers 1987, 26, 1087.

32. Koehler, J. E. H.; Saenger, W.; van Gunsteren, W. F.; Eur. Biophys. J. 1988, 16, 153.

33. Koehler, J. E. H.; Saenger, W.; van Gunsteren, W. F.; J. Mol. Biol. 1988, 203, 241.

34. Köhler, J. E. H.; Saenger, W.; Van Gunsteren, W. F.; J. Biomol. Struct. Dyn. 1988, 6, 181.

35. Pozuelo, J.; Madrid, J. M.; Mendicuti, F.; Mattice, W. L.; Comput. Theor. Polym. Sci. 1996, 6, 125.

36. Kozár, T.; Venanzi, C. A.; THEOCHEM 1997, 395-396, 451.

37. Manunza, B.; Deiana, S.; Pintore, M.; Gessa, C.; THEOCHEM 1997, 419, 133.

38. Winkler, R. G.; Fioravanti, S.; Ciccotti, G.; Margheritis, C.; Villa, M.; J. Comput.-Aided Mol. Des. 2000, 14, 659.

39. Braesicke, K.; Steiner, T.; Saenger, W.; Knapp, E. W.; J. Mol. Graphics Modell. 2000, 18, 143.

40. Starikov, E. B.; Bräsicke, K.; Knapp, E. W.; Saenger, W.; Chem. Phys. Lett. 2001, 336, 504.

41. Bonnet, P.; Jaime, C.; Morin-Allory, L.; J. Org. Chem. 2002, 67, 8602.

42. Lawtrakul, L.; Viernstein, H.; Wolschann, P.; Int. J. Pharm. 2003, 256, 33.

43. Naidoo, K. J.; Chen, J. Y. J.; Jansson, J. L. M.; Widmalm, G.; Maliniak, A.; J. Phys. Chem. B 2004, 108, 4236.

44. Mark, A. E.; van Helden, S. P.; Smith, P. E.; Janssen, L. H. M.; van Gunsteren, W. F.; J. Am. Chem. Soc. 1994, 116, 6293.

45. Köhler, J. E. H.; Hohla, M.; Richters, M.; Köning, W. A.; Angew. Chem., Int. Ed. 1992, 31, 319; Köhler, J. E. H.; Hohla, M.; Richters, M.; Köning, W. A.; Chem. Ber. 1994, 127, 119.

46. Lipkowitz, K. B.; Pearl, G.; Coner, B.; Peterson, M. A.; J. Am. Chem. Soc. 1997, 119, 11269.

47. Lipkowitz, K. B.; Coner, R.; Peterson, M. A.; Morreale, A.; Shackelford, J.; J. Org. Chem. 1998, 63, 732.

48. Grigera, J. R.; Caffarena, E. R.; de Rosa, S.; Carbohydr. Res. 1998, 310, 253.

49. Dodziuk, H.; Lukin, O.; Chem. Phys. Lett. 2000, 327, 18.

50. Choi, Y. H.; Yang, C. H.; Kim, H. W.; Jung, S. H.; J. Inclusion Phenom. Macrocyclic Chem. 2001, 39, 71.

51. Varady, J.; Wu. X. W.; Wang, S. M.; J. Phys. Chem. B 2002, 106, 4863.

52. Fermeglia, M.; Ferrone, M.; Lodi, A.; Pricl, S.; Carbohydr. Polym. 2003, 53, 15.

53. Jara, F.; Mascayano, C.; Rezende, M. C.; Tirapegui, C.; Urzua, A.; J. Inclusion Phenom. Macrocyclic Chem. 2006, 54, 95.

54. Cornell, W. D.; Cieplak, P.; Bayly, C. I.; Gould, I. R.; Merz, K.M.; Ferguson, D.M. Spellmeyer, D. C.; Fox, T.; Caldwell, J. W.; Kollman, P. A.; J. Am. Chem. Soc. 1995, 117, 5179.

55. Woods, R. J.; Chappelle, R.; THEOCHEM 2000, 527, 149.

56. Ponder, J.W. TINKER: Software Tools for Molecular Design, version 3.8; Washington University School of Medicine, USA, 2001.

57. Verlet, L.; Phys. Rev. 1967, 159, 98.

58. Andersen, H. C.; J. Comput. Phys. 1983, 52, 24.

59. Berendsen, H. J. C.; Postma, J. P. M.; van Gunsteren, W. F.; di Nola, A.; Haak, J. R.; J. Chem. Phys. 1984, 81, 3684.

60. Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H. and Pedersen, L. G.; J. Chem. Phys. 1995, 103, 8577.

61. Pérez, S.; Kouwijzer, M.; Mazeau, K.; Engelsen, S. B.; J. Mol. Graph. 1996, 14, 307, and references therein.

62. Damm, W.; Frontera, A.; Tirado-Rives, J.; Jorgensen, W. L.; J.Comput. Chem. 1997, 18, 1955.

63. Fabricius, J.; Engelsen, S. B.; Rasmussen, K. J.; Carbohydr. Chem. 1997, 16, 751.

64. Momany, F. A.; Willett, J. L.; Carbohydr. Res. 2000, 326, 194.

65. Kony, D.; Damm, W.; Stoll, S.; van Gunsteren, W. F.; J. Comp. Chem. 2002, 23, 1416.

66. Pérez, S.; Imberty, A.; Engelsen, S. B.; Gruza, J.; Mazeau, K.; Jimenez-Barbero, J.; Poveda, A.; Espinosa, J. -F.; van Eyck, B. P.; Johnson, G.; French, A. D.; Kouwijzer, M. L. C. E.; Grootenuis, P. D. J.; Bernardi, A.; Raimondi, L.; Senderowitz, H.; Durier, V.; Vergoten, G.; Rasmussen, K.; Carbohydr. Res. 1998, 314, 141.

67. Kuttel, M., Brady, J.W., Naidoo, K.J.; J. Comput. Chem. 2002, 23, 1236.

68. Lins, R. D.; Hünenberger P. H.; J. Comput. Chem. 2005, 26, 1400.

69. Bayly, C. I.; Cieplak, P.; Cornell, W. D.; Kollman, P. A.; J. Phys. Chem. 1993, 97, 10269.

70. Harrison, R. J.; Nichols, J. A.; Straatsma, T. P.; Dupuis, M.; Bylaska, E. J.; Fann, G. I.; Windus, T. L.; Apra, E.; Anchell, J.; Bernholdt, D.; Borowski, P.; Clark, T.; Clerc, D.; Dachsel, H.; de Jong, B.; Deegan, M.; Dyall, K.; Elwood, D.; Fruchtl, H.; Glendenning, E.; Gutowski, M.; Hess, A.; Jaffe, J.; Johnson, B.; Ju, J.; Kendall, R.; Kobayashi, R.; Kutteh, R.; Lin, Z.; Littlefield, R.; Long, X.; Meng, B.; Nieplocha, J.; Niu, S.; Rosing, M.; Sandrone, G.; Stave, M.; Taylor, H.; Thomas, G.; van Lenthe, J.; Wolinski, K.; Wong, A.; Zhang, Z.; NWChem, A Computational Chemistry Package for Parallel Computers, Version 4.0; Pacific Northwest National Laboratory: USA, 2000.

71. Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Zakrzewski, V. G.; Montgomery Jr., J. A.; Stratmann, R. E.; Burant, J. C.; Dapprich, S.; Millam, J. M.; Daniels, A. D.; Kudin, K. N.; Strain, M. C.; Farkas, O.; Tomasi, J.; Barone, V.; Cossi, M.; Cammi, R.; Mennucci, B.; Pomelli, C.; Adamo, C.; Clifford, S.; Ochterski, J.; Petersson, G. A.; Ayala, P. Y.; Cui, Q.; Morokuma, K.; Malick, D. K.; Rabuck, A. D.; Raghavachari, K.; Foresman, J. B.; Cioslowski, J.; Ortiz, J. V.; Baboul, A. G.;. Stefanov, B. B.; Liu, G.; Liashenko, A.; Piskorz, P.; Komaromi, I.; Gomperts, R.; Martin, R. L.; Fox, D. J.; Keith, T.; Al-Laham, M. A.; Peng, C. Y.; Nanayakkara, A.; Gonzalez, C.; Challacombe, M.; Gill, P. M. W.; Johnson, B.; Chen, W.; Wong, M. W.; Andres, J. L.; Gonzalez, C.; Head-Gordon, M.; Replogle, E. S.; Pople, J. A.; Gaussian 98, Revision A.7; M. J. Gaussian, Inc.: USA, 1998.

72. Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R.W.; Klein, M. L.; J. Chem. Phys. 1983, 79, 926.

73. Dodziuk, H.; J. Mol. Struct. 2002, 614, 33.

74. Dodziuk, H.; Nowinski, K.; THEOCHEM 1994, 110, 61.

75. Saenger, W.; Jacob, J.; Gessler, K.; Steiner, T.; Hoffmann, D.; Sanbe, H.; Koizumi, K.; Smith, S. M.; Takaha, T.; Chem. Rev. 1998, 98, 1787.

76. Lichtenthaler, F. W.; Immel, S.; Liebigs Ann. 1996, 1, 27.

77. Betzel, C.; Saenger, W.; Hingerty, B. E.; Brown, G.M.; J. Am. Chem. Soc. 1984, 106, 7545.

78. Zabel, V; Saenger, W.; Mason, S. A.; J. Am. Chem. Soc. 1986, 108, 3664.

79. Impey, R. W.; Madden, P. A.; MacDonald, I. R.; J. Phys. Chem. 1983, 87, 5071.

80. Sabadini, E.; Cosgrove, T.; Egidio, F.C.; Carbohydr. Res. 2006, 341, 270.

Received: June 2, 2006

Web Release Date: August 6, 2007

FAPESP helped in meeting the publication costs of this article.

Publication Dates

• Publication in this collection
  25 Oct 2007
• Date of issue
  2007

History

• Accepted
  06 Aug 2007
• Received
  02 June 2006