Docking (chimica)

Diagramma che illustra il docking di una piccola molecola ligando (marrone) ad una proteina recettore (verde) per produrre un complesso.
Una piccola molecola legata a una proteina dopo docking.

Il docking molecolare, nel campo della modellistica molecolare, è un metodo che predice l'orientamento preferito di una molecola verso una seconda quando queste si legano fra di loro per formare un complesso stabile.[1]

La conoscenza dell'orientamento preferito può essere utilizzata per predire la forza di un'associazione o di un legame proteina-ligando tra due molecole utilizzando per esempio le funzioni di score.

Le associazioni tra molecole biologicamente rilevanti come proteine, acidi nucleici, carboidrati e lipidi giocano un ruolo centrale nella trasduzione del segnale. Il relativo orientamento delle due molecole che interagiscono può influenzare il tipo di segnale prodotto (ad esempio agonista verso antagonista). Il docking è quindi utile per predire sia la forza che il tipo di segnale prodotto.

Il docking è usato frequentemente per predire l'orientamento del legame di una piccola molecola farmacologicamente attiva alla sua proteina bersaglio, così da poter prevedere l'affinità e l'attività di questa molecola. Il docking quindi gioca un ruolo importante nella progettazione razionale del farmaco.[2][3] Dato il significato biologico e farmaceutico del docking molecolare, sono stati fatti considerevoli sforzi per migliorare i metodi utilizzati per predire il docking.

Definizione del problema

[modifica | modifica wikitesto]

