Neste trabalho, a docagem e a dinâmica molecular foram utilizadas para investigar a interação entre oligonucleotídeos e duas bases de Tröger: uma base simétrica - contendo duas proflavinas – e uma base assimétrica – contendo uma proflavina e uma fenantrolina. Nas docagens que utilizaram DNA canônico com gap como receptor, os resultados reproduziram corretamente a estereosseletividade da interação, que favorece os isômeros levorrotatórios em relação aos dextrorrotatórios, e apontou para os modos de interação sugeridos pela literatura: ligação ao sulco menor para a Tröger simétrica e intercalação com ligação ao sulco menor (modo misto) para a Tröger assimétrica. Para as simulações, foram escolhidos quatro complexos contendo os enantiômeros levorrotatórios das bases de Tröger e oligonucleotídeos com e sem gap. Embora os complexos representassem pontos de partida diferentes, houve convergência entre aqueles que continham o mesmo oligonucleotídeo como receptor. Os complexos com DNA com gap (complexos de intercalação) convergiram para um modo misto de interação, com intercalação de uma proflavina, enquanto os complexos com DNA sem gap (complexos de ligação ao sulco menor) convergiram para um modo de interação novo, em que as bases de Tröger se ligam ao sulco menor unicamente através da ponte diazocina, com os substituintes projetados para fora. Os complexos de intercalação apresentaram tempos de residência altos (20 ns) e diminuíram o ângulo de torção da dupla hélice na região do gap, levando os ângulos x e y da cadeia principal para uma região não-canônica, associada à forma A do DNA. Já nos complexos de ligação ao sulco menor, houve mudança tanto do modo como do sítio de interação e, embora tenham ocorrido eventuais curvaturas no DNA, não houve alterações significativas nos ângulos da cadeia principal, de modo que a estrutura global do DNA se manteve dentro da região considerada B-canônica. Dessa forma, os resultados sugerem que a intercalação (com adicionais contatos no sulco menor) seja o modo preferencial de interação dessas bases de Tröger com o DNA. O modo de ligação ao sulco menor, por depender quase que exclusivamente da ponte diazocina, seria independente da quiralidade e, portanto, não seria capaz de explicar enantiosselevidade das bases de Tröger. / In this work, molecular docking and molecular dynamics simulations were applied to investigate the interaction between DNA oligomers and two Tröger bases: a symmetric Tröger base derived from proflavine and a mixed Tröger base containing proflavine and phenanthroline. The dockings that used canonical DNA with an artificial intercalation gap as receptor correctly reproduced the stronger DNA binding affinity of the levorotatory isomers when compared to the dextrorotatory ones. Moreover, the results pointed out to the binding modes suggested, yet not established, by experimental data: minor groove binding for symmetric Tröger base and intercalation with additional minor groove binding for the asymmetric one. For the simulations, four complexes were selected, containing the (-)-isomers of the Tröger bases and DNA oligomers with or without intercalation gap. Interestingly, the systems containing the same DNA as receptor converged to similar binding modes: those containing DNA with gap (intercalation complexes) resulted in a mixed binding – with one proflavine intercalated and the other moiety occupying the minor groove – while those containing DNA without gap (minor groove complexes) lead to a new binding mode, in which the Tröger bases interact with the minor groove mainly through the diazocin bridge. In the intercalation complexes, the ligands presented long residence times (20 ns) and also unwound the double helix in the region of the intercalation gap, leading some backbone angles (x and y) to assume non-canonical values typical of the A form of DNA. In minor groove complexes, the ligands displayed enhanced mobility – changing both the binding mode and the binding site – and, although they promoted some bending of the double helix, the main chain torsion angles of DNA remained in the B-canonical region. Altogether, these results suggest that intercalation of proflavine (with additional contacts in the minor groove) might be the stronger DNA binding mode for these Tröger bases. Since the minor groove binding depends almost exclusively of the diazocin bridge, it appears to perform independently of the ligand chirality and, therefore, could not explain the enantioselectivity of the interaction.

