UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE CENTRO DE BIOCIÊNCIAS PROGRAMA DE PÓS-GRADUAÇÃO EM CIÊNCIAS BIOLÓGIAS MESTRADO EM CIÊNCIAS BIOLÓGICAS JÉSSICA DE FÁTIMA VIANNA BIOQUÍMICA QUÂNTICA DA CAPREOMICINA E DA ESTREPTOMICINA EM COMPLEXO COM O RIBOSSOMO BACTERIANO NATAL-RN 2017 JÉSSICA DE FÁTIMA VIANNA BIOQUÍMICA QUÂNTICA DA CAPREOMICINA E DA ESTREPTOMICINA EM COMPLEXO COM O RIBOSSOMO BACTERIANO Orientador: Prof. Dr. Umberto Laino Fulco NATAL-RN 2017 Dissertação apresentada ao programa de Pós-Graduação em Ciências Biológicas da Universidade Federal do Rio Grande do Norte, como requisito para obtenção do título de Mestre em Ciências Biológicas. Universidade Federal do Rio Grande do Norte - UFRN Sistema de Bibliotecas – SISBI Catalogação de Publicação na Fonte. UFRN - Biblioteca Setorial do Centro de Biociências- CB Vianna, Jéssica de Fátima. Bioquímica quântica da capreomicina e da estreptomicina em complexo com o ribossomo bacteriano / Jéssica de Fátima Vianna. - Natal, 2017. 68 f.: il. Dissertação (Mestrado) - Universidade Federal do Rio Grande do Norte. Centro de Biociências. Programa de Pós-Graduação em Ciências Biológicas. Orientador: Prof. Dr. Umberto Laino Fulco. 1. Tuberculose - Dissertação. 2. DFT - Dissertação. 3. Antibióticos - Dissertação. I. Fulco, Umberto Laino. II. Universidade Federal do Rio Grande do Norte. III. Título. RN/UF/BSE-CB CDU 616.24-002.5 AGRADECIMENTOS Agradeço à minha mãe pelo grande apoio, dedicação e carinho em todos os momentos desde sempre, e sem dúvidas estará ao meu lado nos próximos desafios. A ela toda a minha gratidão. Ao prof. Dr. Umberto Fulco, por ter me acolhido em seu laboratório, e por toda a confiança depositada em mim. Muito obrigada! Ao prof. Dr. Jonas Oliveira, por ter aberto as primeiras portas para me inserir no mundo científico, e desde então me acompanhar em todos os projetos desenvolvidos. Ao Prof. Dr. Eudenilson Lins de Albuquerque pela contribuição à minha formação, ao compartilhar seus valiosos conhecimentos e experiências. Ao Vinícius Gonçalves pelo imensurável apoio, sempre disposto a ajudar em todas as etapas e por ter estado presente em todos os momentos. Assim como a Délio Delfino por ter me auxiliado a superar as principais limitações que me impediam de ter uma vida plena. A todos os professores da Pós-Graduação em Ciências Biológicas pela transmissão imprescindível de seus conhecimentos, em especial à Profª. Drª. Vanessa Soares Rachetti, por ser sempre um grande exemplo de competência e simpatia. Aos meus colegas do laboratório de Biofísica e Simulação Computacional, por toda a parceria, em especial a Katyanna, Dalila e Gabriela, por terem se tornado grandes amigas e companheiras dentro e fora do ambiente de trabalho, e a José Xavier por dar o suporte ao trabalho de todos. À CAPES pelo apoio financeiro. E a todos que passaram pelos meus caminhos e contribuíram de alguma forma para que esse trabalho pudesse ser desenvolvido, muito obrigada! “A tarefa não é tanto ver aquilo que ninguém viu, mas pensar o que ninguém ainda pensou sobre aquilo que todo mundo vê.” Arthur Schopenhauer RESUMO A tuberculose é uma doença bacteriana provocada pelo Mycobacterium tuberculosis, e de acordo com a Organização Mundial de Saúde, apenas em 2015 foram 10,4 milhões de novos casos relatados e 1,4 milhão de mortes. Cresce o número de casos de pacientes infectados com cepas resistentes aos antimicrobianos mais comumente utilizados, fazendo-se necessário uso de drogas de segunda-linha. A capreomicina e a estreptomicina encaixam-se nesse grupo, e são antibióticos que possuem como mecanismo de atuação a inibição da síntese proteica. Entretanto, seus mecanismos de ligação em seus sítios são distintos: a capreomicina é capaz de se ligar a ambas subunidades ribossomais (30S e 50S), enquanto que a estreptomicina liga-se à subunidade ribossomal menor (30S), e interage com alguns pontos da proteína S12. Através de dados cristalográficos e simulações computacionais, foi calculada a energia de interação da capreomicina e da estreptomicina com cada um dos resíduos constituintes de seus sítios utilizando a Teoria Funcional da Densidade (DFT) e do Método de Fracionamento Molecular com Capas Conjugadas (MFCC). Os resultados revelaram valores energéticos de cada nucleotídeo pertencente ao sítio de ligação desses dois medicamentos, como também dos aminoácidos da proteína S12 com os quais a estreptomicina interage. Assim, para a capreomicina na subunidade 30S, foram avaliados resíduos presentes em um raio de até 14 Å distantes do fármaco, totalizando 44 resíduos; e na subunidade 50S, 30 nucleotídeos foram analisados, e estavam distribuídos até o raio de 30 Å de distância. Com a estreptomicina foram levados em consideração 60 nucleotídeos distribuídos até 12,5 Å de distância da droga na subunidade 30S, e 25 aminoácidos da proteína S12 com até 15 Å de distância. Identificamos também as contribuições das ligações de hidrogênio e das interações hidrofóbicas nas interações fármaco-receptor; as regiões dos fármacos que mais contribuíram para as fixações desses em seus sítios de ligação; como também a identificação dos resíduos que são mais associados às mutações e consequente resistência. Palavras-Chave: Capreomicina. Estreptomicina. Tuberculose. Modelagem Molecular. Mecânica Quântica. MFCC. DFT. ABSTRACT Tuberculosis is a disease caused by Mycobacterium tuberculosis, and according to the World Health Organization, only in 2015 occurred 10.4 million new cases reported and 1.4 million deaths. The number of cases of patients infected with antimicrobial resistant strains most used is increasing, requiring the use of second- line drugs. Capreomycin and streptomycin are part of the group, and are antibiotics whose mechanism of action is the inhibition of protein synthesis. However, its binding mechanisms in their sites are distinct: capreomycin is able to bind to both ribosomal (30S and 50S) subunits, whereas streptomycin binds to the smaller ribosomal subunit (30S), and interacts with some points of S12 protein. Through crystallographic data and computational simulations, we calculated the interaction energy of capreomycin and streptomycin with each of the residues component of their sites using the Density Functional Theory (DFT) and Molecular Fractionation with Conjugated Caps (MFCC). The results showed energy values of each nucleotide belonging to binding site of these two drugs, as well as the amino acids of the S12 protein with which streptomycin interacts. Thus, for capreomycin in the 30S subunit, residues present in a radius of up to 14 Å distant from the drug, totaling 44 residues; and in the 50S subunit, 30 nucleotides were analyzed, and were distributed up to the 30Å radius distance. Regarding streptomycin, 60 nucleotides distributed up to 12.5 Å away from the drug in the 30S subunit, and 25 amino acids of the S12 protein with up to 15 Å were taken into account. We also identify the contributions of hydrogen bonds and hydrophobic interactions in drug-receptor interactions; the regions of the drugs that most contributed to the anchorages of these in their binding sites; as well as the identification of residues that are most associated with mutations and consequent resistance. Palavras-Chave: Capreomycin. Streptomicin. Tuberculosis, Molecular Modeling, Quantum Mechanics, MFCC, DFT LISTA DE FIGURAS Figura 1 Estrutura da parece celular bacteriana. Os componentes incluem (A) membrana plasmática, (B) peptidoglicano, (C) arabinogalactana, (D) manosídeo lipoarabinomanana, (E) proteínas associadas ao plasma e à parece celular, (F) ácidos micólicos, (G) glicolipideos de superfície associados a ácidos micólicos. Adaptado de KARAKOUSIS, Petros C.; BISHAI, William R.; DORMAN, Susan E.. Mycobacterium tuberculosiscell envelope lipids and the host immune response. Cellular Microbiology, v. 6, n. 2, p.105-116, 2004 ..................................................... 15 Figura 2 Duas vias levando à resistência. Adaptado de WHO guidelines for the programmatic management of drug-resistant tuberculosis, 2014 .............................. 20 Figura 3 O fluxo de informação genética do DNA ao RNA (transcrição) e do RNA à proteína (tradução). Fonte: ALBERTS, Bruce et al. Molecular biology of the cell. 5. ed. Newyork: Garland Science, 2008 ........................................................................ 21 Figura 4 Transcrição pela enzima RNA polimerase. Adaptado de ALBERTS, Bruce et al. Molecular biology of the cell. 5. ed. Newyork: Garland Science, 2008. ............ 22 Figura 5 Nucleotídeos. (a) Estrutura geral dos nucleotídeos. (b) O composto ancestral das bases pirimídicas e púricas. Fonte: NELSON, David L.; COX, Michael M.. Princípios de Bioquímica de Lehninger. 6. ed. Artmed, 2014 .............................. 23 Figura 6 Representação esquemática de ambas subunidades Ribossomais (30S e 50S), e os sítios de ligação para o RNA mensageiro e para os Transportadores com os sítios E, P e A. Adaptado de ALBERTS, Bruce et al. Molecular biology of the cell. 5. ed. Newyork: Garland Science, 2008 .................................................................... 24 Figura 7 Tradução de uma molécula de RNAm. Adaptado de ALBERTS, Bruce et al. Molecular biology of the cell. 5. ed. Newyork: Garland Science, 2008. ..................... 25 Figura 8 (A) Subunidades 30S e 50S em complexo com a Capreomicina (CAP), e (B), Subunidade 30S em complexo com a Estreptomicina (SRY), com destaque em amarelo para a proteína ribossomal S12. ................................................................. 33 Figura 9 Estruturas moleculares da Estreptomicina (A) e da Capreomicina (B), com a localização das cargas presentes em pH fisiológico e enumeração de suas regiões. ..................................................................................................................... 34 Figura 10 Esquematização dos subsistemas formados para aplicação do MFCC. C*2C*1: Nucleotídeos antecessores ao resíduo em estudo; Nn: Nucleotídeo em estudo; C*1C*2: Nucleotídeos sucessores ao Nn; .................................................... 36 Figura 11 Gráficos com os valores energéticos dos 20 resíduos mais importantes para capreomicina (30S e 50S) e estreptomicina (30S), e os aminoácidos da proteína S12 com os maiores e menores resultados. As distâncias de cada resíduo em relação à droga são apresentadas em angstrons (Å), com o detalhamento dos grupos funcionais dos resíduos e dos ligantes mais próximos. ................................. 43 Figura 12 Resíduos capazes de interagir através de ligações de hidrogênio nas subunidades 30S e 50S com a capreomicina (Cap). ................................................ 46 Figura 13 Resíduos capazes de interagir através de ligações de hidrogênio nas subunidades 30S e proteína S12 com a estreptomicina (SRY). ................................ 47 Figura 14 Regiões da capreomicina (Cap) com interações na subunidade 30S e 50S: (i) Região I, (ii) Região II, (iii) Região III, (iv) Região IV e (v) Região V. ............ 49 Figura 15 Regiões da estreptomicina (SRY) com interações na subunidade 30S e proteína S12: (i) Região I, (ii) Região II e (iii) Região III. ........................................... 50 Figura 16 Resíduos comuns entre a capreomicina (CAP) e estreptomicina (SRY), e suas respectivas energias. ........................................................................................ 51 Figura 17 Estrutura química da família de antibióticos tuberactimicinas (A) e estrutura da Capreomicina (B). Fonte: Adaptado de WANK, H., ROGERS, J., DAVIES, J, SCHROEDER, R. Peptide Antibiotics of the Tuberactinomycin Family as Inhibitors of Group I Intron RNA Splicing. J. Mol. Biol, 236, 1001- 1010, 1994. ........ 52 LISTA DE TABELAS Tabela 1 Energias de interação da Capreomicina em complexo com as subunidades 30S e 50S; raio de distância, e as regiões e grupos funcionais com os quais há as menores distâncias fármaco-resíduo. ....................................................................... 40 Tabela 2 Energias de interação da Estreptomicina com a subunidade 30S e proteína S12; raio de distância, e as regiões e grupos funcionais com os quais há as menores distâncias fármaco-resíduo. ...................................................................................... 41 SUMÁRIO 1 INTRODUÇÃO ................................................................................................................ 14 1.1 CAPREOMICINA .................................................................................................... 26 1.2 ESTREPTOMICINA ............................................................................................... 27 1.3 MODELAGEM MOLECULAR ............................................................................... 27 2 OBJETIVOS .................................................................................................................... 30 2.1 OBJETIVOS GERAIS ............................................................................................ 30 2.2 OBJETIVOS ESPECÍFICOS ................................................................................ 30 3 METODOLOGIA ............................................................................................................. 32 3.1 DADOS CRISTALOGRÁFICOS ........................................................................... 32 3.2 OTIMIZAÇÃO DE GEOMETRIA .......................................................................... 33 3.3 MFCC ....................................................................................................................... 35 4 RESULTADOS E DISCUSSÕES ................................................................................ 39 5 CONCLUSÕES ............................................................................................................... 57 6 PERSPECTIVAS ............................................................................................................ 60 7 REFERÊNCIAS .............................................................................................................. 62 INTRODUÇÃO 14 1 INTRODUÇÃO A tuberculose é uma doença infecciosa causada pelo Mycobacterium tuberculosis, e continua sendo um problema de saúde global apesar de todos os esforços envolvendo entidades governamentais e não-governamentais para reduzir sua incidência. Apenas em 2015 foram 10,4 milhões de novos casos relatados e 1,4 milhão de mortes (WHO, 2016). O Brasil encontra-se entre os 20 países com as maiores taxas de incidência, com uma estimativa de 90.000 casos em 2015. O gênero Mycobacterium engloba microrganismos em forma de bastão, imóveis, e não formadores de esporos. Possuem paredes ricas em lipídeos, que formam uma superfície hidrofóbica e produzem resistência a vários desinfetantes, detergentes e antibióticos comuns; além de ser um dos determinantes para o lento crescimento bacteriano. Assim, a parede celular dessas bactérias tem uma grande complexidade, e em sua composição é possível encontrar: uma membrana plasmática interna, revestida por uma espessa camada de peptidoglicano; proteínas; manosídeos; fosfatidilinositol; lipoarabinomanana (LAM) e outros lipídeos; glicolipídeos e peptidoglicolipídeos (MURRAY; ROSENTHAL; PFALLER, 2013) (Figura 1). Esses componentes são responsáveis, por dentre outras funções, manter a integridade bacteriana e auxiliar no metabolismo, além de estarem envolvidos com o despertar a resposta imune quando o microrganismo se instala em um hospedeiro. O agente etiológico da tuberculose é transmitido por via aérea, no momento em que um indivíduo com a tuberculose pulmonar ou faríngea tosse, espirra ou fala, por exemplo. Indivíduos saudáveis que estejam próximos podem inalar as bactérias lançadas no ar, e serem infectados (CHAPARAS; GOOD, 1982). Os fatores que mais contribuem para uma maior probabilidade de transmissão do M. tuberculosis, são: número de organismos expelidos no ar, pois pequenas partículas contendo de um a três bacilos que conseguem chegar aos espaços alveolares são capazes de estabelecer a infecção; a concentração desses microrganismos no ar, determinada pelo volume do espaço e sua ventilação; o tempo de exposição de uma pessoa em um ambiente com ar contaminado; e o estado do sistema imune do indivíduo exposto (ATS; CDC; IDSA, 2000). 15 M. tuberculosis pode ingressar nas vias aéreas e chegar aos alvéolos, onde é fagocitado por macrófagos alveolares, e consegue escapar da sua eliminação ao inibir a fusão do fagossomo com os lisossomos. Ao mesmo tempo, é capaz de induzir o fogossomo a fusionar-se a outras vesículas intracelulares, fornecendo acesso a nutrientes e facilitando a sua replicação. (MURRAY; ROSENTHAL; PFALLER, 2013). O estabelecimento uma infecção, dependerá da virulência das cepas e da resposta imune do hospedeiro (EDWARDS; KIRKPATRICK, 1986; DANNENBERG, 1989). Figura 1 Estrutura da parede celular bacteriana. Os componentes incluem (A) membrana plasmática, (B) peptidoglicano, (C) arabinogalactana, (D) manosídeo lipoarabinomanana, (E) proteínas associadas ao plasma e à parece celular, (F) ácidos micólicos, (G) glicolipideos de superfície associados a ácidos micólicos. Fonte: Adaptado de KARAKOUSIS, Petros C.; BISHAI, William R.; DORMAN, Susan E.. Mycobacterium tuberculosiscell envelope lipids and the host immune response. Cellular Microbiology, v. 6, n. 2, p.105-116, 2004. Após a infecção primária, a tuberculose pode se desenvolver na forma ativa (pulmonar ou extrapulmonar), ou pode ficar em estado latente durante uma parte ou em toda a vida do paciente. O bacilo consegue se difundir via sistemas linfático e 16 sanguíneo, e alcançar sítios mais distantes e se instalar em rins, ossos e cérebro, gerando a forma extrapulmonar (SWEANY, 1955). A destruição tissular provocada pela tuberculose está mais associada com a ação do sistema imune sobre as áreas infectadas (CRUZ-KNIGHT; BLAKE-GUMBS, 2013). O sistema imune do hospedeiro responde a essa infecção liberando citosinas e toxinas, e as células humanas infectadas com essas bactérias em seu interior (macrófagos alveolares, células epiteliais e células gigantes de Langhans), formam o centro uma massa necrótica circundada por uma densa parede de macrófagos e células do tipo T, CD4, CD8 e células Natural Killer (NK). Essa estrutura é chamada de granuloma e é uma tentativa de prevenir a disseminação da bactéria, originando uma necrose caseosa. Com isso, se houver uma infecção com grande quantidade da bactéria, um grande granuloma pode ser formado (MURRAY; ROSENTHAL; PFALLER, 2013). Pessoas imunocompetentes que foram expostas a esse bacilo, conseguem controlar sua replicação e infecção, porém poucos conseguem eliminar completamente todos os agentes invasores (AHMAD; MOKADDAS, 2014). Após a instalação, a bactéria pode permanecer num estado latente por um longo período e posteriormente ser reativada, quando a capacidade de resposta do sistema imune cai como resultado de terapia imunossupressiva, ou mesmo avanço da idade, levando à tuberculose ativa (AHMAD; MOKADDAS, 2014; CRUZ-KNIGHT; BLAKE-GUMBS, 2013). Outros fatores que aumentam o risco de reativação incluem: infecção por HIV, contato com indivíduo infectado, paciente em tratamento por diálise, recepção de órgãos transplantados, silicose, situação carcerária e uso de drogas, por exemplo (WHO, 2015). A tuberculose latente é definida como um estado de resposta imune persistente a M. tuberculosis, sem manifestação clínica (MACK et al., 2009). Os grupos de risco nessa situação englobam crianças com menos de cinco anos que possuem contato com casos de tuberculose na forma ativa e pessoas com HIV (WHO, 2015). Estima-se que um terço da população mundial tenha essa forma de tuberculose (DYE, 1999). O diagnóstico da forma latente pode ser feito através do teste da Tuberculina na pele ou pelos ensaios de liberação de interferon-Gamma (IGRAs), uma vez que os exames bacteriológicos costumam ser negativos. Entretanto, esses exames não conseguem distinguir a infecção latente 17 com microrganismos viáveis, das infecções já curadas ou tratadas. Além disso, não fornece uma estimativa de potencial avanço para a forma ativa da doença, motivo pelo qual a tuberculose latente deve ser tratada a fim de evitar essa progressão (WHO, 2015). Os sintomas produzidos pela tuberculose são sistêmicos, e dentre os mais frequentes estão febre, diminuição de apetite, perda de peso, fraqueza, suores noturnos e mal-estar. Nos casos de tuberculose pulmonar, tosse é o sintoma mais comum, e com o desenvolvimento da doença há produção de escarro, que é comumente utilizado para diagnóstico. A tuberculose extrapulmonar é mais difícil de ser diagnosticada, já que os sítios geralmente são inacessíveis e as concentrações dos bacilos nas áreas são baixas. A vacina Bacillus Calmete-Guerin (BCG) é a única licenciada para a tuberculose atualmente, consistindo em cepas atenuadas de M. bovis com eficácia variando entre 0 - 80%. É recomendada para ser administrada em crianças residentes em áreas de risco, com tuberculose endêmica ou para crianças expostas à tuberculose em ambiente familiar. Ela é capaz de fornecer proteção contra severas formas extrapulmonares em crianças, mas não oferece tal proteção contra a forma pulmonar, que é a forma mais comum e com maiores índices de mortalidade (COLDITZ, 1994). Dessa forma, há necessidade de novas vacinas que sejam seguras, efetivas e acessíveis, o que representaria um grande avanço no controle da doença (MCSHANE et al., 2005). Além disso, a revacinação durante a adolescência, em uma população já vacinada durante a fase infantil, não aumenta a eficácia da vacina (RODRIGUES et al., 2005). As micobactérias podem ser isoladas a partir de uma grande variedade de amostras clínicas, incluindo amostras das vias respiratórias superiores (escarro, lavado brônquico, lavado broncoalvelar e biópsias brônquicas), urina, fezes, sangue, líquor, biópsias de tecidos e aspirados de qualquer órgão ou tecido (KRASNOW; WAYNE, 1969). Segundo as III Diretrizes para Tuberculose da Sociedade Brasileira de Pneumologia e Tisiologia, o diagnóstico da tuberculose é feito inicialmente por meio da radiografia de tórax e a pesquisa de Bacilos Álcool-Ácido Resistentes (BAAR) no escarro. As amostras podem ser semeadas em meios microbiológicos enriquecidos que favorecem o crescimento bacteriano, o que permite realizar o 18 isolamento e a identificação do M. tuberculosis, e as análises de sensibilidade medicamentosa. Inovações como testes moleculares possibilitam uma maior celeridade ao processo de identificação bacteriana, e baseiam-se na avaliação do material genético do microrganismo. A tuberculose em sua forma latente pode ser diagnosticada quando ainda não há estabelecimento da doença, mas o teste da tuberculina apresenta resultado positivo. Além de objetivar a redução dos riscos ao paciente com tuberculose, o tratamento dessa doença é necessário também para minimizar a transmissão do agente patológico para outras pessoas, e assim trazer benefícios tanto para o paciente quanto para a comunidade. Por ser um microrganismo com crescimento lento e grande capacidade de entrar em estado de latência, o tratamento para essa doença requer o uso combinado de várias drogas por período predominantemente longo. O esquema mais recomendado inicialmente inclui uso de medicamentos enquadrados como de primeira linha, que são isoniazida, rifampicina, pirazinamida e etambutol por dois meses, constituindo a fase inicial, seguidos por uma fase de continuação por quatro meses usando isoniazida e rifampicina (ATS; CDC; IDSA, 2003). Esses medicamentos possuem diversos mecanismos de ação. A isoniazida e etambutol interferem na síntese da parede bacteriana, afetando a síntese de ácidos micólicos, e a de arabinogalactano, respectivamente. A rifampicina liga-se à enzima RNA-polimerase e inibe a síntese de RNA, enquanto a pirazinamida, ainda não foi totalmente esclarecido seu mecanismo de ação (MURRAY; ROSENTHAL; PFALLER, 2013). O uso de múltiplos agentes é necessário para prevenir a emergência de populações bacterianas resistentes. A resistência do M. tuberculosis a alguns medicamentos aparece como um grande obstáculo para os programas de controle da tuberculose, e uma vez identificada, o esquema de tratamento é alterado passando a incluir drogas de segunda linha, que são mais tóxicas e menos efetivas que as de primeira-linha (ATS; CDC; IDSA, 2003). Os agentes de segunda-linha podem ser divididos em três grupos (de 2 a 4), de acordo com seus modos de ação, rotas de entrada e potência. O grupo 2 inclui aminoglicosídeos injetáveis (estreptomicina, canamicina e amicacina) e os polipeptídeos cíclicos (capreomicina e 19 viomicina). No grupo 3, estão fluoroquinolonas (ofloxacina, levofloxacina, moxifloxacina e gatifloxacina), e o grupo 4 reúne bacteriostáticos menos eficazes, mais caros e com mais efeitos tóxicos, que são representados por etionamida, protionamida, cicloserina, terizidona e ácido para-amino salicílico, que são usados de acordo com a susceptibilidade das cepas de M. tuberculosis (ATS; CDC; IDSA, 2003; MITCHISON; DAVIES, 2012) A capacidade que as micobactérias possuem de se multiplicar, mesmo na presença de isoniazida e rifampicina caracteriza a multirresistência. No ano de 2014, foram 480.000 de novos casos (WHO, 2015), enquanto que no ano de 2016 foram estimados 580.000, e desses, 111000 pacientes iniciaram tratamento com drogas de segunda-linha (WHO, 2016). Dentre essas estimativas, 9,7% incluíram a forma extensivamente resistente, definida quando há concomitantemente, bloqueio da atuação dos grupos: (1) rifampicina e isoniazida; (2) qualquer uma das fluoroquinolonas; (3) uma das drogas de segunda-linha injetáveis: capreomicina, canamicina ou amicacina (CAMINERO, 2006; CHAN et al., 2004; EKER et al., 2008). Resistência a drogas além das listadas também têm sido relatadas, e foram definidas como tuberculose totalmente resistente (KLOPPER et al., 2013). Há duas vias principais que levam ao desenvolvimento de resistência aos medicamentos: (i) via primária; e (ii) via secundária, que é adquirida. Essas vias estão interconectadas e possuem vários fatores contribuintes (WHO, 2014) (Figura 2). A resistência adquirida é resultado de um tratamento inadequado, incompleto ou improdutivo, que acaba permitindo a seleção de cepas com mutações que podem perder a susceptibilidade aos medicamentos disponíveis. Assim, um tratamento apropriado requer uma combinação de vários medicamentos com qualidade garantida. A resistência primária caracteriza os casos em que o paciente é infectado por uma cepa já resistente, pela mesma forma de transmissão das cepas susceptíveis (WHO, 2014). 20 Figura 2 Duas vias levando à resistência. Fonte: Adaptado de WHO guidelines for the programmatic management of drug-resistant tuberculosis, 2014. O esquema de tratamento da tuberculose multirresistente também consiste em uma poliquimioterapia, e a recomendação da WHO inclui uma fase intensiva (com duração de pelo menos oito meses) que pode ser integrada por pelo menos quatro drogas de segunda-linha, que podem ser escolhidas entre: pirazinamida; uma das fluoroquinolonas; uma das injetáveis (capreomicina, canamicina ou amicacina); etionamida; cicloserina ou PAS (ácido paraminosalicilico). O tempo total de 21 tratamento para essa condição pode chegar a vinte meses, variando conforme a resposta do paciente ao tratamento (WHO, 2011). A estreptomicina é recomendada em casos em que há resistência a todas as drogas de segunda-linha (WHO, 2014), o que reduz a amplitude de sua utilização como primeira opção. Estreptomicina e capreomicina enquadram-se dentre os medicamentos de segunda linha utilizados nos casos de multirresistência, e estão listados como uns dos medicamentos essenciais pela Organização Mundial de Saúde (OMS, 2015). Ambos são injetáveis e atuam inibindo a síntese proteica ao interagir com o ribossomo bacteriano. Os ribossomos são formados pela associação de RNAr e algumas proteínas. Nas bactérias, essas organelas possuem uma subunidade maior (50S) e uma menor (30S), onde a 50S possui os RNAr 23S e 5S, além de 31 proteínas; e a 30S contém o RNAr 16S e 20 proteínas (PUGLISI; GREEN; BLANCHARD, 2000). A síntese proteica é a maneira pela qual as células expressam suas informações genéticas. Para isso, as sequências de nucleotídeos presentes no DNA são copiadas na forma de moléculas de RNA (também compostas por nucleotídeos) pelo processo de transcrição, e posteriormente esse produto é traduzido em cadeias de aminoácido, que geram proteínas (ALBERTS et al., 2008) (Figura 3). Figura 3 O fluxo de informação genética do DNA ao RNA (transcrição) e do RNA à proteína (tradução). Fonte: ALBERTS, Bruce et al. Molecular biology of the cell. 5. ed. Newyork: Garland Science, 2008. 22 Durante a transcrição, uma das fitas do DNA age como molde para a síntese da molécula de RNA, e assim, a sequencia de nucleotídeos desta é determinada pela complementariedade com a fita-molde. Os nucleotídeos componentes são covalentemente ligados por meio da enzima RNA-polimerase, que catalisa a formação das ligações fosfodiéster à medida que os nucleotídeos são adicionados na direção 5' - 3' (ALBERTS et al., 2008; NELSON; COX, 2014) (Figura 4). Figura 4 Transcrição pela enzima RNA polimerase. Fonte: Adaptado de ALBERTS, Bruce et al. Molecular biology of the cell. 5. ed. Newyork: Garland Science, 2008. O RNA é um polímero linear de fita simples, constituído por nucleotídeos, que podem ser derivados de componentes purinas e pirimidinas (adenina, guanina, citosina e uracila), possuem em sua estrutura três componentes característicos: base nitrogenada, pentose e fosfato (Figura 5). Esses subcomponentes são ligados entre si pelo grupo 5’-fosfato de uma unidade nucleotídica com o grupo 3’-hidroxil do próximo, formando ligações fosfodiéster (ALBERTS et al., 2008). Os grupos fosfatos com um pKa próximo a zero, são completamente ionizados e carregados negativamente em pH neutro (NELSON; COX, 2014). 23 Figura 5 Nucleotídeos. (a) Estrutura geral dos nucleotídeos. (b) O composto ancestral das bases pirimídicas e púricas. Fonte: NELSON, David L.; COX, Michael M.. Princípios de Bioquímica de Lehninger. 6. ed. Artmed, 2014 O RNA originado da transcrição dos genes é chamado RNA mensageiro (RNAm), e carrega informações sobre a sequencia de aminoácidos a ser organizada. Na tradução, há o auxílio do RNA transportador (RNAt), que funciona como um adaptador capaz de reconhecer e se ligar ao códon do RNAm e, em um outro local da sua superfície, ao aminoácido correspondente. (ALBERTS et al., 2008). A partir dessa fase o processo ocorre nos ribossomos, que são grandes complexos catalíticos, compostos por proteínas e RNA ribossomais (RNAr) que se organizam em duas subunidades: uma menor, responsável por formar o complexo com o RNAt, e intermediando o pareamento com os códons do RNAm; e uma subunidade maior, que catalisa a formação de ligações peptídicas que unem os aminoácidos, para originar uma cadeia polipeptídica. Além disso, contém quatro sítios de ligação para as moléculas de RNA: um para o RNAm, e três para os RNAt, chamados de sítios A, P e E (Figura 6) (ALBERTS et al., 2008). 24 Figura 6 Representação esquemática de ambas subunidades Ribossomais (30S e 50S), e os sítios de ligação para o RNA mensageiro e para os Transportadores com os sítios E, P e A. Fonte: Adaptado de ALBERTS, Bruce et al. Molecular biology of the cell. 5. ed. Newyork: Garland Science, 2008 Uma vez iniciada a tradução, cada aminoácido é adicionado à cadeia de alongamento, em um ciclo de reações contendo de um modo geral, três grandes passos: (1) ligação do RNAt; (2) formação das ligações peptídicas, e (3) translocações nas subunidades maior e menor. A figura 7 ilustra o processo: no passo 1, um RNAt trouxe um aminoácido e ligou-se ao sítio A, enquanto que o sítio E se esvazia; passo 2, há catalisação da ligação peptídica no aminoácido trazido; passo 3, a subunidade maior transloca-se para frente, transferindo os RNAt para os sítios adjacentes; passo 4, a subunidade menor também se transloca para frente, resultando no esvaziamento do sítio A e transferindo os RNAt para os sítios de saída. Assim, o ciclo pode recomeçar e novos aminoácidos serem incorporados à cadeia crescente, até que todo o RNAm seja traduzido. 25 Figura 7 Tradução de uma molécula de RNAm. Fonte: Adaptado de ALBERTS, Bruce et al. Molecular biology of the cell. 5. ed. Newyork: Garland Science, 2008. Portanto, a síntese proteica é um processo vital, e medicamentos que atuem em alguma de suas etapas conseguem retardar ou impedir a proliferação de microrganismos como bactérias. Capreomicina e estreptomicina são alguns dos fármacos com esse mecanismo de ação, e são utilizados em esquemas de tratamento para cepas bacterianas resistentes. Assim, o estudo dos mecanismos moleculares envolvidos na atuação deles em seus sítios de ligação, é fundamental para uma melhor compreensão de suas farmacodinâmicas, como também para produzir dados necessários ao desenvolvimento de novos fármacos. A capreomicina possui um mecanismo de atuação ímpar, capaz de inibir o crescimento bacteriano ao bloquear sua síntese proteica, ligando-se à ambas 26 subunidades ribossomais (STANLEY et al., 2010), interferindo na formação do complexo de iniciação na subunidade 30S, e consequentemente bloqueando a translocação do RNAt do sítio A para o sítio P (CAMMARATA, 1973; STANLEY et al., 2010; MODOLFLL; VAZQUEZ, 1977), enquanto a estreptomicina liga-se predominantemente à subunidade 30S, aumentando a afinidade do RNAt pelo sítio A, estabilizando essas moléculas no ribossomo e evitando a progressão do processo, ou ainda provocando erros de leitura (NOLLER, 1991; PESKE et al., 2004). 1.1 CAPREOMICINA A capreomicina pertence ao grupo das tuberactinomicinas (FARMER; KIM, 1998; MARRERO; CABAÑAS; MODOLELL, 1980), e é um medicamento que possui menos efeitos ototóxicos que aminoglicosídeos (STURDY et al., 2011), e é considerada a droga de escolha nos casos do isolado ser resistente simultaneamente a estreptomicina e canamicina, ou à canamicina e amicacina (WHO, 2014). Reações adversas incluem nefrotoxicidade (20-25%), ototoxicidade, dor no local da injeção e anormalidades nos eletrólitos, podendo ferar hipocalemia, hipocalcemia e hipomagnesemia (WHO, 2014). Seu mecanismo de ação consiste em inibir a síntese proteica ao interagir com ambas subunidades ribossomais, ligando-se especificamente com a hélice 69 do RNAr 23S e a hélice 44 do RNAr 16S (STANLEY et al., 2010), gerando o bloqueio da formação do complexo de iniciação na subunidade 30S (JOHNSTON, 1969), além de prejudicar o movimento de translocação do RNAt do sítio A para o sítio P (MODOLFLL; VAZQUEZ, 1977; PESKE et al., 2004). Os resíduos A1913 e C1914 do RNAr 23S formam um bolsão de interação juntamente com os resíduos A1493 e G1494, enquanto que os resíduos A1492, A1493 e A1913 parecem sofrer modificação nos seus posicionamentos por esse fármaco, o que contribui para seu efeito inibitório, uma vez que suas conformações ideais são críticas para o processo de decodificação do RNAm (STANLEY et al., 2010). Além de sua efetividade contra a MDR-TB, a capreomicina é uma das poucas drogas ativas contra o M. tuberculosis latente, que consiste na condição em que o 27 indivíduo é portador da bactéria, mas seu sistema imune consegue inibir o crescimento bacteriano, e por sua vez não desenvolve os sintomas da doença e não é capaz de transmiti-la. Entretanto, esses indivíduos podem um dia desenvolver a doença, por isso se faz necessário o tratamento ainda na forma latente (FU; SHINNICK, 2007). 1.2 ESTREPTOMICINA A estreptomicina foi o primeiro medicamento a mostrar efeito terapêutico contra tuberculose, mas o aumento da prevalência de resistência, fez com que seu uso passasse a ser reduzido (ATS; CDC; IDSA, 2003), sendo recomendada a escolha para o tratamento em casos em que há resistência à todas as drogas de segunda-linha (WHO, 2014). As reações adversas podem incluir nefrotoxicidade, ototoxicidade, toxicidade vestibular, dor no local da injeção e anormalidades eletrolíticas como hipocalemia, hipocalcemia e hipomagnesemia (WHO, 2014). Esse medicamento pertence à classe dos aminoglicosídeos, que são bactericidas que atuam prejudicando a síntese proteica ao interagir com o RNAr 16S da subunidade ribossomal menor, mais especificamente no sítio-A no interior da hélice 44 (DAVIES; GORINI; GILBERT, 1964; PESKE et al., 2004; WONG et al., 1998). Assim, esse antimicrobiano provoca erros de decodificação, resultando na incorporação de aminoácidos errados no polipeptídeo em formação, que por sua vez se dissocia rapidamente do ribossomo (RODNINA; PAPE; WINTERMEYER, 2000; RUUSALA; KURLAND, 1984). A estreptomicina está fortemente ligada ao RNAr 16S, com os resíduos 13, C526, 915 e C1490 (RAMAKRISHNAN et al., 2000), que estão relacionados às bases de proteção (MOAZED; NOLLER, 1987) e dados de mutagênese. Além disso, também interage com a proteína ribossomal S12, através do resíduo K45 (RAMAKRISHNAN et al., 2000). 1.3 MODELAGEM MOLECULAR A modelagem molecular se tornou uma área bem estabelecida durante a última década graças aos avanços em hardware e software, e tem fornecido 28 poderosas ferramentas para construção, visualização, análise e armazenamento de modelos de sistemas complexos, que podem auxiliar na interpretação das relações estrutura-atividade. Assim, é capaz de fornecer informações que guiam o design de novas estruturas (COHEN et al., 1990). Uma das aplicações da modelagem é a bioquímica computacional, que permite diversos estudos em sistemas biológicos com a utilização de softwares que processam cálculos físicos, gerando dados relacionados a representação das estruturas moleculares, bem como a simulação de seu comportamento em ambiente fisiológico. Uma importante abordagem da bioquímica computacional é o estudo da interação entre fármacos e seus receptores, tanto no desenvolvimento de novas drogas, quanto no melhoramento das existentes. Assim, as inovações permitem identificar possíveis sítios e quantificar a energia de interação entre fármacos e suas proteínas-alvo, além de determinar orientações e conformações de moléculas de interesse farmacológico (VERLI e BARREIRO, 2005), o que pode fornecer dados para descoberta de novas drogas e melhoramento de regimes de tratamento (KITANO, 2002). Essas análises podem ser realizadas através da utilização de estruturas tridimensionais obtidas por técnicas de difração de raios-X, RMN ou geradas por modelagem molecular (BLUNDELL et al., 2006; RAHA, 2007; MATTA, 2010; ZHOU, 2010; HENNIG e SATTLER, 2014; JASKOLSKI et al., 2014). O método a ser utilizado dependerá da complexidade do sistema (quantidade e tipo de átomos, tipo de ligações envolvidas, etc.), do tipo de informação que se deseja obter, além da avaliação do custo e poder computacional. 29 OBJETIVOS 30 2 OBJETIVOS 2.1 OBJETIVOS GERAIS Este trabalho pretende, como objetivo geral, apresentar uma investigação quântica, a fim de mensurar a interação dos medicamentos capreomicina e estreptomicina com o ribossomo bacteriano, identificando os resíduos componentes de seus sítios de ligação, bem com avaliar as peculiaridades que envolvem os complexos fármaco-receptor. 2.2 OBJETIVOS ESPECÍFICOS  Identificar os resíduos dos RNAs e proteínas ribossomais que são alvos dos fármacos estudados, analisando os tipos de ligações químicas que ocorrem entre eles e suas respectivas energias;  Estabelecer uma ordem de importância para os resíduos que mais contribuem para a ligação da capreomicina e da estreptomicina aos seus sítios, dimensionando as distâncias resíduo-fármaco;  Analisar as contribuições das regiões dos fármacos em estudo para a fixação dos medicamentos em seus sítios, e assim identificar as áreas que possuem grupos funcionais mais importantes para as contribuições energéticas, como também as que menos contribuem;  Evidenciar as semelhanças e diferenças entre os mecanismos de atuação dos dois fármacos, a fim de produzir dados que contribuam para estudos futuros que almejem os desenhos de novos fármacos, cada vez mais seletivos . 31 METODOLOGIA 32 3 METODOLOGIA 3.1 DADOS CRISTALOGRÁFICOS Para análise das interações fármaco-receptor, foram realizados cálculos energéticos a partir de estruturas cristalográficas disponíveis no Research Collaboratory for Structural Bioinformatics - Protein Data Bank (RCSB - PDB), que se trata de um repositório de dados de estruturas macromoleculares como proteínas e ácidos nucleicos. Esses dados estão identificados pelos códigos 4V7M (STANLEY et al., 2010) e 1FJG (RAMAKRISHNAN et al., 2000), contendo a capreomicina em complexo com as duas subunidades ribossomais 30S e 50S, e a estreptomicina com a subunidade 30S, respectivamente (Figura 8). Os cristais utilizados foram produzidos por meio da técnica de difração de raio-X, possuindo resoluções de 3,45 Å para o 4V7M e 3,0 Å para o 1FJG. A difração de raio-X é capaz de descrever as posições atômicas nas estruturas cristalinas, e assim predizer a estrutura tridimensional de macromoléculas a partir de seus cristais. As descrições produzidas nos fornecem informações sobre o conteúdo da célula unitária do cristal, e com o auxílio das transformações de Fourier, é possível distribui-los no espaço molecular real, originando a imagem da molécula cristalizada de forma tridimensional. (BÜTTNER; RENNER-SCHNECK; STEHLE, 2015). Assim, a determinação da estrutura de um cristal ou estrutura cristalina envolve a descrição da disposição no espaço de todas as entidades químicas existentes na amostra. Dentre as possibilidades do uso de dados de difração de raios-X para a reconstrução do arranjo triperiódico de unidades fundamentais, destaca-se o uso das técnicas denominadas ab initio. Classicamente este termo é utilizado dentro do mundo científico para procedimentos e análises baseadas em “primeiros princípios” físico-matemáticos. No caso da cristalografia via difração de raios-X, os dados experimentais e o modelo físico teórico, necessários para as análises, estão disponíveis. 33 Figura 8 (A) Subunidades 30S e 50S em complexo com a Capreomicina (CAP), e (B), Subunidade 30S em complexo com a Estreptomicina (SRY), com destaque em amarelo para a proteína ribossomal S12. Fonte: Elaborado pela autora 3.2 OTIMIZAÇÃO DE GEOMETRIA Os átomos de hidrogênio apresentam densidade eletrônica mínima, o que faz com que eles acabem não sendo identificados, ou ainda, posicionados 34 incorretamente nos arquivos de coordenadas cristalográficas. Portanto, foi necessário adicioná-los e em seguida proceder com a otimização de suas geometrias a fim de encontrar conformações de baixa energia. Para isso, foi utilizado o módulo Forcite do pacote de softwares Acelerys Materials Studio (DELLEY, 1990), que consiste em uma coleção de ferramentas de mecânica molecular reunindo diferentes campos de força com abordagens clássicas dentre eles o CVFF (Consistent-Valence ForceField), que foi aplicado nessa etapa uma vez que é adaptado para sistemas com biomoléculas. Com auxílio do software Marvin Sketch® (versão 5.4.1.0), identificamos o estado de protonação assumido pelos fármacos de estudo em pH fisiológico, o qual foi considerado 7,4 (NOGRADY; WEAVER, 2005). A análise identificou a presença de quatro cargas positivas na capreomicina, e três cargas positivas na estreptomicina, todas em grupamentos de aminas e amidas. Além disso, as estruturas moleculares das drogas foram divididas esquematicamente em regiões, para que uma melhor análise das interações fármaco-resíduo pudesse ser realizada. Assim, a estreptomicina passou a contar com três regiões: i, ii e iii; e a capreomicina com cinco: i, ii, iii, iv e v (Figura 9). Figura 9 Estruturas moleculares da Estreptomicina (A) e da Capreomicina (B), com a localização das cargas presentes em pH fisiológico e enumeração de suas regiões. Fonte: Elaborado pela autora 35 3.3 MFCC A acurácia dos cálculos energéticos com abordagem da mecânica quântica em moléculas biológicas como proteínas e ácidos nucleicos, constitui-se como um dos maiores desafios em química computacional e biologia. O método de Fracionamento Molecular com Capas Conjugadas (MFCC) permite análises quantitativas e eficientes para cálculos ab initio das energias de interação entre duas moléculas, como fármacos e seus receptores. A ideia básica dessa metodologia consiste em particionar a análise da interação energética entre uma molécula relativamente pequena e uma macromolécula, de modo que as somas individuais de interação que possam ser facilmente calculadas de forma ab initio (ZHANG; ZHANG, 2003), e assim serem capazes de representar a energia do sistema maior. Nesse trabalho, os RNAs ribossômicos foram decompostos em nucleotídeos e as energias de interação entre cada um deles com os fármacos em estudo foram calculadas. O mesmo foi realizado para análise da estreptomicina com a proteína S12, onde foram avaliados os aminoácidos componentes dessa macromolécula. Assim, para cada resíduo foram criados quatro subsistemas que assumiram as seguintes configurações, de acordo com a figura 10: (E I) dois resíduos antecessores (C*2C*1) + resíduo em estudo (Nn) + dois resíduos sucessores (C *1C*2) + fármaco (L); (E II) dois resíduos antecessores + dois resíduos sucessores + fármaco; (E III) dois resíduos antecessores + resíduo em estudo + dois resíduos sucessores; e (E IV) dois resíduos antecessores + dois resíduos sucessores. A energia total (L - Nn) para cada um desses subsistemas foi calculada, para então estimarmos a energia final da interação entre o nucleotídeo ou aminoácido em estudo com o fármaco, através da seguinte expressão: E (L - Nn) = E (L - C*2C*1 Nn C *1 C *2 ) – E (L - C*2C*1 C *1 C *2 ) – E (C*2C*1 Nn C *1 C *2 ) + E (C*2C*1 C *1 C *2 ) (1) Uma importante característica do MFCC é a introdução das capas moleculares conjugadas em cada quebra de ligação química, que são representadas pelos resíduos que antecedem e sucedem o resíduo em estudo. Essas capas não apenas executam a tarefa de fechar as valências abertas durante as quebras na 36 estrutura, mas também são capazes de representar a influência das áreas dos cortes nos orbitais eletrônicos das partes remanescentes (ZHANG; CHEN; ZHANG, 2003). Assim, nesse estudo foram utilizadas duas capas em cada área de corte da sequencia de aminoácidos e nucleotídeos estudados. Figura 10 Esquematização dos subsistemas formados para aplicação do MFCC. C*2C*1: Nucleotídeos antecessores ao resíduo em estudo; Nn: Nucleotídeo em estudo; C*1C*2: Nucleotídeos sucessores ao Nn; Fonte: Elaborado pela autora Após a criação de todos os subsistemas, eles foram otimizados e os cálculos das energias de interação foram realizados dentro do formalismo da Teoria do Funcional da Densidade (DFT), utilizando o módulo DMol3 (DELLEY, 1990; DELLEY, 2000). O DFT permite a realização de cálculos energéticos a partir da densidade eletrônica do sistema (HOHENBERG; KOHN, 1964), ao invés de utilizar os movimentos de cada elétron individualmente. Essa teoria possibilita em alguns casos, melhores acordos com dados experimentais disponíveis. Isso se deve ao fato 37 de que sistemas de tamanho moderado a grande possam ser estudados com uma precisão química aceitável, a um custo computacional que algumas vezes corresponde a uma fração daquele obtido utilizando-se métodos correlacionados tradicionais. O uso do DFT deve-se à abordagem de Kohn e Sham, que permite uma exata descrição da interação de sistemas com muitas partículas em termos de um sistema inerte, no qual seu potencial pode ser demonstrado como sendo completamente determinado pela densidade eletrônica do sistema que realiza interações. Uma das componentes da equação da energia total do estado fundamental de Kohn-Sham, inclui a correção da energia cinética no termo de troca e correlação, a qual não é conhecida (KOHN; SHAM, 1965). Assim, para que seja possível a utilização dessa abordagem, é necessário determinar uma boa aproximação para o termo de troca e correlação. Nesse estudo, o nível do potencial de troca-correlação utilizado foi a aproximação de densidade local (LDA), que providencia cálculos mais rápidos, em conjunto com o funcional de densidade PWC (PERDEW; WANG, 1992). Assim como no estudo desenvolvido por Mota et al. (2016), foi utilizado um raio r a fim limitar a quantidade de resíduos analisados, sem perder interações importantes. A partir desse raio r, esferas imaginárias (considerando r = n/2; n = 1, 2, 3, 4, ...) foram traçadas delimitando distâncias crescentes a partir do fármaco, a fim de se analisar as somas das energias individuais dos resíduos presentes em cada uma delas, e assim ser possível determinar a partir de qual distância dentro do sítio de ligação, o fármaco passa a ter pouca influência energética. Esse estado é observado quando a variação de energia entre os raios subsequentes é menor do que 10% (WU et al., 2007), e foi considerado como o limite para os resíduos a serem considerados nas análises de interação fármaco-receptor realizadas. Desse modo, para a capreomicina na subunidade 30S, foram avaliados resíduos presentes em um raio de até 14 Å distantes do fármaco, totalizando 44 resíduos; e na subunidade 50S, 30 nucleotídeos foram analisados, e estavam distribuídos até o raio de 30 Å de distância. Com a estreptomicina foram levados em consideração 60 nucleotídeos distribuídos até 12,5 Å de distância da droga na subunidade 30S, e 25 aminoácidos da proteína S12 com até 15 Å de distância. 38 RESULTADOS E DISCUSSÕES 39 4 RESULTADOS E DISCUSSÕES Os cálculos foram realizados a partir das estruturas cristalográficas do ribossomo em complexo com a capreomicina e com a estreptomicina. As Tabelas 1 e 2 apresentam todos os valores energéticos obtidos. Sabendo-se que a interação fármaco-receptor trata-se de um processo espontâneo, os valores energéticos são calculados em termos de energia livre de Gibbs. Considera-se que quanto mais negativa for a energia, maior é a atração entre as moléculas estudadas. A metodologia do MFCC permite o estudo das estruturas dos complexos fármaco-receptor, fragmentando-as. Assim, por meio dessa análise identificamos os resíduos que mais e os que menos contribuem energeticamente para a fixação desses antibióticos em seus sítios, como também aqueles que atuam promovendo repulsões devido às suas propriedades químicas. A capreomicina e a estreptomicina são antibióticos que possuem como mecanismo de atuação a inibição da síntese proteica bacteriana. Entretanto, seus mecanismos de ligação em seus sítios são distintos: a capreomicina é capaz de se ligar às duas subunidades ribossomais (30S e 50S), enquanto que a estreptomicina liga-se, mormente, à subunidade ribossomal menor (30S), e interage com alguns pontos da proteína S12. Assim, foram calculados os valores energéticos desses fármacos com cada resíduo componente dos seus sítios, e dentre eles, destacam-se como os principais, os possuidores das maiores energias de atração. Para a interação capreomicina-30S, foram os nucleotídeos: A1493 (-29,69 Kcal/mol); G1494 (-26,94 Kcal/mol); G1491 (-19,27 Kcal/mol); A1492 (-14,08 Kcal/mol) e A1408 (- 12,83 Kcal/mol). Enquanto que para capreomicina-50S, destacam-se como principais: A1913 (-25,26 Kcal/mol); C1914 (-9,01 Kcal/mol); U1915 (-4,58 Kcal/mol); G1964 (-4,46 Kcal/mol) e G1921 (-3,92 Kcal/mol). Na interação estreptomicina-30S, os principais são os resíduos: A914 (-38,71 Kcal/mol); A913 (-19,84 Kcal/mol); C1490 (-19,68 Kcal/mol); C912 (-16,07 Kcal/mol); G527 (-10,01 Kcal/mol). Já na interação estreptomicina-proteínaS12, há os aminoácidos PRO94 (-69,20 Kcal/mol); LEU93 (-54,92 Kcal/mol); ARG89 (-38,63 Kcal/mol); ASN49 (-24,36 Kcal/mol); ASP92 (-19,08 Kcal/mol). 40 Tabela 1 Energias de interação da Capreomicina em complexo com as subunidades 30S e 50S; raio de distância, e as regiões e grupos funcionais com os quais há as menores distâncias fármaco- resíduo. Capreomicina - 30S Capreomicina - 50S Resíduo Grupo Atômico r (Å) Energia (Kcal/mol) Resíduo Grupo Atômico r (Å) Energia (Kcal/mol) C26 i(C)O 2 -11,713 A1913 ii(NH)OH 2 -25,259 A1408 iii(NH₄⁺)NH₃ 2,5 -12,828 C1914 v(NH₂)C 2,5 -9,011 G1491 iii(NH₄⁺)NH₃ 2,5 -19,272 A1912 ii(C)OH 6 -3,800 A1492 iv(C)C 2,5 -14,755 U1915 v(NH₂)C 6,5 -4,585 A1493 ii(OH)H; iv (NH)O 2,5 -29,694 A1916 iii(O)NH₂ 9 -3,243 G1494 ii(NH₄⁺)N 2,5 -26,943 U1917 iii(O)NH₂ 12 -2,474 A36 ii(NH₄⁺)N 3 -9,669 U1911 ii(C)OH; i(NH₄⁺)O 12,5 -2,565 C1409 iii(NH)O 3,5 -5,968 A1918 ii( C)C 12,5 -2,988 A25 i(C)O 3,5 -5,951 A1919 iii(NH₂)C 13 -2,985 G1410 v(O) C 5 -4,846 C1920 ii(C)C; iii(C)C 16 -2,767 C27 v(O) C 5,5 -3,899 G1910 i(NH₂) OH 17,5 -3,639 G35 i(C)OH 5,5 -4,089 G1921 ii(C)C 20 -3,924 U38 i(NH₄⁺)C 5,5 -5,905 C1962 iii(NH₂)P 23 -2,958 C1407 iii(NH₄⁺)NH₂ 6 -1,946 C1909 i(NH₄⁺)O 23,5 -2,487 U1495 iii(NH₄⁺)C 6 -4,104 G1922 ii(C)C 24 -3,333 C24 i(C)O 6 -4,059 C1961 iii(C)O 24,5 -2,174 C1490 iii(C)NH₂ 7 -2,342 U1963 iii(NH₂)O 24,5 -3,649 C39 i(NH₄⁺)NH₂ 7 -2,432 U1946 iii(O)OH 26 -1,707 A30 i(NH₄⁺)NH₂ 8 -2,159 C1947 iii(O)C 26 -2,327 C519 iv(C)O 8,5 -4,345 C1908 iii(O)NH₂ 27 -2,126 C1411 v(O)O 8,5 -2,149 G1948 iii(O)C 27,5 -3,480 G530 iv(NH₄⁺)NH₂ 9 -3,094 G1945 v(O)NH₂ 28 -1,301 U1406 iii(NH₄⁺)O 9 -1,721 A1932 v(NH₂)O 28,5 -1,834 C1496 iii(NH₄⁺)NH₂ 9 -2,209 G1933 iii(NH₂)O 28,5 -3,329 G28 i(NH₄⁺)O 9 -3,840 G1935 iii(NH₂)O 28,5 -2,934 C40 I(NH₄⁺)NH₂ 9 -2,506 U1923 iii(O)O 28,5 -1,875 C518 iv(NH₄⁺)OH 9,5 -3,087 U1931 iii(O)C 29 -1,972 G1489 iii(C)O 10 -2,363 C1934 iii(NH₂)NH₂ 29 -2,383 G10 i(NH₄⁺)C 10 -2,292 G1964 iii(NH₂)O; v(O)O 30 -4,462 G529 iv(C)NH 10,5 -2,147 G1907 iii(O)O 30 -1,264 G1405 iii(NH₄⁺)O 10,5 -1,879 U31 i(NH₄⁺)NH 10,5 -2,615 G42 i(C)Nh 10,5 -2,416 A520 iv(C)O 11 -2,500 C1404 iii(NH₄⁺)NH₂ 11 -1,729 C11 v(NH2)H 11 -2,208 U34 i(C)OH 11 -2,591 G23 i(NH₄⁺)OH 11,5 -2,656 G1497 iii(NH₄⁺)O 12 -3,269 C1412 v(O)C 12,5 -3,873 G29 i(NH₄⁺)O 12,5 -3,475 U32 i(NH₄⁺)O 12,5 -2,637 G521 iv(C)N 13 -2,879 C1403 iii(NH₄⁺)NH₂ 14 -2,545 Fonte: Elaborado pela autora 41 Tabela 2 Energias de interação da Estreptomicina com a subunidade 30S e proteína S12; raio de distância, e as regiões e grupos funcionais com os quais há as menores distâncias fármaco-resíduo. Estreptomicina - 30S Estreptomicina - S12 Resíduo Grupo Atômico r (Å) Energia (Kcal/mol) Resíduo Grupo Atômico r (Å) Energia (Kcal/mol) U14 iii(O)O 2 -4,250 LYS47 i(OH)C 2 1,804 G527 iii(N)O 2 -10,677 LYS46 i(OH)NH₄⁺ 2,5 4,130 C912 i(O)O 2 -16,071 PRO48 ii(C)C 2,5 -12,581 A913 i(NH₂)O 2 -19,844 LYS91 ii(O)C 2,5 8,627 A914 i(NH₂)O; iii(N)O 2 -38,717 ASP92 ii(C)OH 5,5 -19,089 C1490 i(NH₂)O 2 -19,687 PRO94 i(NH2)C 5,5 -69,204 A915 i(NH₂)P 2,5 -9,009 ASN49 ii(C)NH₃ 6 -24,364 G1491 i(NH₂)O 2,5 -10,018 ARG89 i(C)NH₃; ii(O)NH₄⁺ 7 -38,631 U12 iii(N)O 3 -2,769 PRO45 ii(C)C 8 -15,540 U13 iii(O)O 3 -7,541 LEU93 ii(C)NH₃ 8,5 -54,921 C526 iii(N)O 3 -2,963 ARG53 ii(C)NH₃ 9 1,787 U911 i(NH)O 4 -2,826 VAL90 ii(C)O 9,5 -0,011 G886 i(NH₂)H 4,5 -2,526 GLY95 i(OH)NH₃ 9,5 0,010 G1489 i(NH₂)H 4,5 -2,678 SER50 ii(C)NH₃ 10 -0,070 G15 iii(O)O 5 -3,291 THR44 ii(C)O 11 -0,037 C525 ii(O)H 5,5 -2,441 VAL96 i(OH)NH₃ 11 -0,062 A523 ii(O)N 6 -1,717 ALA51 ii(C)NH₃ 12,5 -0,071 G885 i(NH₂)H 6,5 -1,834 ARG97 i(OH)C 13 1,711 G887 i(NH₂)H 6,5 -3,002 VAL43 i(OH)C; ii(C)C 13 0,024 G916 i(NH2)P 6,5 -3,967 LYS23 i(C)NH₄⁺; ii(O)NH₄⁺ 13,5 1,825 A1492 i(H)O 6,5 -3,952 HIS99 ii(C)C 14 -0,060 A16 iii(O)O 7 -2,614 TYR69 ii(C)C 14,5 -0,049 G22 ii(H)N 7 -1,865 LEU52 ii(C)NH₃ 15 0,087 C528 ii(O)O 7 -2,471 THR67 i(OH)C 15 -0,035 G11 iii(C)H 8 -1,928 GLY88 ii(O)C 15 0,058 A1396 iii(O)O 8 -2,670 C1403 iii(O)P 8 -3,138 G21 iii(C)N 8,5 -1,523 C522 ii(C)N 8,5 -1,572 C1397 iii(O)C 8,5 -2,518 C1395 i(NH₂)O 9 -2,614 G1410 i(C)N 9 -1,360 C23 ii(C)O 10 -2,398 G505 iii(C)NH₂ 10 -1,530 G524 ii(C)O 10 -1,864 G529 ii(C)O 10 -2,147 G888 i(NH₂)C 10 -2,440 C910 i(C)H 10 -1,987 C1402 iii(OH)C 10 -2,420 C1404 i(N)P 10 -2,913 G1488 iii(N)P 10 -1,620 U24 ii(O)C 10,5 -2,068 G506 ii(O)NH₂ 10,5 -1,530 C518 ii(O)NH₂ 10,5 -1,892 A535 ii(O)C 10,5 -2,045 U17 iii(OH)O 11 -2,169 G521 ii(C)O 11 -1,490 A1394 i(NH2)O 11 -2,086 C18 iii(OH)O 11,5 -1,685 C19 iii(OH)NH₂ 11,5 -1,635 G1401 iii(OH)O 11,5 -1,753 G1405 iii(OH)P 11,5 -2,600 C1411 i(N)NH₂ 11,5 -1,736 C1412 i(N)C 11,5 -1,733 A1493 i(N)C 11,5 -2,429 U20 iii(C) N 12 -1,751 42 A909 i(N)C 12 -1,778 A1408 i(N)NH₂ 12 -1,458 C1409 i(C)NH₂ 12 -1,481 G530 iii(OH)O 12,5 -1,842 Fonte: Elaborado pela autora Além disso, foram identificados resíduos dentro do sítio de ligação com os menores valores atrativos. No complexo capreomicina-30S, foram: U1406 (-1,72 Kcal/mol); C1404 (-1,73 Kcal/mol); G1405 (-1,88 Kcal/mol); C1407 (-1,94 Kcal/mol); G529 (-2,14 Kcal/mol). Enquanto que para a capreomicina-50S: G1907 (-1,26 Kcal/mol); G1945 (-1,30 Kcal/mol); U1946 (-1,71 Kcal/mol); A1932 (-1,83 Kcal/mol); U1923 (-1,87 Kcal/mol). Entre as menores interações da estreptomicina com a subunidade 30S, há os resíduos: G1410 (-1,36 Kcal/mol); A1408 (-1,45 Kcal/mol); C1409 (-1,48 Kcal/mol); G521 (-1,49 Kcal/mol); G21 (-1,52 Kcal/mol). Com a proteína S12, houve aminoácidos com forças repulsivas: LYS91 (8,62 Kcal/mol); LYS46 (4,13 Kcal/mol); LYS23 (1,82 Kcal/mol); LYS47 (1,80 Kcal/mol) e ARG53 (1,78 Kcal/mol). Nota-se que os nucleotídeos com menores valores energéticos possuem como característica comum estarem a longas distâncias das drogas, de modo que os listados para a capreomicina com a subunidade menor, encontram-se entre 6 a 11 Å distantes da droga, enquanto que com a subunidade maior, estão entre 26 e 30 Å. Para a estreptomicina, estão entre 8,5 e 12 Å. Esse padrão não ocorre nas interações proteicas, uma vez que as localizações dos aminoácidos repulsivos mostraram-se diversas, ocorrendo desde os raios mais próximos a estreptomicina (2 Å de distância), quanto os mais afastados. Assim, na análise dos nucleotídeos, como todos eles possuem as mesmas propriedades químicas, a distância atua como um fator de grande importância para os mecanismos de influência da droga. O que revela o quão importante é o design molecular dos fármacos que possuem como alvo os RNAs bacterianos, vez que seus grupos funcionais deverão atuar ativamente na busca por pontos de interações com os nucleotídeos, e serem capazes de formarem sítios de ligação com longos alcances, ou seja, com grandes raios de distância. Enquanto que nas relações com aminoácidos, além da estrutura da droga, os diferentes aminoácidos presentes no sítio devem ser levados em consideração, vez que na natureza há cerca de 20 que são essenciais aos organismos vivos, cada um com suas particularidades atômicas (há neutros, carregados positiva e 43 negativamente), o que gera uma maior diversidade nas interações químicas. Desse modo, a construção de um fármaco que possua proteínas como alvo, pode requerer mais atenção no que concerne ao posicionamento estratégico de grupos funcionais, e assim ser capaz de interagir de maneira efetiva. Os gráficos da figura 11 trazem os valores energéticos encontrados para alguns dos resíduos considerados no sítio de ligação dos fármacos em estudo, bem como as distâncias em angstroms (Å) entre cada resíduo-fármaco, relacionadas com as regiões das drogas e os grupos funcionais mais próximos. Figura 11 Gráficos com os valores energéticos dos 20 resíduos mais importantes para capreomicina (30S e 50S) e estreptomicina (30S), e os aminoácidos da proteína S12 com os maiores e menores resultados. As distâncias de cada resíduo em relação à droga são apresentadas em angstrons (Å), com o detalhamento dos grupos funcionais dos resíduos e dos ligantes mais próximos. Fonte: Elaborado pela autora 44 A interação de um fármaco com seu sítio de ação no sistema biológico é determinada por forças intermoleculares, dentre as quais se podem listar as interações hidrofóbicas, eletrostáticas e estéricas (PRATT; TAYLOR, 1990). As forças eletrostáticas são aquelas resultantes da interação entre dipolos e/ou íons de cargas opostas, cuja magnitude depende diretamente da constante dielétrica do meio e da distância entre as cargas (BARREIRO; FRAGA, 2008). Esse tipo de interação é uma das mais presentes nos sistemas em estudo, pelo fato do ribossomo ser uma molécula com grande potencial eletrostático, uma vez que possui inúmeras cargas negativas devido às ligações fosfodiéster que conectam os nucleotídeos formadores da cadeia de RNA. Além disso, a proteína S12 também possui alguns aminoácidos que se apresentam ionizados em pH fisiológico, e no sítio estudado, houve incidência de argininas e lisinas que possuem cargas positivas. A capreomicina possui quatro grupos NH4 + presentes nas regiões i, iii e iv, totalizando quatro cargas positivas, o que confere um grande potencial para interação com as cargas negativas presentes nos nucleotídeos que formam o seu sítio receptor. Na superfície de interação com a subunidade 30S, essas três regiões estão envolvidas energeticamente com alguns dos resíduos definidos como principais no sítio de ligação, enquanto que a interação com a subunidade 50S, apenas a região iv não está realizando interação com alguns dos importantes resíduos. As áreas carregadas na estreptomicina são i e iii, também em grupamentos NH4+. Analisando os pontos de contato com o RNA 16S, é possível notar que essas cargas trazem para essas regiões um maior nível de interações, vez que há um maior número de resíduos mais próximos a elas quando comparadas com a que não possui tamanha polaridade. A região ii mostra-se mais ativa na análise do complexo com a proteína S12, por se aproximar a um maior número de resíduos, ainda que dentre eles incluem-se aminoácidos carregados positivamente. Mas a presença de carbonilas e hidroxilas nessa área podem estar evitando a origem de forças repulsivas, que acabam ocorrendo na região i, por concentrar duas cargas positivas e ao seu redor haver aminoácidos com cargas também positivas. Moléculas nas quais há separação parcial de cargas entre átomos ou grupos funcionais adjacentes podem formar interações dipolo-dipolo. Ambos fármacos 45 objetos desse estudo possuem em sua estrutura grupos amônios (NH4 +) que favorecem interações iônicas, vários grupos amina (NH3) e carbonilas (C=O) envolvidos em interações dipolo-dipolo (NOGRADY; WEAVER, 2005). Os grupos funcionais carbonil (C=O) formam um dipolo, devido à eletropositividade do carbono e eletronegatividade do oxigênio. Eles estão distribuídos em todas as regiões da capreomicina, contribuindo igualmente para todas as interações, fazendo isso tanto com as porções dos grupos fosfatos dos nucleotídeos (grupos negativos), quanto com os grupos aminos normalmente presentes também nos nucleotídeos, que também são grupos com potencial para originar dipolos, enquanto que na estreptomicina está presente apenas na região ii. As ligações de hidrogênio são as mais importantes interações não-covalentes existentes nos sistemas biológicos, essas interações são formadas entre heteroátomos eletronegativos (como oxigênio, nitrogênio ou flúor), e hidrogênio (BARREIRO; FRAGA, 2008). A figura 12 mostra que capreomicina realiza esse tipo de interação com os resíduos U37, A1408, A1493 e G1494 pertencentes à subunidade 30S, em parcial concordância com outros estudos que relatam esse tipo de interação apenas com os resíduos A1408, A1493 e incluem ainda G1491 (AKBERGENOV et al., 2011). Além disso, há ligação de hidrogênio também com os resíduos A1914 e A1913 da subunidade 50S. A estreptomicina realiza esse tipo de interação com os aminoácidos LYS91 e LYS46, pertencentes à proteína S12; e com a subunidade 30S, faz ligações de hidrogênio com os resíduos G1491, C1490, G527, A914 e U14 (Figura 13). 46 Figura 12 Resíduos capazes de interagir através de ligações de hidrogênio nas subunidades 30S e 50S com a capreomicina (Cap). Fonte: Elaborado pela autora 47 Figura 13 Resíduos capazes de interagir através de ligações de hidrogênio nas subunidades 30S e proteína S12 com a estreptomicina (SRY). Fonte: Elaborado pela autora Forças mais fracas como as interações de van der Walls, também podem estar ocorrendo nesses sistemas, caracterizam-se pela aproximação de moléculas apolares com múltiplos dipolos induzidos. Mas mesmo sendo interações de baixa energia, são importantes no processo de reconhecimento molecular do sítio pelo fármaco, e como costumam se revelar em vários pontos, o somatório energético acaba originando contribuições significativas (BARREIRO; FRAGA, 2008). A capreomicina e a estreptomicina caracterizam-se por serem fármacos estruturalmente específicos, uma vez que exercem seus efeitos por meio de 48 interações seletivas com seus alvos ribossomais. Assim, o reconhecimento molecular da droga é dependente do arranjo espacial dos grupamentos funcionais e das propriedades estruturais, que devem ser complementares ao sítio de ligação. Foram destacadas as regiões dos fármacos com as quais os 20 resíduos mais importantes estavam mais próximos. Na capreomicina interagindo com a subunidade 30S (Figura 14), onde a região i mostrou-se envolvida com dez nucleotídeos, diferentemente das regiões ii e v, com ligações com dois, e as iii e iv, com quatro, cada uma. As energias totais de interação por região foram de: -54,76 Kcal/mol, para a região i, -56,64 Kcal/mol para a região ii, -22,90 Kcal/mol para a região iii, -38,37 Kcal/mol para a região iv e -8,72 Kcal/mol para a v. Na subunidade 50S, apenas a região iv da capreomicina não apresentou resíduos mais próximos dentre os 20 principais. Assim, na região i houve ocorrência de três nucleotídeos e energia total de -8,70 Kcal/mol; na região ii, seis resíduos e - 42,07 Kcal/mol; região iii, nove resíduos e -29,51 Kcal/mol e na região v: dois resíduos e -13,60 Kcal/mol. Com a estreptomicina, as regiões e os respectivos resíduos mais próximos também foram destacados (Figura 15). Na subunidade 30S, também foram esquematizados os 20 com maiores valores de forças atrativas: para a região i, houve ocorrência de 10 nucleotídeos, com um somatório energético de -90 Kcal/mol; a região ii, mostrou-se mais próxima a apenas um resíduo, sendo este possuidor de uma energia de -3,97 Kcal/mol; e a região iii, contou com oito resíduos e a energia total foi de -37,30 Kcal/mol. Em seus pontos de menores distâncias com a proteína S12, analisando os 13 resíduos com maiores e menores forças atrativas e repulsivas, apenas a região iii não se mostrou com potencial de interação com os aminoácidos considerados. As regiões i e ii atuaram de forma mais ativa, sendo a primeira com cinco resíduos, e a segunda com oito. As energias totais foram -59,80 Kcal/mol e -144,25 Kcal/mol, respectivamente. Realizando um somatório das energias de todos os nucleotídeos considerados na capreomicina para subunidade 30S, foi obtida -275,18 Kcal/mol de energia, e para a 50S, -112,83 Kcal/mol. Com a estreptomicina, a energia total do sistema desse fármaco em complexo com a subunidade menor foi de -223,24 Kcal/mol, e com a proteína S12, -214,66 Kcal/mol. 49 Figura 14 Regiões da capreomicina (Cap) com interações na subunidade 30S e 50S: (i) Região I, (ii) Região II, (iii) Região III, (iv) Região IV e (v) Região V. Fonte: Elaborado pela autora 50 Figura 15 Regiões da estreptomicina (SRY) com interações na subunidade 30S e proteína S12: (i) Região I, (ii) Região II e (iii) Região III. Fonte: Elaborado pela autora 51 Os sítios de ligação dos fármacos estudados possuem uma grande semelhança pelo fato de compartilharem uma localização muito próxima dentro da subunidade 30S, visto que o RNA ribossômico 16S é comum aos dois fármacos. Assim, foi possível identificar dezessete resíduos presentes em ambos sítios de ligação dessas drogas, evidenciados na figura 16, onde há comparação suas respectivas distribuições e energias. Observa-se que em onze dos dezessete resíduos, a capreomicina apresentou valores atrativos maiores que as interações da estreptomicina, e o contrário ocorreu com os resíduos G1405, C1404, C1403, G1489 e C1490. Esses dados não são suficientes para esclarecerem se esse evento é uma das justificativas comprobatórias do fato da capreomicina ser considerado um medicamento melhor que a estreptomicina, dado que esse primeiro é mais cotado para uso terapêutico (OMS, 2011). Figura 16 Resíduos comuns entre a capreomicina (CAP) e estreptomicina (SRY), e suas respectivas energias. Fonte: Elaborado pela autora Estudos da seletividade medicamentosa mostrou que os resíduos A1408 e G1491 auxiliam na especificidade do sítio para a capreomicina no ribossomo 52 bacteriano (AKBERGENOV et al., 2011), distinguindo-o assim do ribossomo de células eucarióticas. Esse estudo revelou que esses nucleotídeos configuram-se como uns dos mais importantes do ponto de vista energético, justificando assim a grande relevância desses nucleotídeos. A família das tuberactinomicinas são caracterizadas por possuírem em comum um core pentapeptídico, composto por L-serine e resíduos não- proteinogênicos 2,3-diaminoproprionato, L-capreomycidine e β-ureidodehydroalanine (STANLEY et al., 2010), e diferenciam-se pelas modificações em seus radicais, como a presença ou ausência da cadeia β-lysine. Análises cromatográficas mostram que a capreomicina consiste em quatro componentes microbiologicamente ativos (Figura 17), que são as formas IA, IB, IIA e IIB (HERR; REDSTONE, 1966). A cadeia L-capreomicidina possui importante papel no efeito inibitório dessas drogas (NOMOTO; SHIBA, 1977), e pode estar envolvida na restrição estérica da orientação de A1492 e A1493 (STANLEY et al., 2010). Figura 17 Estrutura química da família de antibióticos tuberactimicinas (A) e estrutura da Capreomicina (B). Fonte: Adaptado de WANK, H., ROGERS, J., DAVIES, J, SCHROEDER, R. Peptide Antibiotics of the Tuberactinomycin Family as Inhibitors of Group I Intron RNA Splicing. J. Mol. Biol, 236, 1001- 1010, 1994. A B 53 A estreptomicina é antibiótico pertencente ao grupo dos aminoglicosídeos, que são estruturados a partir de um aminociclitol conectados através de ligações glicosídicas a diferentes açúcares aminados (SPELMAN; MCDONALD; SPICER, 1989). Além disso, as estruturas desses medicamentos contêm muitos grupos amino e hidroxil, que conferem alta polaridade e basicidade, como também cargas positivas em pH fisiológico, que contribuem para a atividade antimicrobiana (JANA; DEB, 2006). Dessa forma, a estreptomicina é um pseudotrissacarídeo, composto por um aminociclitol do tipo estreptidina, aderido a um dissacarídeo (MAGNET; BLANCHARD, 2005). Análises estruturais revelam que o sítio de ligação desses fármacos compreende as bases A1408, A1492 e A1493, que são deslocalizadas durante a fixação da droga, e consequentemente promove alterações conformacionais na hélice 44 do RNA 16S (SILVA; CARVALHO, 2007). A resistência à capreomicina tem sido observada quando há inativação do gene tlyA methyltransferase, pois a enzima TlyA methyltransferase é responsável por metilar os resíduos C1920 e C1409, e a ausência desse processo reduz a susceptibilidade a esse medicamento (JOHANSEN et al., 2006; MAUS; PLIKAYTIS; SHINNICK, 2005). A interação da capreomicina com os resíduos C1920 e C1409, foi de -4,43 Kcal/mol e -5,97 Kcal/mol, respectivamente. Outra mutação associada à resistência a esse medicamento é a que afeta o gene rrs, que quando modificado gera polimorfismos, com as substituições A1401G, C1402T e G1484T, impedindo a fixação da capreomicina em seu sítio (CAMPBELL et al., 2011). Mutações nos resíduos C1409 e G1491 afetam significativamente a susceptibilidade do M. tuberculosis à capreomicina (AKBERGENOV et al., 2011). Trata-se de resíduos que possuem grande importância no sítio de ligação, uma vez que o G1491 possui a terceira maior energia de interação (-19,27 Kcal/mol), e o C1409, a oitava (-5,97 Kcal/mol). Com isso, apesar da capreomicina realizar fortes interações com resíduos do RNAr 23S, polimorfismos que afetem apenas o RNAr 16S são suficientes para reduzir a seletividade da droga (AKBERGENOV et al., 2011). Isso pode ser justificado pelo fato da subunidade menor ter maior representatividade dentre os nucleotídeos que interagem com essa droga, quando comparada à subunidade maior. Assim, a maior quantidade de pontos de ligação da capreomicina com o RNAr 16S pode ser a principal responsável pela especificidade da droga em estudo. 54 As mutações associadas com resistência à estreptomicina envolvem genes que codificam componentes do RNAr 16S e da proteína S12. No RNA 16S, as mutações relatadas são associadas a alterações no gene rrs, e pode-se citar a que inclui o resíduo A866, que pode sofrer uma troca A-G (DOUGLASS; STEYN, 1993). Em nossa análise, é um resíduo que se encontra a 4,5 Å de distância da droga, com uma energia de -2,52 Kcal/mol. Outro gene que comumente sofre mutações gerando resistência a esse fármaco, é o rpsL, que codifica a proteína ribossomal S12, fazendo trocas nos aminoácidos das posições 43 ou 88 (FINKEN et al., 1993; SREEVATSAN et al., 1996). Encontramos que esses aminoácidos correspondem aos resíduos VAL43 e GLY88, que possuem interações com energias de 0,02 Kcal/mol e 0,06 Kcal/mol respectivamente. As trocas ocorridas nessa mutação provavelmente resultam na inserção de resíduos capazes de produzir repulsão da droga em seu sítio, uma vez que os aminoácidos da proteína selvagem não possuem influencia de forma atrativa. As simulações computacionais e análises estruturais confirmaram que o sítio de ligação da capreomicina situa-se entre as duas subunidades ribossomais, e a estreptomicina, à subunidade menor e com aminoácidos da proteína S12. Os nucleotídeos A1493 e C1494, pertencentes ao RNA ribossômico 16S, possuem as maiores energias de interação com capreomicina (-27,23 Kcal/mol e - 21,95 Kcal/mol, respectivamente), e a literatura os apresenta como sendo os principais responsáveis pela manutenção do sítio de ligação desse fármaco, juntamente com os resíduos A1913 (-30,41 Kcal/mol) e A1914 (-21,46 Kcal/mol), presentes no RNA 23S da subunidade 50S. A grande proximidade desses nucleotídeos (aproximadamente 2,5 Å) favorece a formação de ligações de hidrogênio. Esses dados confirmam a literatura (STANLEY et al., 2010), e mutações nesses pontos podem conferir resistência a essa droga. Algumas cepas de Escherichia coli produzem a citotoxina Colicina E3, que uma vez dentro da célula, um de seus domínios (E3 rRNAse) cliva a ligação fosfodiéster entre A1493 e G1494, prejudicando a síntese proteica (SENIOR; HOLLAND, 1971), e a interação de tuberactinomicinas nesse sítio, uma vez que esses nucleotídeos são essenciais para suas ligações, e essas alterações podem levar à resistência (LANCASTER et al., 2008). 55 A restrição estérica na orientação de A1492 e A1493 contribui para o efeito inibitório da capreomicina no processo de síntese proteica, pois na presença desse fármaco, eles não podem retornar à posição que ocupam no apo ribossomo (STANLEY et al., 2010). O A1492 é o resíduo que possui a quarta maior energia de interação com a capreomicina (-12,33 Kcal/mol) na subunidade menor, por possuir uma grande proximidade e estabelecer ligações de hidrogênio com a droga. As posições de A1913 (RNAr, 23S) e A1492 e A1493 (RNAr, 16S) são críticas na decodificação do RNAm, e esses posicionamentos são alterados pela capreomicina (STANLEY et al., 2010), que interage fortemente com todos eles. Os resultados aqui apresentados corroboram com a literatura vigente, no que se refere a consagrados mecanismos de atuação desses dois medicamentos, como também trazemos outras perspectivas e novos dados, revelando nucleotídeos e outros aminoácidos capazes de realizar fortes interações não previstas na literatura, mostrando que essa análise pode auxiliar futuros estudos que objetivem o desenvolvimento de novos medicamentos para a tuberculose, o que se mostra urgente dado os crescentes casos de multirresistência. 56 CONCLUSÕES 57 5 CONCLUSÕES As técnicas computacionais possuem grande relevância para a pesquisa bioquímica e farmacêutica. Os perfis de interação energética da ligação da estreptomicina e da capreomicina ao ribossomo bacteriano valeram-se das estruturas cristalográficas desses complexos, juntamente com os cálculos em um nível quântico, que foram fundamentais para o entendimento da interação dessas drogas que são capazes de eliminar cepas de M. tuberculosis que desenvolvem resistência aos medicamentos convencionais. As difrações de raios-X utilizadas permitiram a identificação dos resíduos que mais contribuem para fixação desses fármacos em seus sítios de ligação, como também o que menos contribuem. Assim, por meio de interações eletrostáticas, ligações de hidrogênio, e interações de van der Walls, por exemplo, esses ligantes são capazes de produzir seus efeitos inibitórios, e os posicionamentos dos grupos funcionais como hidroxilas, aminas, amônios e carboxilas, amplamente presentes nas duas drogas, são cruciais para essas ligações, já que permitem a complementariedade necessária para os fármacos identificarem seus sítios, bem como possuem um grande potencial energético, devido às suas propriedades atômicas. Os resultados obtidos confirmaram a literatura mostrando as ligações de hidrogênio entre a capreomicina e os resíduos A1408 e A1493, e revelou ainda a presença desse tipo de interação com o U37 e G1494, como também os não previstos na literatura com a estreptomicina, G1491, C1490, G527, A914 e U14. Estudos trazem ainda informações de que os nucleotídeos A1408 e G1491 auxiliam na especificidade da capreomicina no sítio de ligação, o que pode ser resultado da grande proximidade, pois nas análises realizadas nesse trabalho, foi identificado que estão à 2,5 Å de distância da capreomicina e configuram-se entre os cinco resíduos com as maiores energias. Pesquisadores que estudaram os mecanismos de resistência a esses medicamentos, apontam como principais as que afetam os resíduos C1920, C1409, A1401, A1402, G1491 na ligação com a capreomicina, e todos eles configuraram em nossas avaliações com altos valores energéticos, e mudanças provocadas em suas estruturas, são capazes de reduzir a sensibilidade da bactéria ao medicamento, assim como no A886 com a estreptomicina. 58 A estreptomicina consegue englobar em seu sítio um maior número de nucleotídeos em um menor raio de distância, quando comparada com a capreomicina, e nas interações com a proteína S12, identificou-se aminoácidos com altos valores energéticos, não relatados pelos estudos experimentais, como os dos resíduos LEU93, ARG89, PRO94 e ASP92. Ainda segundo a literatura, os resíduos A1492 e A1493 sofrem mudanças conformacionais pelos dois fármacos estudados, e inclui ainda o A1408 afetado pela estreptomicina. Os resultados obtidos confirmaram a presença deles no sítio, como também suas energias que foram significativas. Inclusive foi demonstrado que capreomicina realiza ligação de hidrogênio com o A1493. O que nos conduz a suposição de que uma possível ineficiência dessa alteração estrutural pelas drogas podem também caracterizar um mecanismo de resistência. A divisão esquemática das estruturas moleculares dos fármacos em estudo possibilitou uma visão geral sobre os pontos das drogas mais envolvidos com os resíduos que mais contribuem para suas fixações em seus sítios. Assim, na análise da energia total dos resíduos mais próximos, na ligação da capreomicina à subunidade 30S, a região ii apresentou-se como a mais importante (-56,64 Kcal/mol), mas com valor próximo ao da região i (-54,76 Kcal/mol), enquanto que com a subunidade 50S, a região ii também se configurou como a mais importante, apresentando um total de -42,07 Kcal/mol. A estreptomicina, que foi dividida em três porções, revelou como principal a região i, com os nucleotídeos da subunidade ribossomal, e a região ii para as interações com a proteína S12. Por fim, apesar de serem medicamentos com mecanismos de ações semelhantes e com sítios próximos, as estruturas moleculares influenciam de forma substancial na produção dos seus efeitos. O que pode ser notado por meio da quantidade de resíduos atingidos por cada um e na intensidade energética dos nucleotídeos que se apresentaram como sendo comuns aos sítios dos dois. 59 PERSPECTIVAS 60 6 PERSPECTIVAS Este trabalho trouxe a caracterização energética e os componentes do sítio de interação de importantes medicamentos utilizados no tratamento da tuberculose multirresistente, e os dados obtidos são capazes de guiar uma melhor compreensão dos mecanismos de ligação desses fármacos, como também estudos para o desenvolvimento de novos fármacos a partir do conhecimento sobre o comportamento dos nucleotídeos do RNA bacteriano, e da influencia do design do fármaco na complementariedade do sítio de ligação e na extensão da influência energética. Dessa forma, as regiões apontadas como menos influentes nas interações podem ser alvo de remodelagens, para que tenhamos medicamentos cada vez mais eficazes. Além disso, levando-se em consideração as modificações conformacionais promovidas por ambos os medicamentos, posteriores análises estruturais serão necessárias para determinar como essas alterações ocorrem, bem como realizar uma comparação entre a capreomicina e a estreptomicina, visto que modificam os mesmos resíduos. Uma melhor análise do sítio da estreptomicina se faz necessária, dado que há poucos dados na literatura que especifiquem os nucleotídeos participantes do seu sítio de ligação, bem como da real influência da proteína S12, pois diversos aminoácidos não previstos na literatura apresentaram altos valores energéticos, e podem ter grande papel nos mecanismos de resistência facilmente adquiridos pelas bactérias a esse medicamento. Levando-se em consideração que a estreptomicina acaba sendo a última saída nos casos de multirresistência, esses dados auxiliam o desenvolvimento de novos fármacos que tomem sua estrutura molecular como base. 61 REFERÊNCIAS 62 7 REFERÊNCIAS AHMAD, Suhail; MOKADDAS, Eiman. Current status and future trends in the diagnosis and treatment of drug-susceptible and multidrug-resistant tuberculosis. Journal Of Infection And Public Health, v. 7, n. 2, p.75-91, 2014. ALBERTS, Bruce et al. Molecular biology of the cell. 5. ed. Newyork: Garland Science, 2008. American Thoracic Society (ATS); Centers For Disease Control And Prevention (CDC); Infectious Diseases Society Of America (IDSA). Treatment of tuberculosis. American Journal Of Respiratory And Critical Care Medicine, v.167, n. 4, p.603-662, 2003. American Thoracic Society (ATS), Centers For Disease Control And Prevention (CDC), Infectious Disease Society of America (IDSA). Diagnostic standards and classification of tuberculosis in adults and children. Am J Respir Crit Care Med, v. 161, p.1376–95, 2000. AKBERGENOV, R. et al. Molecular basis for the selectivity of antituberculosis compounds capreomycin and viomycin. Antimicrobial agents and chemotherapy. v. 55, n. 10, p.4712-4717, 2011. BARREIRO, E. J; F. Carlos A. M. Química medicinal: As bases moleculares da ação dos fármacos. Artmed. 2ª ed., 2008. BÜTTNER, Felix Michael; RENNER-SCHNECK, Michaela; STEHLE, Thilo. X- ray crystallography and its impact on understanding bacterial cell wall remodeling processes. International Journal Of Medical Microbiology, v. 305, n. 2, p.209-216, 2015. CAMINERO, J.A. Treatment of multidrug-resistant tuberculosis: evidence and controversies. The International Journal of Tuberculosis and Lung Disease, v.10, n. 8, p.829-837, 2006. CAMMARATA, Arthur. The Molecular Basis of Antibiotic Action. Journal Of Pharmaceutical Sciences, v. 62, n. 9, p.1577-1578, 1973. 63 CAMPBELL, P. J. et al. Molecular detection of mutations associated with first- and second-line drug resistance compared with conventional drug susceptibility testing of Mycobacterium tuberculosis. Antimicrobial agents and chemotherapy. v. 55, n. 5, p.2032-2041, 2011. CHAN, Edward D. et al. Treatment and Outcome Analysis of 205 Patients with Multidrug-resistant Tuberculosis. American Journal Of Respiratory And Critical Care Medicine, v. 169, n. 10, p.1103-1109, 2004. CHAPARAS, Sotiros D.; GOOD, Robert C.. The Immunology of Mycobacterial Infections. Crc Critical Reviews In Microbiology, v. 9, n. 2, p.139- 197, 1982. COHEN, N. Claude et al. Molecular modeling software and methods for medicinal chemistry. J. Med. Chem., v. 33, n. 3, p.883-894, 1990. COLDITZ, Graham A. Efficacy of BCG Vaccine in the Prevention of Tuberculosis. Jama, v. 271, n. 9, p.698-702, 1994. CRUZ-KNIGHT, Wanda; BLAKE-GUMBS, Lyla. Tuberculosis. Primary Care: Clinics in Office Practice, v. 40, n. 3, p.743-756, 2013. DANNENBERG, A. M.. Immune Mechanisms in the Pathogenesis of Pulmonary Tuberculosis. Clinical Infectious Diseases, v. 11, n. 2, p.369-378, 1989. DAVIES, Julian; GORINI, Luigi; GILBERT, Walter. Streptomycin, suppression, and the code. Proceedings Of The National Academy Of Sciences Of The United States Of America, v. 5, n. 51, p.883-890, 1964. DELLEY, B.. An all-electron numerical method for solving the local density functional for polyatomic molecules. The Journal Of Chemical Physics, v. 92, n. 1, p.508-517, 1990. DELLEY, B.. From molecules to solids with the DMol [sup 3] approach. The Journal Of Chemical Physics, v. 113, n. 18, p.7756-7764, 2000. DOUGLASS, J.; STEYN, L. M. A Ribosomal Gene Mutation in Streptomycin-Resistant Mycobacterium tuberculosis Isolates. Journal Of Infectious Diseases. v. 167, n. 6, p.1505-1506, 1993. 64 DYE, Christopher. Global Burden of TuberculosisEstimated Incidence, Prevalence, and Mortality by Country. Jama, v. 282, n. 7, p.677-686, 1999. EDWARDS, Dearborn; KIRKPATRICK, Charles H.. The Immunology of Mycobacterial Diseases. Am Rev Respir Dis, v. 134, n. 5, p.1062-1071, 1986. EKER, Barbara et al. Multidrug and Extensively Drug-Resistant Tuberculosis, Germany. Emerging Infectious Diseases, v. 14, n. 11, p.1700-1706, 2008. FARMER, Paul; KIM, Jim Yong. Community based approaches to the control of multidrug resistant tuberculosis: introducing “DOTS-plus”. Bmj : British Medical Journal, v. 317, n. 7159, p.671-674, 1998. FU, Li M.; SHINNICK, Thomas M.. Genome-wide exploration of the drug action of capreomycin on Mycobacterium tuberculosis using Affymetrix oligonucleotide GeneChips. Journal Of Infection, v. 54, n. 3, p.277-284, 2007. HERR, E. B.; REDSTONE, M. O. Chemical and physical characterization of capreomycin. Annals of the new york academy of sciences. v. 135, n. 2, p.940- 946, 1966. HOHENBERG, P.; KOHN, W. The Inhomogeneous Electron Gas. Phys. Rev, v.136, n.3B, p.864-871, 1964. JANA, S.; DEB, J. K.. Molecular understanding of aminoglycoside action and resistance. Applied Microbiology And Biotechnology, v. 70, n. 2, p.140-150, 2006. JOHANSEN, Shanna K. et al. Capreomycin binds across the ribosomal subunit interface using tlyA-Encoded 2′-O-Methylations in 16S and 23S rRNAs. Molecular Cell. v. 23, n. 2, p.173-182, 2006. JOHNSTON, Robert F.. Chemotherapy of Pulmonary Tuberculosis. Annals Of Internal Medicine, v. 70, n. 2, p.359-367, 1969. KARAKOUSIS, Petros C.; BISHAI, William R.; DORMAN, Susan E.. Mycobacterium tuberculosiscell envelope lipids and the host immune response. Cellular Microbiology, v. 6, n. 2, p.105-116, 2004. 65 KIRST, H.; ALLEN, N. Aminoglycoside Antibiotics. Reference Module In Chemistry, Molecular Sciences And Chemical Engineering, p.1-23, 2013. KITANO, Hiroaki. Computational systems biology. Nature, v. 420, n. 6912, p.206-210, 2002. KLOPPER, Marisa et al. Emergence and Spread of Extensively and Totally Drug-Resistant Tuberculosis, South Africa. Emerging Infectious Diseases, v. 19, n. 3, p.449-455, 2013. KOHN, W.; SHAM, L. J. Self-consistent equations including exchange and correlation effects. Physical Review, v. 140, n.4, p.A1133-A1138, 1965. KRASNOW, I.; WAYNE, Lawrence G.. Comparison of Methods for Tuberculosis Bacteriology. Applied Microbiology, v. 18, n. 5, p.915-917, 1969. LANCASTER, Lorna E. et al. Colicin E3 cleavage of 16S rRNA impairs decoding and accelerates tRNA translocation onEscherichia coliribosomes. Molecular Microbiology. v. 69, n. 2, p.390-401, 2008. MACK, U. et al. LTBI: latent tuberculosis infection or lasting immune responses to M. tuberculosis? A TBNET consensus statement. European Respiratory Journal, v. 33, n. 5, p.956-973, 2009. MAGNET, S.; BLANCHARD, J. S. Molecular Insights into Aminoglycoside Action and Resistance. Chemical Reviews, v. 105, n. 2, p.477-498, 2005. MASSA, W. Crystal Structure Determination. Berlin, Germany. Springer. 4° ed, 2004. MAUS, Courtney E.; PLIKAYTIS, Bonnie B.; SHINNICK, Thomas M. Mutation of tlyA confers capreomycin resistance in Mycobacterium tuberculosis. Antimicrobial Agents And Chemotherapy. v. 49, n. 2, p.571-577, 2005. MCSHANE, Helen et al. Boosting BCG with MVA85A: the first candidate subunit vaccine for tuberculosis in clinical trials. Tuberculosis, v. 85, n. 1-2, p.47- 52. 2005. 66 MITCHISON, D.; DAVIES, G.. The chemotherapy of tuberculosis: past, present and future. The International Journal Of Tuberculosis And Lung Disease, v. 16, n. 6, p.724-732, 2012. MOAZED, Danesh; NOLLER, Harry F.. Interaction of antibiotics with functional sites in 16S ribosomal RNA. Nature, v. 327, n. 327, p.389-394, 1987. MODOLFLL, Juan; VAZQUEZ, David. The Inhibition of Ribosomal Translocation by Viomycin. Eur J Biochem, v. 81, n. 3, p.491-497, 1977. MOTA, K.b. et al. A quantum biochemistry model of the interaction between the estrogen receptor and the two antagonists used in breast cancer treatment. Computational And Theoretical Chemistry, v. 1089, p.21-27, 2016. MURRAY, Patrick R.; ROSENTHAL, Ken S.; PFALLER, Michael A.. Medical microbiology. 7. ed. Elsevier, 2013. NELSON, David L.; COX, Michael M.. Princípios de Bioquímica de Lehninger. 6. ed. Artmed, 2014. NOGRADY, Thomas; WEAVER, Donald F. Medicinal Chemistry: A Molecular and Biochemical Approach. Oxford University Press. 3ª ed., 2005. NOLLER, Harry F.. Ribosomal RNA and Translation. Annual Review Of Biochemistry, v. 60, n. 1, p.191-227, 1991. PERDEW, John P.; WANG, Yue. Accurate and simple analytic representation of the electron-gas correlation energy. Physical Review B, v. 45, n. 23, p.13244-13249, 1992. PESKE, Frank et al. Conformational Changes of the Small Ribosomal Subunit During Elongation Factor G-dependent tRNA–mRNA Translocation. Journal Of Molecular Biology, v. 343, n. 5, p.1183-1194, 2004. PRATT, William B; TAYLOR, Palmer. Principles of drug action: the basis of pharmacology. Churchill Livingstone. 3ª ed., 1990. 67 PUGLISI, Joseph D.; GREEN, Rachel; BLANCHARD, Scott C.. Approaching translation at atomic resolution. Nature Structural Biology, v. 7, n. 10, p.855-861, 2000. RAMAKRISHNAN, V. et al. Functional insights from the structure of the 30S ribosomal subunit and its interactions with antibiotics. Nature, v. 407, n. 6802, p.340-348, 2000. RODRIGUES, Laura C et al. Effect of BCG revaccination on incidence of tuberculosis in school-aged children in Brazil: the BCG-REVAC cluster- randomised trial. The Lancet, v. 366, n. 9493, p.1290-1295, 2005. SENIOR, B. W.; HOLLAND, I. B. Effect of Colicin E3 upon the 30S ribosomal subunit of Escherichia coli. Proceedings Of The National Academy Of Sciences Of The United States Of America. v. 68, n. 5, p.959-963, 1971. SILVA, J. G.; CARVALHO, I. New insights into aminoglycoside antibiotics and derivatives. Current Medicinal Chemistry, v.14, p.1101-1119, 2007. Sociedade Brasileira de Pneumologia e Tisiologia. III Diretrizes para Tuberculose da Sociedade Brasileira de Pneumologia e Tisiologia. Jornal Brasileiro de Pneumologia, v. 35, n.10, p.1018-1048, 2009. SPELMAN, D. W.; MCDONALD, M.; SPICER, W. J.. Aminoglycoside antibiotic agents: a review. The Medical Journal of Australia, v. 6, n. 151, p.9-346, 1989. SREEVATSAN, Srinand et al. Characterization of rpsL and rrs mutations in Streptomycin-Resistant Mycobacterium tuberculosis isolates from diverse geographic localities. Antimicrobial Agents And Chemotherapy. v. 40, n. 4, p.1024- 1026, 1996. STANLEY, Robin et al. The structures of the anti-tuberculosis antibiotics viomycin and capreomycin bound to the 70S ribosome. Nature, v. 17, n. 3, p.289-293, 2010. 68 STURDY, A. et al. Multidrug-resistant tuberculosis (MDR-TB) treatment in the UK: a study of injectable use and toxicity in practice. Journal Of Antimicrobial Chemotherapy, v. 66, n. 8, p.1815-1820, 2011. SWEANY, Henry C.. The tubercle bacillus in the pulmonary lesion of man: histobacteriology and its bearing on the therapy of pulmonary tuberculosis. Chest Journal, v. 28, n. 6, p.699-701, 1955. WORLD HEALTH ORGANIZATION. Global tuberculosis report, 2016. WORLD HEALTH ORGANIZATION. Guidelines for the programmatic management of drug-resistant tuberculosis, Update. 1ª ed, 2011. WORLD HEALTH ORGANIZATION. Model List of Essential Medicines. 19 ed., 2015. WORLD HEALTH ORGANIZATION. Guidelines on the management of latent tuberculosis infection, 2015 WORLD HEALTH ORGANIZATION. Companion handbook to the WHO guidelines for the programmatic management of drug-resistant tuberculosis, 2014. WONG, Chi-huey et al. Specificity of aminoglycoside antibiotics for the A- site of the decoding region of ribosomal RNA. Chemistry & Biology, v. 5, n. 7, 397-406, 1998. WU, Emilia L. et al. Quantum and Molecular Dynamics Study for Binding of Macrocyclic Inhibitors to Human α-Thrombin. Biophysical Journal, v. 92, n. 12, 4244-4253, 2007. ZHANG, D. W.; CHEN, X. H.; ZHANG, J. Z. H. Molecular caps for full quantum mechanical computation of peptide-water interaction energy. Journal Of Computational Chemistry, v. 24, n. 15, 1846-1852, 2003. ZHANG, da W.; ZHANG, J. Z. H.. Molecular fractionation with conjugate caps for full quantum mechanical calculation of protein–molecule interaction energy. The Journal Of Chemical Physics, v. 119, n. 7, 3599-3605, 2003.