Si può eseguire la minimizzazione
- con due metodi: Steepest Descent e Conjugate Gradient
- a partire da varie configurazioni NMR
- Portarsi sulla directory di lavoro (es.
~/orac/data
).
- Creare un input modificato (es minim.in) che esegua una minimizzazione
dell'energia della molecola:
&SETUP
[...]
READ_PDB ../pdb/chignolin.pdb
&END
&SIMULATION
MINIMIZE
CG 0.00001
# oppure SD 0.00001
WRITE_GRADIENT
END
&END
&POTENTIAL
[...]
CUTOFF 100.
STRETCHING
&END
&INOUT
ASCII 100.0 OPEN minim.pdb
&END
&RUN
TIME 3000.0
PRINT 100.0
PROPERTY 100.0
&END
- lanciare il programma usando la prima struttura NMR (es. chignolin.pdb)
- Visualizzare i risultati con vmd
- Misurare gli RMSD con vmd
- Analizzare i diversi contributi all'energia:
post-out P="TotPot NonBond Bonded" chigno.out
> ene
plot -1,2-4 ene
- Osservare
- se vi è correlazione tra variazione nei RMSD e nell'energia
- quale modifica conformazionale determina la massima variazione di
energia
- quale parte del potenziale guida la transizione verso il minimo
- Estrarre una struttura NMR dal file del PDB
- Modificare i nomi di atomo come per la struttura precedente
(vedi lo scorso lab.)
-