Archive for abril, 2010

Codi per posar un camp d’un sdf al títol de la molècula

Sovint, mentre manipulem fitxers sdf ens trobem que volem canviar el títol de la molècula. Moltes vegades, com a títol hi volem posar algun camp de les dades addicionals de la molècula en qüestió. M’he fet un petit script per fer aquesta tasca de manera automàtica i senzilla.

El programet llegeix l’sdf i, molècula a molècula, en llegeix el camp que volguem, l’assigna com a títol i desa la molècula en un nou fitxer. A més a més, agafa el títol original i el desa en un camp “old_title”.

?View Code PYTHON
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
"""
Script to move a field to an id in a sdf
Alfons Nonell-Canals - April 2010
"""
 
import optparse
import sys
from pybel import *
 
p = optparse.OptionParser()
p.add_option('--sdf', '-i', help='SDF file to modify')
p.add_option('--field', '-f', help='Field to you want to set as molecule title')
options, arguments = p.parse_args()
 
if not options.sdf  or not options.field:
        p.error('An input sdf file and a field are required')
 
molfile = options.sdf
field = options.field
 
outFile = molfile.replace('.sdf','')+'.Mod.sdf'
out = Outputfile('sdf',outFile)
 
for mol in readfile("sdf", molfile):
        oldTitle = mol.title
        mol.data['old_title'] = oldTitle
        newtitle = mol.data[field]
        mol.OBMol.SetTitle(newtitle)
        out.write(mol)

És hora de tirar-se a la piscina?

Fa temps que tinc una idea al cap, una idea que ja té un domini .com, que barreja xarxes socials amb química i biologia (i algunes coses més que de moment no diré). Hi ha qui m’ha animat a engegar-la i transformar-la en negoci.

Avui he llegit un article a Nature Jobs que diu que hi ha vida més enllà del món acadèmic. Paral·lelament a Expansión he llegit un altre article que parla de capital risc i empreses 2.0 (sembla que ara és el moment).

Fa unes setmanes, en un curs que vaig fer al PRBB ens animaven a engegar idees. Ens van explicar que comporta un risc però també ens van fer veure que el món científico-publico-acadèmic també té el seu risc (quasi cap dels presents al curs tenim contracte més allà del 31 de desembre vinent).

És hora de llençar-se a la piscina? per si de cas, jo la vaig omplint d’aigua i, a estones lliures, treballant en una primera versió de la “idea”.