As funções de correlação temporal para o momento de dipolo induzido do metano CH4 líquido, foram simuladas pelo método de dinâmica molecular em quatro estados termodinâmicos, a densidades de 28,10; 27,29; 26,41 e 25,34 mol/dm3 e temperaturas de 91,2; 100,7; 110,7 e 122,2 K respectivamente. Foram avaliadas as funções de 2- 3- e 4-corpos, mostraram sua importância no efeito de cancelamento e exibiram no tempo longo decaimentos aproximadamente exponenciais. Incluiu-se os mecanismos de indução emanados pelos campos elétricos dos momentos de octopolo e hexadecapolo, seus campos gradientes, assim como a superposição eletrônica. Con- seguiu-se reproduzir com concordância muito boa os espectros experimentais, no intervalo de freqüências de 0-600 cm-1. O dipolo induzido por octopolo representa a mais importante contribuição ao espectro de absorção no infravermelho longínquo, seguido do dipolo induzido por hexadecapolo. Os momentos espectrais zero estão perto dos valores experimentais, tendo o primeiro estado termodinâmico o melhor valor, havendo uma porcentagem média de erro dos quatro estados termodinâmicos igual ao 6 % a mais do valor experimental. / Time correlations functions for the induced dipole moment for liquid methane were simulated by molecular dynamics in four states of density 28,10, 27,29, 26,41, 25,34 mol/dm3 and temperature 91,2, 100,2, 110,7, 122,2 K respectively. The calculation of n-body terms show their importance on cancellation effects and exhibit in long time aproximately exponential decay. The induction mechanisms included were octupole and hexadecapole moments emanated by the electric fields and their field gradients, as well as electronic overlap. The simulated spectra described agrees very good with the experimental spectra in the frequencies of 0-600 cm−1, the octupole induced dipole represents the most important contribution to the far infrared absorption spectrum, followed by the hexadecapole induced dipole. The zero spectral moments were nearly the experimental values, having the first thermodynamic state the best value, resulting in an average percentual of errors in the four thermodynamics states of 6% over the experimental value.

Fosfolipases A2 constituem uma família de enzimas interfaciais que catalisam a hidrólise de fosfolipídios de membrana na posição sn-2. Nas fosfolipases cataliticamente ativas (FLA2 D49), o resíduo Asp49 possui um importante papel na coordenação do íon cálcio divalente que auxilia na estabilização do intermediário tetraédrico. Nesse trabalho, foram realizadas simulações por dinâmica molecular (DM) de uma fosfolipase secretada de veneno de Agkistrodon halys pallas utilizando como solventes: água, metanol e octanol. Devido à baixa afinidade do íon cálcio em ambientes aquosos, diferentes forças de restrição foram aplicadas nesse íon e mediu-se o comportamento desse, através da variação da sua distância em relação ao resíduo Asp49. Todas as simulações foram realizadas por um período de tempo de 5 ns, somando um total de 45 ns de simulação, utilizando o pacote GROMACS implementado com o campo de força GROMOS96 43a2. Os resultados mostraram uma série de mudanças conformacionais importantes na fosfolipase durante a sua aproximação à superfície da membrana, como o aumento da área acessível ao solvente e a reorientação de cadeias laterais de resíduos com importância catalítica. Os dados obtidos com as simulações foram utilizados para propor um modelo interfacial de ativação para FLA2sv. Desta forma, os principais eventos associados com a ativação de FLA2sv envolvem: a) a reorganização da alça de coordenação de cálcio na proximidade da membrana, seguida pela captura do íon Ca2+; b) a reorientação da cadeia lateral de Trp31, que irá definir um novo sítio de especificidade para a cadeia fosfolipídica; c) a redução da distância entre os resíduos His48 e Asp49, para aumentar a nucleofilicidade da His48; d) a reorientação da cadeia lateral do resíduo Lys69 concomitante com a projeção da alça-69 para o solvente, resultando no aumento da acessibilidade ao sítio ativo. Assim, esses resultados, em adição ao conhecimento bioquímico e a dados cristalográficos, abrem acesso a uma compreensão mais detalhada dos mecanismos moleculares associados com a ação de fosfolipases A2. / Phospholipases A2 (PLA2) are interfacial active enzymes that catalyze hydrolysis at the sn-2 position of membrane phospholipids. In the catalytic active PLA2 (D49 PLA2), the Asp49 residue plays a key role in the coordination of the calcium ion that assists in the stabilization of the transition state intermediary. Despite the catalytic importance of Ca2+, it has low affinity for the sPLA2 active site. Molecular dynamics simulations of a D49 svPLA2 from Agkistrodon halys pallas were performed in aqueous environment, methanol and octanol. Because of the svPLA2 low affinity for the calcium ion in aqueous systems, different restrain forces were applied on Ca2 + and evaluated with respect to ion coordination with key residues. All simulations were performed for a time period of 5 ns using the GROMACS package with GROMOS96 43a2 force field. The results have showed a series of important conformation changes of sPLA2v during its approaches to membrane surface, as the increase of the solvent accessible surface area and the reorientation of side chains of residues of catalytic importance. Analyses of the simulations data were utilized to propose an interfacial activation model for svPLA2. According, the principal events associated with svPLA2 activation involve: a) the reorganization of calcium binding loop at membrane proximity followed by the Ca2+ up take; b) side chain reorientation of Trp31, which will define a new specificity pocket for the phospholipid chain; c) reduction of the distance between His48 and Asp49, in order to increase the nucleophilicity of His48; d) side chain reorientation of Lys69 concomitant with the projection of the 69-loop to solvent resulting in an increase of the active site accessibility. These results, in addition with biochemical and crystallographic data knowledge, open access to a more detailed understanding for the molecular mechanism associated with the interfacial activation of phospholipases A2.

