The nature of flavonoids in polar organic solvents solutions is studied using classical molecular dynamics simulations considering quercetin as an archetypical flavonoid and acetone, dimethylformamide and dimethyl sulfoxide as representatives of solvents with different polarity. The solvation, intermolecular forces (hydrogen bonding) and interactions of the flavonoid with the solvents are analyzed. Likewise, the role of quercetin on changing the solvent properties and the possibility of acting as a solubility enhancer for fullerenes (C60) are studied by considering the properties of C60 fullerene in quercetin plus polar solvents solutions. The reported results provide information on the nature of the considered complex liquid mixtures and analyze the possibility of using flavonoids as natural, non-toxic, modifiers of traditional polar organic solvents and to improve the solubility of complex solutes such as fullerene nanoparticles. [link]