Protein pocket and ligand shape comparison and its application in virtual screening.

Understanding molecular recognition is one major requirement for drug discovery and design. Physicochemical and shape complementarity between two binding partners is the driving force during complex formation. In this study, the impact of shape within this process is analyzed. Protein binding pockets and co-crystallized ligands are represented by normalized principal moments of inertia ratios (NPRs). The corresponding descriptor space is triangular, with its corners occupied by spherical, discoid, and elongated shapes. An analysis of a selected set of sc-PDB complexes suggests that pockets and bound ligands avoid spherical shapes, which are, however, prevalent in small unoccupied pockets. Furthermore, a direct shape comparison confirms previous studies that on average only one third of a pocket is filled by its bound ligand, supplemented by a 50 % subpocket coverage. In this study, we found that shape complementary is expressed by low pairwise shape distances in NPR space, short distances between the centers-of-mass, and small deviations in the angle between the first principal ellipsoid axes. Furthermore, it is assessed how different binding pocket parameters are related to bioactivity and binding efficiency of the co-crystallized ligand. In addition, the performance of different shape and size parameters of pockets and ligands is evaluated in a virtual screening scenario performed on four representative targets.

Comparing assays for the measurement of AhR ligand activity

In this thesis (TFG) the results of the comparison of three assays for the measurement of AhR ligand activity are exposed. This study was part of a collaborative project aiming at the characterization of the AhR signaling activities of known naturally occurring compounds to explore the potential of using non-toxic compounds to treat inflammatory diseases via oral administration. The first goal of this project was to find an assay able to measure AhR-activity, so the comparison of different assays has been done in order to find the most convenient one according to the efficiency, sensitivity and precision. Moreover, other elements with operational nature such as price, toxicity of components or ease of use has been considered. From the use of compounds known from the literature to be AhR ligands, three assays have been tested: (1) P450-GloTM CYP1A2 Induction/Inhibition assay, (2) quantitative Polymerase Chain Reaction (qPCR) and (3) DR. CALUX® Bioassay. Moreover, a different experiment using the last assay was performed for the study in vivo of the transport of the compounds tested. The results of the TFG suggested the DR. CALUX® Bioassay as the most promising assay to be used for the screening of samples as AhR-ligands because it is quicker, easier to handle and less expensive than qPCR and more reproducible than the CYP1A2 Induction/Inhibition assay. Moreover, the use of this assay allowed having a first idea of which compounds are uptaken by the epithelial barrier and in with direction the transport happens.

Estimating the value of the metal-ligand bond dissociation enthalpy<D> (M-L) for adducts using empirical equations supported by TG data

In this work is presented and tested (for 106 adducts, mainly of the zinc group halides) two empirical equations supported in TG data to estimate the value of the metal-ligand bond dissociation enthalpy for adducts: <D> (M-O) = t i / g if t i < 420 K and <D> (M-O) = (t i / g ) - 7,75 . 10-2 . t i if t i > 420 K. In this empirical equations, t i is the thermodynamic temperature of the beginning of the thermal decomposition of the adduct, as determined by thermogravimetry, andg is a constant factor that is function of the metal halide considered and of the number of ligands, but is not dependant of the ligand itself. To half of the tested adducts the difference between experimental and calculated values was less than 5%. To about 80% of the tested adducts, the difference between the experimental (calorimetric) and the calculated (using the proposed equations) values are less than 15%.

Role of PheE15 gate in ligand entry and nitric oxide detoxification function of Mycobacterium tuberculosis truncated hemoglobin N

The truncated hemoglobin N, HbN, of Mycobacterium tuberculosis is endowed with a potent nitric oxide dioxygenase (NOD) activity that allows it to relieve nitrosative stress and enhance in vivo survival of its host. Despite its small size, the protein matrix of HbN hosts a two-branched tunnel, consisting of orthogonal short and long channels, that connects the heme active site to the protein surface. A novel dual-path mechanism has been suggested to drive migration of O(2) and NO to the distal heme cavity. While oxygen migrates mainly by the short path, a ligand-induced conformational change regulates opening of the long tunnel branch for NO, via a phenylalanine (PheE15) residue that acts as a gate. Site-directed mutagenesis and molecular simulations have been used to examine the gating role played by PheE15 in modulating the NOD function of HbN. Mutants carrying replacement of PheE15 with alanine, isoleucine, tyrosine and tryptophan have similar O(2)/CO association kinetics, but display significant reduction in their NOD function. Molecular simulations substantiated that mutation at the PheE15 gate confers significant changes in the long tunnel, and therefore may affect the migration of ligands. These results support the pivotal role of PheE15 gate in modulating the diffusion of NO via the long tunnel branch in the oxygenated protein, and hence the NOD function of HbN.