Il docking molecolare si può pensare come un problema di “serratura e chiave”, dove bisogna capire il corretto orientamento della “chiave” che aprirà la “serratura” (in quale superficie della serratura c'è il buco, in quale direzione bisogna girare la chiave dopo che è inserita, ecc. ). La proteina può essere immaginata come la "serratura" e il ligando come "chiave". Il docking molecolare può essere definito come un problema di ottimizzazione nel quale bisogna trovare l'orientamento migliore di un legando che si lega ad una particolare proteina. Tuttavia, siccome sia il ligando che la proteina sono flessibili, l'analogia migliore è quella della "mano nel guanto".[4] Durante il corso del processo, il ligando e la proteina modificano la loro conformazione per ottenere l'orientamento migliore e questo tipo di aggiustamento di conformazione è chiamato “adattamento indotto”.[5]

Lo scopo del docking molecolare è di simulare tramite computer il processo di riconoscimento molecolare. Il fine è quindi quello di ottenere una conformazione ottimizzata contemporaneamente per la proteina, il ligando e il loro orientamento relativo tale che l'energia libera del sistema sia minimizzata.

Approcci al docking

[modifica | modifica wikitesto]

Due approcci sono particolarmente popolari nella comunità del docking molecolare. Un approccio usa una tecnica che descrive la proteina e il ligando come superfici complementari.[6][7] Il secondo approccio simula il processo di docking nel quale sono calcolate le energie di interazione della coppia ligando-proteina.[8] Entrambi gli approcci hanno vantaggi significativi, ma anche delle limitazioni. Entrambi gli aspetti vengono sottolineati qui di seguito.

Complementarità della forma

[modifica | modifica wikitesto]

I metodi di corrispondenza geometrica/complementarità della forma, descrivono una serie di proprietà che rendono possibile il legame tra proteina e ligando.[9] Queste proprietà possono includere i descrittori superficie molecolare/superficie complementare. In questo caso la superficie molecolare del recettore è descritta in termini di area superficiale accessibile al solvente e la superficie molecolare del ligando è descritta in termini della sua superficie di legame. La complementarità tra le due superfici risulta dalla descrizione della corrispondenza della forma che può aiutare a trovare la pose complementare del docking tra le molecole del recettore e del ligando. Un ulteriore approccio è descrivere le proprietà idrofobiche della proteina utilizzando turns negli atomi della catena principale. Un ulteriore approccio è utilizzare una tecnica Fourier shape descriptor.[10][11][12] Se gli approcci basati sulla complementarità della forma sono tipicamente veloci e robusti, non possono di solito modellare accuratamente i movimenti o i cambiamenti dinamici nella conformazione ligando/proteina, nonostante recenti sviluppi permettano a questi metodi di investigare la flessibilità del ligando.

I metodi di complementarità della forma possono analizzare in pochi secondi varie migliaia di ligandi e verificare se possono legarsi al sito attivo della proteina; inoltre normalmente sono scalabili alle interazioni proteina-proteina.

La simulazione del processo di docking è un processo molto più complicato. In questo approccio, la proteina e il ligando sono separati da una distanza fisica. Il ligando trova la sua posizione nel sito attivo della proteina dopo un certo numero di "mosse" nel suo spazio conformazionale. Le mosse incorporano le trasformazioni del corpo rigido come traslazioni e rotazioni, ma anche cambiamenti interni della struttura del ligando come rotazione e torsione. Ogni mossa nello spazio conformazionale del ligando induce un costo energetico totale del sistema e quindi dopo ogni mossa viene ricalcolata l'energia totale del sistema.

Meccanismi del docking

[modifica | modifica wikitesto]

Per effettuare un docking il primo requisito è la struttura della proteina. Normalmente la struttura è stata determinata tramite una tecnica chiamata cristallografia a raggi X o più raramente viene utilizzata la NMR di proteine. La struttura della proteina e un database di potenziali ligandi servono come input al programma. Il successo del docking dipende da due componenti: l'algoritmo di ricerca e la funzione di score.

Algoritmo di ricerca

[modifica | modifica wikitesto]

Lo spazio di ricerca in teoria consiste in ogni possibile orientamento e conformazione della proteina accoppiata al ligando. Tuttavia in pratica con le correnti risorse computazionali è impossibile esplorare esaustivamente tutto lo spazio di ricerca comprendente tutte le possibili distorsioni di una molecole (le molecole sono dinamiche possono avere un insieme di stati conformazionali) e tutte le possibili orientazioni del ligando relativo alla proteina ad un dato livello di granularità. La maggior parte dei programmi di docking in uso prevedono un ligando flessibile e alcuni tentativi per agganciarsi ad un recettore flessibile nella proteina. Ogni "snapshot" della coppia è chiamato pose.

Ci sono una serie di strategie di ricerca che vengono applicate al ligando e al recettore, queste includono:

Flessibilità del ligando

[modifica | modifica wikitesto]

Le conformazioni del ligando possono essere generate in assenza del recettore a cui successivamente verrà applicato il docking[13] o conformazioni che possono essere generate al volo in presenza del recettore,[14] o con una flessibilità rotazionale di ogni angolo diedro utilizzando il docking basato sui frammenti.[15] La valutazione dell'energia del campo di forza è utilizzata molto spesso per scegliere le conformazioni ragionevoli dal punto di vista energetico,[16] ma possono essere utilizzati anche metodi basati sulla conoscenza.[17]

Flessibilità del recettore

[modifica | modifica wikitesto]

La capacità computazionale si è incrementata notevolmente negli ultimi anni, rendendo possibile l'uso di tecniche più sofisticate e metodi computazionalmente intensivi. Tuttavia la flessibilità del recettore rimane un problema complesso. La ragione principale per questa difficoltà è il grande numero di gradi di libertà che bisogna prendere in considerazione in questo tipo di calcoli. Non considerarli, tuttavia, conduce a risultati poco significativi.[18]

Per emulare la flessibilità del recettore vengono utilizzate molte strutture statiche sperimentalmente determinate per la stessa proteina in diverse conformazioni.[19] In alternativa, per generare una conformazione energetica proteica ragionevole vengono utilizzate le conoscenze sulla struttura proteica delle catene di amminoacidi che circondano la cavità di legame.[20][21]

Funzione di score

[modifica | modifica wikitesto]

La funzione di score prende in input una "pose" e restituisce un numero che indica la possibilità che la pose rappresenti un'interazione di legame favorevole.

La maggior parte delle funzioni di score sono basate sul force field della meccanica molecolare che stima l'energia di una pose; una bassa (negativa) energia indica un sistema stabile e quindi una possibile interazione di legame. Un approccio alternativo è di creare un potenziale statistico per le interazioni da un grosso database di complessi proteina-ligando come il Protein Data Bank e valutarne il risultato.

Grazie alla cristallografia a raggi X sono state risolte un gran numero di strutture di proteine e ligandi con alta affinità, ma i ligandi con bassa affinità tendono ad essere sfuggenti a causa del fatto che sono meno stabili e quindi più difficili da cristallizzare. Le funzioni di score che si ottengono utilizzando questi dati possono simulare correttamente ligandi ad alta affinità, ma possono anche fornire plausibili dock per ligandi che non legano. Questo porta ad un gran numero di falsi positivi, come ad esempio un ligando che si lega ad una proteina, anche se questo non si verifica nella fase sperimentale.

Un modo per limitare il numero di falsi positivi è quello di ricalcolare l'energia delle migliori pose utilizzando delle tecniche (potenzialmente) più accurate, ma computazionalmente più intense come i metodi Generalized Born o Poisson-Boltzmann.[8]

Un'interazione di legame tra una piccola molecola ligando e una proteina enzima può attivare o inibire l'enzima. Se la proteina è un recettore, il legame col ligando può risultare agonista o antagonista. Il docking è comunemente usato nel campo della progettazione del farmaco — la maggior parte dei farmaci sono delle piccole molecole organiche a cui il docking può essere applicato, soprattutto ai campi indicati qui sotto:

  • hit identification – il docking combinato con una funzione di score può essere usato per monitorare velocemente grandi database di potenziali farmaci in silico, per identificare le molecole che più probabilmente si leghino al bersaglio di interesse.
  • lead optimization – il docking può essere usato per predire dove e in quale orientamento relativo un ligando si lega ad una proteina. Quest'informazione può essere utilizzata per creare analoghi più potenti e selettivi.
  • Biorisanamento – Il docking proteina ligando può essere utilizzato per predire gli inquinanti che possono essere degradati dagli enzimi.[22]
  1. ^ (EN) Lengauer T, Rarey M, Computational methods for biomolecular docking, in Curr. Opin. Struct. Biol., vol. 6, n. 3, 1996, p. 402–6, DOI:10.1016/S0959-440X(96)80061-3, PMID 8804827.
  2. ^ (EN) Kitchen DB, Decornez H, Furr JR, Bajorath J, Docking and scoring in virtual screening for drug discovery: methods and applications, in Nature reviews. Drug discovery, vol. 3, n. 11, 2004, pp. 935–49, DOI:10.1038/nrd1549, PMID 15520816.
  3. ^ Docking molecolare, su pianetachimica.it. URL consultato il 12 giugno 2012 (archiviato dall'url originale il 21 maggio 2012).
  4. ^ (EN) Jorgensen WL, Rusting of the lock and key model for protein-ligand binding, in Science, vol. 254, n. 5034, 1991, pp. 954–5, DOI:10.1126/science.1719636, PMID 1719636.
  5. ^ (EN) Wei BQ, Weaver LH, Ferrari AM, Matthews BW, Shoichet BK, Testing a flexible-receptor docking algorithm in a model binding site, in J. Mol. Biol., vol. 337, n. 5, 2004, pp. 1161–82, DOI:10.1016/j.jmb.2004.02.015, PMID 15046985.
  6. ^ (EN) Meng EC, Shoichet BK, Kuntz ID, Automated docking with grid-based energy evaluation, in Journal of Computational Chemistry, vol. 13, n. 4, 2004, pp. 505–524, DOI:10.1002/jcc.540130412.
  7. ^ (EN) Morris GM, Goodsell DS, Halliday RS, Huey R, Hart WE, Belew RK, Olson AJ, Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function, in Journal of Computational Chemistry, vol. 19, n. 14, 1998, pp. 1639–1662, DOI:10.1002/(SICI)1096-987X(19981115)19:14<1639::AID-JCC10>3.0.CO;2-B.
  8. ^ a b Feig M, Onufriev A, Lee MS, Im W, Case DA, Brooks CL, Performance comparison of generalized born and Poisson methods in the calculation of electrostatic solvation energies for protein structures, in Journal of Computational Chemistry, vol. 25, 2004, pp. 265–84, DOI:10.1002/jcc.10378, PMID 14648625.
  9. ^ (EN) Shoichet BK, Kuntz ID, Bodian DL, Molecular docking using shape descriptors, in Journal of Computational Chemistry, vol. 13, n. 3, 2004, pp. 380–397, DOI:10.1002/jcc.540130311.
  10. ^ (EN) Cai W, Shao X, Maigret B, Protein-ligand recognition using spherical harmonic molecular surfaces: towards a fast and efficient filter for large virtual throughput screening, in J. Mol. Graph. Model., vol. 20, gennaio 2002, pp. 313–28, DOI:10.1016/S1093-3263(01)00134-6, PMID 11858640.
  11. ^ (EN) Morris RJ, Najmanovich RJ, Kahraman A, Thornton JM, Real spherical harmonic expansion coefficients as 3D shape descriptors for protein binding pocket and ligand comparisons, in Bioinformatics, vol. 21, n. 10, maggio 2005, pp. 2347–55, DOI:10.1093/bioinformatics/bti337, PMID 15728116.
  12. ^ (EN) Kahraman A, Morris RJ, Laskowski RA, Thornton JM, Shape variation in protein binding pockets and their ligands, in J. Mol. Biol., vol. 368, n. 1, aprile 2007, pp. 283–301, DOI:10.1016/j.jmb.2007.01.086, PMID 17337005.
  13. ^ (EN) Kearsley SK, Underwood DJ, Sheridan RP, Miller MD, Flexibases: a way to enhance the use of molecular docking methods, in J. Comput. Aided Mol. Des., vol. 8, n. 5, ottobre 1994, pp. 565–82, DOI:10.1007/BF00123666, PMID 7876901.
  14. ^ (EN) Friesner RA, Banks JL, Murphy RB, Halgren TA, Klicic JJ, Mainz DT, Repasky MP, Knoll EH, Shelley M, Perry JK, Shaw DE, Francis P, Shenkin PS, Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy, in J. Med. Chem., vol. 47, n. 7, marzo 2004, pp. 1739–49, DOI:10.1021/jm0306430, PMID 15027865.
  15. ^ (EN) Zsoldos Z, Reid D, Simon A, Sadjad SB, Johnson AP, eHiTS: A new fast, exhaustive flexible ligand docking system, in Journal of Molecular Graphics and Modelling, vol. 26, n. 1, luglio 2007, pp. 198–212, DOI:10.1016/j.jmgm.2006.06.002, PMID 16860582.
  16. ^ (EN) Wang Q, Pang YP, Preference of small molecules for local minimum conformations when binding to proteins, in Floyd Romesberg (a cura di), PLoS ONE, vol. 2, n. 9, 2007, pp. e820, DOI:10.1371/journal.pone.0000820, PMC 1959118, PMID 17786192.
  17. ^ (EN) Klebe G, Mietzner T, A fast and efficient method to generate biologically relevant conformations, in J. Comput. Aided Mol. Des., vol. 8, n. 5, ottobre 1994, pp. 583–606, DOI:10.1007/BF00123667, PMID 7876902.
  18. ^ (EN) Cerqueira NM, Bras NF, Fernandes PA, Ramos MJ, MADAMM: a multistaged docking with an automated molecular modeling protocol, in Proteins, vol. 74, n. 1, gennaio 2009, pp. 192–206, DOI:10.1002/prot.22146, PMID 18618708.
  19. ^ (EN) Totrov M, Abagyan R, Flexible ligand docking to multiple receptor conformations: a practical alternative, in Curr. Opin. Struct. Biol., vol. 18, n. 2, aprile 2008, pp. 178–84, DOI:10.1016/j.sbi.2008.01.004, PMC 2396190, PMID 18302984.
  20. ^ (EN) Hartmann C, Antes I, Lengauer T, Docking and scoring with alternative side-chain conformations, in Proteins, vol. 74, n. 3, febbraio 2009, pp. 712–26, DOI:10.1002/prot.22189, 18704939.
  21. ^ (EN) Taylor RD, Jewsbury PJ, Essex JW, FDS: flexible ligand and receptor docking with a continuum solvent model and soft-core energy function, in J Comput Chem, vol. 24, n. 13, ottobre 2003, pp. 1637–56, DOI:10.1002/jcc.10295, PMID 12926007.
  22. ^ Suresh PS, Kumar A, Kumar R, Singh VP, An in silico [correction of insilico] approach to bioremediation: laccase as a case study, in J. Mol. Graph. Model., vol. 26, n. 5, gennaio 2008, pp. 845–9, DOI:10.1016/j.jmgm.2007.05.005, PMID 17606396.

Voci correlate

[modifica | modifica wikitesto]

Altri progetti

[modifica | modifica wikitesto]