Next Article in Journal
Improving the Performances of Perovskite Solar Cells via Modification of Electron Transport Layer
Next Article in Special Issue
Influence of the Soluble–Insoluble Ratios of Cyclodextrins Polymers on the Viscoelastic Properties of Injectable Chitosan–Based Hydrogels for Biomedical Application
Previous Article in Journal
Anatase TiO2-Decorated Graphitic Carbon Nitride for Photocatalytic Conversion of Carbon Dioxide
Previous Article in Special Issue
Preparation of β-Cyclodextrin Inclusion Complex and Its Application as an Intumescent Flame Retardant for Epoxy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cavity Closure of 2-Hydroxypropyl-β-Cyclodextrin: Replica Exchange Molecular Dynamics Simulations

by
Khanittha Kerdpol
1,
Jintawee Kicuntod
2,
Peter Wolschann
2,3,4,
Seiji Mori
5,
Chompoonut Rungnim
6,
Manaschai Kunaseth
6,
Hisashi Okumura
7,
Nawee Kungwan
1,8,* and
Thanyada Rungrotmongkol
2,9,10,*
1
Department of Chemistry, Faculty of Science, Chiang Mai University, Chiang Mai 50200, Thailand
2
Structural and Computational Biology Research Unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand
3
Department of Pharmaceutical Chemistry, University of Vienna, Vienna 1090, Austria
4
Institute of Theoretical Chemistry, University of Vienna, Vienna 1090, Austria
5
Institute of Quantum Beam Science, Graduate School of Science and Engineering, Ibaraki University, 2-1-1 Bunkyo, Mito, Ibaraki 310-8512, Japan
6
National Nanotechnology Center (NANOTEC), National Science and Technology Development Agency (NSTDA), Pathum Thani 12120, Thailand
7
Research Center for Computational Science, Institute for Molecular Science, Okazaki, Aichi 444-8585, Japan
8
Center of Excellence in Materials Science and Technology, Chiang Mai University, Chiang Mai 50200, Thailand
9
Ph.D. Program in Bioinformatics and Computational Biology, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand
10
Molecular Sensory Science Center, Faculty of Science, Chulalongkorn University, 254 Phayathai Road, Patumwan, Bangkok 10330, Thailand
*
Authors to whom correspondence should be addressed.
Polymers 2019, 11(1), 145; https://doi.org/10.3390/polym11010145
Submission received: 13 December 2018 / Revised: 9 January 2019 / Accepted: 11 January 2019 / Published: 16 January 2019
(This article belongs to the Special Issue Cyclodextrin-Containing Polymers)

Abstract

:
2-Hydroxypropyl-β-cyclodextrin (HPβCD) has unique properties to enhance the stability and the solubility of low water-soluble compounds by inclusion complexation. An understanding of the structural properties of HPβCD and its derivatives, based on the number of 2-hydroxypropyl (HP) substituents at the α-d-glucopyranose subunits is rather important. In this work, replica exchange molecular dynamics simulations were performed to investigate the conformational changes of single- and double-sided HP-substitution, called 6-HPβCDs and 2,6-HPβCDs, respectively. The results show that the glucose subunits in both 6-HPβCDs and 2,6-HPβCDs have a lower chance of flipping than in βCD. Also, HP groups occasionally block the hydrophobic cavity of HPβCDs, thus hindering drug inclusion. We found that HPβCDs with a high number of HP-substitutions are more likely to be blocked, while HPβCDs with double-sided HP-substitutions have an even higher probability of being blocked. Overall, 6-HPβCDs with three and four HP-substitutions are highlighted as the most suitable structures for guest encapsulation, based on our conformational analyses, such as structural distortion, the radius of gyration, circularity, and cavity self-closure of the HPβCDs.

Graphical Abstract

1. Introduction

Cyclodextrins (CDs) are cyclic oligosaccharides, which received much attention in various technical applications, due to their unique properties. In the pharmaceutical industry, they are widely used to improve the stability and the solubility of insoluble drugs in water or organic solvent–water mixtures by molecular encapsulation [1,2,3,4,5,6,7,8,9]. The ability of the encapsulation of CDs with drugs strongly depends on the structural nature of the applied CDs. There are several different kinds of CDs that are defined by number of glucose units, and the most common CDs are α-, β-, and γ-CD, with different number of α-d-glucopyranose units (α = 6, β = 7 and γ = 8). Among these CDs, a derivative of β-cyclodextrin (βCD) named as 2-hydroxypropyl-β-cyclodextrin (HPβCD), shown in Figure 1, has been found to be more soluble and to have a lower toxicity than βCD [10,11,12,13].
HPβCD, a partially substituted poly(hydroxypropyl) ether of βCD, is commercially available as a mixture with a certain range of degrees of substitution [11,13,14,15,16]. The separation of the individual derivatives and isomers from each other is rather difficult, in particular, at the industrial scale. Also, a selective synthesis of the individual isomer is not easy; therefore, a characterization of various derivatives with regard to their inclusion ability is not possible. However, theoretical studies, particularly molecular dynamics simulations, have become a popular tool for providing this important information. By investigating the molecular structural behavior at the atomistic level of each substituted derivative, compounds with a certain degree of substitution of HPβCD can be evaluated. Note that in the technical product, the degree of substitution exists randomly, with different numbers of HP at other positions (e.g., at different glucose units at the 2, 3, and 6 positions), which depends on the concentration of the alkaline solution during the synthesis. Normally, βCD alkylations are observed at the O2 and O6 positions [13,17,18]. Therefore, HPβCD with substitutions at O2 (the most acidic position) and O6 (the most accessible position) with different degrees of substitution have been intensively investigated [11,13,14,15,16,19,20,21,22].
Generally, it is difficult to explore the conformations of biomolecules with complicated free energy surfaces and a larger number of local minima by single conventional simulation at a low temperature condition. To overcome this problem, the replica-exchange molecular dynamics simulations (REMD) [23,24], which is one of the most effective methods with generalized-ensemble algorithms, is applied by non-interacting replica exchange at various temperatures [23,24,25,26,27]. This could allow us to obtain the information of the temperature dependence. In the present study, we applied the REMD method on βCD and HPβCDs models to study the conformational change affected by the different numbers of HP-substitutions on the O2 and/or O6 atoms. The structural behaviors of all models were analyzed based on structural distortion analysis, the radius of gyration, circularity, and the cavity self-closure of the HPβCDs. By doing so, the best candidate of the HPβCDs with a suitable amount of HP-substitution, representing the perfect conical shape for molecular encapsulation, will be provided as a useful guideline for the suitable substitution degree of HPβCD.

