Analisi di dati di DNA microarray: reti di regolazione

I processi cellulari coinvolgono milioni di molecole che svolgono un ruolo coerente al fine di scambiare materia, energia e informazione con l’ambiente. Questi processi vengono regolati dai geni, la cui espressione è a loro volta regolata da una rete di interazioni fra altri geni, proteine e piccole molecole, detta rete di regolazione genica. Le reti di regolazione riguardano una svariata serie di segnali sia intracellulari, come le concentrazioni di mRNA, sia intercellulari, come gli ormoni, sia extracellulari, quali le variazioni delle condizioni ambientali (come gli stress termici). In maggior dettaglio, i geni vengono tradotti in proteine tramite un sistema di trascrizione (da DNA a mRNA) e traduzione (da mRNA a proteine) che è controllato da vari meccanismi (Nelson e Cox, 2002, Kohane et al., 2003). La trascrizione ha inizio in una particolare regione del DNA nota come regione del promotore, che è solitamente caratterizzata dalla presenza di segnali di riconoscimento per fattori di trascrizione, come la "TATA box" (TATAAA), "GC box", o "CAAT box" (CCAAT). Alla TATA box si lega una proteina, detta TBP (“TATA-binding” protein”), che, insieme ad altre proteine, consente all’enzima che permette la trascrizione, la RNA-polimerasi, di posizionarsi correttamente e venire attivata per iniziare la sintesi di RNA. I meccanismi di regolazione della trascrizione coinvolgono una serie di proteine, dette fattori di trascrizione, che possono bloccare il legame al promotore, o viceversa aumentare la probabilità che la RNA-polimerasi vi si leghi. In aggiunta a questi meccanismi, vi sono altri meccanismi di controllo che possono dar luogo, ad esempio, a splicing alternativi e modifiche post-trascrizionali. Tutti questi meccanismi influenzano la traduzione del messaggero in proteine. Per altro, proteine diverse possono essere formate a seguito dell’attivazione di un medesimo gene anche per effetto di modifiche post-traduzionali, anch’esse governate dalla presenza di proteine specifiche. La concentrazione di ogni proteina inoltre, è fortemente dipendente dalla presenza di altre molecole sia all’interno sia al di fuori della cellula, come i metaboliti o altre proteine; il ruolo svolto da queste molecole dipende da un complesso sistema noto come “signalling network”, che comprende recettori di membrana, proteine di “segnalazione” intra e intercellulari e i destinatari dei messaggi, ad esempio enzimi, altre proteine o complessi proteici (Roberts et al., 2000). E’ evidente quindi come lo studio delle reti di regolazione genica richieda di tener conto di aspetti genomici, proteomici e di regolazione metabolica. Non è perciò difficile apprezzare l’enorme complessità che queste reti possono raggiungere né la difficoltà di una loro modellizzazione matematica. La disponibilità dei DNA microarrays ha reso oggi possibile effettuare delle misure di espressione a livello dell’intero genoma e, grazie ai dati raccolti nel tempo, di osservarne la dinamica in diversi contesti sperimentali (Kohane et al., 2003). Grazie a queste osservazioni, è possibile inferire delle relazioni temporali fra le espressioni dei geni e, quindi, formulare delle ipotesi sulla struttura delle reti di regolazione. Questa attività è chiamata “reverse-engineering”, in quanto affronta il problema inverso di formulare delle ipotesi sul sistema dinamico che ha generato i dati sulla base dell’osservazione dei dati stessi (de Jong, 2002, Brazhnik et al., 2002).Nel corso degli ultimi anni sono stati proposti molti metodi per la generazione di reti di regolazione sulla base di dati di DNA microarrays (de Jong, 2002). Gli “ingredienti” della maggior parte di questi metodi sono una rete di interazioni, che esprime quali sono i legami fra i gen