{"id":64,"date":"2015-04-02T10:38:18","date_gmt":"2015-04-02T08:38:18","guid":{"rendered":"https:\/\/sperimentando.com\/?p=64"},"modified":"2015-04-10T20:23:50","modified_gmt":"2015-04-10T18:23:50","slug":"introduzione-alla-bioinformatica-le-sequenze-di-dna-sono-stringhe-di-caratteri","status":"publish","type":"post","link":"https:\/\/sperimentando.com\/?p=64","title":{"rendered":"Introduzione alla bioinformatica: le sequenze di DNA sono stringhe di caratteri – versione per studenti con preparazione in informatica"},"content":{"rendered":"
\"Immagine<\/a>

Immagine realizzata da Michael Str\u00f6ck e rilasciata con licenza libera http:\/\/en.wikipedia.org\/wiki\/File:DNA_Overview.png<\/p><\/div>\n

Il DNA \u00e8 un polimero organico costituito da monomeri chiamati nucleotidi. Tutti i nucleotidi sono costituiti da tre componenti fondamentali: un gruppo fosfato, un deossiribosio\u00a0e una base azotata. Il gruppo fosfato ed il ribosio si legano a formare lo scheletro della doppia elica, mentre i pioli sono costituiti da coppie di basi azotate. Le basi azotate che possono essere utilizzate nella formazione dei nucleotidi da incorporare nella molecola di DNA sono quattro: adenina (A),\u00a0guanina(G), citosina (C) e timina (T).<\/p>\n

L’ordine nella disposizione sequenziale dei nucleotidi costituisce l’informazione genetica.<\/p>\n

Possiamo in effetti rappresentare questa informazione come una sequenza di lettere, o stringa di caratteri. Le stringhe di caratteri che otteniamo possono essere analizzate in vari modi e restituire informazioni con un significato biologico.<\/p>\n

Il terminale<\/h2>\n

Iniziamo a familiarizzare con questa idea eseguendo dei semplici comandi su un file che contiene le sequenze nucleotidiche di alcuni geni.<\/p>\n

Per prima cosa stampate a video il nome della directory in cui vi trovate: il comando\u00a0pwd <\/strong>(print working directory) vi consente di farlo<\/p>\n

Quindi spostatevi nella vostra home con il comando cd ~ <\/strong><\/p>\n

Controllate utilizzando di nuovo pwd <\/strong>di avere eseguito il comando con successo e di essere in effetti nella home<\/p>\n

A\u00a0questo punto stampate a video la lista dei files contenuti nella vostra home con il comando\u00a0ls<\/strong>\u00a0(list)<\/p>\n

Nella home dovreste aver trovato la cartella\u00a0“esercitazione_bio”, stampate a video la lista dei files presenti in quella cartella digitando il comando\u00a0ls esercitazione_bio (ricordavi che il TAB \u00e8 un potentissimo strumento<\/a> di aiuto!!!)<\/strong><\/p>\n

Ora spostatevi nella cartella dell’esercitazione utilizzando\u00a0scrivendo\u00a0cd eser <\/strong>e poi premendo TAB<\/strong><\/p>\n

Iniziamo a scoprire come sono fatti i files contenenti le sequenze di DNA stampandone a video il contenuto: cat sequenze.fasta<\/a><\/strong><\/p>\n

Se il file \u00e8 molto lungo e quindi facciamo fatica a leggerlo possiamo iniziare stampando a video solo le prime righe, ad esempio digitando: head -n 20\u00a0sequenze.fasta<\/strong><\/strong><\/p>\n

Un comando simile ci consente di stampare solo le ultime righe del file: tail -n 10 sequenze.fasta<\/strong><\/strong><\/p>\n

Il file contenente le sequenze di DNA<\/h2>\n

Come avrete notato osservando con attenzione\u00a0le vostre\u00a0stampe a video, il\u00a0file “sequenze” ha\u00a0la seguente\u00a0struttura: all’inizio della riga c’\u00e8 un simbolo di maggiore seguito dal nome della sequenza;\u00a0poi c’\u00e8 un accapo; a partire dalla nuova riga c’\u00e8 una lunga fila di lettere a, c, g, t senza spazi disposta su pi\u00f9 righe\u00a0che rappresenta\u00a0appunto la sequenza di DNA.\u00a0Sapreste dire quante sequenze di DNA distinte identificate da un nome ci sono in quel file?<\/p>\n

Con il comando\u00a0grep “>”\u00a0sequenze.fasta<\/strong> <\/strong>stampate a video tutte le righe in cui \u00e8\u00a0presente\u00a0il simbolo “>” (grep \u00e8 un comando molto utile con le sequenze di DNA, vale la pena approfondire)<\/p>\n

E ora per sapere quante sono cosa fate, le contate con il dito dal video? mah no, potete lanciare due comandi in sequenza in questo modo:\u00a0grep “>”\u00a0sequenze.fasta<\/strong> | wc -l\u00a0<\/strong><\/p>\n

Per salvare il risultato che avete ottenuto in un file, basta direzionare il risultato del vostro comando in un file\u00a0grep “>”\u00a0sequenze.fasta<\/strong> | wc -l > numero_sequenze.txt<\/strong><\/p>\n

Il nostro esercizio di oggi: utilizziamo le sequenze di DNA per capire cosa \u00e8 successo ai nostri pazienti!<\/h2>\n

Cosa \u00e8 successo<\/h3>\n

 <\/p>\n