2. Computational Methods

The optimized structures of native βCD and Hep6-HPβCD were taken from our previous studies [28,29]. The other HPβCD derivatives were prepared by different numbers of 2-hydroxypropyl (HP) substitutions at the O2 or O6 positions on α-d-glucopyranose units with substitution degrees of around 0.14–1.14 from one to eight HP-substitutions on βCD, as shown in Figure 2. The HP-substitutions on βCD in this work were divided into two groups, single- and double-sided substitutions. For the single-sided HP-substitution (Figure 2a), the HP groups were substituted with from one up to seven HP, only on O6 atoms of the primary rim (called as 6-HPβCDs), because the O6 is more reactive compared to O2 atoms of the secondary rim. In the case of double-sided HP-substitution (Figure 2b), the structures were generated by introducing from one up to four HP groups at the O2 atoms of the non-substituted glucose units (called as 2,6-HPβCD), as defined in Table 1.
Overall, 10 CD structures were generated, to study how different numbers of HP influence the structural behavior, using REMD simulations. A detailed information of the REMD method is given elsewhere [23,24]. The REMD simulations were performed by the Amber 14 package [30]. The parameters of βCD and HPβCDs were taken from the Glycam06 carbohydrate force field [31,32], with the solvation model based on the generalized Born (GB) implicit solvent model, Igb5, which gives a suitable description of cyclodextrin εCD, as reported by Khuntawee et al. [33], and the smaller sizes of CDs (αCD-δCD, unpublished data) relative to the available crystal structures and MD studies in the explicit solvent model. The initial structures of βCD and all HPβCDs were fully minimized with 2000 steps of the steepest descent method, followed by 1000 steps of the conjugated gradient method to relax the structures before simulation. The REMD simulations were performed for 30 ns per replica, including an equilibration step for 5 ns, and the conformations at all temperatures were collected every 1 ps for 25 ns. The temperature distribution and the number of replicas were tested to obtain a reasonable replica exchange simulation. The overlapping between the potential energy distributions of eight replicas is shown in Figure S1 of the Supplementary Materials. The results confirm that the temperature from 269.5 K to 570.9 K, with interval steps of around 30–60 K, was proper for the present case.
The structural distortions of all HPβCDs were analyzed by distances analysis, as defined in Figure 3. Firstly, the distances of adjacent glucopyranose units, d1(i)[O2(i)–O3(i+1)] was defined as the distance between the secondary hydroxyl groups related to the intramolecular hydrogen bonds of the wider CD rim, as labeled in Figure 3a. Secondly, the distance between the glycosidic oxygen atoms, d2(i)[O4(i)–O4(i+1)] is calculated for monitoring the ellipticity of the CDs. The probability distributions of the d1(i) and d2(i) were calculated using Equation (1) in terms of free energy:
F ( x , y ) = k B T   log [ P ( x , y ) ]
where k B is the Boltzmann constant, T is the absolute temperature, and P ( x , y ) is the probability of x for d1(i)[O2(i)–O3(i+1)] and y for d2(i)[O4(i)–O4(i+1)] distances.
Next, the flips of glucose subunits for various conformations were monitored by angle between a pair of adjacent glucose units, θ(i)[C6(i)–C2(i+1)–C6(i+1)], as defined in Figure 3a. Without the flipping, the C6 atoms are in the same site with the adjacent unit. Thus, the θ(i)[C6(i)–C2(i+1)–C6(i+1)] values in no-flip structures are less than 90 degrees. In contrast, when the unit was flipped, the C6 atoms were opposite to those of the adjacent unit, and its θ(i)[C6(i)–C2(i+1)–C6(i+1)] was higher than 90 degrees. In addition, the radius of gyration ( R g ) was calculated. The R g represents the mass-weighted scalar length of each atom from the center of gravity of the molecule, calculated using Equation (2).
R g = i = 1 N m i ( r i r c m ) 2 i = 1 N m i
where N is the number of atoms, m i is the mass of atom i , r i is the Cartesian position vector of atom i , and r c m is the center of mass of CD.
The next parameter is a modified circularity ( C ), which is a dimensionless shape factor based on our previous work [34]. Here, we measured the diameter d3(i)[Cg(Glu(i))–O4(i+3)], which is a distance between the center of mass of the i-th glucose unit (Cg(Glu(i))) and the glycosidic oxygen atom of the opposite glucose unit (O4(i+3)), as depicted in Figure 3b. There are seven d3(i) parameters for each HPβCDs snapshot, since there are seven glucose subunits in βCDs. Hence, C is defined as:
C = min i { 1 7 } ( d 3 ( i ) ) max i { 1 7 } ( d 3 ( i ) )
Note that the value of C is in a range of 0 < C ≤ 1. If C is equal to 1.0, where min(d3(i)) = max(d3(i)), it means that the cavity has a perfectly circular shape. However, when C becomes deviated from 1.0, this is an indication for conformational changes to an elliptical shape. Finally, in order to study the cavity self-closure of HPβCDs, we define the d4(i)[Cg(βCD)–Cg(HP(i))] parameter, which is the distance of the center of mass between the βCD ring and each HP group, as demonstrated in Figure 3c. All parameters were analyzed from 25,000 snapshots taken from REMD simulations. The analyzed results at 300 K will be presented and discussed in the next section, while those of other temperatures are provided in the Supplementary Materials.

3. Results and Discussion

3.1. Structural Analysis

3.1.1. Structural Distortion of Glucose Units in the HPβCDs

