Aggregation of polyglutamine peptides with β‐sheet structures is related to some important neurodegenerative diseases such as Huntington's disease. However, it is not clear how polyglutamine peptides form the β‐sheets and aggregate. To understand this problem, we performed all‐atom replica‐exchange molecular dynamics simulations of one and two polyglutamine peptides with 10 glutamine residues in explicit water molecules. Our results show that two polyglutamine peptides mainly formed helix or coil structures when they are separated, as in the system with one‐polyglutamine peptide. As the interpeptide distance decreases, the intrapeptide β‐sheet structure sometimes appear as an intermediate state, and finally the interpeptide β‐sheets are formed. We also find that the polyglutamine dimer tends to form the antiparallel β‐sheet conformations rather than the parallel β‐sheet, which is consistent with previous experiments and a coarse‐grained molecular dynamics simulation. © 2014 Wiley Periodicals, Inc.