The saccharinate anion: a versatile and fascinating ligand in coordination chemistry

The saccharinate anion, obtained by deprotonation of the N-H moiety of saccharin (o-sulfobenzimide) is a very versatile and polyfunctional ligand in coordination chemistry. In this review the different forms of metal-to-ligand interactions involving this anion and some other coordination peculiarities are briefly discussed on the basis of some selected examples.

Building complexity in O2-binding copper complexes : site-selective metalation and intermolecular O2-binding at dicopper and heterometallic complexes derived from an unsymmetric ligand

A novel unsymmetric dinucleating ligand (LN3N4) combining a tridentate and a tetradentate binding sites linked through a m-xylyl spacer was synthesized as ligand scaffold for preparing homo- and dimetallic complexes, where the two metal ions are bound in two different coordination environments. Site-selective binding of different metal ions is demonstrated. LN3N4 is able to discriminate between CuI and a complementary metal (M′ = CuI, ZnII, FeII, CuII, or GaIII) so that pure heterodimetallic complexes with a general formula [CuIM′(LN3N4)]n+ are synthesized. Reaction of the dicopper(I) complex [CuI 2(LN3N4)]2+ with O2 leads to the formation of two different copper-dioxygen (Cu2O2) intermolecular species (O and TP) between two copper atoms located in the same site from different complex molecules. Taking advantage of this feature, reaction of the heterodimetallic complexes [CuM′(LN3N4)]n+ with O2 at low temperature is used as a tool to determine the final position of the CuI center in the system because only one of the two Cu2O2 species is formed

Electrocatalytic oxidation of methanol by the [Ru3O(OAc)6(py)2(CH3OH)]3+cluster: improving the metal-ligand electron transfer by accessing the higher oxidation states of a multicentered system

The [Ru3O(Ac)6(py)2(CH3OH)]+ cluster provides an effective electrocatalytic species for the oxidation of methanol under mild conditions. This complex exhibits characteristic electrochemical waves at -1.02, 0.15 and 1.18 V, associated with the Ru3III,II,II/Ru3III,III,II/Ru 3III,III,III /Ru3IV,III,III successive redox couples, respectively. Above 1.7 V, formation of two RuIV centers enhances the 2-electron oxidation of the methanol ligand yielding formaldehyde, in agreement with the theoretical evolution of the HOMO levels as a function of the oxidation states. This work illustrates an important strategy to improve the efficiency of the oxidation catalysis, by using a multicentered redox catalyst and accessing its multiple higher oxidation states.

Preconcentration of copper from natural water samples using ligand-less in situ surfactant-based solid phase extraction prior to FAAS determination

In the present work, a simple and rapid ligand-less, in situ, surfactant-based solid phase extraction for the preconcentration of copper in water samples was developed. In this method, a cationic surfactant (n-dodecyltrimethylammonium bromide) was dissolved in an aqueous sample followed by the addition of an appropriate ion-pairing agent (ClO4-). Due to the interaction between the surfactant and ion-pairing agent, solid particles were formed and subsequently used for the adsorption of Cu(OH)2 and CuI. After centrifugation, the sediment was dissolved in 1.0 mL of 1 mol L-1 HNO3 in ethanol and aspirated directly into the flame atomic absorption spectrometer. In order to obtain the optimum conditions, several parameters affecting the performance of the LL-ISS-SPE, including the volumes of DTAB, KClO4, and KI, pH, and potentially interfering ions, were optimized. It was found that KI and phosphate buffer solution (pH = 9) could extract more than 95% of copper ions. The amount of copper ions in the water samples varied from 3.2 to 4.8 ng mL-1, with relative standard deviations of 98.5%-103%. The determination of copper in water samples was linear over a concentration range of 0.5-200.0 ng mL-1. The limit of detection (3Sb/m) was 0.1 ng mL-1 with an enrichment factor of 38.7. The accuracy of the developed method was verified by the determination of copper in two certified reference materials, producing satisfactory results.

Study of ligand-based virtual screening tools in computer-aided drug design