To investigate the conformational changes of βCD and all selected HPβCDs, the probability distributions of d1(i)[O2(i)–O3(i+1)] and d2(i)[O4(i)–O4(i+1)] distances for all models at 300 K were calculated in terms of free energy, using Equation (1), and these are plotted in Figure 4a. The contour graphs give the values of the distance probability, where the denoted darkest blue color (ranked from dark red to blue) represents the lowest free energy.
In the case of βCD, three conformational minima (M1, M2, and M3) were detected as shown in Figure 4(a1), which is in agreement with previous results from MD in aqueous solution [29]. The most likely distribution at M1 was found with the d1(i) and d2(i) distances at around 3.5 and 4.5 Å, respectively, describing the almost perfect conical shape of the βCD ring. The second and third main populations were M2 and M3, with d1(i)/d2(i) of around 5.5 Å/4.5 Å and 5.0 Å/5.5 Å, respectively. When comparing βCD with a large-ring CD (εCD) by the REMD method [33], three important minima were found to be similar to the βCD results, but the M3 population of εCD was higher than that of βCD, due to the large size and the high flexibility of the macroscopic ring, having a higher probability of flipping the glucopyranose units in the large-ring system.
From all of the free energy plots of HPβCDs with single- and double-sided HP-substitutions, the d2(i)[O4(i)–O4(i+1)] values fluctuated in the range of 4.0–5.0 Å, implying that the backbone structures of HPβCD were not greatly affected by the HP-substitution groups. In contrast, it was noticeable that d1(i)[O2(i)–O3(i+1)] was in a wide range around 2.5–6.5 Å. Therefore, we plotted the probability of d1(i)[O2(i)–O3(i+1)] and d2(i)[O4(i)–O4(i+1)], as depicted in Figure 4b,c, respectively, in order to compare the conformational minima among the models. For HPβCDs with single-sided HP-substitutions, the most dominant population of d1(i)[O2(i)–O3(i+1)] for the Tri6- and Tet6-HPβCDs was around 3.5 Å, indicating the narrow behavior of their secondary rim. Additionally, the conformational changes of other models were quite similar to the native βCD, except for Hep6-HPβCD with d1(i)[O2(i)–O3(i+1)] at around 4.0-5.0 Å. Consequently, M3 was found in Hep6-HPβCD higher than βCDs and the other 6-HPβCDs, resulting from the steric hindrance of seven HP groups at the primary rim.
In the case of HPβCDs with double-sided HP-substitution, the conformational changes were quite similar to the βCD, in which three main probability distributions were found, following this order: M1 > M2 > M3. When increasing the number of HP at the O2 positions, the M1 population decreased, while the M2 population increased, as depicted in Figure 4(a7–a10), which were related to the distribution of d1(i)[O2(i)–O3(i+1)] (Figure 4(b2)). In addition, the new probability distribution with d1(i)/d2(i) around 4.3 Å/3.5 Å was observed when adding a HP group on the secondary rim, and this probability distribution was enhanced when increasing the number of HP at the O2 positions.
The results indicate that the substitution of the HP groups on narrow and/or wider rims directly affected the distortion in the βCD ring. It is a worth noting that the M1 population was dramatically increased, whilst M2 was decreased and M3 completely disappeared due to complexation with low-water-soluble compounds [35,36,37].

3.1.2. Flipping of the HPβCDs

It is noticeable that the flips of glucose subunits were related directly to the distortion of the CD structure. As mention earlier, a flip of the glucose units was counted when θ(i)[C6(i)–C2(i+1)–C6(i+1)] was higher than 90 degrees. According to the θ(i)[C6(i)–C2(i+1)–C6(i+1)] analysis, the results can be divided into three main different flip-conformations: which are (1) no flip; (2) one flip; (3) two flips. Example snapshots of flipped βCD conformations are shown in Figure 5.
The probabilities of different numbers of flip glucose subunits for all models at 300 K are summarized in Table 2. The no-flip conformation was the main population in all models, while other higher flips (one- and two-flip) were found to be of minor importance during the simulation. In all HPβCDs, more than 58% of no-flip conformations were observed, whereas βCD showed higher values of one- and two-flip conformations. The highest percentages of one- and two-flip conformations in βCD corresponded to the long length of d1(i)[O2(i)–O3(i+1)] detected around 5.0–6.0 Å in the M2 population, as shown in Figure 4(a1). With an increasing number of HP only at O6 positions of glucose subunits, the no-flip population increased by up to 75–78% in Tri6- and Tet6-HPβCDs. Meanwhile, the one-flip population decreased from 35% to 21%. These results were related to the large population of d1(i)[O2(i)–O3(i+1)] at 3.5 Å for Tri6- and Tet6-HPβCDs. For HPβCDs with double-sided HP-substitution, the no-flip population also increased when compared with those of native βCD. However, the one- and two-flip populations of these of double-sided HP-substitutions were higher than those of Tet6-HPβCDs. In addition, the trends of these flips were not significantly different at other temperatures, as shown in Table S1. The percentage of no-flip angles correlates with M1, which is the most populated state during simulation, and the percentages of flip angles relate with M2 and M3, in which the glucopyranose subunits flip and the intramolecular hydrogen bonds of the wider CD rim disappear. Thus, glucose subunits in both 6-HPβCDs and 2,6-HPβCDs have a lower chance to flip (22–31% of the flip angle) than in βCD (42% of flip angle), which is in agreement with their inclusion efficiency reported in our previous study [22].

3.1.3. Radius of Gyration

The flipping of the glucose subunits directly influences the distortion of the macrocyclic ring, as shown in Figure 5 (top view of the different flipped conformations). Therefore, the shapes of all models were investigated in terms of the radius of gyration ( R g ), and the results are shown in Figure 6. The average of R g of native βCD, Mon6-, Di6-, Tri6-, Tet6-, and Hep6-HPβCDs are 6.18, 6.38, 6.39, 6.45, 6.55, and 6.85 Å, respectively, and those of Mon2Tet6-, Di2Tet6-, Tri2Tet6-, and Tet2Tet6-HPβCDs are 6.45, 6.46, 6.65, and 6.85 Å. The results show that the average R g trends to increase when adding more HP groups, which are similar to the R g results from a molecular dynamics simulation, as reported by Yong et al. [17]. However, for Hep6-HPβCD and Tri2Tet6-HPβCD with the same degree of substitution (1.00), the average R g of Hep6-HPβCD was higher than that of Tri2Tet6-HPβCD, because of the higher steric hindrance of the seven HP groups at the primary rim of Hep6-HPβCD. Overall, the R g in the substituted models was observed to be higher than the native βCD, because of the fluctuation and steric hindrance of the HP groups.

