RECURSOS DE BIOINFORMÁTICA APLICADOS ÀS CIÊNCIAS ÔMICAS COMO GENÔMICA, TRANSCRIPTÔMICA, PROTEÔMICA, INTERATÔMICA E METABOLÔMICA - PDF
Download "RECURSOS DE BIOINFORMÁTICA APLICADOS ÀS CIÊNCIAS ÔMICAS COMO GENÔMICA, TRANSCRIPTÔMICA, PROTEÔMICA, INTERATÔMICA E METABOLÔMICA"
Luciana Maria do Pilar Campos Lencastre
1 RECURSOS DE BIOINFORMÁTICA APLICADOS ÀS CIÊNCIAS ÔMICAS COMO GENÔMICA, TRANSCRIPTÔMICA, PROTEÔMICA, INTERATÔMICA E METABOLÔMICA BIOINFORMATIC RESOURCES APPLIED ON THE OMIC SCIENCES AS GENOMIC, TRANSCRIPTOMIC, PROTEOMIC, INTERATOMIC AND METABOLOMIC Foued Salmen ESPINDOLA 1 ; Luciana Karen CALÁBRIA 1, ; Alexandre Azenha Alves de REZENDE 1, ; Boscolli Barbosa PEREIRA 1, ; Flávia Assumpção SANTANA 1, ; Isabel Marques Rodrigues AMARAL 1, ; Janaina LOBATO 1, ; Juliana Luzia FRANÇA 1, ; Justino Luiz MARIO 1, ; Leonardo Bruno FIGUEIREDO 1, ; Luana Pereira dos SANTOS- LOPES 1, ; Neire Moura de GOUVEIA 1, ; Rafael NASCIMENTO 1, ; Renata Roland TEIXEIRA 1, ; Taís Alves dos REIS ; Thaise Gonçalves de ARAÚJO 1, 1. Instituto de Genética e Bioquímica, Universidade Federal de Uberlândia, Campus Umuarama, Uberlândia, MG, Brasil. Programa de Pós-Graduação em Genética e Bioquímica da Universidade Federal de Uberlândia. Programa de Pós-Graduação em Odontologia da Universidade Federal de Uberlândia RESUMO: As ciências ômicas tratam da análise global dos sistemas biológicos, integrando diferentes áreas do conhecimento, como a bioquímica, genética, fisiologia e computação, com o objetivo de isolar e caracterizar genes, proteínas e metabólitos, assim como estudar as interações entre eles, com base em técnicas experimentais, softwares e bancos de dados. A bioinformática por sua vez, propõe novas formas de ciência baseada na experimentação in silico, sendo muito dinâmica na sua atualização e fornecendo a base para geração de novos dados e conhecimentos
2 que podem ser aplicados na pesquisa básica e na aplicada com o desenvolvimento de novos produtos e soluções. Este processo está intimamente relacionado à inovação tecnológica, que é conseguida unindo-se a biotecnologia e a bioinformática. Contudo, o objetivo desta revisão é apresentar uma pequena abordagem dos recursos de bioinformática aplicados às ciências ômicas, como genômica, transcriptômica, proteômica, interatômica, metabolômica, farmacogenômica, dentre outras. PALAVRAS-CHAVE: Ômica. Bioinformática. Biotecnologia. Bancos de dados ARTIGO DE REVISÃO Os dados gerados pelo sequenciamento dos genomas de diferentes organismos transformaram a biologia. A integração de várias áreas do conhecimento permitiu avançar os estudos em relação à genômica, os processos de transcrição das informações contidas nos genes, a transcriptômica, bem como a compreensão do conjunto dos produtos destes genes pela proteômica. No início desta década com o advento do genoma humano também se iniciava as discussões e as ações para uma nova era da biologia, a era pós-genômica. Neste contexto, promoveu-se o desenvolvimento e o aperfeiçoamento das técnicas que permitiram os avanços destas novas ciências ômicas (Figura 1), como a transcriptômica, proteômica e metabolômica, com o objetivo de isolar e caracterizar o RNA, as proteínas e os metabólitos, respectivamente; sendo possível devido também ao desenvolvimento da bioinformática. O termo ômicos refere-se à análise global dos sistemas biológicos. Além das citadas anteriormente, uma variedade de subdisciplinas ômicas têm surgido, cada uma com seu próprio conjunto de instrumentos, técnicas, softwares e base de dados. Entre as tecnologias ômicas que impulsionam estas novas áreas de investigação, mencionam-se as tecnologias de DNA e microarrays, a espectrometria de massas e uma série de outras tecnologias e instrumentação que permitiram uma alta capacidade de análise (WINGENDER et al., 00).
3 O domínio da bioinformática cresceu em paralelo e com a internet, em que a rápida análise de dados e a troca de informações sobre os códigos biológicos e computacionais estão em convivência harmônica, por meio de suas múltiplas ramificações, gerenciando e integrando bancos de dados aplicáveis, e construindo sistemas in silico para simulação de formas naturais e modificadas de produtos específicos. Todos os projetos de sequenciamento genômico realizados e em andamento, tanto de procariotos como de eucariotos, continuam a nos lembrar que o nosso conhecimento sobre o funcionamento de um organismo ou célula, a nível molecular, é realmente muito limitado. Desta forma, o aumento substancial de sequências e de informações produzidas pelo rápido avanço das ciências ômicas está ajudando a prover novos caminhos da exploração de textos pela bioinformática (YANDELL; MAJOROS, 00) Figura 1: Representação esquemática das principais ciências ômicas. 1 Tantos são os temas a serem discutidos sobre os vários aspectos da bioinformática, todos eles na fronteira do conhecimento. Por isso, apresentamos aqui uma pequena abordagem baseada nos seminários desenvolvidos no Programa de Pós-Graduação em Genética e Bioquímica da Universidade Federal de Uberlândia na disciplina de Bioinformática. O quadro 1 apresenta páginas da web com recursos de bioinformática.
4 Mineração de textos e dados (text mining e data mining) O PubMed contém mais de 1 milhões de artigos científicos, e em 00 aproximadamente artigos foram adicionados ao Medline, sendo que o número de artigos adicionados por ano aumentou mais de 0.000/ano entre 000 e 00. Além de dados publicados no formato de artigos, bancos on-line são criados sobre muitos aspectos da biologia, contendo informações sobre interações moleculares, vias metabólicas, caracterização de genes, compostos químicos, estrutura protéica, doenças e organismos. Esta produção de dados científicos em larga escala impossibilita a extração ou análise destes dados por especialistas através de métodos manuais tradicionais, e para solucionar este problema, na década passada, métodos de text mining foram desenvolvidos. Quadro 1: Lista de páginas da Web com recursos de bioinformática e que são abordadas neste artigo. NOME PFAM SANGER Blast CaM Target CAP CAS Clustal Cytoscape Drug DataBase Easy Align Entrez Protein GenBank GENE D Gene Ontology Google GoogleScholar HiMAP HoGenom INSDC Interpare InterPro KEGG NCBI Osprey PANTHER Pfam Phrap PHYLIP PIRSF ENDEREÇO calcium.uhnres.utoronto.ca/ctdb/ctdb/home.html
5 PRF PRINTS ProDom PROSITE PubChem PubMed RNAMOTIF SBBiotec SCOP SetupX SMART String SUPERFAMILY SwissProt TIGRFAMs Tree View UniProt UniProtKB O text mining é uma ferramenta de obtenção de dados não-estruturados, escritos em linguagem natural, extraídos a partir de um banco de dados estruturado, com o auxílio de algoritmos para análise de textos não-estruturados. O processo de text mining envolve três subáreas: recuperação, extração, questões e respostas da informação, permitindo a identificação de identidades biológicas e suas interações, facilitando a análise de dados. A primeira sub-área (information retrieval), e mais comum em biologia molecular, consiste na extração de documentos a partir de uma grande coleção. Neste caso, há dois tipos de estratégia de busca: uma é a partir da combinação de palavras-chave e busca usando documentos como referência, para selecionar outros documentos semelhantes. Como muitas palavras são encontradas com grande freqüência e levam ao encontro de documentos com pouca informação, elas são excluídas durante o processo de busca. Uma ferramenta de busca de dados amplamente usada é o sistema de Information Retrieval do EntrezPubMed fornecido pelo NCBI (NATIONAL CENTER FOR BIOTECHNOLOGY INFORMATION). O popular Google foi recentemente incorporado como ferramenta de busca específica para literatura acadêmica com o GoogleScholar, que recupera artigos científicos, livros e reportagens. Na figura estão esquematizados outros bancos de dados separados de acordo com a aplicação biológica.
6 Por outro lado, o data mining é uma das novidades da Ciência da Computação, que utiliza vastos repertórios para tentar descobrir se há algum conhecimento escondido neles. A definição mais importante de data mining foi elaborada em 1 por Fayyad et al. (1):... Data mining é um processo não-trivial de identificar, em dados, padrões válidos, novos, potencialmente úteis e ultimamente compreensíveis. Este processo vale-se de diversos algoritmos desenvolvidos recentemente que processam os dados e encontram padrões válidos e novos. Embora os algoritmos atuais sejam capazes de descobrir esses padrões, os analistas humanos são os principais responsáveis por essa determinação Figura : Uma visão do processamento da linguagem natural biológica (BioNLP) e aplicações de text mining na biologia. O tópico central é envolvido por sete círculos com suas aplicações correspondentes dadas por outras caixas posteriores: NER (nome da identidade reconhecida), SDI (informação seletiva disseminada), Information extration (obtenção de informações), Análise por microarrays, Recuperação de artigos, Interações protéicas e relações e BioNLP. Modificado de: Krallinger e Valencia (00). A evolução da informática conta um pouco sobre o surgimento do processo de data mining. Nos anos 0, os computadores tinham capacidade precária não dispondo de mecanismos eficientes para armazenamento de grandes volumes de dados. No início da década de 0 até 0, um grande avanço marcou os meios físicos de armazenamento de dados, o desenvolvimento de
7 1 1 softwares para o gerenciamento de dados, denominados de Sistemas Gerenciadores de Bancos de Dados, e logo após o surgimento do modelo relacional, permitindo rápida recuperação de dados dirigidos. Nos anos 0, surgem os Bancos de Dados Multidimensionais ou Data Warehouses que propiciam o processo analítico on-line (OLAP). A diferença entre o OLAP e o data mining é que no primeiro o analista gera as hipóteses que podem ser validadas ou negadas, e no segundo o próprio sistema gera as hipóteses. O processo geral de descoberta de conhecimento em banco de dados é composto por diversas etapas. As principais tarefas são associação, agrupamento e descoberta de regras de classificação. A tarefa de classificação pode ser realizada por algoritmos convencionais ou por métodos de inteligência artificial como, por exemplo, as redes neurais, algoritmos evolucionários, dentre outros. Além disso, trabalhos têm demonstrado a importância desta ferramenta para estudos científicos. Baseando-se na literatura e análises in silico, pesquisadores utilizaram o data mining e fizeram uma seleção de 1 candidatos à vacina contra o M. tuberculosis. Este repertório foi ranqueado para gerar uma lista com os melhores antígenos, selecionando genes que abrangem todos os estágios da infecção, sendo incorporados no rbcg ou vacinas baseadas em subunidades (ZVI et al., 00) Sequenciamento do DNA e a Genômica Os mecanismos envolvidos na expressão e interação dos genes, assim como a compreensão das redes funcionais estabelecidas pelas proteínas, fazem com que, no cenário científico atual, a genômica e a proteômica estejam cada vez mais em evidência. Quatro são os principais bancos de dados utilizados para as diferentes análises de nucleotídeos. Um deles é o INSDC (INTERNATIONAL NUCLEOTIDE SEQUENCE DATABASE) que disponibiliza um repertório de sequências e é resultado da associação de três bancos de dados parceiros, o DDBJ (DATA BANK OF JAPAN), o EMBL (EMBL NUCLEOTIDE SEQUENCE DATABASE) e o GenBank. Os registros da associação
8 EMBL/GenBank/DDBJ incluem genes individuais, genomas completos, RNAs, anotações, sequências expressas, cdnas e sequências sintéticas. Devido à sua designação como sendo um provedor de dados primários, o banco EMBL/DDBJ/GenBank é a fonte inicial de muitos bancos de dados em biologia molecular (TATENO et al., 00; KANZ et al., 00 e BENSON et al., 00). Um exemplo de banco de dados de sequências genômicas secundárias de nucleotídeos é o Ensembl, uma fonte compreensível de anotações estáveis, em que genes são anotados por evidências derivadas de proteínas conhecidas, cdnas e sequências expressas. Novos genes são determinados pelo sistema de construção de genes, incorporando uma variedade de métodos, incluindo homologia e predição pela aplicação do HMM (HIDDEN MARKOV MODEL) (HUBBARD et al., 00). O RefSeq (REFERENCE SEQUENCE) é um banco que disponibiliza sequências compreensíveis, integradas e não-redundantes, incluindo DNA genômico, transcritos e proteínas para diversos organismos (PRUITT et al., 00). Por outro lado, o Genome Review representa uma versão da sequência original de um cromossomo ou plasmídeo, com informações importadas de fontes que incluem o UniProt (UNIVERSAL PROTEIN RESOURCE), Gene Ontology (GO), projeto GOA (GO ANNOTATION), InterPro e HoGenom, além de serem disponibilizadas referências cruzadas com 1 bancos de dados (KERSEY et al., 00). Milhares de sequências são obtidas através de técnicas, como o shotgun que foi usada, entre outros, no genoma de Apis mellifera (THE HONEYBEE GENOME SEQUENCING CONSORTIUM, 00). As sequências obtidas são analisadas e reunidas por software, como CAP. Após a organização das sequências é verificada a existência de redundância, a identificação das regiões codificadoras e, a seguir, identificação de funções. Os níveis de redundância são comparados usando o software Phrap, que foi usado no genoma da cana de açúcar (VETTORE et al., 00). Além disso, para a identificação do gene, em seres eucariotos, dois softwares foram muito utilizados, o Glimmer (DELCHER et al., 1) e o Genemark (BORODOVSKY; MCININCH, ). Por meio deles se analisa as
9 janelas de leitura na sequência (ORFs), sendo cada uma alinhada e comparada com as de outras espécies conhecidas, depositadas em bancos de dados como o GenBank e o SwissProt. A anotação funcional se dá através da comparação das sequências obtidas com as depositadas em bancos de dados como o GenBank e o Blast, sendo o último uma ferramenta mais amplamente utilizada para esse tipo de comparação (ALTSCHUL et al., ). A anotação é considerada completa quando o genoma está decodificado e minimamente anotado, com seus genes identificados e conferidos. Em 001, por junção da iniciativa pública (INTERNATIONAL HUMAN GENOME SEQUENCING CONSORTIUM, 001) e privada (CELERA) (VENTER et al., 001) realizou-se o projeto Genoma Humano. Os objetivos deste projeto foram identificar todos os genes estimados do DNA humano, determinar as sequências de bases, armazenar as informações em banco de dados e desenvolver ferramentas para a análise dos dados. Entretanto, de todos os genes que já foram seqüenciados, em média, apenas 0% codificam proteínas de função conhecida. Sobre o genoma humano em particular, o banco Genew como parte do HUGO (HUMAN GENOME ORGANIZATION), mantém um depósito de nomes e símbolos de genes, para se definir uma nomenclatura de dados submetidos por este Genoma (WAIN et al., 00). Outros projetos também foram desenvolvidos, como o Projeto EST (EXPRESSED SEQUENCE TAG), que ao invés de sequenciar todo o genoma de um organismo e depois tentar descobrir quais são seus genes, apenas os genes expressos pelo organismo são capturados e seqüenciados. Contudo, surgiu um importante problema computacional, o clustering. Ou seja, qual o modo de agrupar todos os ESTs que correspondem ao mesmo gene em um único grupo (cluster)? Assim, vários métodos foram propostos e a partir do agrupamento de tais dados, foram criados bancos de dados EST (dbest), que agrupam a informação de milhões de ESTs. Mesmo com o grande avanço provocado pela bioinformática, ainda persistem vários desafios, como por exemplo, o Gene Ontology que ainda não possui condições consistentes para diferenciar todos os processos biológicos, e os dados de microarray que ainda não possuem
10 uma grande reprodução de dados, resultando em um baixo poder estatístico. Muitos métodos estatísticos padrões falham por causa de problemas com o tamanho das amostragens podendo levar a uma desatualização dos bancos de dados comumente utilizados Filogenia Uma árvore filogenética é uma representação gráfica em forma de árvore, que retrata a história de parentesco entre as espécies, sendo que cada aresta da árvore representa uma mutação; cada vértice, um ponto na escala evolutiva; e cada folha uma espécie. Além disso, os indivíduos partem de um ancestral comum, a raiz (PROSDOCIMI et al., 00). A análise das relações filogenéticas pode ser dividida em duas etapas principais. A primeira consiste na construção das listas de caracteres para a análise, ou seja, o estabelecimento das homologias, sua codificação e ordenação, resultando na matriz de dados. Já na segunda etapa, se escolhe a árvore ótima por meio da análise numérica dessa matriz. A construção da matriz de dados é um passo muito importante na obtenção de filogenias, pois a qualidade da árvore filogenética final depende diretamente de sua construção. Com esse propósito, alguns programas estão disponíveis gratuitamente para auxiliar nas análises filogenéticas. O PHYLIP (THE PHYLOGENY INFERENCE PACKAGE) é um pacote para a análise filogenética criado por Joseph Felsenstein, na Universidade de Washington. Esse pacote pode resolver a maioria das análises filogenéticas existentes na literatura atual, incluindo métodos como o da parcimônia, o da matriz de distância e o likelihood. Como também, pode trabalhar com os seguintes tipos de entrada de dados: sequências moleculares, freqüência de gene, matriz de distância e características discretas. Outro programa de acesso livre é o Easy Align que foi desenvolvido em plataforma Windows, e possibilita o alinhamento de sequências de DNA e proteínas pelo critério de máxima parcimônia. Por outro lado, o Clustal X (THOMPSON et al., 1) é o algoritmo mais
11 utilizado em análises filogenéticas de dados moleculares, e para finalizar a análise, a construção da árvore é feita pelo software Tree View, muito utilizado nesta área RNAs e a transcriptômica Atualmente, estudos moleculares requerem a interação entre análises genômicas, celulares e dados de bioinformática, a qual apresenta, gradativamente, um papel essencial na geração de resultados aliados a alta tecnologia. A evolução da bioinformática, iniciada com análises de sequenciamento, tem oferecido avanços nas ciências ômicas, principalmente nas anotações dos transcriptomas, permitindo a interrelação entre o genoma funcional e a informação codificada. Contudo, o sequenciamento do genoma humano demonstrou que aproximadamente % de todos os produtos transcritos em humanos correspondem a RNAs não-codificantes (ncrna) (MATTICK; GARDEN, 001). Essa descoberta levantou questões referentes à diferenciação e desenvolvimento, tanto de espécies quanto de mecanismos moleculares individuais (MATTICK, 001). Embora o repertório de isoformas protéicas expressas nos organismos complexos seja significativamente incrementado pelo processamento pós-transcricional alternativo (splicing alternativo), a combinação desses polipeptídeos com os sinais ambientais fornecem informações insuficientes sobre os processos biológicos, de maneira que a maioria dos sistemas regulatórios é controlada por moléculas de RNA (MATTICK, 00). A predição da estrutura do RNA é normalmente baseada nas características termodinâmicas do folding da molécula ou na conservação filogenética das regiões de pareamento de bases e, nesse sentido, existem duas estratégias básicas na predição de ncrnas. A primeira é baseada na homologia genômica que existe ao longo da evolução e um dos métodos computacionais é o Rfam, capaz de alinhar mais de 00 famílias de ncrnas e determinar a estrutura predominante que resulta desse alinhamento. Outros softwares desenvolvidos para detectar ncrnas é o T RNA SCAN -SE, especializado na busca por trnas com elevada
12 sensibilidade e baixa proporção de falso positivos; e os PROMIR, MIR-ABELA, MIR-SCAN, MIRSEEKER e HARVESTER, dedicados na busca por mirnas em humanos, mamíferos, vertebrados, D. melanogaster e plantas, respectivamente. Esses programas utilizam tanto métodos heurísticos quanto modelos probabilísticos para capturar as sequências e características estruturais dos mirnas (MEYER, 00). Se a função do ncrna depende de elementos estruturais bem definidos e há uma sequência ortóloga de um organismo correlato, pode-se usar dois outros programas: o INFERNAL ou CMFINDER (YAO et al., 00). Caso exista apenas um único ncrna cuja estrutura funcional é determinada, mas nenhuma sequência equivalente é conhecida em outro genoma relacionado evolutivamente, pode-se buscar uma sequência alvo para regiões que são similares ao RNA conhecido. Nesse caso, encontra-se disponível o programa RSEARCH (KLEIN; EDDY, 00). Como o método compara apenas a sequência alvo a uma sequência única de RNA e não a um conjunto de sequências relacionadas filogeneticamente, sua sensibilidade é tipicamente menor. Portanto, como método alternativo, encontra-se o programa RNAMOTIF, o qual permite uma definição manual de motivos, capturando a sequência e as estruturas secundárias e terciárias (MACKE et al., 001). Alguns ncrnas não dependem de uma estrutura bem definida, permitindo-se optar por uma busca baseada essencialmente na similaridade de sequência conservada ao longo da evolução, sendo suficiente para identificá-los no genoma. Nesse sentido, existem programas especializados na busca por apenas homologia de sequência (perfil-hmms) e aqueles disponíveis para realizar o alinhamento e folding ao mesmo tempo. Dentre esses se encontram: o FOLDALIGN, capaz de detectar estruturas locais ao invés de identificar estruturas globais com vários loops (GORODKIN et al., 001), e o DYNALIGN, que reduz a complexidade computacional limitando o espaço de busca e o tamanho dos loops presentes nas estruturas internas do RNA (MATHEWS; TURNER, 00). As desvantagens desses dois últimos residem
13 na sua incapacidade de explicitar regiões que não adquirem uma determinada estrutura e na ineficiência do alinhamento baseado na estrutura da molécula. Em casos excepcionais, em que as sequências correspondem a determinados transcritos e sabe-se que a estrutura global desempenha um papel essencial na sua função, pode-se predizer essa estrutura para cada sequência individual através de softwares específicos. Nesse sentido, as duas ferramentas mais conhecidas são MFOLD e RNAFOLD, que predizem a estrutura mais estável a partir de uma sequência pré-estabelecida (ZUKER, 1; SCHUSTER et al., 1). A segunda estratégia inclui a predição ab-initio de ncrnas, o que constitui o maior desfio na busca dessas moléculas não-codificantes. Também existem algoritmos computacionais eficientes que objetivam predizer estruturas estáveis de RNA em grande escala genômica, como o RNAPLFOLD. Contudo, o potencial desse programa na detecção de ncrnas ainda não foi sistematicamente investigado (MEYER, 00) Esquema 1: Experimento de microarray. Uma das ferramentas mais utilizadas na análise transcriptômica é a tecnologia de microarrays que constitui uma das principais ferramentas para estudos de expressão gênica (SCHENA et al., 1), sendo muito aproveitada na avaliação de aspectos da biologia de sistemas e o estudo dos perfis de interação entre diversas biomoléculas (KITANO, 00).
14 O primeiro microarray surgiu em meados da década de e possuía sondas de cdna (SCHENA et al., 1). Com os aprimoramentos tecnológicos, no ano seguinte à sua publicação, pesquisadores apresentaram trabalhos com cerca de 00 sondas de arrays (SCHENA et al., 1; SHALON et al., 1), sendo que atualmente é comum encontrar trabalhos que utilizem dezenas de milhares de sondas. A AFFYMETRIX foi a empresa pioneira em microarrays, trabalhando com a metodologia de apenas um canal (uma cor). Na tecnologia de microarrays com lâminas de vidro, várias sequências de DNA conhecidas (sondas), são impressas em uma mesma lâmina. Nos arrays de duas cores, os mrnas são extraídos de células pertencentes às duas condições distintas e por meio de transcrição reversa, utilizando oligonucleotídeos marcados, o cdna é obtido. Os oligonucleotídeos são marcados com corantes fluorescentes (cianinas), sendo o corante Cy verde, e o Cy, vermelho como mostra o esquema 1. Após toda a experimentação biológica, as marcações são interpretadas por um software específico e os dados são analisados por ferramentas estatísticas Proteômica e estrutura de proteínas Como, onde, quando e por que são produzidas centenas de milhares de proteínas individuais em um organismo vivo? Como elas interagem entre si e com outras moléculas para construir uma célula? Como elas funcionam e conduzem o desenvolvimento e crescimento programado e interagem com os ambientes biótico e abiótico? Responder todas essas questões é o objetivo da proteômica, que como uma metodologia, deve ser considerada parte de uma análise integrativa e multidisciplinar em diferentes níveis, estendendo desde os genes até o fenótipo expresso nas proteínas. Estas análises devem envolver as tecnologias ômicas (genômica, transcriptômica, proteômica e metabolômica) bem como as técnicas de bioquímica clássica e biologia celular. No estudo completo das proteínas, integrando estrutura e função, os pesquisadores utilizam bancos de dados diversos que possam atender os diferentes ramos da proteômica. Um
15 dos mais usados é o banco de dados Entrez Protein, um depósito de sequências disponibilizado pelo NCBI e compilado através de uma variedade de fontes. O banco contêm as sequências de proteínas submetidas aos bancos PIR (PROTEIN INFORMATION RESOURCE) (WU et al., 00), UniProtKB/Swiss-Prot, PRF (PROTEIN RESEARCH FOUNDATION) e PDB. Outro, também muito utilizado é o UniProt, um catálogo de dados de sequências e funções de proteínas, mantido pelo consórcio UniProt. O consórcio é uma colaboração entre o SIB (SWISS INSTITUTE OF BIOINFORMATICS), o EBI (EUROPEAN BIOINFORMATICS INSTITUTE) e o PIR. O banco UniProt é compreendido por três componentes, o acurado UniProtKB (UNIPROT KNOWLEDGEBASE), que continuou o trabalho do UniProtKB/Swiss- Prot; o UniProtKB/TrEMBL (BOECKMANN et al., 00) e o PIR. O UniProtKB/Swiss- Prot é um banco anotado manualmente com informações extraídas da literatura e análises computacionais, contendo níveis mínimos de redundância e alto nível de integração com outros bancos de dados (BAIROCh et al., 00). Na análise de dados obtidos utilizando a eletroforese bidimensional, o banco de dados SWISS-DPAGE (HOOGLAND et al., 00) é o mais útil, pois armazena resultados experimentais que utilizam esta metodologia e acrescenta uma variedade de referências cruzadas com outros bancos de dados semelhantes, além do UniProtKB/Swiss-Prot. No entanto, se o objetivo é descrever a função molecular, o contexto biológico e a localização celular do produto gênico, o Gene Ontology é o mais indicado (CAMON et al., 00). O grande desafio enfrentado por estudiosos e bioinformatas é descobrir qual a estrutura tridimensional adotada pelas proteínas a partir da estrutura primária. No entanto, as ferramentas in silico disponíveis atualmente ainda não são totalmente confiáveis. Os métodos experimentais utilizados para obtenção da estrutura tridimensional são cristalografia por difração de raio-x e ressonância magnética nuclear. Entretanto, esses métodos podem ser onerosos e de difícil execução, além de apresentarem limitações técnicas. Estas e outras dificuldades fazem com que
16 a quantidade de estruturas de proteínas decifradas ainda compõe uma pequena fração do total de proteínas existentes (PROSDOCIMI et al., 00). Um método alternativo e não-experimental é a modelagem molecular, baseada em conhecimentos estereoquímicos dos aminoácidos. Uma das maneiras de se fazer a modelagem molecular é através da homologia entre sequências, em que uma delas já possui forma tridimensional definida. O primeiro passo é a pesquisa de proteínas homólogas em bancos de dados de estruturas de proteínas como o PDB (PROTEIN DATABASE BANK) (HULO et al., 00), que é uma colaboração entre o RCSB (RESEARCH COLLABORATORY FOR STRUCTURAL BIOINFORMATICS), o MSD-EBI (MACROMOLECULAR STRUCTURAL DATABASE) e o PDBj (PROTEIN DATA BANK OF JAPAN) (BERMAN et al., 000). A seguir, deve-se realizar o alinhamento das sequências de aminoácido das proteínas homólogas e a proteína-alvo, através do Clustal, por exemplo. A modelagem é realizada através de softwares como o Modeller, SWISS-MODEL, D-PSSM, dentre outros. Esses programas normalmente procuram encontrar a estrutura terciária que melhor se aproxime da disposição dos átomos das proteínas utilizadas como modelo, e ao mesmo tempo atenda às restrições físico-químicas (FORSLUND et al., 00). Outro tipo de modelagem é o Threading, que compara estrutura de uma proteína teste com a estrutura de outra proteína conhecida com uma pequena similaridade de sequência. Neste modelo é levada em consideração a distância entre os resíduos de aminoácidos, a estrutura secundária e as características físico-químicas (RATTEI et al., 00). 1 Interatômica O termo interatoma é menos bem definido do que genoma, proteoma ou transcriptoma, porque diferentes comunidades usam o termo interação de proteína para se referir de qualquer interação física às interações funcionais. No entanto, Sharan e Ideker (00) caracterizaram interatoma como o estudo do conjunto de interações macromoleculares, físicas e genéticas, e
17 uma das chaves da análise em larga escala é o alinhamento ou comparação global de duas ou mais redes (alinhamento múltiplo) para identificar regiões similares. No estudo interatômico, uma forma de avaliar a qualidade das redes de interação proteínaproteína é comparando as interações sugeridas com a localização subcelular ou as classes funcionais da proteína, tais como o Gene Ontology (BADER; HOGUE, 00). A suposição de tal análise é que os integrantes à interação devem pertencer à mesma categoria e a validez depende fortemente da escolha das classes. Além disso, a co-expressão dos genes correspondentes também é usada como um critério de avaliação (KEMMEREN et al., 00). Outra forma de validar interações é associar as proteínas dentro de uma via metabólica. Para isso existe um banco de dados, ou também nomeada de enciclopédia metabólica, chamada KEGG. Desde 1 têm sido desenvolvidos métodos bioinformáticos para descobrir comportamentos sistêmicos de informações bioquímicas e/ou genéticas. Os resultados são armazenados nesse banco de dados que possibilita pesquisa básica e aplicada das vias descobertas, como também interações com drogas (AOKI; KANEHISA, 00). As redes de interação proteína-proteína em larga escala têm sido avaliadas para diferentes modelos de organismos. Os primeiros estudos foram realizados por Sanchez et al. (1) em Drosophila melanogaster, por Ito et al. (001) e Ge et al. (001) em S. cerevisiae, e em 00 foi publicado por Rual et al. o interatoma humano. Além disso, interatomas específicos como o de desordens neurodegenerativas (LIMVIPHUVADH et al., 00) e de esquizofrenia (CAMARGO et al., 00) têm sido revelados. Na busca por domínios e famílias gerais o banco de dados mais utilizado é o Pfam. Nele, cada família é manualmente refinada e representada por dois alinhamentos múltiplos de sequência, dois perfis HMMs e um arquivo de anotação. Outro recurso muito utilizado, mas baseado na similaridade da sequência ou estrutura é o SCOP (STRUCTURAL CLASSIFICATION OF PROTEINS). Se a intenção é a busca por domínios específicos existem
18 bancos de dados disponíveis para os diferentes tipos como, por exemplo, para motif ligante de calmodulina, Calmodulin Target Database. Por outro lado, o Interpro (MULDER et al., 00) realiza buscas contra diferentes bancos de dados de domínios e famílias de proteínas, integrando os serviços oferecidos pelo Pfam, Uniprot, PROSITE, SMAR, PANTHER, PIRSF, SUPERFAMILY PRINTS, ProDom, GENE D e TIGRFAMs. Este banco de dados combina os diferentes métodos de reconhecimento de proteínas e na ausência da caracterização bioquímica, a predição de domínios pode ser um bom guia em direção à sua função (QUEVILLON et al., 00). O Interpare, da mesma forma que o Interpro, é um banco de dados para busca de domínios em conjunto com o PDB, SCOP, Uniprot e Swiss-Prot. Contudo, este banco de dados também utiliza um método computacional para identificar sítios de interação e moléculas ligantes, e classifica as proteínas pelos alvos de interação com drogas. Na construção das redes interatômicas, os softwares mais utilizados são: String, Cytoscape, Osprey e HiMAP. A escolha depende do organismo estudado, como também do banco de dados no qual está depositada a sequência estudada, se NCBI ou Swiss-Prot. Além disso, é importante que as redes interatômicas sejam feitas por diferentes softwares e depois comparadas, para confiabilizar os dados finais. O objetivo final do interatoma é unir as informações do genoma, proteoma e metaboloma, gerando informações que auxiliam no entendimento de funções e ações direcionadas a fármacos e moléculas biologicamente ativas. 0 1 Metabolômica O termo metabolômica foi criado e tem sido usado na última década para abranger o estudo do metabolismo sob perturbações ambientais e genéticas. No entanto, os primeiros trabalhos que envolviam técnicas relacionadas a metabólitos foram publicados há mais de 0 anos, a fim diagnósticos médicos (HORNING; HORNING, 11).
19 1 Os resultados em metabolômica são geralmente ricos em dados, sendo necessário o uso de ferramentas estatísticas e de bioinformática para avaliação e sistematização dos dados, em que propriedades bioquímicas e relações celulares podem ser mapeadas em plataformas de software que podem reforçar a interpretabilidade dos dados como, por exemplo, o SetupX que organiza e armazena os resultados de várias pesquisas em metabolômica. Em 00, uma série de relatos destacaram a importância de se fornecer informações. Entre elas, a base de dados ArMet, que descreve a arquitetura geral para metabolômica (JENKINS et al., 00) e MIAMet, que demonstra considerações sobre o mínimo de informações de um experimento em metabolômica (Bino et al., 00). Estas considerações têm sido concretizadas apenas parcialmente em bases de dados disponíveis sobre metabolômica de plantas (KOPKA et al., 00). Para uma série de compostos vegetais, várias empresas de agro-biotecnologia têm publicado dados dos metabólitos referentes ao valor nutricional das culturas. O mais abrangente é o CAS (CHEMICAL ABSTRACTS), que inclui informações sobre milhões de compostos, entre eles, metabólitos biogênicos. No entanto, este serviço vem com elevados encargos e não contém links para bases de dados genômicos Farmacogenômica Um objetivo nos estágios iniciais do desenvolvimento de fármacos é a identificação de um ou mais compostos bioativos. Um composto bioativo é qualquer substância que apresenta a atividade biológica que se procura (BUCHWALD; BODOR, 1). Qualquer composto com atividade farmacológica ou compostos similares normalmente possuem atividades parecidas, mas variam em sua potência e especificidade. Baseados em um composto bioativo, os cientistas investigam um grande número de moléculas parecidas de forma a otimizar as propriedades farmacológicas desejadas. Para uma busca sistemática, seria muito importante o entendimento de como as variações nas características estruturais e físico-químicas da família de moléculas estão relacionadas com
20 suas propriedades farmacológicas. O problema é que existem muitos descritores diferentes para caracterizar as moléculas. Eles incluem características estruturais, como a natureza e distribuição dos substituintes; características experimentais, como solubilidade em solventes aquosos e orgânicos, ou momentos dipolo; e características calculadas computacionalmente, como cargas parciais dos átomos. Estes fatores sejam eles de caráter eletrônico, hidrofóbico ou estérico, influenciam na interação do fármaco com a biofase, e na sua distribuição nos compartimentos que compõem o sistema biológico. Assim, dois fármacos com estruturas químicas semelhantes, diferenciando-se apenas por um átomo ou posição que este ocupa na molécula, podem apresentar diferenças quanto às suas propriedades físico-químicas e, conseqüentemente, quanto à atividade biológica, tanto do ponto de vista quantitativo como qualitativo (ESTRADA, 00). Os bancos de dados mais utilizados na análise da interação de fármacos ou compostos ativos e outra molécula biologicamente ativa são: KEGG, Drug DataBase e PubChem. A farmacogenômica surgiu em 1, da união da farmacogenética com a genômica e a biotecnologia (NEBERT; VESELL, 00), sendo definida como o estudo da expressão de genes individuais relevantes na susceptibilidade a doenças, bem como resposta a fármacos em níveis celular, tecidual, individual ou populacional (PIRAZZOLI; RECCHIA, 00). Como muitos outros ramos das ciências biomédicas, foi impulsionada pelos avanços da genômica, que conduziram às expectativas de que a segurança e a eficácia dos medicamentos seriam melhoradas pela personalização da terapêutica, com base nos dados genéticos (FONTANA et al., 00). Para o seu estudo, a farmacogenômica utiliza técnicas genômicas, como o sequenciamento de DNA, mapeamento genético e a bioinformática para facilitar as pesquisas na identificação das bases genéticas da variação inter-individual e inter-racial na eficácia, metabolismo e transporte com fármacos (MANCINELLI et al., 000). A genômica combinada com as ferramentas da bioinformática permite dissecar as bases genéticas das doenças multifatoriais e
Review Article 463 RECURSOS DE BIOINFORMÁTICA APLICADOS ÀS CIÊNCIAS ÔMICAS COMO GENÔMICA, TRANSCRIPTÔMICA, PROTEÔMICA, INTERATÔMICA E METABOLÔMICA BIOINFORMATIC RESOURCES APPLIED ON THE OMIC SCIENCES AS