Este estudo tem como finalidade principal a análise termodinâmica e estatística de proteínas através de modelos mínimos. Uma proteína é um polímero de aminoácidos, cuja função está essencialmente relacionada às conformações espaciais que ela adota em solução aquosa. Na forma funcional (dita nativa), essas conformações flutuam levemente em torno de um mínimo de energia-livre. O processo pelo qual uma cadeia protéica transita de estados não-nativos para a estrutura nativa é chamado de enovelamento, ou dobramento. Uma questão em aberto no campo de estudo de proteínas consiste justamente em entender a fundo o processo de enovelamento, cujo avanço tem um vasto potencial de aplicação, desde a predição de estruturas a partir de sequências de aminoácidos até o planejamento de fármacos e moléculas bioativas. Nossa investigação teórica procura abordar aspectos do enovelamento expressos através de grandezas termodinâmicas (energia média, calor específico, número de ligações de hidrogênio, entre outras) derivadas de modelos estatísticos na rede. Assim, num primeiro momento, analisamos o chamado modelo HP, ora por meio de enumeração exata, para cadeias curtas, ora por simulações de Monte Carlo, para cadeias maiores. No primeiro caso, propusemos a existência de uma relação entre a ocorrência de um segundo pico no calor específico associado na literatura à transição de congelamento com uma drástica redução no número de configurações entre os primeiros estados excitados e aqueles de menor energia. Observamos, também, que esse pico pode aparecer tanto para homopolímeros quanto para heteropolímeros, em ambas as redes quadrada e triangular. Num segundo momento, nosso enfoque se voltou para a inclusão de um solvente aquoso (dado pelo modelo de Bell-Lavis) ao sistema inicial. Isso nos possibilitou verificar, usando exclusivamente simulações de Monte Carlo e o algoritmo de Metropolis, o comportamento e a competição das ligações de hidrogênio água-água, água-proteína, proteína-proteína e na primeira camada de solvatação. O modelo acoplado exibiu algumas características do enovelamento, como o colapso hidrofóbico e a separação de monômeros (apolares no núcleo e polares na superfície), embora não capture a desnaturação fria. No apêndice, adicionamos algumas propostas para realização do cálculo numérico da pressão no ensemble canônico, desenvolvidas em paralelo ao projeto principal desta dissertação, mas que, numa primeira análise, verificamos serem consistentes e passíveis de futuros desdobramentos. / The finality of this study is to analyse proteins thermodynamics and statistics through minimal models. A protein is a polymer of amino acids, whose spatial conformations in aqueous solution determine its function. In the functional form (said native), those conformations fluctuates slightly around a free-energy minimum. The process by which a protein chain passes from non-native states to a stable native structure is called protein folding. An open question in the field of protein studies is to understand more deeply the folding process, whose advance can find a wide range of potential applications, since ab initio structure prediction from the amino acids sequence to biomolecules design. The theoretical approaches used here focus on aspects of protein folding given by some thermodynamic quantities (as mean energy, specific heat, number of hydrogen bonds and so on) obtained from statistical lattice models. Initially, we analyse the so-called HP model, at first using exact enumeration for short chains, then by Monte Carlo simulations for longer chains. In the first case, we propose a correlation between the occurrence of a second peak in the specific heat associated in the literature with a freezing transition and a sharp reduction on the number of configurations from the first excited states to the lowest energy states. In addition, we observe that this peak may appear to both homopolymers and heteropolymers on square and triangular lattices. At a second moment, our focus turned to the introduction of a water-like solvent (Bell-Lavis model) to the initial system. This allowed us to verify, exclusively by means of Monte Carlo simulations with Metropolis algorithm, the behavior and competition of hydrogen bonds between water-water molecules, water-protein, and protein-protein monomers and at the first hydration layer. The combined model showed some classical folding properties, as hydrophobic collapse and monomers segregation (apolar residues at the core and polar residues at the surface), although it did not capture cold denaturation. We have included in the appendix some proposals to perform numerical calculations of the canonical pressure, which were developed alongside the main subject of this thesis and a first analysis has proved to be consistent and susceptible to further developments.
Identifer | oai:union.ndltd.org:IBICT/oai:teses.usp.br:tde-20022018-190857 |
Date | 22 September 2017 |
Creators | Fernando Takeshi Tanouye |
Contributors | Vera Bohomoletz Henriques, Carlos Eduardo Fiore dos Santos, Mário Noboru Tamashiro |
Publisher | Universidade de São Paulo, Física, USP, BR |
Source Sets | IBICT Brazilian ETDs |
Language | Portuguese |
Detected Language | Portuguese |
Type | info:eu-repo/semantics/publishedVersion, info:eu-repo/semantics/masterThesis |
Source | reponame:Biblioteca Digital de Teses e Dissertações da USP, instname:Universidade de São Paulo, instacron:USP |
Rights | info:eu-repo/semantics/openAccess |
Page generated in 1.0135 seconds