Amyloidogenic diseases, such as Alzheimer's disease (AD) and Type II Diabetes (T2D), are characterized by the accumulation of amyloid aggregates. Despite having very different amino-acid sequences, the underlying amyloidogenic proteins form similar supramolecular fibril structures that are highly stable and resistant to physical and chemical denaturation. AD is characterized by two toxic lesions: extracellular amyloid β-peptide (Aβ) plaques and intracellular neurofibrillary tangles composed of microtubule-associated protein tau. Similarly, a feature of T2D is the deposition of islet amyloid polypeptide (IAPP) aggregates in and around the pancreas. The mechanisms by which Aβ, tau, and IAPP aggregate, and cause cell death is unknown; thus, gaining greater insight into the stabilizing forces and initial unfolding events is crucial to our understanding of these amyloidogenic diseases. This work uses molecular dynamics (MD) simulations to study the secondary, tertiary, and quaternary structure of Aβ, tau, and IAPP. Specifically, this work used the Drude polarizable force field (FF), which explicitly represents electronic polarization allowing charge distributions to change in response to perturbations in local electric fields. This model allows us to describe the role charge plays on protein folding and stability and how perturbations to the charge state drive pathology. Studies were conducted to address the following questions: 1) What are the stabilizing forces of fibril and oligomeric structures? 2) How do charge-altering mutations modulate the conformational ensemble and thermodynamic properties of Aβ? 3) How do charge-altering post-translational modifications of Aβ and tau modulate changes in the conformational ensembles? These studies establish that shifts in local microenvironments play a role in fibril and oligomer stability. Furthermore, these studies found that changes in protein sequence and charge are sufficient to disrupt and change the secondary and tertiary structure of these amyloidogenic proteins. Overall, this dissertation describes how charge modulates protein unfolding and characterizes the mechanism of those changes. In the long term, this work will help in the development of therapeutics that can target these changes to prevent protein aggregation that leads to cell death. / Doctor of Philosophy / Protein aggregation is the hallmark of many chronic diseases, such as Alzheimer's disease (AD) and Type II Diabetes (T2D). The formation of two toxic aggregates: amyloid β-peptide (Aβ) plaques and neurofibrillary tangles composed of microtubule-associated protein tau are some of the key characteristics of AD. In addition, the formation of islet amyloid polypeptide (IAPP) aggregates in the pancreas is thought to play a role in the development of T2D. The pathways by which the proteins Aβ, tau, and IAPP aggregate are unknown; thus, gaining a greater insight into the properties that may cause these diseases is necessary to develop treatments. By studying these proteins at the atomistic level, we can understand how small changes to these proteins alter how they misfold in a way that promotes toxicity. Herein, we used a computational technique called molecular dynamics (MD) simulations to gain new insights into how protein structure changes. We explored the dynamics of these proteins and investigated the role that charge plays in protein folding and described how charge modulates protein folding and characterized the mechanism of those changes. This work serves as a characterization of protein folding and sets the ground for future structural studies and drug development.

Doctor of Philosophy / Department of Chemistry / Paul E. Smith / The Kirkwood-Buff (KB) theory of solutions, published in 1951, established a route from integrals over radial (pair) distribution functions (RDFs) in the grand canonical ensemble to a set of thermodynamic quantities in an equivalent closed ensemble. These “KB integrals” (KBIs) can also be expressed in terms of the particle-particle (i.e., concentration or density) fluctuations within grand canonical ensemble regions. Contributions by Ben-Naim in 1977 provided the means to obtain the KBIs if one already knew the set of thermodynamic quantities for the mixture of interest; that is, he provided the inversion procedure. Thus, KB theory provides a two-way bridge between local (microscopic) and global (bulk/thermodynamic) properties. Due to its lack of approximations, its wide ranging applicability, and the absence of a competitive theory for rigorously understanding liquid mixtures, it has been used to understand solution microheterogeneity, solute solubility, cosolvent effects on biomolecules, preferential solvation, etc. Here, after using KB theory to test the accuracy of pair potentials, we present and illustrate two extensions of the theory, resulting in a general Fluctuation Solution Theory (FST). First, we generalize KB theory to include two-way relationships between the grand canonical ensemble’s particle-energy and energy-energy fluctuations and additional thermodynamic quantities. This extension allows for non-isothermal conditions to be considered, unlike traditional KB theory. We illustrate these new relationships using analyses of experimental data and molecular dynamics (MD) simulations for pure liquids and binary mixtures. Furthermore, we use it to obtain conformation-specific infinitely dilute partial molar volumes and compressibilities for proteins (other properties will follow) from MD simulations and compare the method to a non-FST method for obtaining the same properties. The second extension of KB theory involves moving beyond doublet particle fluctuations to additionally consider triplet and quadruplet particle fluctuations, which are related to derivatives of the thermodynamic properties involved in regular KB theory. We present these higher order fluctuations obtained from experiment and simulation for pure liquids and binary mixtures. Using the newfound experimental third and fourth cumulants of the distribution of particles in solution, which can be extracted from bulk thermodynamic data using this extension, we also probe particle distributions’ non-Gaussian nature.

