def binding_site_in_pymol(binding_site):
    """
    Get PyMol commands to visualize binding site.
    """
    
    pymol_command = (
        f'fetch 6LU7\n'
        f'remove solvent\n'
        f'select pocket, 6LU7 and resi {"+".join([str(i) for i in binding_site.sort_values("residue_id").residue_id.to_list()])}\n'
        f'show cartoon\n'
        f'hide lines\n'
        f'color blue, pocket\n'
        f'show lines, pocket'
    )
    
    return pymol_command