Acidentes envolvendo abelhas africanizadas são freqüentemente relatados, particularmente na América do Sul. As picadas de abelhas causam reações localizadas e sistêmicas e os sintomas do envenenamento incluem náuseas, vômitos, hemólise, falência renal e coagulação intravascular disseminada. Durante muito tempo todas as reações tóxicas eram atribuídas à presença de uma fosfolipase A2, mesmo sendo o veneno uma mistura complexa de substâncias. A melitina, o componente mais abundante e tóxico do veneno de abelha, é um peptídeo de 26 aminoácidos com a habilidade de interagir e danificar membranas celulares. A melitina é também capaz de modular muitas proteínas, aumentando a diversidade de atividades biológicas do peptídeo. Até recentemente, acreditava-se que a estrutura tri-dimensional biologicamente ativa do peptídeo, que possui conformação em hélice, era um tetrâmero. Neste trabalho avaliamos a conformação da melitina e seus estados oligoméricos em solução por dinâmica molecular e a interferência da melitina na função plaquetária. Aqui está demonstrado que a melitina possui uma conformação randômica em condições fisiológicas e que sua estrutura tri-dimensional sofre alterações de acordo com as condições ambientais. Ainda, foi demonstrada uma nova atividade biológica do peptídeo melitina. O peptídeo é capaz de induzir a 8 agregação plaquetária de forma dose-dependente e de interagir diretamente com a superfície de plaquetas. A correlação entre a conformação da melitina e suas atividades biológicas é discutida. Os resultados aqui apresentados podem ser valiosos no entendimento do papel da melitina nas coagulopatias induzidas por veneno de abelha. O estudo estrutural mostrado aqui pode ser aplicado para explicar as diferentes atividades do peptídeo. / Accidents involving africanized bees are frequently reported, particularly in South America. Bee stings cause localized and systemic reactions and the symptoms of envenomation include nausea, vomiting, hemolysis, kidney failure and disseminated intravascular coagulation. For a long time, all toxic reactions were ascribe to the presence of a phospholipase A2, despite of being the venom a complex mixture of substances. Melittin, the most abundant and the major toxic component of bee venom, is a 26 amino acid peptide with the ability to interact and disrupt cell membranes. Melittin is also able to modulate many proteins, enhancing the wide range of the peptide biological activities. The biologically active tridimensional structure of the peptide, which has a helical conformation, has been described until now as a tetramer. In this work we evaluated the conformation of melittin and its oligomeric states in solution by molecular dynamics simulations and performed studies of melittin effect on platelet function. Here we demonstrate that melittin has a random conformation under physiological conditions and its tridimensional structure changes under different environmental conditions. Moreover, here we describe a new biological activity of melittin. The peptide is able to induce platelet aggregation in a dose10 dependent manner and can interact directly with the platelet surface. The correlation between melittin conformation and biological activity is discussed. Our results might contribute to elucidate the role of melittin in bee venom induced coagulopathies. The structural data gathered in this work may explain the different activities of the peptide.

