An increasingly comprehension of the folding intermediate states of DNA G-quadruplexes (G4s) iscurrently an important scientific challenge, especially for the human telomeric (h-tel) G4s-forming sequences, characterized by a highly polymorphic nature. Despite the G-triplex conformation was proposed as one of the possible folding intermediates for the antiparallel and hybrid h-tel G4s, for the parallel h-tel topology with an all-anti guanine orientation, a vertical strand-slippage involving the G-triplets was proposed in previous works through microseconds-long standard molecular dynamics simulations (MDs). Here, in order to get further insights into the vertical strand-slippage and the folding intermediate states of the parallel h-tel G4s, we have carried out a Well-Tempered Metadynamics simulation (WT-MetaD), which allowed us to retrieve an ensemble of six G4s having two/G-tetrad conformations derived by the G-triplets vertical slippage. The insights highlighted in this work are aimed at rationalizing the mechanistic characterisation of the parallel h-tel G4 folding process.G-quadruplexes 1,2 (G4s) are non canonical DNA structures present at the telomere ends and involved in the genome stability and in the inhibition of the telomerase enzyme, a reverse transcriptase responsible for the cellular immortalization 3,4 and up-regulated in ~85% of human tumors 3 . The evidence that G4s formation at the 3′-end of telomere DNA can effectively hamper telomerase from adding further repeats 3 suggests G4s as a promising target for anticancer therapy 5-7 . In order to identify new compounds as potential anticancer drugs, the structure-based drug design targeting G4s could be a crucial step in gaining accurate information about the structural features of G4s 8-13 . Unfortunately, the rational design of novel G4 binders is hampered both by the structural polymorphism of these non-canonical conformations and by the lack of a well-defined ligand binding site. Therefore, enhanced sampling methods, like Metadynamics, have been required in order to deeply define the ligand binding mode 14,15 . Monomolecular G4s of the human telomeric sequence d[AG 3 (TTAGGG) 3 ] (h-tel) can adopt different folding topologies according to the relative orientations of the four strands 16 . In fact, while the arrangement of the h-tel monomolecular G4s has been well characterized by NMR in Na + to adopt a basket-type conformation with the single strands assuming an antiparallel orientation 17 , in presence of K + ions the folding topology of the h-tel monomolecular G4s includes both parallel-stranded and hybrid (mixed parallel/antiparallel) conformations, raising the question of which conformation is biologically more favoured. In 2006, the NMR hybrid syn/anti conformation was proposed as the most favoured one under K + physiological conditions 18 , whose folding organization is different from that reported in K + ions by crystallography consisting, instead, in the parallel-stranded topology with three G-tetrads and three symmetrical external d(TTA) l...