3.1.4. Circularity

To determine the effect of the HP-substitutions on the geometry of the CD cavity, we used the circularity ( C ) calculated following Equation (3). The C values of all models at 300 K are listed in Table 3. βCD had the lowest C (0.727 ± 0.087), indicating that its conformational change is higher than those of other models. Moreover, the C values of Tet2Tet6-HPβCDs are also very low (0.742 ± 0.080). The reason for this more flexible conformation may be the result from the higher steric hindrance of HP groups on both the O2 and O6 positions of the glucose subunits. In addition, taking the standard deviation into account, the C values of Tri6- and Tet6-HPβCDs were found to be 0.814 ± 0.075 and 0.815 ± 0.076, respectively. The data suggests that the cavity becomes more circular in shape, when increasing the number of substituting the HP at the O6 primary hydroxyl groups of up to three or four residues, corresponding to a high population of no-flip conformations. Although the population of no-flip conformations of Tri6-, Tet6- and Mon2Tet6-HPβCDs, as well as that of Di2Tet6-HPβCD, are quite similar (higher than 75%), the C values of Mon2Tet6- and Di2Tet6-HPβCDs are lower than the Tri6- and Tet6-HPβCDs, due to the steric hindrance of the HP groups, leading to a distortion of some glucose subunits in the no-flip conformation. For that reason, Tri6- and Tet6-HPβCDs are highlighted as proper structures for forming inclusion complexes with guest molecules, as indicated from their high circularities with a lower possibility of flipping.

3.2. Cavity Self-Closure

The arrangement of the HP groups during simulations influences the CD cavity accessibility directly. Some HP groups can point toward the CD interior, leading to self-closure of the CD cavity. The arrangements of the HP substituents at the glucose subunits were monitored via the distance of the center of mass between the βCD ring and each HP group, called d4(i)[Cg(βCD)–Cg(HP(i))]. The arrangement of HP for Mon6-HPβCD is plotted in Figure 7, while those of the other HPβCDs are shown in Figure S2.
During the REMD simulation of Mon6-HPβCD, almost 98% of HP groups pointed toward the CD exterior, with d4(i)[Cg(βCD)–Cg(HP(i))] around 7–10 Å. For Mon6-HPβCD with d4(i)[Cg(βCD)–Cg(HP(i))] in the rank from 7 to 12 Å, the shape of Mon6-HPβCD seemed to be like a bowl (Figure 7a). When the d4(i)[Cg(βCD)–Cg(HP(i))] was less than 3 Å, the HP group was hindered and rotated into the cavity of the βCD ring, leading to self-closure of the CD cavity (Figure 7c). We also found that the cavity self-closure was related to flipping of the HP-substituted glucose subunits, which triggers HP entrance into the CD cavity (d4(i)[Cg(βCD)–Cg(HP(i))] ~ 1 Å). Some snapshots show that the substituents enter the cavity at the narrow rim with d4(i)[Cg(βCD)–Cg(HP(i))] ~ 2 Å. For that reason, we defined the self-closure of CD cavity or HP occupied in the CD cavity when d4(i)[Cg(βCD)–Cg(HP(i))] < 3 Å. The numbers of d4(i)[Cg(βCD)–Cg(HP(i))] < 3 Å were called n(HPinserted), as illustrated in Equation (4). For all REMD simulations of HPβCDs at 300 K, n(HPinserted) are plotted in Figure 8 with the probability in terms of the percentage of 25,000 snapshots. The results at the various temperatures are summarized in Table S2.
n(HPinserted) = n(d4(i)[Cg(βCD)–Cg(HP(i))] < 3 Å)
For HPβCDs with single-sided HP-substitutions, it was noticeable that when n(HPinserted) > 0, the cavity self-closure occurred. Example snapshots of HPβCDs with different numbers of n(HPinserted) are given in Figure 9. As shown in the figure, the flipped HP-substituted glucopyranose subunits and the HP groups inserted into the hydrophobic cavity led to cavity self-closure in the HPβCDs. When the number of HP groups was increased, the percentage of cavity self-closure significantly increased from 2.41 to 57.59% in Hep6-HPβCD. This indicates that the cavity of the HPβCD structure was almost blocked when the number of HP-substitutions was increased. Furthermore, the maximum of n(HPinserted) was equal to 2, even in the HPβCD with a high degree of HP-substitutions, such as Hep6-HPβCD and Tet2Tet6-HPβCD.
For HPβCDs with double-sided HP-substitutions, the HP groups, both on the narrow and wider rim of HPβCD could point into the cavity. For that reason, the percentage of cavity self-closure was higher than those of single-sided HP-substitutions, because of the large fluctuation of their structures. When increasing the number of HP groups at the O2 positions on Tet6-HPβCD from one to four, the probability of cavity self-closure rose dramatically from 38.34% (Tet6-HPβCD) to 55.72, 66.60, 73.54, and 61.42% for Mon2Tet6-, Di2Tet6-, Tri2Tet6-, and Tet2Tet6-HPβCDs, respectively. With increasing temperatures, the glucose subunits distorted, and the HP groups fluctuated far from the center of mass of the βCD ring, leading to a lower probability of cavity self-closure at higher temperatures for all models (Table S2).
In summary, all data showed that the HP-substitution strongly affects the cavity self-closure of HPβCD, especially when the number of the HP substituents increases. To prevent cavity self-closure, which might block the inclusion of the guest molecules, we suggest that the HP-substitution should be a single-sided substitution, with the degree of substitution being less than 0.57 such as Di-, Tri-, and Tet6-HPβCDs. However, it should be mentioned that the above results were only valid in cases where there was no interaction between the cyclodextrins (at higher concentrations). Moreover, the change of the dielectric properties, e.g., by using solvent mixtures, might influence the conformations of the considered molecules, as well as their conformational equilibria.

