Script para separar estados conformacionales de un PDB.
El siguiente script .pml para ejecutar desde PyMol separa los estados de un objeto cargado en una escena, por ejemplo, las 3 subunidades del trímero de 7SSQ.
pythonfrom pymol import cmdobj = "7SSQ" # Cambiar si tu objeto tiene otro nombren_states = cmd.count_states(obj)for i in range(1, n_states + 1): new_obj = f"{obj}_state{i}" cmd.create(new_obj, obj, i, 1) if cmd.count_atoms(new_obj) == 0: print(f"[AVISO] Estado {i} está vacío. Se elimina {new_obj}.") cmd.delete(new_obj) else: print(f"[OK] Se creó {new_obj} desde el estado {i}.")python end
¿Cómo usarlo?
Guarda ese contenido como separar_estados.pml y desde la consola de PyMol, ejecuta: run ruta/del/archivo/separar_estados.pml

Última actualización: 20/05/2025, 07:30:40