Human telomerase that activates within cancer cells has a telomeric sequence at the 3' end. Each factor that stabilizes the G-quadruplex in guanine-rich telomeric sequences can inhibit the regular telomerase activity. Therefore, the telomeric G-quadruplex is known as a promising target in cancer treatment. In this work, we studied the binding of positively charged distamycin A and its uncharged derivative to the G-quadruplex in a solution environment by Molecular Dynamics (MD) simulation. The binding mechanism and subtle conformational changes were investigated as a result of the ligand attachment. Moreover, binding free energy and clustering analysis describe the stability and flexibility of G-quadruplexes upon ligand binding. Structural analyses displayed that the favorable binding of both ligands imposes significant stability and rigidity in G-quadruplex conformation compared to free G-quadruplex, especially charged distamycin. Hydration pattern and ion distribution were different for free G-quadruplex and both of the ligand complexes. Energy decomposition reveals the electrostatic effect on the stability of G-quadruplex. The radial distribution function displayed the solvent shell and ion moving away from the groove. The hydrogen bond played an essential role in the binding of both ligands, especially for the charged derivative. van der Waals interaction is the only factor that is more important in binding uncharged distamycin into G-quadruplex than the charged one. The calculated ΔG showed the stability of both ligands within grooves and good agreement with the experimental binding free energy data. Finally, the results suggest that ligand modification improves the binding mode toward stabilizing G-quadruplexes.