4. Conclusions

REMD simulations have been performed on HPβCDs with various degrees of substitution (DS = 0.14–1.14) to study the structural behaviors and the effects of HP substituents. Several parameters that influence such changes have been identified. The circularity and the radius of gyration explain the size and shape of the cavity of the ring, and the flip angle and important distances describe the conformational changes and the flexibilities of the HP groups. The results show that HPβCDs have a more pronounced conical shape than βCD; however, cavity self-closure occurs because some glucopyranoses with HP groups flip or distort, followed by one or two HP groups coming close to the CD cavity, thus hindering drug inclusion. HPβCDs with high DS are more likely to be blocked, while HPβCDs with double-sided HP-substitutions are even more likely to be blocked. Among the nine HPβCDs, all analysis parameters point out that Tri- and Tet6-HPβCDs with three and four HP-substitutions on the primary rim have a distinctive conformation, being mostly circular with a low possibility of flipping and cavity enclosure. Thus, these HPβCDs could serve as more proper hosts for the encapsulation of low-water soluble compounds.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4360/11/1/145/s1, Figure S1: The overlapping between the potential energy distributions of each replica temperature, ranging from 269.5 K to 570.9 K. Figure S2: The distance of the centers of mass between the βCD ring and the HP group, d4(i)[Cg(βCD)–Cg(HP(i))] for (a) Di6-HPβCD, (b) Tri6-HPβCD, (c) Tet6-HPβCD, (d) Hep6-HPβCD, (e) Mon2Tet6-HPβCD, (f) Di2Tet6-HPβCD, (g) Tri2Tet6-HPβCD, and (h) Tet2Tet6-HPβCD at 300 K. Figures S3–S12: Contour graphs of the native βCD (S3), Mon6-HPβCD (S4), Di6-HPβCD (S5), Tri6-HPβCD (S6), Tet6-HPβCD (S7), Hep6-HPβCD (S8), Mon2Tet6-HPβCD (S9), Di2Tet6-HPβCD (S10), Tri2Tet6-HPβCD (S11), and Tet2Tet6-HPβCD (S12), probability distributions of 25,000 snapshots with the glycam06 force field at various temperatures. Figure S13: The radius of gyration at various temperatures for (a) βCD, (b) Mon6-HPβCD, (c) Di6-HPβCD, (d) Tri6-HPβCD, (e) Tet6-HPβCD, (f) Hep6-HPBCD, (g) Mon2Tet6-HPBCD, (h) Di2Tet6-HPBCD, (i) Tri2Tet6-HPBCD, and (j) Tet2Tet6-HPBCD. Table S1: The probabilities of different numbers of flip glucose subunits in βCD and all HPβCDs, using the flip angle parameter, θ(i)[C6(i)–C2(i+1)–C6(i+1)] at various temperatures (criteria: a value higher than 90 degree) compared with the classical MD simulation in the parenthesis. Table S2: The probability of the overall snapshots with different numbers of HP occupied in the CD cavity (criteria: d4(i)[Cg(βCD)–Cg(HP(i))] < 3 Å) for all HPβCDs at various temperatures.

Author Contributions

K.K. analyzed all data and wrote the original draft. Some results of the classical MD simulation were analyzed by J.K., T.R. and N.K. conceived the conceptualization. P.W., S.M., C.R., M.K., H.O., N.K., and T.R. supervised the project. All authors discussed the results and contributed to the final manuscript.

Acknowledgments