In most applications concerning a fluid flowing over a solid surface, the no-slip velocity condition was widely used because it is simple and produces the results in agreement with experiments. However, this dynamical boundary condition is not appropriate when the flow under consideration is at a micro or nano length scale.In order to model this effect at the macroscopic scale, the Navier boundary conditions have been introduced, with the slip length as a parameter. When the fluid is a gas, this length is related to the tangential momentum accommodation coefficient (TMAC) and the mean free path, according to the Maxwell model. The aim of this work is to systematically address this model using a multi-scale approach and to extend it by incorporating both the morphology and the anisotropy of a surface. The thesis consists of five chapters. In Chapter 1, the basics of the kinetic theory of gases, the Boltzmann equation and related solutions (Navier-Stokes-Fourier, Burnett, Grad, Direct Simulation Monte Carlo ...) are briefly presented. The models of gas-wall interaction and slip models introduced in the fluid mechanics are also recalled. The chapter ends with a description of the computational method used for the molecular dynamics simulations performed in this work. Chapter 2 is dedicated to the development of a simple technique to simulate the pressure driven flows. The principle is to rely on the atomistic formulas of the stress tensor (Irving Kirkwood, Method of Plane, Virial Stress) and to modify the periodic conditions by maintaining the difference between the kinetic energy of the ingoing and outgoing particles of the simulation domain. Several types of channels are studied with this technique. The results (temperature, velocity ...) are discussed and compared. Chapter 3 deals with the study of the gas-wall interaction potential by the ab-initio method. The code CRYSTAL 09 is used to obtain the potential between an atom of argon (Ar) and a surface of platinum (Pt) <111> as a function of distance. Then the gas-wall potential is decomposed into binary potential and approached by an analytic function. This function is then implemented in a MD code to simulate the gas-wall collisions and determine the TMAC coefficient. In Chapter 4, the effect of morphology is studied. The multi-body Quantum Sutton Chen (QSC) potential is used for Pt <100> solid and the binary potential proposed in the previous chapter for the Ar-Pt couple is employed. The QSC potential is needed to reproduce the surface effects that affect the final results. Different surfaces are treated : smooth, nanostructured surface and, random surface obtained by Chemical vapor deposition (CVD). The TMAC is determined using a generalized approach, i.e. depending on the angle of incident flux of gas atoms on the surface. The surface anisotropy and the scattering kernel are also examined. In Chapter 5, we propose a model of anisotropic slip for fluids based on accommodation tensor. The model is obtained by the analytical approximate calculations developed in the framework of the kinetic theory. We thus generalize Maxwell's equation by showing that the slip length tensor is directly related to the accommodation tensor. The model is in good agreement with the MD results. Thanks to our MD simulations, we develop a suitable technique for reproducing the anisotropy of the accommodation tensor. The thesis ends with a conclusion section in which we suggest some perspectives for a continuation of this work