O presente trabalho relata a parametrização do modelo de contínuo IEF-PCM para o estudo da solvatação de ânions monovalentes em acetonitrila e N,N-dimetilformamida. Nós propomos duas metodologias de parametrização do modelo de contínuo polarizável para este estudo teórico. A primeira proposta é baseada unicamente no cálculo do termo eletrostático da energia livre e na equiparação deste termo com a energia livre de solvatação experimental. Na segunda proposta de metodologia, a parametrização do modelo de contínuo polarizável foi realizada com o auxílio das simulações de dinâmica molecular e dentro desta segunda proposta nós calculamos todos os termos que contribuem para a energia livre de solvatação, a energia de cavitação. a energia eletrostática e a energia de van der Waals. Com a aplicação da primeira metodologia, é obtido um erro médio absoluto de 2,1 kcal/mol nas energias de solvatação dos ânions em acetonitrila enquanto que o erro médio absoluto calculado para energias de solvatação dos ânions em N,N-dimetilformamida é 2,8 kcal/mol. A aplicação da segunda metodologia resulta num erro médio absoluto de 1,4 kcal/mol nas energias de solvatação dos ânions em acetonitrila e um erro médio absoluto de 1,2 kcal/mol nas energias de solvatação calculadas para os ânions em N,N-dimetilformamida. Estes resultados indicam que o modelo de contínuo polarizável IEF-PCM é adequado para o estudo dos efeitos termodinâmicos nas soluções de ânions nos solventes acetonitrila e N,N-dimetilformamida.

