def METHOD_NAME(self, molecule):
    moltype = molecule.meta.get('moltype')
    if not moltype:
        raise ValueError('The molecule does not have a moltype name.')
    add_virtual_sites(molecule, prefix=moltype)
    includes = molecule.meta.get('post_section_lines', {})
    for section in self.sections:
        section_includes = includes.get(section, [])
        section_includes.append('#include "{moltype}_{section}_VirtGoSites.itp"'
                                .format(moltype=moltype, section=section))
        includes[section] = section_includes
    molecule.meta['post_section_lines'] = includes
    return molecule