Virtual screening is a central technique in drug discovery today. Millions of molecules can be tested in silico with the aim to only select the most promising and test them experimentally. The topic of this thesis is ligand-based virtual screening tools which take existing active molecules as starting point for finding new drug candidates. One goal of this thesis was to build a model that gives the probability that two molecules are biologically similar as function of one or more chemical similarity scores. Another important goal was to evaluate how well different ligand-based virtual screening tools are able to distinguish active molecules from inactives. One more criterion set for the virtual screening tools was their applicability in scaffold-hopping, i.e. finding new active chemotypes. In the first part of the work, a link was defined between the abstract chemical similarity score given by a screening tool and the probability that the two molecules are biologically similar. These results help to decide objectively which virtual screening hits to test experimentally. The work also resulted in a new type of data fusion method when using two or more tools. In the second part, five ligand-based virtual screening tools were evaluated and their performance was found to be generally poor. Three reasons for this were proposed: false negatives in the benchmark sets, active molecules that do not share the binding mode, and activity cliffs. In the third part of the study, a novel visualization and quantification method is presented for evaluation of the scaffold-hopping ability of virtual screening tools.

Conformational Regulation of alpha2beta1 Integrin in Ligand Binding

Integrins are heterodimeric cell adhesion receptors involved in cell-cell and cell-extracellular matrix (ECM) interactions. They transmit bidirectional signals across the cell membrane. This results in a wide range of biological events from cell differentiation to apoptosis. alpha2beta1 integrin is an abundant collagen receptor expressed on the surface of several cell types. In addition to ECM ligands, alpha2beta1 integrins are bound by echovirus 1 (EV1) which uses alpha2beta1 as a receptor to initiate its life cycle in the infected cell. The aim of this thesis project was to provide further insight into the mechanisms of alpha2beta1 integrin ligand recognition and receptor activation. Collagen fibrils are the principal tensile elements of the ECM. Yet, the interaction of alpha2beta1 integrin with the fibrillar form of collagen I has received relatively little attention. This research focused on the ability of alpha2beta1 integrin to act as a receptor for type I collagen fibrils. Also the molecular requirements of the EV1 interaction with alpha2beta1 were studied. Conventionally, ligand binding has been suggested to require integrin activation and the binding may further trigger integrin signalling. Another main objective of this study was to elucidate both the inside-out and outside-in signalling mechanisms of alpha2beta1 integrin in adherent cells. The results indicated that alpha2beta1 integrin is the principal integrin-type collagen receptor for type I collagen fibrils, and alpha2beta1 may participate in the regulation of pericellular collagen fibrillogenesis. Furthermore, alpha2beta1 integrin inside-out activation appeared to be synergistically regulated by integrin clustering and conformational activation. The triggering of alpha2beta1 integrin outside-in signalling, however, was shown to require both conformational changes and clustering. In contrast to ECM ligands, EV1 appeared to take advantage of the bent, inactive form of alpha2beta1 integrin in initiating its life cycle in the cell. This research together with other recent studies, has shed light on the molecular mechanisms of integrin activation. It is becoming evident that large ligands are able to bind to the bent form of integrin, which has been previously considered to be physiologically inactive. Consequently, our understanding of the conformational modulation of integrins upon activation is changing.

Electrochemical properties of Cu4[PhN3C6H4N3(H)Ph]4(µ-O)2, a tetranuclear Copper(II) complex with 1-phenyltriazenido-2-phenyltriazene-benzene as ligand