Dans leurs membranes natives, les protéines membranaires peuvent former des assemblages multimériques impliqués dans un grande variété de processus biologiques, de la transduction du signal à la structuration d'organelles. Leur agrégation est influencée par relations qu'elles entretiennent avec la membrane. Au cours de cette thèse, j'étudie l'influence mutuelle des protéines membranaires et de leur environnement au moyen de simulations numériques. L'étude d'un un complexe peptide-détergent modèle a révélé la présence d'une compétition entre deux forces : la cohésion entre les détergents maintenant la structure micellaire, et les interactions adhésives à certaines surfaces peptidiques. L'équilibre entre ces forces conduit à une frustration du système et à l'exposition de surfaces hydrophobes au solvant, expliquant les phénomènes d'agrégations observées expérimentalement. Ces résultats montrent l'importance de la topologie des surfaces du peptide dans l'organisation du complexe. Pour caractériser l'influence de la composante protéique sur les amphiphiles, un estimateur d'entropie configurationelle des chaînes grasses a été développé. Les mesures réalisées sur 15 systèmes membranaires gros grains montrent que la capacité des protéines à apparier l'épaisseur hydrophobe membranaire à la leur est le facteur principal responsable des variations entropiques. Cependant, au delà des simples contributions de l'épaisseur hydrophobe protéique, la diversité des comportements observés dans la première couche de lipides au contact avec la protéine montre bien un rôle de la topologie dans la modulation des interactions. Les variations observées à plus longues distances semblent alors être un facteur clé dans les interactions entre protéines membranaires. Dans le cadre d'une collaboration, la description d'un membranaire par microscopie à force atomique haute vitesse a permis de caractériser le paysage énergétique des interactions entre protéines membranaires. En combinant ces résultats avec ceux obtenus par simulations de dynamiques moléculaires, des chemins basses énergies on pu être identifiés et soulignent de nouveau l'importance des lipides dans l'organisation et la dynamique des agrégats protéiques. L'ensemble de ces données permettent alors de mieux comprendre l'influence mutuelle entre les protéines et les amphiphiles constituant leur environnement et les phénomènes d'agrégation protéiques dans un contexte biologique et expérimental. / In their native membranes, membrane proteins can form multimeric assemblies implicated in a wide range of biological processes, from signal transduction to organelle structure. Their aggregation is influenced by the relations they have with the membrane. Through this thesis, I study the mutual influence between membrane proteins and their environment using numerical simulations. The study of a model protein detergent complex has revealed the presence of two competing forces : cohesion between detergents maintaining the micellar structure, and adhesive interactions to certain peptide surfaces. The balance between these forces leads to a certain frustration of the system, and to the exposure of hydrophobic surfaces to the solvent, explaining phenomena observed experimentally. These results show the importance of peptide surface topology on the organization of the complex. To characterize the influence of the protein component on amphiphiles, a configuartional entropy estimator has been developed. Measures realized on 15 coarse grained membrane systems show that the protein ability to match the membrane hydrophobic thickness to its own is the main factor responsible for entropic variations. However, beyond the simple contributions of protein hydrophobic width, the variety of behaviors in the first shell of lipids show a role of proteins in interaction modulation. Variations observed at longer ranges then seem to be a key factor in the interactions between membrane proteins. In the context of a collaboration, the description of a membrane plane by high speed atomic force microscopy allowed to characterize the energetic landscape of interactions between membrane proteins. By combining these results with those obtained by coarse grained molecular dynamics simulations, low energy paths have been identified and underly the importance of lipids in the organization and dynamics of proteins aggregates. All these data allow to further understand the mutual influence between proteins and amphiphiles constituting their environment and the aggregations phenomena in a biological and experimental context.

