Submitted by TABATA CRUZ MANOEL null (tabatacm@hotmail.com) on 2016-03-18T13:04:54Z No. of bitstreams: 1 ultima versao defesa.pdf: 3979971 bytes, checksum: 9e1ee9e6a385c2ba53c2341fe29ccecc (MD5) / Approved for entry into archive by Ana Paula Grisoto (grisotoana@reitoria.unesp.br) on 2016-03-21T19:29:22Z (GMT) No. of bitstreams: 1 manoel_tc_me_sjrp.pdf: 3979971 bytes, checksum: 9e1ee9e6a385c2ba53c2341fe29ccecc (MD5) / Made available in DSpace on 2016-03-21T19:29:22Z (GMT). No. of bitstreams: 1 manoel_tc_me_sjrp.pdf: 3979971 bytes, checksum: 9e1ee9e6a385c2ba53c2341fe29ccecc (MD5) Previous issue date: 2016-02-23 / Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) / As ciclodextrinas são um grupo de oligossacarídeos cíclicos, que apresentam um formato de cone truncado com cavidade hidrofóbica e exterior hidrofílico. Isto faz com que as ciclodextrinas possam encapsular fármacos, o que propicia melhoras na biodisponibilidade, estabilidade e proteção das moléculas encapsulada. Por esta característica, as ciclodextrinas constituem uma nova classe de excipientes farmacêuticos de grande utilidade, por formar complexos de inclusão reversíveis com moléculas apolares e atuarem como moléculas carregadoras de substancias de interesse biológico e baixa solubilidade. Dentre os compostos com ação farmacológica encontra-se os flavonoides, em destaque a molécula de quercetina, um fármaco natural com diversas propriedades biológicas, como anticâncer, antioxidante, antialérgica, porém devido sua baixa solubilidade torna-se necessária a formação de complexos de inclusão com moléculas carregadoras como a β-ciclodextrina. O presente trabalho busca investigar, por meio de simulações computacionais, as diferentes formas de interação dos flavonoides com a ciclodextrina e do complexo formado por bicamadas lipídicas, a fim de contribuir para um melhor entendimento do processo de formação com complexo de inclusão ciclodextrina/flavonoide e de entrega do fármaco à célula. / The cyclodextrins are cyclic oligosaccharides that present a structure in form of a truncated cone, with a hydrophobic inside surfaces and hydrophilic outside. This causes the cyclodextrin can encapsulate drugs provides improvements in bioavailability, stability and protection of encapsulated molecules. By this feature, the cyclodextrins are a new class of pharmaceutical excipients useful for forming inclusion complexes with reversible nonpolar molecules and act as carrier of molecules with high biological interest and low solubility. Among the compounds with pharmacological action it is the flavonoid, featured the quercetin molecule; a drug with natural biological proprieties such as anticancer, antioxidant, antialergic, however, due to its low solubility, it becomes necessary the formation of inclusion complexes with carrier molecules. This study aims to investigate, through computer simulation, different forms of interaction of flavonoids with cyclodextrin and the complex formed by lipid bilayers, in order to contribute to a better understanding of the formation process with inclusion complex cyclodextrin/flavonoids and the drug delivery to cell.

Submitted by Rafaela Magalhães Tavares null (rafaelarmt@yahoo.com.br) on 2016-04-19T11:40:30Z No. of bitstreams: 1 RAFAELA MAGALHAES TAVARES - DISSERTAÇÃO.pdf: 2642568 bytes, checksum: d735adb9c5a571ebaf0a5f7274658a8b (MD5) / Approved for entry into archive by Juliano Benedito Ferreira (julianoferreira@reitoria.unesp.br) on 2016-04-26T12:34:50Z (GMT) No. of bitstreams: 1 tavares_rm_me_sjrp.pdf: 2642568 bytes, checksum: d735adb9c5a571ebaf0a5f7274658a8b (MD5) / Made available in DSpace on 2016-04-26T12:34:50Z (GMT). No. of bitstreams: 1 tavares_rm_me_sjrp.pdf: 2642568 bytes, checksum: d735adb9c5a571ebaf0a5f7274658a8b (MD5) Previous issue date: 2016-03-01 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Nos últimos anos, o interesse por estudar peptídeos com atividade antimicrobiana e anticâncer tem aumentado consideravelmente. Neste trabalho, direcionamos nossos estudos para os peptídeos antimicrobianos da classe dos matoparanos extraídos do veneno da vespa social Polybia paulista: o Polybia-MP III (MP-III), Polybia-MP I (MP-I) e seu sintético análogo, o Asn2-Polybia-MP I (NMP-I). Com o objetivo de estudar a interação desses peptídeos com um modelo de membrana composta por lipídeos do tipo POPC (Palmitoil-Oleil-Fosfatidil-Colina) em solução aquosa, realizamos dois tipos de simulações de Dinâmica Molecular (DM). Primeiramente, foram realizadas simulações no equilíbrio com a finalidade de amostrar o comportamento geral do sistema. Em seguida, para investigarmos especificamente o processo de adsorção do peptídeo à bicamada, realizamos simulações de DM com a utilização do método Adaptive Biasing Force (ABF), o que nos permitiu calcular o perfil de energia livre desse processo. Dentre os três peptídeos estudados, o que mais se destacou com relação a sua interação com o modelo de membrana, na simulação no equilíbrio, foi o peptídeo MP-I, por ser o único a se adsorver na bicamada com 200ns de simulação. Com a utilização do método ABF, verificamos que a posição mais estável para cada peptídeo é a posição na qual estes estão paralelos à face da bicamada, com a face hidrofóbica de cada peptídeo voltada para o interior da membrana, e a face hidrofílica voltada para o meio aquoso. Além disso, os resíduos hidrofílicos estão em contato com o grupo polar dos fosfolipídeos e com a água, e o resíduo de triptofano encontra-se posicionado na interface hidrofóbica/hidrofílica. Os resultados obtidos na simulação no equilíbrio para o peptídeo MP-I, que foi o único a se adsorver na bicamada, estão de acordo com os resultados obtidos pelo método ABF, e ambos concordam com resultados experimentais da literatura. / In recent years, interest in studying peptides with antimicrobial and anti-cancer activity has increased considerably. In this study, we focus our studies on antimicrobial peptides of the mastoparans class extracted from the venom of the social wasp Polybia paulista: the Polybia-MP III (MP III), Polybia-MP I (MP-I) and its synthetic analogue, the Asn2-Polybia-MP I (NMP-I). In order to study the interaction of these peptides with a model of membrane composed of POPC (Palmitoyl-Oleoyl-Phosphatidyl-Choline) lipids type in aqueous solution, we conducted two types of molecular dynamics simulations (MD). At first, a balance system simulation was performed in order to get a sampling of the general behavior of the system. Then, to investigate the adsorption of the peptide to the bilayer, MD simulations using Adaptive Biasing Force (ABF) method was performed, which allowed us to calculate the free energy profile of this process. Among the three studied peptides, the MP-I peptide was the one that stood out related to its interaction with the membrane model, in the simulation on balance, for being the one to adsorb the bilayer with simulation of 200ns. Using the ABF method, we verified that the most stable position for each peptide is the position in which they are parallel to the surface of the bilayer with the hydrophobic side of each peptide facing into the membrane, and the hydrophilic side facing the aqueous solution. In addition, the hydrophilic residues are in contact with the polar group of the phospholipid and water, and the tryptophan residue is positioned on the hydrophobic/hydrophilic interface. The results of the simulation in balance for the MP-I peptide, which was the only adsorbed in the bilayer, are in agreement with the results obtained by the method ABF, and both agree with experimental results in the literature.

This work is concerned with sampling and computation of rare events in molecular systems. In particular, we present new methods for sampling the canonical ensemble corresponding to the Boltzmann-Gibbs probability measure. We combine an equation for controlling the kinetic energy of the system with a random noise to derive a highly degenerate diffusion (i.e. a diffusion equation where diffusion happens only along one or few degrees of freedom of the system). Next the concept of hypoellipticity is used to show that the corresponding Fokker-Planck equation of the highly degenerate diffusion is well-posed, hence we prove that the solution of the highly degenerate diffusion is ergodic with respect to the Boltzmann-Gibbs measure. We find that the new method is more efficient for computation of dynamical averages such as autocorrelation functions than the commonly used Langevin dynamics, especially in systems with many degrees of freedom. Finally we study the computation of free energy using an adaptive method which is based on the adaptive biasing force technique.

Currently the treatment of alcohol use disorder is very difficult and often requires the combination of therapy and medications, with many who undertake treatment experiencing relapse over time. There is also no treatment in use to prevent the development of alcohol use disorder. It is the aim of this work to provide information that may be useful for the development of a preventative treatment for developing alcohol use disorder by elucidating more of the acute effects of alcohol use. It is known that these effects originate in the brain. Within the brain are circuits made up of neurons that communicate with each other through chemical synapses. These chemical synapses involve the release of neurotransmitters from one neuron that are detected by another neuron, which initiates its own response. It is known that ethanol can change how much neurotransmitter is released from a neuron, depending on the specific neuron tested, and many researchers have implicated the "release machinery" as a target. It is also known that alcohol can affect lipid membrane properties that are important for the fusion of the vesicle membrane, encapsulating the neurotransmitter, with the cell membrane for release of the neurotransmitter outside of the neuron. It is not known if alcohol directly affects the SNARE proteins ("release machinery") or the lipid membranes to initiate the change in neurotransmitter release previously observed. Within this work you will find a discussion of the steps of neurotransmitter release and the known effects of anesthetics on components of this process, as an introduction to the topic (Chapters 1 and 2). In Chapters 3-5 you will find studies that successively dive deeper and deeper into the effects of alcohol on the SNARE proteins and lipid membranes. We show that ethanol is effective at a dose of 0.4% v/v or 64 mM at increasing fusion probability in a model of neurotransmitter release that uses the 3 SNARE proteins to drive fusion of a vesicle with a supported membrane. We also show that alcohol has little direct effect on the SNARE proteins themselves. In addition, we provide evidence that alcohol alters fusion oppositely, depending on which membrane leaflet it has most direct access to. In Chapter 5 we show that alcohol increases the probability of lipid tail protrusion in silico. Previously it has been shown that protrusion of one fatty acid tail of one lipid can initiate fusion of that membrane with an apposing membrane. These data provide further insight into the effects of alcohol on a neuron and we would argue are valuable to research pursuing treatment and prevention of alcohol use disorder.

No description available.

En physique statistique computationnelle, de bonnes techniques d'échantillonnage sont nécessaires pour obtenir des propriétés macroscopiques à travers des moyennes sur les états microscopiques. La principale difficulté est que ces états microscopiques sont généralement regroupés autour de configurations typiques, et un échantillonnage complet de l'espace configurationnel est donc typiquement très complexe à réaliser. Des techniques ont été proposées pour échantillonner efficacement les états microscopiques dans l'ensemble canonique. Un exemple important de quantités d'intérêt dans un tel cas est l'énergie libre. Le calcul d'énergie libre est très important dans les calculs de dynamique moléculaire, afin d'obtenir une description réduite d'un système physique complexe de grande dimension. La première partie de cette thèse est consacrée à une extension de la méthode adaptative de force biaisante classique (ABF), qui est utilisée pour calculer l'énergie libre associée à la mesure de Boltzmann-Gibbs et une coordonnée de réaction. Le problème de cette méthode est que le gradient approché de l'énergie libre, dit force moyenne, n'est pas un gradient en général. La contribution à ce domaine, présentée dans le chapitre 2, est de projeter la force moyenne estimée sur un gradient en utilisant la décomposition de Helmholtz. Dans la pratique, la nouvelle force gradient est obtenue à partir de la solution d'un problème de Poisson. En utilisant des techniques d'entropie, on étudie le comportement à la limite de l'équation de Fokker-Planck non linéaire associée au processus stochastique. On montre la convergence exponentielle vers l'équilibre de l'énergie libre estimée, avec un taux précis de convergence en fonction des constantes de l'inégalité de Sobolev logarithmiques des mesures canoniques conditionnelles à la coordonnée de réaction. L'intérêt de la méthode d'ABF projetée par rapport à l'approche originale ABF est que la variance de la nouvelle force moyenne est plus petite. On observe que cela implique une convergence plus rapide vers l'équilibre. En outre, la méthode permet d'avoir accès à une estimation de l'énergie libre en tout temps. La deuxième partie (voir le chapitre 3) est consacrée à étudier l'existence locale et globale, l'unicité et la régularité des solutions d'une équation non linéaire de Fokker-Planck associée à la méthode adaptative de force biaisante. Il s'agit d'un problème parabolique semilinéaire avec une non-linéarité non locale. L'équation de Fokker-Planck décrit l'évolution de la densité d'un processus stochastique associé à la méthode adaptative de force biaisante. Le terme non linéaire est non local et est utilisé lors de la simulation afin d'éliminer les caractéristiques métastables de la dynamique. Il est lié à une espérance conditionnelle qui définit la force biaisante. La preuve est basée sur des techniques de semi-groupe pour l'existence locale en temps, ainsi que sur une estimée a priori utilisant une sursolution pour montrer l'existence globale / In computational statistical physics, good sampling techniques are required to obtain macroscopic properties through averages over microscopic states. The main difficulty is that these microscopic states are typically clustered around typical configurations, and a complete sampling of the configurational space is thus in general very complex to achieve. Techniques have been proposed to efficiently sample the microscopic states in the canonical ensemble. An important example of quantities of interest in such a case is the free energy. Free energy computation techniques are very important in molecular dynamics computations, in order to obtain a coarse-grained description of a high-dimensional complex physical system. The first part of this thesis is dedicated to explore an extension of the classical adaptive biasing force (ABF) technique, which is used to compute the free energy associated to the Boltzmann-Gibbs measure and a reaction coordinate function. The problem of this method is that the approximated gradient of the free energy, called biasing force, is not a gradient. The contribution to this field, presented in Chapter 2, is to project the estimated biasing force on a gradient using the Helmholtz decomposition. In practice, the new gradient force is obtained by solving Poisson problem. Using entropy techniques, we study the longtime behavior of the nonlinear Fokker-Planck equation associated with the ABF process. We prove exponential convergence to equilibrium of the estimated free energy, with a precise rate of convergence in terms of the Logarithmic Sobolev inequality constants of the canonical measure conditioned to fixed values of the reaction coordinate. The interest of this projected ABF method compared to the original ABF approach is that the variance of the new biasing force is smaller, which yields quicker convergence to equilibrium. The second part, presented in Chapter 3, is dedicated to study local and global existence, uniqueness and regularity of the mild, Lp and classical solution of a nonlinear Fokker-Planck equation, arising in an adaptive biasing force method for molecular dynamics calculations. The partial differential equation is a semilinear parabolic initial boundary value problem with a nonlocal nonlinearity and periodic boundary conditions on the torus of dimension n, as presented in Chapter 3. The Fokker-Planck equation rules the evolution of the density of a given stochastic process that is a solution to Adaptive biasing force method. The nonlinear term is non local and is used during the simulation in order to remove the metastable features of the dynamics

Cette thèse comporte quatre fichiers vidéo. This thesis comes with four video files. / L'industrie chimique mondiale est en pleine mutation, cherchant des solutions pour rendre la synthèse organique classique plus durable. Une telle solution consiste à passer de la catalyse chimique classique à la biocatalyse. Bien que les avantages des enzymes incluent leur stéréo, régio et chimiosélectivité, cette sélectivité réduit souvent leur promiscuité. Les efforts requis pour adapter la fonction enzymatique aux réactions désirées se sont révélés d'une efficacité modérée, de sorte que des méthodes rapides et rentables sont nécessaires pour générer des biocatalyseurs qui rendront la production chimique plus efficace. Dans l’ère de la bioinformatique et des outils de calcul pour soutenir l'ingénierie des enzymes, le développement rapide de nouvelles fonctions enzymatiques devient une réalité. Cette thèse commence par un examen des développements récents sur les outils de calcul pour l’ingénierie des enzymes. Ceci est suivi par un exemple de l’ingénierie des enzymes purement expérimental ainsi que de l’évolution des protéines. Nous avons exploré l’espace mutationnel d'une enzyme primitive, la dihydrofolate réductase R67 (DHFR R67), en utilisant l’ingénierie semi-rationnelle des protéines. La conception rationnelle d’une librarie de mutants, ou «Smart library design», impliquait l’association covalente de monomères de l’homotétramère DHFR R67 en dimères afin d’augmenter la diversité de la librairie d’enzymes mutées. Le criblage par activité enzymatique a révélé un fort biais pour le maintien de la séquence native dans un des protomères tout en tolérant une variation de séquence élevée pour le deuxième. Il est plausible que les protomères natifs procurent l’activité observée, de sorte que nos efforts pour modifier le site actif de la DHFR R67 peuvent n’avoir été que modérément fructueux. Les limites des méthodes expérimentales sont ensuite abordées par le développement d’outils qui facilitent la prédiction des points chauds mutationnels, c’est-à-dire les sites privilégiés à muter afin de moduler la fonction. Le développement de ces techniques est intensif en termes de calcul, car les protéines sont de grandes molécules complexes dans un environnement à base d’eau, l’un des solvants les plus difficiles à modéliser. Nous présentons l’identification rapide des points chauds mutationnels spécifiques au substrat en utilisant l'exemple d’une enzyme cytochrome P450 industriellement pertinente, la CYP102A1. En appliquant la technique de simulation de la dynamique moléculaire par la force de polarisation adaptative, ou «ABF», nous confirmons les points chauds mutationnels connus pour l’hydroxylation des acides gras tout en identifiant de nouveaux points chauds mutationnels. Nous prédisons également la conformation du substrat naturel, l’acide palmitique, dans le site actif et nous appliquons ces connaissances pour effectuer un criblage virtuel d'autres substrats de cette enzyme. Nous effectuons ensuite des simulations de dynamique moléculaire pour traiter l’impact potentiel de la dynamique des protéines sur la catalyse enzymatique, qui est le sujet de discussions animées entre les experts du domaine. Avec la disponibilité accrue de structures cristallines dans la banque de données de protéines (PDB), il devient clair qu’une seule structure de protéine n’est pas suffisante pour élucider la fonction enzymatique. Nous le démontrons en analysant quatre structures cristallines que nous avons obtenues d’une enzyme β-lactamase, parmi lesquelles un réarrangement important des résidus clés du site actif est observable. Nous avons réalisé de longues simulations de dynamique moléculaire pour générer un ensemble de structures suggérant que les structures cristallines ne reflètent pas nécessairement la conformation de plus basse énergie. Enfin, nous étudions la nécessité de compléter de manière informatisée un hémisphère où l’expérimental n’est actuellement pas possible, à savoir la prédiction de la migration des gaz dans les enzymes. À titre d'exemple, la réactivité des enzymes cytochrome P450 dépend de la disponibilité des molécules d’oxygène envers l’hème du site actif. Par le biais de simulations de la dynamique moléculaire de type Simulation Implicite du Ligand (ILS), nous dérivons le paysage de l’énergie libre de petites molécules neutres de gaz pour cartographier les canaux potentiels empruntés par les gaz dans les cytochromes P450 : CYP102A1 et CYP102A5. La comparaison pour les gaz CO, N2 et O2 suggère que ces enzymes évoluent vers l’exclusion du CO inhibiteur. De plus, nous prédisons que les canaux empruntés par les gaz sont distincts des canaux empruntés par le substrat connu et que ces canaux peuvent donc être modifiés indépendamment les uns des autres. / The chemical industry worldwide is at a turning point, seeking solutions to make classical organic synthesis more sustainable. One such solution is to shift from classical catalysis to biocatalysis. Although the advantages of enzymes include their stereo-, regio-, and chemoselectivity, their selectivity often reduces versatility. Past efforts to tailor enzymatic function towards desired reactions have met with moderate effectiveness, such that fast and cost-effective methods are in demand to generate biocatalysts that will render the production of fine and bulk chemical production more benign. In the wake of bioinformatics and computational tools to support enzyme engineering, the fast development of new enzyme functions is becoming a reality. This thesis begins with a review of recent developments on computational tools for enzyme engineering. This is followed by an example of purely experimental enzyme engineering and protein evolution. We explored the mutational space of a primitive enzyme, the R67 dihydrofolate reductase (DHFR), using semi-rational protein engineering. ‘Smart library design’ involved fusing monomers of the homotetrameric R67 DHFR into dimers, to increase the diversity in the resulting mutated enzyme libraries. Activity-based screening revealed a strong bias for maintenance of the native sequence in one protomer with tolerance for high sequence variation in the second. It is plausible that the native protomers procure the observed activity, such that our efforts to modify the enzyme active site may have been only moderately fruitful. The limitations of experimental methods are then addressed by developing tools that facilitate computational mutational hotspot prediction. Developing these techniques is computationally intensive, as proteins are large molecular objects and work in aqueous media, one of the most complex solvents to model. We present the rapid, substrate-specific identification of mutational hotspots using the example of the industrially relevant P450 cytochrome CYP102A1. Applying the adaptive biasing force (ABF) molecular dynamics simulation technique, we confirm the known mutational hotspots for fatty acid hydroxylation and identify a new one. We also predict a catalytic binding pose for the natural substrate, palmitic acid, and apply that knowledge to perform virtual screening for further substrates for this enzyme. We then perform molecular dynamics simulations to address the potential impact of protein dynamics on enzyme catalysis, which is the topic of heated discussions among experts in the field. With the availability of more crystal structures in the Protein Data Bank, it is becoming clear that a single protein structure is not sufficient to elucidate enzyme function. We demonstrate this by analyzing four crystal structures we obtained of a β-lactamase enzyme, among which a striking rearrangement of key active site residues was observed. We performed long molecular dynamics simulations to generate a structural ensemble that suggests that crystal structures do not necessarily reflect the conformation of lowest energy. Finally, we address the need to computationally complement an area where experimentation is not currently possible, namely the prediction of gas migration into enzymes. As an example, the reactivity of P450 cytochrome enzymes depends on the availability of molecular oxygen at the active-site heme. Using the Implicit Ligand Sampling (ILS) molecular dynamics simulation technique, we derive the free energy landscape of small neutral gas molecules to map potential gas channels in cytochrome P450 CYP102A1 and CYP102A5. Comparison of CO, N2 and O2 suggests that those enzymes evolved towards exclusion of the inhibiting CO. In addition, we predict that gas channels are distinct from known substrate channels and therefore can be engineered independently from one another.