G-quadruplexes (GQs) are noncanonical nucleic acid structures that form in guanine-rich sequences, importantly gene promoters and telomeres. The guanines within GQs arrange in square, planar tetrads through Hoogsteen hydrogen bonding. These tetrads stack on top of one another to make a highly ordered structure with an electronegative core. Coordinated cations are required to counteract the core's electronegativity and stabilize the GQ. GQs play important roles in gene expression and genomic stability and are functionally relevant in several human proto-oncogenes. One such proto-oncogene, c-myc, is linked to several human cancers such as Burkitt's lymphoma and multiples types of myeloma due to GQ formation within the promoter region of the gene acting as a gene silencer. As such, the c-myc GQ is a novel drug target for cancer therapy. Two solved GQ structures of interest in the c-myc promoter are the 1:2:1 and 1:6:1 GQs, named for differences in central loop length. Loop regions of GQs are known to affect topology, stability, and overall structure. Moreover, studies have shown nucleolin preferentially binds to the 1:6:1 GQ, prompting our present study on the dynamics of c-myc GQs. To analyze the importance of loop sequence on structural dynamics and gain information on this potential drug target, we performed molecular dynamics simulations utilizing the CHARMM36 non-polarizable and the Drude-2017 polarizable force fields (FF). The Drude FF accounts for electronic degrees of freedom via negatively charged particles placed on every heavy atom in the systems and has provided new insights on GQ dynamics and DNA-ion interactions. Our simulations of the 1:2:1 and 1:6:1 GQs demonstrate that polarization is important to model the stability, ion-ion, and ion-DNA interactions that underly GQ dynamics, emphasizing the importance of polarizable molecular dynamics in studying GQs.