As metaloproteinases de matriz (MMP) são enzimas superexpressas em quase todos os tumores humanos, sendo que os subtipos MMP-2 e MMP-9 têm sido associados ao potencial metastático e prognóstico desfavorável em neoplasias malignas como, por exemplo, melanoma metastático e glioma. Compostos capazes de inibir a atividade destas enzimas podem representar potenciais agentes terapêuticos. O composto 4-nerolidilcatecol (4-NC), isolado de plantas do gênero Pothomorphe, apresentou resultados promissores para o tratamento do melanoma e glioma e foi capaz de atuar em várias etapas bioquímicas importantes envolvidas na progressão dessas patologias, inclusive inibindo MMP-2 e MMP-9. No entanto, o mecanismo de ação do 4-NC não está completamente elucidado. O presente estudo envolveu a aplicação de métodos de modelagem molecular e de formalismos do planejamento de novas moléculas auxiliado por computador, CAMD (Computer-Aided Molecular Design) a fim de explorar a interação entre esta molécula e as enzimas MMP-2 e MMP-9, além de planejar novos inibidores para estes alvos. Análise exploratória de dados, que compreende a análise de agrupamentos hierárquicos e de componentes principais. foi desenvolvida para um conjunto de hidroxamatos (N=64) descritos como inibidores de MMP-2 e MMP-9, a fim de identificar as propriedades moleculares que mais influenciavam o processo de discriminação dos compostos. As propriedades termodinâmicas, eletrônicas e estéricas foram importantes para descrever os compostos mais ativos no conjunto de dados da MMP-2. Para a MMP-9, o coeficiente de distribuição (ClogD) em pH 1,5 foi relevante no processo de discriminação do conjunto. A presença de substituintes volumosos na porção R3 parece ser crucial para o conjunto de inibidores investigados. Esta região está envolvida em interações moleculares com a cavidade S1 de ambas as enzimas, mas há um limite de volume a ser considerado para estes substituintes. O formalismo QSAR-4D independente do receptor (IR) foi aplicado ao mesmo conjunto de dados e permitiu estabelecer o mapeamento do farmacóforo, além de explorar diferentes alinhamentos para a obtenção da hipótese de conformação bioativa prevista pelo melhor modelo de QSAR. OS modelos QSAR apresentaram boa capacidade de previsão, auxiliaram na proposição de novos inibidores e estimaram a atividade do 4-NC. Com o melhor modelo QSAR para MMP-9 (N=64), a atividade prevista para o 4-NC foi classificada na faixa dos inibidores com atividade moderada. Entretanto, o melhor modelo QSAR obtido para MMP-2 (N=38) não foi capaz de prever, de forma adequada, a atividade de compostos com arcabouço químico diferente daqueles utilizados na construção dos modelos. Estudos de ancoramento molecular foram desenvolvidos para investigar a orientação do 4-NC no sitio catalítico das duas enzimas e as interações que poderiam ser estabelecidas nestes complexos. Duas conformações favoráveis foram encontradas. Simulações computacionais de dinâmica molecular foram desenvolvidas com os complexos mais promissores selecionados nos estudos de ancoramento, a fim de obter informações mais detalhadas e de maior confiabilidade. sobre suas interações intermoleculares. O 4-NC tende a se orientar no sítio de forma a acomodar sua cadeia lateral no bolso S1 adjacente ao sítio catalítico em ambas as enzimas. Ensaios de zimografia também foram realizados com o objetivo de elucidar possíveis contribuições da cadeia lateral e do núcleo catecólico do 4-NC na atividade inibitória frente às enzimas em estudo. O núcleo catecólico parece ser o responsável por sua atividade, pois o composto 1,2dimetoxibenzeno, que possui as hidroxilas bloqueadas por grupos metil, não foi capaz de exercer atividade inibitória significante frente à MMP-2 e MMP-9. Estudos de voltametria reforçaram a hipótese de que o 4-NC tem a capacidade de quelar os íons zinco presentes no tampão de incubação. / Matrix metalloproteinases (MMP) enzymes are overexpressed in almost all human tumors, and MMP-2 and MMP-9 subtypes have been associated with metastatic potential and poor prognosis in malignant tumors, such as metastatic melanoma and glioma. Compounds capable of inhibiting the activity of theses enzymes would be considered as potential therapeutic agents. The 4-nerolidylcatechol compound (4-NC), isolated from plants of genus Pothomorphe, has showed promising results in the treatment of melanoma and glioma, and was able to act in several important biochemical steps involved in the progression of these diseases, as well as inhibiting MMP-2 and MMP-9. However, the 4-NC mechanism of action is not completely understood. This study has involved the application of molecular modeling methods and formalisms of computer-aided molecular design (CAMD) in order to explore the interaction between 4-NC and MMP-2/MMP-9, and to design new inhibitors for these targets. Exploratory data analysis, which comprises hierarchical cluster analysis and principal components analysis, was performed to a set of hydroxamates (N=64). previously reported as MMP-2 and MMP-9 inhibitors, in order lo identify the molecular properties that is most critical for the discrimination process regarding the investigated compounds. The thermodynamic, electronic, and steric properties were: quite important to describe the highly active compounds in the data set of MMP-2, whereas the apparent partition coefficient (ClogD) at pH 1.5 was the property more relevant for MMP-9 data set. The presence of bulky substituents on the R3 moiety seems to be crucial for this set of inhibitors due to the molecular interaction with the S1 subsite of both enzymes. However, there is a limit regarding the substituents volume in this region. Receptor independent (RI) 4D-QSAR analysis was applied lo the same data set and it was possible to establish the pharmacophore mapping, besides to explore different alignments in order to generate the hypothesized bioactive conformation through the best QSAR model. The QSAR models have presented good predictability, assisted in proposing new inhibitors, and estimated the activity of 4-NC. Regarding the best QSAR model for MMP-9 (N=64), the 4-NC predicted activity was classified in the range of the moderate active inhibitors. The best QSAR model obtained for MMP-2 (N=38), however was not able to properly predict the activity for compounds with different chemical scaffold from those used to build up the QSAR model. Molecular docking studies have been developed to investigate the 4-NC binding mode into the catalytic site of the two enzymes and the interactions that could be established in those complexes. The results have shown two favorable conformers regarding the MMP inhibition. Molecular dynamics computational simulation were combined to molecular docking studies in order to obtain more detailed and reliable information regarding the intermolecular interactions of each complex. The 4-NC molecule tends to accommodate the side chain in the S1 pocket adjacent to the catalytic site in both enzymes. Experimental zymography assays were also performed to elucidate the possible contribution of the side chain and the catechol core in the 4-NC inhibitory activity against the MMP-2 and MMP-9 enzymes. The catechol core seems to be responsible for its activity, since the 1,2 dimethoxybenzene compound, which has the hydroxyl blocked by a methyl group, was not able to exert any significant inhibition on enzymes. Voltametric assays confirmed the hypothesis that 4-NC chelates zinc ions present in the incubation buffer.