Fosfolipases A2 constituem uma família de enzimas interfaciais que catalisam a hidrólise de fosfolipídios de membrana na posição sn-2. Nas fosfolipases cataliticamente ativas (FLA2 D49), o resíduo Asp49 possui um importante papel na coordenação do íon cálcio divalente que auxilia na estabilização do intermediário tetraédrico. Nesse trabalho, foram realizadas simulações por dinâmica molecular (DM) de uma fosfolipase secretada de veneno de Agkistrodon halys pallas utilizando como solventes: água, metanol e octanol. Devido à baixa afinidade do íon cálcio em ambientes aquosos, diferentes forças de restrição foram aplicadas nesse íon e mediu-se o comportamento desse, através da variação da sua distância em relação ao resíduo Asp49. Todas as simulações foram realizadas por um período de tempo de 5 ns, somando um total de 45 ns de simulação, utilizando o pacote GROMACS implementado com o campo de força GROMOS96 43a2. Os resultados mostraram uma série de mudanças conformacionais importantes na fosfolipase durante a sua aproximação à superfície da membrana, como o aumento da área acessível ao solvente e a reorientação de cadeias laterais de resíduos com importância catalítica. Os dados obtidos com as simulações foram utilizados para propor um modelo interfacial de ativação para FLA2sv. Desta forma, os principais eventos associados com a ativação de FLA2sv envolvem: a) a reorganização da alça de coordenação de cálcio na proximidade da membrana, seguida pela captura do íon Ca2+; b) a reorientação da cadeia lateral de Trp31, que irá definir um novo sítio de especificidade para a cadeia fosfolipídica; c) a redução da distância entre os resíduos His48 e Asp49, para aumentar a nucleofilicidade da His48; d) a reorientação da cadeia lateral do resíduo Lys69 concomitante com a projeção da alça-69 para o solvente, resultando no aumento da acessibilidade ao sítio ativo. Assim, esses resultados, em adição ao conhecimento bioquímico e a dados cristalográficos, abrem acesso a uma compreensão mais detalhada dos mecanismos moleculares associados com a ação de fosfolipases A2. / Phospholipases A2 (PLA2) are interfacial active enzymes that catalyze hydrolysis at the sn-2 position of membrane phospholipids. In the catalytic active PLA2 (D49 PLA2), the Asp49 residue plays a key role in the coordination of the calcium ion that assists in the stabilization of the transition state intermediary. Despite the catalytic importance of Ca2+, it has low affinity for the sPLA2 active site. Molecular dynamics simulations of a D49 svPLA2 from Agkistrodon halys pallas were performed in aqueous environment, methanol and octanol. Because of the svPLA2 low affinity for the calcium ion in aqueous systems, different restrain forces were applied on Ca2 + and evaluated with respect to ion coordination with key residues. All simulations were performed for a time period of 5 ns using the GROMACS package with GROMOS96 43a2 force field. The results have showed a series of important conformation changes of sPLA2v during its approaches to membrane surface, as the increase of the solvent accessible surface area and the reorientation of side chains of residues of catalytic importance. Analyses of the simulations data were utilized to propose an interfacial activation model for svPLA2. According, the principal events associated with svPLA2 activation involve: a) the reorganization of calcium binding loop at membrane proximity followed by the Ca2+ up take; b) side chain reorientation of Trp31, which will define a new specificity pocket for the phospholipid chain; c) reduction of the distance between His48 and Asp49, in order to increase the nucleophilicity of His48; d) side chain reorientation of Lys69 concomitant with the projection of the 69-loop to solvent resulting in an increase of the active site accessibility. These results, in addition with biochemical and crystallographic data knowledge, open access to a more detailed understanding for the molecular mechanism associated with the interfacial activation of phospholipases A2.