Bis-(µ2-oxo)-tetrakis{[1-feniltriazene-1,3-diil)-2-(phenyltriazenil)benzene copper(II) is a tetranuclear complex which shows four Cu(II) ions coordinated by four 1,2-bis(phenyltriazene)benzene bridged ligands, with one diazoaminic deprotonated chain, and two O2- ligands. The complex reduces at E1/2 = -0.95 V vs Fc+/Fc, a two electrons process. Cyclic voltammetric and spectroelectrochemical studies showed a reversible process. When immobilized on carbon paste electrode, the complex electrocatalyses the reduction of O2 dissolved on aqueous solution at -0.3 V vs SCE potential. The obtained current shows linearity with O2 concentration.

Alpha2-adrenoceptors: Structure and ligand binding properties at the molecular level

Alpha₂-Adrenoceptors: structure and ligand binding properties at the molecular level The mouse is the most frequently used animal model in biomedical research, but the use of zebrafish as a model organism to mimic human diseases is on the increase. Therefore it is considered important to understand their pharmacological differences from humans also at the molecular level. The zebrafish Alpha_2-adrenoceptors were expressed in mammalian cells and the binding affinities of 20 diverse ligands were determined and compared to the corresponding human receptors. The pharmacological properties of the human and zebrafish Alpha_2--adrenoceptors were found to be quite well conserved. Receptor models based on the crystal structures of bovine rhodopsin and the human Beta₂-adrenoceptor revealed that most structural differences between the paralogous and orthologous Alpha_2--adrenoceptors were located within the second extracellular loop (XL2). Reciprocal mutations were generated in the mouse and human Alpha_2--adrenoceptors. Ligand binding experiments revealed that substitutions in XL2 reversed the binding profiles of the human and mouse Alpha_2--adrenoceptors for yohimbine, rauwolscine and RS-79948-197, evidence for a role for XL2 in the determination of species-specific ligand binding. Previous mutagenesis studies had not been able to explain the subtype preference of several large Alpha_2--adrenoceptor antagonists. We prepared chimaeric Alpha_2--adrenoceptors where the first transmembrane (TM1) domain was exchanged between the three human Alpha_2--adrenoceptor subtypes. The binding affinities of spiperone, spiroxatrine and chlorpromazine were observed to be significantly improved by TM1 substitutions of the Alpha_2a--adrenoceptor. Docking simulations indicated that indirect effects, such as allosteric modulation, are more likely to be involved in this phenomenon rather than specific side-chain interactions between ligands and receptors.

Prediction of 3D structure and ligand-binding properties : ABC transporters and flavodiiron proteins from Synechocystis sp. strain PCC 6803

Cyanobacteria are unicellular, non-nitrogen-fixing prokaryotes, which perform photosynthesis similarly as higher plants. The cyanobacterium Synechocystis sp. strain PCC 6803 is used as a model organism in photosynthesis research. My research described herein aims at understanding the function of the photosynthetic machinery and how it responds to changes in the environment. Detailed knowledge of the regulation of photosynthesis in cyanobacteria can be utilized for biotechnological purposes, for example in the harnessing of solar energy for biofuel production. In photosynthesis, iron participates in electron transfer. Here, we focused on iron transport in Synechocystis sp. strain PCC 6803 and particularly on the environmental regulation of the genes encoding the FutA2BC ferric iron transporter, which belongs to the ABC transporter family. A homology model built for the ATP-binding subunit FutC indicates that it has a functional ATPbinding site as well as conserved interactions with the channel-forming subunit FutB in the transporter complex. Polyamines are important for the cell proliferation, differentiation and apoptosis in prokaryotic and eukaryotic cells. In plants, polyamines have special roles in stress response and in plant survival. The polyamine metabolism in cyanobacteria in response to environmental stress is of interest in research on stress tolerance of higher plants. In this thesis, the potd gene encoding an polyamine transporter subunit from Synechocystis sp. strain PCC 6803 was characterized for the first time. A homology model built for PotD protein indicated that it has capability of binding polyamines, with the preference for spermidine. Furthermore, in order to investigate the structural features of the substrate specificity, polyamines were docked into the binding site. Spermidine was positioned very similarly in Synechocystis PotD as in the template structure and had most favorable interactions of the docked polyamines. Based on the homology model, experimental work was conducted, which confirmed the binding preference. Flavodiiron proteins (Flv) are enzymes, which protect the cell against toxicity of oxygen and/or nitric oxide by reduction. In this thesis, we present a novel type of photoprotection mechanism in cyanobacteria by the heterodimer of Flv2/Flv4. The constructed homology model of Flv2/Flv4 suggests a functional heterodimer capable of rapid electron transfer. The unknown protein sll0218, encoded by the flv2-flv4 operon, is assumed to facilitate the interaction of the Flv2/Flv4 heterodimer and energy transfer between the phycobilisome and PSII. Flv2/Flv4 provides an alternative electron transfer pathway and functions as an electron sink in PSII electron transfer.

Probabilistic contact preferences in protein-ligand and protein-protein complexes

Modeller för intermolekulär växelvärkan utnyttjas brett inom biologin. Analys av kontakter mellan proteiner och läkemedelsforskning representerar typiska tillämpningsområden för dylika modeller. En modell som beskriver sådana molekylära växelverkningar kan utformas med hjälp av biofysisk teori, vilket tenderar att resultera i ytterst tung beräkningsbörda även för enkla tillämpningar. Ett alternativt sätt att formulera modeller är att utnyttja stora databaser som innehåller strukturmätningar gjorda med hjälp av till exempel röntgendiffraktion. Då man använder sig av empiriska mätdata direkt, möjliggör en statistisk modell att osäkerheten och inexaktheten i datat tas till hänsyn på ett adekvat sätt, samtidigt som beräkningsbördan håller sig på en rimligare nivå jämfört med kvantmekaniska metoder som i princip borde ge de optimala resultaten. I avhandlingen utvecklades en 3D modell för numerisk undersökning av intermolekulär växelverkan baserad på Bayesiansk statistik. Modellens syfte är att åstadkomma prognoser för det hurdana eller vilka molekylstrukturer prefereras i en given kontext, d.v.s. är mer sannolika inom ramen för interaktion. Modellen testades i essentiella molekyläromgivningar - en liten molekyl vid sin bindningsplats hos ett protein och en gränsyta mellan proteinerna i ett komplex. De erhållna numeriska resultaten motsvarar väl experimentella resultat som tidigare rapporterats i litteraturen, exempelvis kvalitativa bindningsaffiniteter och kemisk kännedom av vissa aminosyrors rumsliga förmågor att utgöra bindningar. I avhandlingen gjordes ytterligare preliminära tester av den statistiska ansatsen för modellering av den centrala molekylära strukturella anpassningsbarheten. I praktiken är den utvecklade modellen ämnad som ett led i en mer omfattande analysmetod, så som en s.k. farmakofor modell. Molekyylivuorovaikutusten mallintamista hyödynnetään laajasti biologisten kysymysten tarkastelussa. Tyypillisiä esimerkkejä sovelluskohteista ovat proteiinien väliset kontaktit ja lääkesuunnittelu. Vuorovaikutuksia kuvaavan mallin lähtökohta voi olla molekyyleihin liittyvä teoria, jolloin soveltamiseen liittyvä laskenta saattaa olla erityisen raskasta, tai suuri havaintojoukko joka on saatu aikaan esimerkiksi mittaamalla rakenteita röntgendiffraktio menetelmällä. Tilastollinen malli mahdollistaa havaintoaineistossa olevan epätarkkuuden ja epävarmuuden huomioimisen, samalla pitäen laskennallisen kuorman pienempänä verrattuna periaatteessa parhaan tuloksen antavaan kvanttimekaaniseen mallinnukseen. Väitöstyössä kehitettiin bayesiläiseen tilastotieteeseen perustuva 3D malli molekyylien välisten vuorovaikutusten laskennalliseen tarkasteluun. Mallin tehtävä on tuottaa ennusteita sen suhteen, minkä tai millaisten molekyylirakenteiden väliset kompleksit ovat etusijalla, toisin sanoen todennäköisempiä, vuorovaikutustilanteessa. Työssä kehitetyn menetelmän toimivuutta testattiin käyttötarkoituksen suhteen olennaisissa molekyyliympäristöissä - pieni molekyyli sitoutumiskohdassaan proteiinissa sekä rajapinta kahden proteiinin välilllä proteiinikompleksissa. Saadut laskennalliset tulokset vastasivat hyvin vertailuun käytettyjä kirjallisuudesta saatuja kokeellisia tuloksia, kuten laadullisia sitoutumisaffiniteetteja, sekä kemiallista tietoa esimerkiksi tiettyjen aminohappojen avaruudellisesta sidoksenmuodostuksesta. Väitöstyössä myös alustavasti testattiin tilastollista lähestymistapaa tärkeän molekyylien rakenteellisen mukautuvuuden mallintamiseen. Käytännössä malli on tarkoitettu osaksi jotakin laajempaa analyysimenetelmää, kuten farmakoforimallia.

Structure and function of the selectin ligand PSGL-1

P-selectin glycoprotein ligand-1 (PSGL-1) is a dimeric mucin-like 120-kDa glycoprotein on leukocyte surfaces that binds to P- and L-selectin and promotes cell adhesion in the inflammatory response. The extreme amino terminal extracellular domain of PSGL-1 is critical for these interactions, based on site-directed mutagenesis, blocking monoclonal antibodies, and biochemical analyses. The current hypothesis is that for high affinity interactions with P-selectin, PSGL-1 must contain O-glycans with a core-2 branched motif containing the sialyl Lewis x antigen (NeuAca2®3Galß1®4[Fuca1®3]GlcNAcß1®R). In addition, high affinity interactions require the co-expression of tyrosine sulfate on tyrosine residues near the critical O-glycan structure. This review addresses the biochemical evidence for this hypothesis and the evidence that PSGL-1 is an important in vivo ligand for cell adhesion.