Orientador: José Roberto Ruggiero / Banca: Pedro Geraldo Pascutti / Banca: José Maria Pires / Banca: Alexandre Suman de Araújo / Banca: Sabrina Thais Broggio Costa / Resumo: Peptídeos antimicrobianos são moléculas biologicamente ativas que, geralmente, tem as membranas fosfolipídicas como alvo primário. Resultados de diferentes técnicas experimentais têm sugerido que esses peptídeos permeabilizam as membranas pela formação de poros. Parte dos peptídeos caracterizados apresentam especificidade de disrupção para membranas de bactérias, em detrimento das membranas dos hospedeiros. Essa característica tem atraído a atenção da comunidade científica internacional, porque indica que estas moléculas podem ser modelos para o desenvolvimento de novos antibióticos, portanto o entendimento do mecanismo de ação, ou seja, do mecanismo de formação de poro, tem extrema importância. Simulações por Dinâmica Molecular foram produzidas para investigarmos o impacto que peptídeos antimicrobianos da família Mastoparano tem sobre membranas lipídicas modelo. Dois cenários foram explorados: (i) de baixa concentração peptídeo/lipídeo, P/L=1/128, que consistia de simulações fine-grained das interações de um peptídeo com uma bicamada pura de 128 lipídeos aniônicos (POPG) ou zwiteriônicos (POPC); (ii) de alta concentração, P/L=1/21, que abordava as interações de seis peptídeos com uma bicamada mista de 128 lipídeos POPC/POPG (1/1) usando uma modelagem coarse-grained. Tomando o peptídeo MP1 como caso paradigmático, verificamos que em baixo P/L é possível sugerir que sua característica seletiva surge da capacidade de coordenar e perturbar maior número de lipídeos em membrana aniônica comparada à neutra. Essa capacidade fica acentuada nas simulações com membrana mista, onde a atração dos lipídeos aniônicos pelos peptídeos catiônicos guiou a separação local e a formação de domínios de lipídeos aniônicos, o que facilitou o afinamento local da membrana e a formação de poro transmembrânico. Esses achados ajudam a explicar como peptídeos / Abstract: Antimicrobial peptides are biologically active molecules that, usually, have the phospholipid membranes as a primary target. Results from different experimental techniques have suggested these peptides permeabilize membranes by the pore formation. Part of the characterized peptides have specificity of disruption for bacterial membranes, instead of host membrane. This feature has attracted the attention of the international scientific community, because it indicates that these molecules can be models for the development of novel antibiotics, so understanding the mechanism of action, ie, the mechanism of pore formation, is extremely important. Molecular dynamics simulations were performed to investigate the impact of antimicrobial peptides from the Mastoparano family have on model lipid membranes. Two scenarios were explored: (i) of low peptide/lipid concentration, P/L=1/128, which consisted of fine-grained simulations of the interactions of a peptide with a pure bilayer of 128 anionic (POPG) or zwitterionic (POPC) lipids; (ii) of high concentration, P/L=1/21, which addressed the interactions of six peptides with a mixed bilayer of 128 POPC/POPG (1/1) lipids, using a coarse-grained modeling. Taking the MP1 peptide as a paradigmatic case, we found that in low P/L is possible to suggest that its selective feature arises of its ability to coordinate and disturb large number of lipids in the anionic membrane compared to neutral one. This ability is accentuated in simulations with mixed membrane, where the attraction of the anionic lipids by the cationic peptides led to the local segregation and formation of POPG lipid domains, which facilitated the local thinning of the membrane and the formation of transmembrane pore. These findings help to explain how short peptides, such as MP1, are able of forming pores in a membrane whose thickness is larger than the length of the peptide / Doutor