T.R. and N.K. gratefully acknowledge the financial support from the Thailand Research Fund (RSA5980069 to T.R. and RSA6180044 to N.K.). This work was supported by a Grant-in-Aid for Scientific Research on Innovative Areas in “Precisely Designed Catalysts with Customized Scaffolding” (JSPS KAKENHI Grant Number JP18H04233 to S.M.). This research work was partially supported by Chiang Mai University. K.K. thanks W. Khuntawee, C. Hanpiboon, B. Nutho, and T. Saelee for resolving some technical problems. P.W. thanks Chulalongkorn University (CU) for a short-term visit grant. Travel grants for research visits, and for the research reported in this publication, were also supported by the ASEAN-European Academic University Network (ASEA-UNINET). The Center of Excellence Computational Chemistry (CU), and the Vienna Scientific Cluster (VSC-2) are gratefully thanked for providing computational facilities.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Uekama, K.; Hirayama, F.; Irie, T. Cyclodextrin Drug Carrier Systems. Chem. Rev. 1998, 98, 2045–2076. [Google Scholar] [CrossRef] [PubMed]
  2. Hirayama, F.; Uekama, K. Cyclodextrin-based controlled drug release system. Adv. Drug Deliv. Rev. 1999, 36, 125–141. [Google Scholar] [CrossRef]
  3. Aachmann, F.L.; Otzen, D.E.; Larsen, K.L.; Wimmer, R. Structural background of cyclodextrin-protein interactions. Protein Eng. 2003, 16, 905–912. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Davis, M.E.; Brewster, M.E. Cyclodextrin-based pharmaceutics: Past, present and future. Nat. Rev. Drug Discov. 2004, 3, 1023–1035. [Google Scholar] [CrossRef] [PubMed]
  5. Challa, R.; Ahuja, A.; Ali, J.; Khar, R.K. Cyclodextrins in drug delivery: An updated review. AAPS PharmSciTech 2005, 6, E329–E357. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Loftsson, T.; Jarho, P.; Másson, M.; Järvinen, T. Cyclodextrins in drug delivery. Expert Opin. Drug Deliv. 2005, 2, 335–351. [Google Scholar] [CrossRef] [PubMed]
  7. Brewster, M.E.; Loftsson, T. Cyclodextrins as pharmaceutical solubilizers. Adv. Drug Deliv. Rev. 2007, 59, 645–666. [Google Scholar] [CrossRef]
  8. Marques, H.M.C. A review on cyclodextrin encapsulation of essential oils and volatiles. Flavour Frag. J. 2010, 25, 313–326. [Google Scholar] [CrossRef]
  9. Crini, G. Review: A History of Cyclodextrins. Chem. Rev. 2014, 114, 10940–10975. [Google Scholar] [CrossRef] [PubMed]
  10. Gould, S.; Scott, R.C. 2-Hydroxypropyl-β-cyclodextrin (HP-β-CD): A toxicology review. Food Chem. Toxicol. 2005, 43, 1451–1459. [Google Scholar] [CrossRef] [PubMed]
  11. Schönbeck, C.; Westh, P.; Madsen, J.C.; Larsen, K.L.; Städe, L.W.; Holm, R. Hydroxypropyl-Substituted β-Cyclodextrins: Influence of Degree of Substitution on the Thermodynamics of Complexation with Tauroconjugated and Glycoconjugated Bile Salts. Langmuir 2010, 26, 17949–17957. [Google Scholar] [CrossRef] [PubMed]
  12. Reyes-Reyes, M.L.; Roa-Morales, G.; Melgar-Fernández, R.; Reyes-Pérez, H.; Gómez-Oliván, L.; Gonzalez-Rivas, N.; Bautista-Renedo, J.; Balderas-Hernández, P. Chiral recognition of abacavir enantiomers by (2-hydroxy)propyl-β-cyclodextrin: UHPLC, NMR and DFT studies. J. Incl. Phenom. Macrocycl. Chem. 2015, 82, 373–382. [Google Scholar] [CrossRef]
  13. Yuan, C.; Liu, B.; Liu, H. Characterization of hydroxypropyl-β-cyclodextrins with different substitution patterns via FTIR, GC–MS, and TG–DTA. Carbohydr. Polym. 2015, 118, 36–40. [Google Scholar] [CrossRef] [PubMed]
  14. Buvári-Barcza, Á.; Barcza, L. Influence of the guests, the type and degree of substitution on inclusion complex formation of substituted β-cyclodextrins. Talanta 1999, 49, 577–585. [Google Scholar] [CrossRef]
  15. Yuan, C.; Jin, Z.; Li, X. Evaluation of complex forming ability of hydroxypropyl-β-cyclodextrins. Food Chem. 2008, 106, 50–55. [Google Scholar] [CrossRef]
  16. Concha-Santos, S.; Pérez-Casas, S.; Brocos, P.; Piñeiro, Á. Testing the effect of the cavity size and the number of molecular substitutions in host–guest complexes formed by 2-hydroxypropyl-cyclodextrins and n-octyl-β-D-glucopyranoside. J. Chem. Thermodyn. 2013, 67, 112–119. [Google Scholar] [CrossRef]
  17. Yong, C.; Washington, C.; Smith, W. Structural Behaviour of 2-Hydroxypropyl-β-Cyclodextrin in Water: Molecular Dynamics Simulation Studies. Pharm. Res. 2008, 25, 1092–1099. [Google Scholar] [CrossRef] [PubMed]
  18. Pitha, J.; Trinadha Rao, C.; Lindberg, B.; Seffers, P. Distribution of substituents in 2-hydroxypropyl ethers of cyclomaltoheptaose. Carbohydr. Res. 1990, 200, 429–435. [Google Scholar] [CrossRef]
  19. Trinadha Rao, C.; Pitha, J.; Lindberg, B.; Lindberg, J. Distribution of substituents in O-(2-hydroxypropyl) derivatives of cyclomalto-oligosaccharides (cyclodextrins): Influence of increasing substitution, of the base used in the preparation, and of macrocyclic size. Carbohydr. Res. 1992, 223, 99–107. [Google Scholar] [CrossRef]
  20. Buvári-Barcza, Á.; Bodnár-Gyarmathy, D.; Barcza, L. Hydroxypropyl-β-cyclodextrins: Correlation between the stability of their inclusion complexes with phenolphthalein and the degree of substitution. J. Incl. Phenom. Mol. Recognit. Chem. 1994, 18, 301–306. [Google Scholar] [CrossRef]
  21. Malanga, M.; Szemán, J.; Fenyvesi, É.; Puskás, I.; Csabai, K.; Gyémánt, G.; Fenyvesi, F.; Szente, L. “Back to the Future”: A New Look at Hydroxypropyl Beta-Cyclodextrins. J. Pharm. Sci. 2016, 105, 2921–2931. [Google Scholar] [CrossRef] [PubMed]
  22. Nutho, B.; Nunthaboot, N.; Wolschann, P.; Kungwan, N.; Rungrotmongkol, T. Metadynamics supports molecular dynamics simulation-based binding affinities of eucalyptol and beta-cyclodextrin inclusion complexes. RSC Adv. 2017, 7, 50899–50911. [Google Scholar] [CrossRef]
  23. Sugita, Y.; Okamoto, Y. Replica-exchange molecular dynamics method for protein folding. Chem. Phys. Lett. 1999, 314, 141–151. [Google Scholar] [CrossRef]
  24. Hukushima, K.; Nemoto, K. Exchange Monte Carlo Method and Application to Spin Glass Simulations. J. Phys. Soc. Jpn. 1996, 65, 1604–1608. [Google Scholar] [CrossRef] [Green Version]
  25. Hansmann, U.H.E. Parallel tempering algorithm for conformational studies of biological molecules. Chem. Phys. Lett. 1997, 281, 140–150. [Google Scholar] [CrossRef] [Green Version]
  26. Nymeyer, H.; Gnanakaran, S.; García, A.E. Atomic Simulations of Protein Folding, Using the Replica Exchange Algorithm. In Methods in Enzymology; Academic Press: Cambridge, MA, USA, 2004; Volume 383, pp. 119–149. [Google Scholar]
  27. Cheng, X.; Cui, G.; Hornak, V.; Simmerling, C. Modified Replica Exchange Simulation Methods for Local Structure Refinement. J. Phys. Chem. B 2005, 109, 8220–8230. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Snor, W.; Liedl, E.; Weiss-Greiler, P.; Karpfen, A.; Viernstein, H.; Wolschann, P. On the structure of anhydrous β-cyclodextrin. Chem. Phys. Lett. 2007, 441, 159–162. [Google Scholar] [CrossRef]
  29. Kicuntod, J.; Khuntawee, W.; Wolschann, P.; Pongsawasdi, P.; Chavasiri, W.; Kungwan, N.; Rungrotmongkol, T. Inclusion complexation of pinostrobin with various cyclodextrin derivatives. J. Mol. Graph. Model. 2016, 63, 91–98. [Google Scholar] [CrossRef] [PubMed]
  30. Case, D.A.; Babin, V.; Berryman, J.T.; Betz, R.M.; Cai, Q.; Cerutti, D.S.; Cheatham, T.E.; Darden, T.A.; Duke, R.E.; Gohlke, H.; et al. AMBER 14; University of California, San Francisco: San Francisco, CA, USA, 2014. [Google Scholar]
  31. Kirschner, K.N.; Yongye, A.B.; Tschampel, S.M.; González-Outeiriño, J.; Daniels, C.R.; Foley, B.L.; Woods, R.J. GLYCAM06: A generalizable biomolecular force field. Carbohydrates. J. Comput. Chem. 2008, 29, 622–655. [Google Scholar] [CrossRef]
  32. Tessier, M.B.; DeMarco, M.L.; Yongye, A.B.; Woods, R.J. Extension of the GLYCAM06 biomolecular force field to lipids, lipid bilayers and glycolipids. Mol. Simul. 2008, 34, 349–364. [Google Scholar] [CrossRef]
  33. Khuntawee, W.; Rungrotmongkol, T.; Wolschann, P.; Pongsawasdi, P.; Kungwan, N.; Okumura, H.; Hannongbua, S. Conformation study of ε-cyclodextrin: Replica exchange molecular dynamics simulations. Carbohydr. Polym. 2016, 141, 99–105. [Google Scholar] [CrossRef]
  34. Khuntawee, W.; Kunaseth, M.; Rungnim, C.; Intagorn, S.; Wolschann, P.; Kungwan, N.; Rungrotmongkol, T.; Hannongbua, S. Comparison of Implicit and Explicit Solvation Models for Iota-Cyclodextrin Conformation Analysis from Replica Exchange Molecular Dynamics. J. Chem. Inf. Model. 2017, 57, 778–786. [Google Scholar] [CrossRef] [PubMed]
  35. Wongpituk, P.; Nutho, B.; Panman, W.; Kungwan, N.; Wolschann, P.; Rungrotmongkol, T.; Nunthaboot, N. Structural dynamics and binding free energy of neral-cyclodextrins inclusion complexes: Molecular dynamics simulation. Mol. Simul. 2017, 43, 1356–1363. [Google Scholar] [CrossRef]
  36. Mahalapbutr, P.; Nutho, B.; Wolschann, P.; Chavasiri, W.; Kungwan, N.; Rungrotmongkol, T. Molecular insights into inclusion complexes of mansonone E and H enantiomers with various β-cyclodextrins. J. Mol. Graph. Model. 2018, 79, 72–80. [Google Scholar] [CrossRef] [PubMed]
  37. Mahalapbutr, P.; Thitinanthavet, K.; Kedkham, T.; Nguyen, H.; thi ha Theu, L.; Dokmaisrijan, S.; Huynh, L.; Kungwan, N.; Rungrotmongkol, T. A theoretical study on the molecular encapsulation of luteolin and pinocembrin with various derivatized beta-cyclodextrins. J. Mol. Struct. 2018, in press. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of 2-hydroxypropyl-β-cyclodextrin (HPβCD), which comprises seven glucopyranose units.
