We report molecular dynamics simulations of cyclohenicosakis-[(1-->2)-beta-D-gluco-henicosapyranosyl], termed 'cyclosophohenicosamer', a member of a class of cyclic (1-->2)-beta-D-glucans ('cyclosophoraoses'). Our goals were to provide insights into the conformational preferences of these cyclosophoraoses. Simulated annealing and constant-temperature molecular dynamics calculations were performed on the DP 21 cyclosophohenicosamer. The radius of gyration (R(G)) of the molecule and the conformation of glycosidic dihedral angles were used to analyze the result of computational studies. Most glycosidic linkages were concentrated in the lowest-energy region of the phi-psi energy map, and the values of radius of gyration from our simulations were consistent with the reported experimental value. The simulations produced various types of compact and asymmetric conformations within reasonable ranges of the glycosidic linkage conformation and radius of gyration. The results indicate the presence of a high degree of molecular flexibility of cyclosophohenicosamer and suggest the uniqueness of inclusion complexation with other molecules through this molecular flexibility.