Submitted by Rosivalda Pereira (mrs.pereira@ufma.br) on 2017-06-14T16:54:49Z No. of bitstreams: 1 AlbertoLopes.pdf: 7651509 bytes, checksum: 542872f26d9cc8bf6b0b65d8a7151d16 (MD5) / Made available in DSpace on 2017-06-14T16:54:49Z (GMT). No. of bitstreams: 1 AlbertoLopes.pdf: 7651509 bytes, checksum: 542872f26d9cc8bf6b0b65d8a7151d16 (MD5) Previous issue date: 2015-07-07 / Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPQ) / Fundação de Amparo à Pesquisa e ao Desenvolvimento Científico e Tecnológico do Maranhão (FAPEMA) / The tick Rhipicephalus microplus is the major cattle ectoparasite of the world accounting for losses of billions of dollars that directly affect the return of such livestock. Its control is difficult due to the resistance of ticks to all chemical bases commercially available miticides. Acaricidal activity of terpenes has been evaluated in several studies that show satisfactory results, indicating these compounds are potential sources of new acaricidal products. The aim of this work was to select terpenes with potential activity against enzyme acetylcholinesterase (AChE) from R. microplus. Properties of the molecular volume, geometric parameters and vibrational terpenes were obtained from quantum chemical calculations the density functional theory level. Bioinformatic methodologies were applied to study the interaction of terpenes identified in essential oils of Citrus spp. and Lippia spp. with three AChE R. microplus. Since there are no available experimental structures, models of the three AChE were generated by homology modeling and then refined by molecular dynamics simulations. Soon after, were studies of molecular docking to detect best energy conformation of interaction and molecular dynamics simulations of this complex were carried out to study the behavior of this interaction. Our results suggest that the known acaricide activity of carvacrol is associated with its interaction with AChEs, while the acaricide activity of thymol is not associated with inhibition of that enzyme. Also, as expected, showed an excellent interaction coumafos acaricide and reports the first record of interaction of AChE from R. microplus with gammamuuruleno and elemol terpenes, molecules with few studies and that now configure themselves as candidates potential new acaricidal products. / O carrapato Rhipicephalus microplus é o principal ectoparasita da bovinocultura mundial sendo responsável por perdas de bilhões de dólares que afetam diretamente o retorno de tal produção animal. Seu controle é difícil devido à resistência dos carrapatos a todas as bases químicas de acaricidas comercialmente disponíveis. A atividade acaricida de terpenos tem sido avaliada em vários estudos que mostram resultados satisfatórios, indicando que estes compostos são fontes potenciais de novos produtos acaricidas. O objetivo desse trabalho foi selecionar terpenos com potencial atividade sobre a enzima acetilcolinesterase (AChE) de R. microplus. As propriedades do volume molecular, os parâmetros geométricos e vibracionais de terpenos foram obtidos a partir de cálculos de química quântica no nível da teoria do funcional de densidade. Metodologias de bioinformática foram aplicadas para estudar a interação de terpenos identificados em óleos essenciais de Citrus spp. e Lippia spp. com as três AChE de R. microplus. Como não existem estruturas experimentais disponíveis, modelos das três AChE foram gerados por modelagem por homologia e em seguida refinadas por simulações de dinâmica molecular. Logo após, foram realizados estudos de docagem molecular para detectar a melhor conformação energética de interação e simulações de dinâmica molecular desse complexo foram realizadas para estudarmos o comportamento dessa interação. Os resultados sugerem que a conhecida atividade carrapaticida do carvacrol está associada com a sua interação com a AChE, enquanto que a atividade carrapaticida do timol não está associado com a inibição dessa mesma enzima. Além disso, como esperado, mostrou uma excelente interação do carrapaticida cumafós e relata o primeiro registro da interação da AChE de R. microplus com os terpenos gama-muuruleno e elemol, moléculas com poucos estudos e que a partir de agora configuram-se como candidatos potenciais para novos produtos acaricidas.

