A new time-dependent density-functional method for molecular plasmonics: Formalism, implementation, and the Au144(SH)60 case study