Through atomistic molecular dynamic simulations using a GROMOS53a6 force field for the carbohydrate, we studied the lyotropic reverse hexagonal phase HII from a glycolipid, namely the Guerbet branched-chain β-d-glucoside, at 14% and 22% water concentrations. Our simulations showed that at low water concentration (14%) the sugar head group overlapped extensively and protruded into the water channel. In contrast, in the 22% concentration system a water column free from the sugar headgroup ('free' water) was formed as expected for the system close to the limit of maximum hydration. In both concentrations, we found anomalous water diffusion in the xy-plane, i.e. the two-dimensional space confined by the surface of the cylinder. On the other hand, along the z-axis, the water diffusion obeyed the Einstein relation for the 22% system, while for the 14% system it was slightly anomalous. For the 22% system, the diffusion along the z-axis of the 'free' water obeyed the Einstein relation, while that of the 'bound' water is slightly anomalous. The xy-plane displacement of the 'bound' water was higher than that for the 'free' water at times longer than 200 ps, as a consequence of the exchange of water molecules between the two regions. Based on our findings, we proposed an alternative explanation to the observed spatial heterogeneity in the HII phase from probe diffusion by Penaloza et al. (Phys. Chem. Chem. Phys., 2012, 14(15), 5247-5250). We found the extent of contact with water was different at different oxygen atoms within the sugar ring. Generally, a higher probability of hydrogen bonding but a shorter lifetime was found in 22% water compared to the case of 14% water. Finally, we examined the extension and compression of the alkyl chain of a columnar.
A theoretical study of a series of five glucose based glycolipid crown ethers and their complexes with Na(+) and K(+) was performed using the density functional theory with B3LYP/6-31 G* to obtain the optimized geometrical structures and electronic properties. The local nucleophilicity of the five molecules was investigated using Fukui function, while the global nucleophilicity was calculated from the ionization potential and electron affinity. The structures and coordination of the complexes were studied to identify the best match of the glycolipid crown ethers with cations. In general, it was found that the oxygen atoms pairs O2 and O3 (or O4 and O6) on the sugar ring are constrained from moving toward the cation, which results in a weaker O-cation coordination strength for the oxygen pair compared to the other oxygen atoms in the crown ether ring. The thermodynamic properties of the binding of the complexes and the exchange reaction in gas phase were evaluated. The cation selectivity pattern among the five molecules was in good agreement with the experiment.
Fully atomistic molecular dynamics simulation studies of thermotropic bilayers were performed using a set of glycosides namely n-octyl-β-D-glucopyranoside (β-C8Glc), n-octyl-α-D-glucopyranoside (α-C8Glc), n-octyl-β-D-galactopyranoside (β-C8Gal), and n-octyl-α-D-galactopyranoside (α-C8Gal) to investigate the stereochemical relationship of the epimeric/anomeric quartet liner glycolipids with the same octyl chain group. The results showed that, the anomeric stereochemistry or the axial/equatorial orientation of C1-O1 (α/β) is an important factor controlling the area and d-spacing of glycolipid bilayer systems in the thermotropic phase. The head group tilt angle and the chain ordering properties are affected by the anomeric effect. In addition, the L(C) phase of β-C8Gal, is tilting less compared to those in the fluid L(α). The stereochemistry of the C4-epimeric (axial/equatorial) and anomeric (α/β) centers simultaneously influence the inter-molecular hydrogen bond. Thus, the trend in the values of the hydrogen bond for these glycosides is β-C8Gal > α-C8Glc > β-C8Glc > α-C8Gal. The four bilayer systems showed anomalous diffusion behavior with an observed trend for the diffusion coefficients; and this trend is β-C8Gal > β-C8Glc > α-C8Gal > α-C8Glc. The "bent" configuration of the α-anomer results in an increase of the hydrophobic area, chain vibration and chain disorganization. Since thermal energy is dispensed more entropically for the chain region, the overall molecular diffusion decreases.
The molecular dynamics of a synthetic branched chain glycolipid, 2-decyl-tetradecyl-β-d-maltoside (C14-10G2), in the dry assemblage of smectic and columnar liquid crystal phases has been studied by dielectric spectroscopy as a function of frequency and temperature during the cooling process. Strong relaxation modes were observed corresponding to the tilted smectic and columnar phases, respectively. At low frequency (∼900 Hz to 1 kHz) in the smectic phase, Process I* was observed due to the tilted sugar bilayer structure. The process continued in the columnar phase (Process I) with an abrupt dynamic change due to phase transition in the frequency range of ∼1.3 kHz to 22 kHz. An additional process (Process II) was observed in the columnar phase with a broader relaxation in the frequency range of ∼10 Hz to 1 kHz. A bias field dependence study was performed in the columnar phase and we found that the relaxation strength rapidly decreased with increased applied dc bias field. This relaxation originates from a collective motion of polar groups within the columns. The results of dielectric spectroscopy were supported by a molecular dynamics simulation study to identify the origin of the relaxation processes, which could be related to the chirality and hydrogen bonds of the sugar lipid.