Submitted by Rosivalda Pereira (mrs.pereira@ufma.br) on 2017-06-14T17:55:12Z No. of bitstreams: 1 RosaHelenaSilva.pdf: 2716670 bytes, checksum: 07dee1e77d704e91281c33c31e8c4938 (MD5) / Made available in DSpace on 2017-06-14T17:55:12Z (GMT). No. of bitstreams: 1 RosaHelenaSilva.pdf: 2716670 bytes, checksum: 07dee1e77d704e91281c33c31e8c4938 (MD5) Previous issue date: 2016-10-19 / Borreria verticillata (L.) G. Mey species known as broom vassourinha has antibacterial, antimalarial, hepatoprotective, antioxidative, analgesic and antiinflammatory activities; however, its antinociceptive action still demands more thorough investigation. The present study was to assess the antinociceptive activity of B. verticillata crude hydroalcoholic extract (EHBv) and the ethyl acetate fraction (FAc) by means of in vivo and in silico studies. In vivo assessment included the paw edema test, the writhing test, the formalin test and the tail flick test. Wistar rats and Swiss mice were divided into 6 groups and given the following treatments oral: 0.9% NaCl control group (CTL), 10 mg/kg memantine (MEM), 10 mg/kg indomethacin (INDO), 500 mg/kg EHBv (EHBv 500), 25 mg/kg FAC (FAc 25), 50 mg/kg and FAc (FAC 50). EHBv, FAc 25 and 50 treatments exhibited anti-edematous and peripheral antinociceptive effects. For in silico assessment, compounds found in FAc were subjected to molecular docking, and the leading compound was selected for molecular dynamics (MD) simulations. Ursolic acid exhibited better affinity parameters with the enzyme COX-2 and the NMDA receptor subunits GluN1a and GluN2B on molecular docking. In MD simulations, AU exhibited highly frequent interactions with residues Arg120 and Glu524 in the COX-2 active site and NMDA, whereby it might prevent COX-2 and NMDA receptor activation. Treatment with ursolic acid 10mg / Kg (AU) showed peripheral and central antinoceceptivo effect. The antinociceptive effect of B. verticillata might be predominantly attributed to peripheral actions, including the participation of anti-inflammatory components. Ursolic acid is the main active component and seems to be a promising source of COX-2 inhibitors and NMDA receptor antagonists / Borreria verticillata (L.) G. Mey espécie conhecida como vassourinha apresenta atividade antibacteriana, antimalárica, hepatoprotetora, antioxidante, analgésica e anti-inflamatória, entretanto sua atividade antinociceptiva é pouco estudada. O objetivo deste trabalho foi avaliar atividade antinociceptiva do extrato hidroalcoólico bruto (EHBv) e fração acetato de etila (FAc) de B. verticillata realizando estudos in vivo e in silico. Para avaliação in vivo, foram utilizados os testes do edema de pata, contorções abdominais, formalina e tail flick. Ratos Wistar e camundongos Swiss foram tratados via oral e divididos em 6 grupos: controle-NaCl 0.9%(CTL), memantina 10 mg/Kg (MEM), indometacina 10 mg/Kg (INDO), EHBv 500 mg/kg (EHBv 500), FAc 25 mg/Kg (FAc 25), FAc 50 mg/Kg (FAc 50). O tratamento com EHBv 500, FAc 25 e 50 apresentou efeito antiedematogênico e antinociceptivo periférico. Para avaliação in silico os compostos identificados na FAc foram submetidos a docagem molecular, o melhor composto foi selecionado para simulações de dinâmica e testado in vivo molecular. O ácido ursólico apresentou melhores parâmetros de afinidade com COX-2, GluN1a e GluN2B durante a docagem molecular. Nas simulações por dinâmica molecular, o ácido ursólico apresentou alta frequência de contatos com Arg120 e Glu524 do local ativo da COX- 2 e com o domínio LBD da Glun1a e GluN2B podendo com isso, impedir a ativação da COX-2 e do receptor NMDA. O tratamento com ácido ursólico 10mg/Kg (AU) apresentou efeito antinoceceptivo periférico e central. Sugere-se que o efeito antinociceptivo periférico de B. verticillata pode ser atribuído predominantemente à ação de compostos com ação anti-inflamatória. O ácido ursólico é o principal composto ativo, sendo um composto promissor para o desenvolvimento de fármacos inibidores da COX-2 e antagonistas dos receptores NMDA.