Figure 1. Schematic representation of 2-hydroxypropyl-β-cyclodextrin (HPβCD), which comprises seven glucopyranose units.
Polymers 11 00145 g001
Figure 2. Models introducing 2-hydroxypropyl (HP) groups at the O2 (red) and/or O6 (blue) positions on glucose subunits.
Figure 2. Models introducing 2-hydroxypropyl (HP) groups at the O2 (red) and/or O6 (blue) positions on glucose subunits.
Polymers 11 00145 g002
Figure 3. (a) Cyclodextrin (CD) fragment showing the atomic labels and important structural parameters, d1(i)[O2(i)–O3(i+1)], d2(i)[O4(i)–O4(i+1)], and θ(i)[C6(i)–C2(i+1)–C6(i+1)]. (b) Set of diameters, d3(i)[Cg(Glu(i))–O4(i+3)], for circularity. (c) Set of d4(i)[Cg(βCD)–Cg(HP(i))] for studying cavity self-closure.
Figure 3. (a) Cyclodextrin (CD) fragment showing the atomic labels and important structural parameters, d1(i)[O2(i)–O3(i+1)], d2(i)[O4(i)–O4(i+1)], and θ(i)[C6(i)–C2(i+1)–C6(i+1)]. (b) Set of diameters, d3(i)[Cg(Glu(i))–O4(i+3)], for circularity. (c) Set of d4(i)[Cg(βCD)–Cg(HP(i))] for studying cavity self-closure.
Polymers 11 00145 g003
Figure 4. (a) Contour graphs of the probability distributions in terms of free energy, (b) the probability of d1(i), and (c) the probability of d2(i) from a total of 25,000 snapshots at 300 K of REMD simulations with the glycam06 force field for βCD and HPβCDs.
Figure 4. (a) Contour graphs of the probability distributions in terms of free energy, (b) the probability of d1(i), and (c) the probability of d2(i) from a total of 25,000 snapshots at 300 K of REMD simulations with the glycam06 force field for βCD and HPβCDs.
Polymers 11 00145 g004
Figure 5. The flipped conformations of βCD with different numbers of flip glucose subunits, which were selected from the simulation at 300 K.
Figure 5. The flipped conformations of βCD with different numbers of flip glucose subunits, which were selected from the simulation at 300 K.
Polymers 11 00145 g005
Figure 6. The probability of the radius of gyration for βCD and all HPβCDs at 300 K.
Figure 6. The probability of the radius of gyration for βCD and all HPβCDs at 300 K.
Polymers 11 00145 g006
Figure 7. The distance of the centers of mass between the βCD ring and HP group, d4(i)[Cg(βCD)–Cg(HP(i))] for Mon6-HPβCD at 300 K.
Figure 7. The distance of the centers of mass between the βCD ring and HP group, d4(i)[Cg(βCD)–Cg(HP(i))] for Mon6-HPβCD at 300 K.
Polymers 11 00145 g007
Figure 8. The probability of n(HPinserted) from 25,000 snapshots (criteria: d4(i)[Cg(βCD)–Cg(HP(i))] < 3 Å) at 300 K for all HPβCDs.
Figure 8. The probability of n(HPinserted) from 25,000 snapshots (criteria: d4(i)[Cg(βCD)–Cg(HP(i))] < 3 Å) at 300 K for all HPβCDs.
Polymers 11 00145 g008
Figure 9. Example snapshots of different numbers of HP occupied in the CD cavity (criteria: d4(i)[Cg(βCD)–Cg(HP(i))] < 3 Å) for all HPβCDs at 300 K.
Figure 9. Example snapshots of different numbers of HP occupied in the CD cavity (criteria: d4(i)[Cg(βCD)–Cg(HP(i))] < 3 Å) for all HPβCDs at 300 K.
Polymers 11 00145 g009
Table 1. Model summary of introducing HP groups at the O2 and/or O6 positions on glucose subunits.
Table 1. Model summary of introducing HP groups at the O2 and/or O6 positions on glucose subunits.
ModelsDegree of SubstitutionO2 SubstitutionO6 Substitution
βCD0.00NoneNone
Single-sided HP-substitution
Mon6-HPβCD0.14None1 (At glucose unit 1)
Di6-HPβCD0.28None2 (At glucose units 1 and 3)
Tri6-HPβCD0.43None3 (At glucose units 1, 3, and 5)
Tet6-HPβCD0.57None4 (At glucose units 1, 3, 5, and 7)
Hep6-HPβCD1.00None7 (At all glucose units)
Double-sided HP-substitution
Mon2Tet6-HPβCD0.711 (At glucose unit 4)4 (At glucose units 1, 3, 5, and 7)
Di2Tet6-HPβCD0.852 (At glucose units 2 and 6)4 (At glucose units 1, 3, 5, and 7)
Tri2Tet6-HPβCD1.003 (At glucose units 2, 4, and 6)4 (At glucose units 1, 3, 5, and 7)
Tet2Tet6-HPβCD1.144 (At glucose units 2, 4, 6, and 7)4 (At glucose units 1, 3, 5, and 7)
Table 2. The probability of different numbers of flip glucose subunits in βCD and HPβCDs, using the flip angle parameter, θ(i)[C6(i)–C2(i+1)–C6(i+1)] at 300K (criteria: values higher than 90 degrees).
Table 2. The probability of different numbers of flip glucose subunits in βCD and HPβCDs, using the flip angle parameter, θ(i)[C6(i)–C2(i+1)–C6(i+1)] at 300K (criteria: values higher than 90 degrees).
ModelsThe Percentage of the Flip Angle (%)
No FlipOne FlipTwo Flips
βCD58 35 7
Single-sided HP-substitution
Mon6-HPβCD69283
Di6-HPβCD74233
Tri6-HPβCD75241
Tet6-HPβCD78211
Hep6-HPβCD73252
Double-sided HP-substitution
Mon2Tet6-HPβCD77221
Di2Tet6-HPβCD75241
Tri2Tet6-HPβCD70282
Tet2Tet6-HPβCD74242
Table 3. The average and standard deviation of the circularity ( C ) of βCD and HPβCDs at 300 K, using REMD simulations.
Table 3. The average and standard deviation of the circularity ( C ) of βCD and HPβCDs at 300 K, using REMD simulations.
Models C
βCD0.727 ± 0.087
Single-sided HP-substitution
Mon6-HPβCD0.746 ± 0.085
Di6-HPβCD0.751 ± 0.086
Tri6-HPβCD0.814 ± 0.075
Tet6-HPβCD0.815 ± 0.076
Hep6-HPβCD0.773 ± 0.088
Double-sided HP-substitution
Mon2Tet6-HPβCD0.803 ± 0.075
Di2Tet6-HPβCD0.767 ± 0.076
Tri2Tet6-HPβCD0.785 ± 0.074
Tet2Tet6-HPβCD0.742 ± 0.080

