GM1 gangliosides form a microdomain with sphingomyeline (SM) and cholesterol (Chol) and are deeply involved in the aggregation of amyloid beta (A beta) peptides on neural membranes. We performed molecular dynamics simulations on two kinds of lipid bilayers containing GM1 ganglioside: GM1/SM/Chol and GM1/POPC. Both 10 and 100 as simulations and another set of 10 ns simulations with different initial lipid arrangement essentially showed the same computational results. GM1 molecules in the GM1/SM/Chol membrane were condensed, whereas those in GM1/POPC membrane scattered. That is, the formation of GM1 cluster was observed only on the GM1/SM/Chol mixed membrane. There appeared numerous hydrogen bonds among glycan portions of the GM1 clusters due to the condensation. A comparison in distribution of lipid molecules between the two kinds of membranes suggested that cholesterol had important roles to prevent the membrane from interdigitation and to stabilize other lipids for interacting with each other. This property of cholesterol promotes the formation of GM1 clusters.