Theoretical modeling of open-shell molecules in solution: a QM/MM molecular dynamics approach