<p>CYP2C9 is one of our major drug metabolising enzymes and belongs to the cytochrome P450 (CYP) super family. The aim of this thesis was to gain an understanding of the quantitative structure–activity relationships (QSAR) of CYP2C9 substrates and inhibitors. This information will be useful in predicting drug metabolism and the potential for drug–drug interactions. To achieve this, a well characterised data set of structurally diverse, competitive CYP2C9 inhibitors was identified in our laboratory. Several computational methodologies, many based on GRID molecular interaction fields, were applied or developed in order to handle issues such as compound alignment and bioactive conformer selection. First, a traditional 3D QSAR was carried out in GOLPE, generating a predictive model. In this model the selection of a bioactive conformer and alignment was based on docking in a homology model of CYP2C9. Secondly, we introduced the concept of alignment independent descriptors from ALMOND. These descriptors were used to generate quantitatively and qualitatively predictive models. We subsequently derived conformation independent descriptors from molecular interaction fields calculated in FlexGRID. This enabled the derivation of 3D QSAR models without taking into account the selection of an alignment or a bioactive conformer. A subsequent programming effort enabled the conversion of this model back to 3D aligned pharmacophores. Similar alignment independent descriptors were also used in the development of the software MetaSite® that predicts the site of metabolism for CYP2C9 ligands. Finally, as crystal information on this isoform emerged, the performance of molecular dynamics simulations and homology models and the flexibility of the protein were evaluated using statistical analyses.</p><p>These modelling efforts have resulted in detailed knowledge of the structural characteristics in ligand interactions with the cytochrome P450 2C9 isoform.</p>

<p>The aqueous solvation of metal ions is one of the long-standing and complex problems in chemistry, with implications for and applications in a broad range of biochemical and electrochemical systems, where water is the all-pervasive medium.</p><p>This thesis describes computer simulations of Al³⁺(<i>aq</i>), Fe²⁺(<i>aq</i>), Fe³⁺(<i>aq</i>) and Cu²⁺(<i>aq</i>). Various aspects of the solvation of these polyvalent metal ions in water are addressed, at different levels of theory, using Car-Parrinello molecular dynamics, classical molecular dynamics and quantum-mechanical cluster calculations. Polyvalent metal ions are particularly interesting because of their large influence on the solvent structure, dynamics and thermodynamics, as well as on the properties of the individual solvent molecules. Polyvalent metal ions in aqueous solution also constitute a challenging subject for computer simulations since a sophisticated interaction model is needed to incorporate the large many-body effects. </p><p>All the ion-water coordination figures in this thesis are octahedral, except in the Cu²⁺(<i>aq</i>) solution, where the ion is penta-coordinated with four equatorial neighbours in a plane and one axial neighbour located ~0.45 Å further out from the ion. The equatorial ion-water bonds have covalent character, while the axial water molecule is only electrostatically bound. For all the ions, the OD stretching frequencies of the first-shell water molecules are much more downshifted than in liquid water. In the case of Cu²⁺(<i>aq</i>), however, only the OD frequencies of the equatorial water molecules are downshifted with respect to bulk water whereas the OD frequencies of the axial water molecule are slightly upshifted. </p><p>Various limitations of the Car-Parrinello molecular dynamics simulations have been explored and compared, such as finite system-size effects and shortcomings in the electronic structure calculations. The Car-Parrinello simulations are found to give reasonable descriptions of the polyvalent metal ions in aqueous solution.</p>