Prostaglandina Endoperóxido Sintases (PGHSs) são enzimas homodiméricas integrais de membrana, N-glicosiladas, localizadas no lúmen do retículo endoplasmático, onde catalisam o primeiro passo na síntese de prostanóides, produzindo prostaglandina H2 a partir do ácido araquidônico (AA). Constituem-se, assim, em um dos principais alvos terapêuticos no tratamento de doenças inflamatórias. PGHSs apresentam duas isoformas: PGHS-1, expressa constitutivamente em diversos tipos celulares, e PGHS-2, usualmente expressa por indução. Ambas isoformas são glicosiladas nos resíduos Asn68, Asn144 e Asn410, modificações necessárias para a atividade e função adequadas destas enzimas. A reação ciclooxigenase necessita que o AA esteja numa orientação catalítica. Adicionalmente a esta orientação, o banco de dados Protein Data Bank contém estruturas cristalográficas que apresentam uma orientação alternativa para o AA, proposta como não-catalítica. Considerando a importância da glicosilação para a função da enzima, o presente trabalho tem como objetivo a caracterização estrutural e conformacional de PGHSs em suas formas glicosiladas, bem como analisar as diferentes orientações adotadas pelo AA em dados cristalográficos prévios. Assim, modelos de PGHS-1 e PGHS-2 glicosiladas foram construídos combinando diferentes fontes de dados experimentais (bioquímicos, cristalográficos e espectrométricos) e foram submetidos a simulações de dinâmica molecular (DM) usando o pacote de simulação do GROMACS e campo de força GROMOS9643a1. Parâmetros desenvolvidos pelo grupo foram empregados para descrever a porção sacarídica dos sistemas. Os dados obtidos indicam uma diminuição na flexibilidade da proteína na presença de glicosilação no domínio EGF em ambas as isoformas. Além disso, PGHS-1 e PGHS-2 ligadas às duas diferentes orientações do AA foram submetidas a estudos de DM, indicando que a orientação alternativa não é estável, possivelmente sendo uma conformação minoritária em solução. Tal orientação alternativa, contudo, pode sugerir novas abordagens para o planejamento de novos agentes anti-inflamatórios baseado na complementaridade à enzima-alvo. / Prostaglandin Endoperoxide Synthases (PGHSs) are homodimeric integral membrane enzymes, N-glycosylated, located in the endoplasmic reticulum lumen, where they catalyze the first step in the prostanoids synthesis, generating prostaglandin H2 from arachidonic acid (AA), so constituting one of the main therapeutic targets for treating inflammatory diseases. PGHSs present two isoforms: PGHS-1, constitutively expressed in several cell types, and PGHS-2, usually expressed by induction. Both isoforms are glycosylated at residues Asn68, Asn144 and Asn410, modifications necessary for adequate PGHSs enzymatic activities and function. The cyclooxygenase reaction requires AA to be in a catalytic orientation. Additionally to this orientation, the Protein Data Bank contains crystallographic structures presenting an alternative orientation for AA, proposed as non-catalytic. Considering the importance of glycosylation for the enzyme`s function, the current work intents to structurally and conformationally characterize the role of the sacharidic moiety of PGHSs, as well as to analyze the different orientations adopted by AA found in previous crystallographic data. Thus, glycosylated models for PGHS-1 and PGHS-2 were built combining different sources of experimental data (biochemical, crystallographic and spectrometric) and submitted to molecular dynamics (MD) using the GROMACS simulation suite and GROMOS9643a1 force field. Parameters developed by the group were employed in order to describe the carbohydrate moiety. The obtained data indicate a decrease in protein flexibility in the presence of glycosylation on the EGF domain of both isoforms. Also, PGHS-1 and PGHS-2 bound to the two different orientations of AA were submitted to MD studies indicating that the alternative orientation is not stable, possibly being a minority orientation in solution. This orientation, however, may suggest new approaches for structure-based design of new anti-inflammatory agents.

Combinando métodos computacionais da eletrostática e dinâmica molecular, este trabalho teve como objetivo descrever processos de solvatação em sistemas quimicamente importantes. Foram obtidas propriedades termodinâmicas necessárias para o entendimento do processo de solvatação. Estão descritos e testados modelos para descrever a interação soluto-solvente, possibilitando, assim, aprimorar a descrição físico-química dos processos de solvatação. Utilizaram-se programas desenvolvidos em nosso grupo e programas comerciais que permitem os cálculos de dinâmica molecular e química quântica. Uma nova abordagem para o cálculo de energia livre de solvatação foi desenvolvida proporcionando a obtenção acurada e eficiente dessa propriedade, dentro do enfoque da dinâmica molecular. Nessa nova abordagem, novas metodologias para a geração de cavidades moleculares foram propostas e avaliadas. As energias livres de solvatação obtidas estão em boa concordância com os valores experimentais.

Ureases, apesar de sua importância histórica para a Bioquímica, ainda são enigmáticas. Pouco se sabe a respeito de suas propriedades não-catalíticas, por exemplo. Nesse contexto, o presente trabalho submeteu uma urease de Canavalia ensiformis (JBURE-IIB) a simulações de dinâmica molecular, nos níveis de detalhamento atomístico e coarse-grained, buscando obter uma caracterização conformacional dessa urease, tanto como monômero quanto como trímero. Os dados assim obtidos sugerem que os domínios da proteína apresentam diferentes mobilidades, sendo o domínio γ o mais flexível. Esse domínio é também o principal responsável pela compactação do monômero da urease isolado, quando comparado à sua conformação cristalográfica. A mobilidade do domínio γ foi também associada ao processo de oligomerização da enzima. A forma trimérica da urease de C. ensiformis também se tornou mais compacta após a simulação, apontando as ureases como moléculas altamente flexíveis, principalmente devido às regiões de alças que conectam domínios estruturais, no caso de ureases vegetais. Sendo uma propriedade não necessariamente associada à catálise, a análise de tal flexibilidade pode auxiliar tentativas futuras de explicar a participação não-catalítica da urease em suas propriedades multifuncionais. / Ureases, despite their historical importance in biochemistry, are still enigmatic. For instance, little is known about their non-catalytic properties. In this context, the current work submitted a Canavalia ensiformis urease (JBURE-IIB) to molecular dynamics simulations at both atomistic and coarse-grained levels in order to obtain a conformational characterization of such urease in both monomeric and trimeric forms. The data suggest that the domains of the protein present different mobilities, with the γ domain being the most flexible one. This domain is also the main responsible for the compactness of the isolated urease monomer when compared to its crystallographic conformation. The γ domain mobility was also associated with the oligomerization process of the enzyme. The trimeric form of C. ensiformis urease also became more compact after simulation, pointing to ureases as highly flexible molecules, mainly due to coiled regions connecting the structural domains, in the case of plant ureases. As a property not necessarily related to catalysis, the analysis of this flexibility may enlighten further efforts to solve the non-catalytic aspects of urease in its moonlighting properties.

Canais iônicos são nanoporos aquosos formados por proteínas imersas na membrana celular. Sua função biológica básica é permitir o transporte de íons no sentido de seu gradiente eletroquímico. Dada a enorme quantidade de processos biológicos relacionados ao movimento iônico, a importância dos canais iônicos em biologia e medicina dificilmente pode ser exagerada. Neste trabalho, apresentamos uma introdução à vasta área do conhecimento de canais iônicos, discutimos suas perspectivas em pesquisa básica e aplicada e, em especial, a importância do estudo computacional da permeação iônica por canais. Seguimos analisando o estado da arte em relação a modelos computacionais de canais iônicos e discutindo criticamente as vantagens e desvantagens das abordagens existentes, fazendo considerações sobre suas perspectivas. Concluindo que a técnica de Dinâmica Browniana apresenta potencial para o estudo de propriedades básicas da condução iônica por nanoporos, propomos um modelo simples de canal iônico unidimensional através desta técnica para uma geometria que permite solução analítica da equação de Poisson, a qual é apresentada. A simplicidade do modelo proposto lhe confere grande eficiência computacional, a qual é ampliada com a proposta de um critério de “tempo de injeção” dos íons no canal, cuja expressão é demonstrada e que poupa o simulador de considerar explicitamente dois reservatórios eletrolíticos junto ao canal. O simulador desenvolvido baseado neste modelo é estudado e comparado com a literatura, mostrando bons resultados em sua região esperada de validade. Dado o estágio atual do estudo computacional de canais iônicos, é importante dispormos de modelos conceitualmente simples e computacionalmente eficientes. A sua existência possibilita, por exemplo, estudos de Dinâmica Molecular muito eficientes considerando a água explicitamente, algo importante para a análise de novos modelos da água adequados a nanoporos. / Ion channels are nanopores formed through the cell membrane by proteins. Their basic biological function is allowing ionic transport down the electrochemical gradient. Given the huge amount of biological processes that are related to ionic fluxes, the biological and medical importance of ion channels hardly can be overstated. In this work we present a brief introduction to the vast field of ion channels and discuss its perspectives in basic and applied research, emphasizing the importance of computational approaches to study ion permeation through channels. Next, we analyze the state of the art regarding computational models of ion channels, critically discussing the advantages and disadvantages of the current approaches and considering the perspectives of each. We conclude that the Brownian Dynamics approach has potential to be used to study ionic conductance through nanopores and we propose a simple Brownian Dynamics unidimensional channel model for a geometry that allows analytical solution of the Poisson equation, which is presented. The simplicity of the proposed model grants it great computational efficiency, which is further amplified with the proposal of an “injection time” criterion to allow ions to enter the channel. The formula for the “injection time” is demonstrated and its implementation lifts the requirement to simulate two reservoirs next to the channel. The simulator developed based on this model is studied and its results are compared with the literature, showing good agreement in its expected region of validity. Given the current state of computational study of ion channels, it is important to have simple and efficient models, such as ours, available. For instance, their existence allows the realization of very fast Molecular Dynamics studies that consider water explicitly, something very important for the analysis of new explicitly polarizable water models that are required for nanopores.