Fins fa una estona, tenia un problema: estic fent unes dinàmiques moleculars d’un enzim que catal·litza una reacció determinada, en aquest cas, l’oxidació de l’anió formiat a diòxid de carboni. La reacció transfereix un protó a l’enzim que, després s’ha de regenerar per recomençar el cicle. Fer la dinàmica d’aquesta reacció és relativament senzill, he afegit el substrat (formiat) al PDB de la proteïna i he fet la dinàmica. Perfecte, obtinc el CO₂ i l’enzim protonat. Fer les dinàmiques de la regeneració de l’enzim ja no és tant fàcil: vull eliminar producte (el diòxid de carboni) però part del residu introduït al PDB (l’anió formiat) no vull que desapareixi, ja que és el protó que he d’estudiar com deixa l’enzim.

Hi ha opcions complicades: eliminar part del formiat (el CO₂) del PDB, però això implica canviar la llibreria d’aminoàcids. També es pot modificar el PDB per situar el CO₂ llunt del centre actiu, però això fa que el programa no entengui que un residu tingui el seu hidrogen, per exemple, a 100 amstrongs. En fi, ho he provat i no és gens fàcil.

La solució? “silenciar” part del residu formiat. Les dinàmiques que estic fent són una mena de QM/MM i, de fet, en la regeneració de l’enzim ja tinc el diòxid de carboni fora de la part QM mentre el protó sí que hi és. En estar fora de la part QM, es podria dir que la importància de cada àtom rau en la seva càrrega (que s’obté de la llibreria de residus). Doncs, si poso a zero la càrrega del CO₂, problema resolt! es pot fer modificant la llibreria però, per fer-ho de manera senzilla jo ho faig als fitxers d’entrada del càlcul. Com?

1) dic al programa que el residu “formiat” té càrrega zero
2) assigno una càrrega concreta al protó.

En MOLARIS, el programa que estic emprant, es faria així. Abans de començar cada dinàmica cal posar el següent:

enzymix
pre_enz
setcrg0 720
setcrg 11192 0.013
quit
quit