AMBER
(
acronimo
di
Assisted Model Building with Energy Refinement
) e una famiglia di
campi di forze
utilizzati nella
dinamica molecolare
di
biomolecole
, originariamente sviluppati dal gruppo di Peter Kollman all'
Universita della California
,
San Francisco
.
AMBER
e allo stesso tempo anche il nome del
pacchetto di software
con cui si puo simulare il campo di forze.
Con l'espressione "
campo di forze
AMBER" normalmente ci si riferisce alla forma funzionale usata dai campi di forze AMBER. Questa e l'insieme di diversi parametri a cui ogni versione di AMBER assegna un determinato valore.
L'obiettivo del campo di forze e di modellizzare in modo accurato le energie conformazionali di un sistema molecolare e le conseguenti interazioni intermolecolari.
La forma funzionale rappresenta quindi l'energia potenziale del sistema:
[1]
Il campo di forze e la derivata di questa grandezza rispetto alla posizione.
Il significato dei diversi
termini
e il seguente:
- Primo termine (
sommatoria
sui legami): rappresenta l'energia di tutti i legami covalenti della molecola. Questa forza di tipo armonico (una molla ideale) e una buona approssimazione se la distanza tra i due atomi coinvolti nel legame non e troppo lontana dalla distanza di equilibrio, peggiora all'allontanarsi degli atomi.
- Secondo termine (sommatoria sugli angoli di legame): riguarda il contributo energetico dato dalla geometria assunta dagli orbitali elettronici implicati in legami covalenti.
- Terzo termine (sommatoria sugli angoli diedri o di torsione): stima l'energia correlata alla torsione di un legame in base all'ordine di legame e ai legami limitrofi. Un singolo legame puo avere piu termini di questo tipo, per questo motivo l'energia di torsione totale e espressa come una
serie di Fourier
.
- Quarto termine (doppia sommatoria su
e
): rappresenta l'energia di non-legame tra tutti gli atomi, e la somma del contributo energetico delle
forze di Van der Waals
(primo termine della sommatoria) e del contributo energetico delle
forze elettrostatiche
(secondo termine della sommatoria).
Per usare il campo di forze AMBER, e necessario assegnare dei valori ai parametri del campo di forze (e.g. costanti di forza, lunghezza ed angoli di equilibrio di ogni legame, cariche). Esistono numerose collezioni di questi valori, descritte nel dettaglio nel manuale utente del software AMBER. Ognuna di queste collezioni ha un nome e fornisce i parametri per un certo tipo di molecole.
- I dati riguardanti i peptidi, le proteine e gli acidi nucleici sono forniti dai dataset il cui nome e composto da "ff" e l'anno di pubblicazione, come ad esempio "ff99".
- GAFF (campo di forze AMBER generalizzato) fornisce i parametri per piccole molecole organiche per la simulazione di ligandi, farmaci ad esempio, e biomolecole.
- Il campo di forze GLYCAM e stato sviluppato da Rob Woods per la simulazione di carboidrati.
Il pacchetto AMBER fornisce un insieme di programmi con cui applicare il campo di forze alla simulazione di biomolecole. E scritto in
Fortran
90 e
C++
. Lo sviluppo e condotto da un'ampia associazione di laboratori per la maggior parte accademici. Nuove versioni sono distribuite in genere nella primavera di anni pari; AMBER 10 e stato distribuito nell'aprile 2008, mentre AMBER 18 e uscito nell'aprile 2018. Il pacchetto sorgente viene distribuito in formato open secondo le specifiche
GNU GPL
ed e composto dalle due componenti AMBER e AMBERTOOLS, mentre la versione commerciale e a pagamento. Le ultime versioni del pacchetto AMBER supportano l'accelerazione via GPU grazie al supporto del linguaggio
CUDA
.
- LEaP
e usato per preparare i file di input per la simulazione
- Antechamber
si occupa di assegnare valori ai parametri di piccole molecole usando GAFF
- SANDER
(Simulated Annealing with NMR-Derived Energy Restraints) e il programma centrale di simulazione e fornisce algoritmi di minimizzazione dell'energia e di dinamica molecolare
- pmemd
e una re-implementazione limitata di SANDER fatta da Bob Duke. E molto piu performante se utilizzato in parallelo da piu di 8-16 processori
- ptraj
svolge analisi numeriche dei risultati della simulazione. AMBER non permette la visualizzazione dei composti che viene invece effettuata con VMD o Sirius.
- MM-PBSA
consente calcoli con solvente implicito a partire da una struttura di una dinamica molecolare.
- ^
Cornell WD, Cieplak P, Bayly CI, Gould IR, Merz KM Jr, Ferguson DM, Spellmeyer DC, Fox T, Caldwell JW, Kollman PA,
A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules
, in
J. Am. Chem. Soc.
, vol. 117, 1995, p. 5179?5197.