Share and Cite

MDPI and ACS Style

Kerdpol, K.; Kicuntod, J.; Wolschann, P.; Mori, S.; Rungnim, C.; Kunaseth, M.; Okumura, H.; Kungwan, N.; Rungrotmongkol, T. Cavity Closure of 2-Hydroxypropyl-β-Cyclodextrin: Replica Exchange Molecular Dynamics Simulations. Polymers 2019, 11, 145. https://doi.org/10.3390/polym11010145

AMA Style

Kerdpol K, Kicuntod J, Wolschann P, Mori S, Rungnim C, Kunaseth M, Okumura H, Kungwan N, Rungrotmongkol T. Cavity Closure of 2-Hydroxypropyl-β-Cyclodextrin: Replica Exchange Molecular Dynamics Simulations. Polymers. 2019; 11(1):145. https://doi.org/10.3390/polym11010145

Chicago/Turabian Style

Kerdpol, Khanittha, Jintawee Kicuntod, Peter Wolschann, Seiji Mori, Chompoonut Rungnim, Manaschai Kunaseth, Hisashi Okumura, Nawee Kungwan, and Thanyada Rungrotmongkol. 2019. "Cavity Closure of 2-Hydroxypropyl-β-Cyclodextrin: Replica Exchange Molecular Dynamics Simulations" Polymers 11, no. 1: 145. https://doi.org/10.3390/polym11010145

APA Style

Kerdpol, K., Kicuntod, J., Wolschann, P., Mori, S., Rungnim, C., Kunaseth, M., Okumura, H., Kungwan, N., & Rungrotmongkol, T. (2019). Cavity Closure of 2-Hydroxypropyl-β-Cyclodextrin: Replica Exchange Molecular Dynamics Simulations. Polymers, 11(1), 145. https://doi.org/10.3390/polym11010145

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop