Avaliação funcional e estrutural de um novo peptídeo antimicrobiano do escorpião Tityus stigmurus

In Brazil, there is a high incidence of venomous animals. Among them, scorpions are highlighted by their medical importance, and for being their venom a source of several molecules with biological and pharmacological activity not yet fully understood, including several bioactive peptides. Antimicrobial peptides (AMPs) are components of the immune system in prokaryotes and eukaryotes, used in the first line of defense against microorganisms. In the present study, we characterized the first PAM previously identified through transcriptome of the venom gland of the scorpion Tityus stigmurus, named Stigmurin. The characteristics of Stigmurin were investigated by computational modeling and construction of dendrogram. In vitro tests investigated the antibacterial, antifungal, haemolytic and cytotoxic effects of crude venom and Stigmurin. In addition, the structural characteristics of Stigmurin were investigated by circular dochroism in water, 2, 2 , 2- trifluoethanol (TFE) and sodium dodecyl sulfate (SDS) and the models were refined by molecular dynamics simulations. The results showed that the selected sequence encodes a mature protein of 17 amino acid residues and the dendrogram reveals a case of convergent evolution. The crude venom showed no antimicrobial activity, however, Stigmurin exhibited a broad spectrum of antibacterial activity, with minimal inhibitory concentrations (MIC) ranging from 31.25 and 250 µg/mL for different strains, while the hemolytic activity at these concentrations was low. In cytotoxicity studies, the crude venom was unable to reduce cell viability in VERO E6 cells; in contrast, its activity in SiHa cells was significantly higher, corresponding to IC50 of 3.6 µg/mL. For Stigmurin the concentration sable to decrease cell viability of Vero E6 and SiHa cells in 50% were 275.67 µg/mL and 212.54 µg/mL, respectively. The dichroism spectra revealed the conformational flexibility, with predominating extended and β–sheet structures, as well as a remark able renaturation ability. The results suggest that Stigmurin could be considered as a potential antiinfective drug

Modelagem molecular voltada para a construção de um nanobiossensor para detecção do herbicida Imazaquin

In this study, our goal was develop and describe a molecular model of the enzyme-inhibiting interaction which can be used for an optimized projection of a Microscope Force Atomic nanobiosensor to detect pesticides molecules, used in agriculture, to evaluate its accordance with limit levels stipulated in valid legislation for its use. The studied herbicide (imazaquin) is a typical member of imidazolinone family and is an inhibitor of the enzymatic activity of Acetohydroxiacid Synthase (AHAS) enzyme that is responsible for the first step of pathway for the synthesis of side-chains in amino acids. The analysis of this enzyme property in the presence of its cofactors was made to obtain structural information and charge distribution of the molecular surface to evaluate its capacity of became immobilized on the Microscopy Atomic Force tip. The computational simulation of the system, using Molecular Dynamics, was possible with the force-field parameters for the cofactor and the herbicides obtained by the online tool SwissParam and it was implemented in force-field CHARMM27, used by software GROMACS; then appropriated simulations were made to validate the new parameters. The molecular orientation of the AHAS was defined based on electrostatic map and the availability of the herbicide in the active site. Steered Molecular Dynamics (SMD) Simulations, followed by quantum mechanics calculations for more representative frames, according to the sequential QM/MM methodology, in a specific direction of extraction of the herbicide from the active site. Therefore, external harmonic forces were applied with similar force constants of AFM cantilever for to simulate herbicide detection experiments by the proposed nanobiosensor. Force value of 1391 pN and binding energy of -14048.52 kJ mol-1 were calculated.

Stable algorithm for event detection in event-driven particle dynamics : logical states

Acknowledgments The authors gratefully acknowledge the support of the German Research Foundation (DFG) through the Cluster of Excellence ‘Engineering of Advanced Materials’ at the University of Erlangen-Nuremberg and through Grant Po 472/25.

An allosteric role for receptor activity-modifying proteins in defining GPCR pharmacology

G protein-coupled receptors are allosteric proteins that control transmission of external signals to regulate cellular response. Although agonist binding promotes canonical G protein signalling transmitted through conformational changes, G protein-coupled receptors also interact with other proteins. These include other G protein-coupled receptors, other receptors and channels, regulatory proteins and receptor-modifying proteins, notably receptor activity-modifying proteins (RAMPs). RAMPs have at least 11 G protein-coupled receptor partners, including many class B G protein-coupled receptors. Prototypic is the calcitonin receptor, with altered ligand specificity when co-expressed with RAMPs. To gain molecular insight into the consequences of this protein–protein interaction, we combined molecular modelling with mutagenesis of the calcitonin receptor extracellular domain, assessed in ligand binding and functional assays. Although some calcitonin receptor residues are universally important for peptide interactions (calcitonin, amylin and calcitonin gene-related peptide) in calcitonin receptor alone or with receptor activity-modifying protein, others have RAMP-dependent effects, whereby mutations decreased amylin/calcitonin gene-related peptide potency substantially only when RAMP was present. Remarkably, the key residues were completely conserved between calcitonin receptor and AMY receptors, and between subtypes of AMY receptor that have different ligand preferences. Mutations at the interface between calcitonin receptor and RAMP affected ligand pharmacology in a RAMP-dependent manner, suggesting that RAMP may allosterically influence the calcitonin receptor conformation. Supporting this, molecular dynamics simulations suggested that the calcitonin receptor extracellular N-terminal domain is more flexible in the presence of receptor activity-modifying protein 1. Thus, RAMPs may act in an allosteric manner to generate a spectrum of unique calcitonin receptor conformational states, explaining the pharmacological preferences of calcitonin receptor-RAMP complexes. This provides novel insight into our understanding of G protein-coupled receptor-protein interaction that is likely broadly applicable for this receptor class.

The (Un)Folding of Multidomain Proteins Through the Lens of Single-molecule Force-spectroscopy and Computer Simulation

Proteins are specialized molecules that catalyze most of the reactions that can sustain life, and they become functional by folding into a specific 3D structure. Despite their importance, the question, "how do proteins fold?" - first pondered in in the 1930's - is still listed as one of the top unanswered scientific questions as of 2005, according to the journal Science. Answering this question would provide a foundation for understanding protein function and would enable improved drug targeting, efficient biofuel production, and stronger biomaterials. Much of what we currently know about protein folding comes from studies on small, single-domain proteins, which may be quite different from the folding of large, multidomain proteins that predominate the proteomes of all organisms.

In this thesis I will discuss my work to fill this gap in understanding by studying the unfolding and refolding of large, multidomain proteins using the powerful combination of single-molecule force-spectroscopy experiments and molecular dynamic simulations.

The three model proteins studied - Luciferase, Protein S, and Streptavidin - lend insight into the inter-domain dependence for unfolding and the subdomain stabilization of binding ligands, and ultimately provide new insight into atomistic details of the intermediate states along the folding pathway.

How curvature-generating proteins build scaffolds on membrane nanotubes

Bin/Amphiphysin/Rvs (BAR) domain proteins control the curvature of lipid membranes in endocytosis, trafficking, cell motility, the formation of complex sub-cellular structures, and many other cellular phenomena. They form three-dimensional assemblies, which act as molecular scaffolds to reshape the membrane and alter its mechanical properties. It is unknown, however, how a protein scaffold forms and how BAR domains interact in these assemblies at protein densities relevant for a cell. In this work, we employ various experimental, theoretical and simulation approaches to explore how BAR proteins organize to form a scaffold on a membrane nanotube. By combining quantitative microscopy with analytical modeling, we demonstrate that a highly curving BAR protein endophilin nucleates its scaffolds at the ends of a membrane tube, contrary to a weaker curving protein centaurin, which binds evenly along the tube’s length. Our work implies that the nature of local protein-membrane interactions can affect the specific localization of proteins on membrane-remodeling sites. Furthermore, we show that amphipathic helices are dispensable in forming protein scaffolds. Finally, we explore a possible molecular structure of a BAR-domain scaffold using coarse-grained molecular dynamics simulations. Together with fluorescence microscopy, the simulations show that proteins need only to cover 30–40% of a tube’s surface to form a rigid assembly. Our work provides mechanical and structural insights into the way BAR proteins may sculpt the membrane as a high-order cooperative assembly in important biological processes. 

Insight into the mechanism of galactokinase: role of a critical glutamate residue and helix/coil transitions

Galactokinase, the enzyme which catalyses the first committed step in the Leloir pathway, has attracted interest due to its potential as a biocatalyst and as a possible drug target in the treatment of type I galactosemia. The mechanism of the enzyme is not fully elucidated. Molecular dynamics (MD) simulations of galactokinase with the active site residues Arg-37 and Asp-186 altered predicted that two regions (residues 174-179 and 231-240) had different dynamics as a consequence. Interestingly, the same two regions were also affected by alterations in Arg-105, Glu-174 and Arg- 228. These three residues were identified as important in catalysis in previous computational studies on human galactokinase. Alteration of Arg-105 to methionine resulted in a modest reduction in activity with little change in stability. When Arg-228 was changed to methionine, the enzyme’s interaction with both ATP and galactose was affected. This variant was significantly less stable than the wild-type protein. Changing Glu-174 to glutamine (but not to aspartate) resulted in no detectable activity and a less stable enzyme. Overall, these combined in silico and in vitro studies demonstrate the importance of a negative charge at position 174 and highlight the critical role of the dynamics in to key regions of the protein. We postulate that these regions may be critical for mediating the enzyme’s structure and function. 

Coarse Graining to Invesitigate Peptide-Lipid Interactions

Experimental characterization of molecular details is challenging, and although single molecule experiments have gained prominence, oligomer characterization remains largely unexplored. The ability to monitor the time evolution of individual molecules while they self assemble is essential in providing mechanistic insights about biological events. Molecular dynamics (MD) simulations can fill the gap in knowledge between single molecule experiments and ensemble studies like NMR, and are increasingly used to gain a better understanding of microscopic properties. Coarse-grained (CG) models aid in both exploring longer length and time scale molecular phenomena, and narrowing down the key interactions responsible for significant system characteristics. Over the past decade, CG techniques have made a significant impact in understanding physicochemical processes. However, the realm of peptide-lipid interfacial interactions, primarily binding, partitioning and folding of amphipathic peptides, remains largely unexplored compared to peptide folding in solution. The main drawback of existing CG models is the inability to capture environmentally sensitive changes in dipolar interactions, which are indigenous to protein folding, and lipid dynamics. We have used the Drude oscillator approach to incorporate structural polarization and dipolar interactions in CG beads to develop a minimalistic peptide model, WEPPROM (Water Explicit Polarizable PROtein Model), and a lipid model WEPMEM (Water Explicit Polarizable MEmbrane Model). The addition of backbone dipolar interactions in a CG model for peptides enabled us to achieve alpha-beta secondary structure content de novo, without any added bias. As a prelude to studying amphipathic peptide-lipid membrane interactions, the balance between hydrophobicity and backbone dipolar interactions in driving ordered peptide aggregation in water and at a hydrophobic-hydrophilic interface, was explored. We found that backbone dipole interactions play a crucial role in driving ordered peptide aggregation, both in water and at hydrophobic-hydrophilic interfaces; while hydrophobicity is more relevant for aggregation in water. A zwitterionic (POPC: 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and an anionic lipid (POPS: 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine) are used as model lipids for WEPMEM. The addition of head group dipolar interactions in lipids significantly improved structural, dynamic and dielectric properties of the model bilayer. Using WEPMEM and WEPPROM, we studied membrane-induced peptide folding of a cationic antimicrobial peptide with anticancer activity, SVS-1. We found that membrane-induced peptide folding is driven by both (a) cooperativity in peptide self interaction and (b) cooperativity in membrane-peptide interactions. The dipolar interactions between the peptide and the lipid head-groups contribute to stabilizing folded conformations. The role of monovalent ion size and peptide concentration in driving lipid domain formation in anionic/zwitterionic lipid mixtures was also investigated. Our study suggest monovalent ion size to be a crucial determinant of interaction with lipid head groups, and hence domain formation in lipid mixtures. This study reinforces the role of dipole interactions in protein folding, lipid membrane properties, membrane induced peptide folding and lipid domain formation. Therefore, the models developed in this thesis can be used to explore a multitude of biomolecular processes, both at longer time-scales and larger system sizes.

A study of electrochemical transport and diffuse charge dynamics using a Langevin equation

This thesis aims to develop new numerical and computational tools to study electrochemical transport and diffuse charge dynamics at small scales. Previous efforts at modeling electrokinetic phenomena at scales where the noncontinuum effects become significant have included continuum models based on the Poisson-Nernst-Planck equations and atomic simulations using molecular dynamics algorithms. Neither of them is easy to use or conducive to electrokinetic transport modeling in strong confinement or over long time scales. This work introduces a new approach based on a Langevin equation for diffuse charge dynamics in nanofluidic devices, which incorporates features from both continuum and atomistic methods. The model is then extended to include steric effects resulting from finite ion size, and applied to the phenomenon of double layer charging in a symmetric binary electrolyte between parallel-plate blocking electrodes, between which a voltage is applied. Finally, the results of this approach are compared to those of the continuum model based on the Poisson-Nernst-Planck equations.

A influência da incorporação de ácido α-eleosteárico, purificado e isolado do óleo de tungue, na dinâmica molecular de lipossomos

Este trabalho descreve o isolamento e purificação do ácido α-eleosteárico (α-ESA) a partir do óleo de tungue e sua caracterização por espectroscopia de infravermelho com transformada de Fourier (FTIR), cromatografia gasosa acoplada com espectrometria de massas (GC-MS) e espectroscopia de ressonância magnética nuclear (RMN) de 1H e 13C. O α-ESA apresenta atividades biológicas (antitumorais, anti-inflamatórias e antioxidantes), tornando-se importante compreender sua interação com membranas lipídicas. Assim, este trabalho também descreve resultados referentes ao efeito da incorporação de α-ESA na dinâmica molecular de lipossomos compostos por fosfatidilcolina. O sistema lipossomal puro e contendo α-ESA foi caracterizado através do uso de espectroscopia de UV-visível, FTIR, RMN e calorimetria de varredura diferencial (DSC). Como resultados da purificação do α-ESA, obtivemos uma pureza de 95,9% utilizando acetona como solvente de recristalização em detrimento dos 92,2% em solução etanólica. Na incorporação em lipossomos, observou-se uma maior interação do α-ESA com a parte polar, de interface e os primeiros metilenos da região apolar da fosfatidilcolina. Além disso, α-ESA apresentou um efeito de redução da fluidez de lipossomos. Os resultados contribuem para a geração de conhecimento para o desenvolvimento de novos sistemas farmacológicos.

On the temperature stability of extracellular hemoglobin of Glossoscolex paulistus, at different oxidation states: SAXS and DLS studies

Glossoscolex paulistus hemoglobin (HbGp) was studied by dynamic light scattering (DLS) and small angle X-ray scattering (SAXS). DLS melting curves were measured for met-HbGp at different concentrations. SAXS temperature studies were performed for oxy-, cyanomet- and met-HbGp forms, at several pH values. At pH 5.0 and 6.0, the scattering curves are identical from 20 to 60 degrees C, and R-g is 108 angstrom, independent of the oxidation form. At pH 7.0, protein denaturation and aggregation occurs above 55 degrees C and 60 degrees C, for oxy and met-HbGp, respectively. Cyanomet-HbGp, at pH 7.0, is stable up to 60 degrees C. At alkaline pH (8.0-9.0) and higher temperature, an irreversible dissociation process is observed, with a decrease of R-g, D-max and I(0). Analysis by p(r), obtained from GNOM, and OLIGOMER, was used to fit the SAXS experimental scattering curves by a combination of theoretical curves obtained for HbLt fragments from the crystal structure. Our results show clearly the increasing contribution of smaller molecular weight fragments, as a function of increasing pH and temperature, as well as, the order of thermal stabilities: cyanomet-> oxy- > met-HbGp. (C) 2012 Elsevier B.V. All rights reserved.

Fluorinated surfactants in solution: Diffusion coefficients of fluorinated alcohols in water

Intra-diffusion coefficients of three fluorinated alcohols, 2,2,3,3,3-pentafluoropropan-1-ol (PFP), 2,2,3,3,4,4,4-heptafluorobutan-1-ol (HFB) and 2,2,3,3,4,4,5,5,5-nonafluoropentan-1-ol (NFP) in water have been measured by the PFG–NMR spin-echo technique as a function of temperature and composition, focusing on the alcohol dilute region. For comparison, intra-diffusion coefficients of 2,2,2- trifluoroethanol (TFE) and HFB have also been measured in heavy water using the same method and conditions. As far as we know, these are the first experimental measurements of this property for these binary systems. Intra-diffusion coefficients for NFP in water and for TFE and HFB in heavy water have also been obtained by molecular dynamics simulation, complementing those for TFE, PFP and HFB reported in a previous work. The agreement between experimental and simulated results for PFP, HFB and NFP in water is reasonable, although presenting higher deviations than for the TFE/water system. From the dependence of the intra-diffusion coefficients on temperature, diffusion activation energies were estimated for all the solutes in water and heavy water.

Host-guest system of 4-nerolidylcatechol in 2-hydroxypropyl-beta-cyclodextrin: preparation, characterization and molecular modeling

The interaction of 4-nerolidylcatechol (4-NRC), a potent antioxidant agent, and 2-hydroxypropyl-beta-cyclodextrin (HP-beta-CD) was investigated by the solubility method using Fourier transform infrared (FTIR) methods in addition to UV-Vis, (1)H-nuclear magnetic resonance (NMR) spectroscopy and molecular modeling. The inclusion complexes were prepared using grinding, kneading and freeze-drying methods. According to phase solubility studies in water a B(S)-type diagram was found, displaying a stoichiometry complexation of 2:1 (drug:host) and stability constant of 6494 +/- A 837 M(-1). Stoichiometry was established by the UV spectrophotometer using Job's plot method and, also confirmed by molecular modeling. Data from (1)H-NMR, and FTIR, experiments also provided formation evidence of an inclusion complex between 4-NRC and HP-beta-CD. 4-NRC complexation indeed led to higher drug solubility and stability which could probably be useful to improve its biological properties and make it available to oral administration and topical formulations.

Computational and experimental study of the behavior of cyano-based ionic liquids in aqueous solution

The solvation of cyano- (CN-) based ionic liquids (ILs) and their capacity to establish hydrogen bonds (H-bonds) with water was studied by means of experimental and computational approaches. Experimentally, water activity data were measured for aqueous solutions of ILs based on 1-butyl-3-methylimidazolium ([BMIM](+)) cation combined with one of the following anions: thiocyanate ([SCN](-)), dicyanamide ([DCA](-)), or tricyanomethanide ([TCM](-)), and of 1-ethyl-3-methylimidazolium tetracyanoborate ([EMIM][TCB]). From the latter data, water activity coefficients were estimated showing that [BMIM][SCN] and [BMIM][DCA], unlike [BMIM][TCM] and [EMIM][TCB], are able to establish favorable interactions with water. Computationally, the conductor like screening model for real solvents (COSMO-RS) was used to estimate the water activity coefficients which compare well with the experimental ones. From the COSMO-RS results, it is suggested that the polarity of each ion composing the ILs has a strong effect on the solvation phenomena. Furthermore, classical molecular dynamics (MD) simulations were performed for obtaining an atomic level picture of the local molecular neighborhood of the different species. From the experimental and computational data it is showed that increasing the number of CN groups in the ILs' anions does not enhance their ability to establish H-bonds with water but decreases their polarities, being [BMIM][DCA] and [BMIM][SCN] the ones presenting higher propensity to interact.

Ação de anestésicos gerais em canais iônicos

Dissertação (mestrado)—Universidade de Brasília, Instituto de Ciências Biológicas, Departamento de Biologia Molecular, 2016.