Similarity Mapplet: Interactive Visualization of the Directory of Useful Decoys and ChEMBL in High Dimensional Chemical Spaces

An Internet portal accessible at www.gdb.unibe.ch has been set up to automatically generate color-coded similarity maps of the ChEMBL database in relation to up to two sets of active compounds taken from the enhanced Directory of Useful Decoys (eDUD), a random set of molecules, or up to two sets of user-defined reference molecules. These maps visualize the relationships between the selected compounds and ChEMBL in six different high dimensional chemical spaces, namely MQN (42-D molecular quantum numbers), SMIfp (34-D SMILES fingerprint), APfp (20-D shape fingerprint), Xfp (55-D pharmacophore fingerprint), Sfp (1024-bit substructure fingerprint), and ECfp4 (1024-bit extended connectivity fingerprint). The maps are supplied in form of Java based desktop applications called “similarity mapplets” allowing interactive content browsing and linked to a “Multifingerprint Browser for ChEMBL” (also accessible directly at www.gdb.unibe.ch) to perform nearest neighbor searches. One can obtain six similarity mapplets of ChEMBL relative to random reference compounds, 606 similarity mapplets relative to single eDUD active sets, 30 300 similarity mapplets relative to pairs of eDUD active sets, and any number of similarity mapplets relative to user-defined reference sets to help visualize the structural diversity of compound series in drug optimization projects and their relationship to other known bioactive compounds.

Multimodal integration of functional and structural brain connectivity. An application to the characterization of Mild Cognitive Impairment

La investigación para el conocimiento del cerebro es una ciencia joven, su inicio se remonta a Santiago Ramón y Cajal en 1888. Desde esta fecha a nuestro tiempo la neurociencia ha avanzado mucho en el desarrollo de técnicas que permiten su estudio. Desde la neurociencia cognitiva hoy se explican muchos modelos que nos permiten acercar a nuestro entendimiento a capacidades cognitivas complejas. Aun así hablamos de una ciencia casi en pañales que tiene un lago recorrido por delante. Una de las claves del éxito en los estudios de la función cerebral ha sido convertirse en una disciplina que combina conocimientos de diversas áreas: de la física, de las matemáticas, de la estadística y de la psicología. Esta es la razón por la que a lo largo de este trabajo se entremezclan conceptos de diferentes campos con el objetivo de avanzar en el conocimiento de un tema tan complejo como el que nos ocupa: el entendimiento de la mente humana. Concretamente, esta tesis ha estado dirigida a la integración multimodal de la magnetoencefalografía (MEG) y la resonancia magnética ponderada en difusión (dMRI). Estas técnicas son sensibles, respectivamente, a los campos magnéticos emitidos por las corrientes neuronales, y a la microestructura de la materia blanca cerebral. A lo largo de este trabajo hemos visto que la combinación de estas técnicas permiten descubrir sinergias estructurofuncionales en el procesamiento de la información en el cerebro sano y en el curso de patologías neurológicas. Más específicamente en este trabajo se ha estudiado la relación entre la conectividad funcional y estructural y en cómo fusionarlas. Para ello, se ha cuantificado la conectividad funcional mediante el estudio de la sincronización de fase o la correlación de amplitudes entre series temporales, de esta forma se ha conseguido un índice que mide la similitud entre grupos neuronales o regiones cerebrales. Adicionalmente, la cuantificación de la conectividad estructural a partir de imágenes de resonancia magnética ponderadas en difusión, ha permitido hallar índices de la integridad de materia blanca o de la fuerza de las conexiones estructurales entre regiones. Estas medidas fueron combinadas en los capítulos 3, 4 y 5 de este trabajo siguiendo tres aproximaciones que iban desde el nivel más bajo al más alto de integración. Finalmente se utilizó la información fusionada de MEG y dMRI para la caracterización de grupos de sujetos con deterioro cognitivo leve, la detección de esta patología resulta relevante en la identificación precoz de la enfermedad de Alzheimer. Esta tesis está dividida en seis capítulos. En el capítulos 1 se establece un contexto para la introducción de la connectómica dentro de los campos de la neuroimagen y la neurociencia. Posteriormente en este capítulo se describen los objetivos de la tesis, y los objetivos específicos de cada una de las publicaciones científicas que resultaron de este trabajo. En el capítulo 2 se describen los métodos para cada técnica que fue empleada: conectividad estructural, conectividad funcional en resting state, redes cerebrales complejas y teoría de grafos y finalmente se describe la condición de deterioro cognitivo leve y el estado actual en la búsqueda de nuevos biomarcadores diagnósticos. En los capítulos 3, 4 y 5 se han incluido los artículos científicos que fueron producidos a lo largo de esta tesis. Estos han sido incluidos en el formato de la revista en que fueron publicados, estando divididos en introducción, materiales y métodos, resultados y discusión. Todos los métodos que fueron empleados en los artículos están descritos en el capítulo 2 de la tesis. Finalmente, en el capítulo 6 se concluyen los resultados generales de la tesis y se discuten de forma específica los resultados de cada artículo. ABSTRACT In this thesis I apply concepts from mathematics, physics and statistics to the neurosciences. This field benefits from the collaborative work of multidisciplinary teams where physicians, psychologists, engineers and other specialists fight for a common well: the understanding of the brain. Research on this field is still in its early years, being its birth attributed to the neuronal theory of Santiago Ramo´n y Cajal in 1888. In more than one hundred years only a very little percentage of the brain functioning has been discovered, and still much more needs to be explored. Isolated techniques aim at unraveling the system that supports our cognition, nevertheless in order to provide solid evidence in such a field multimodal techniques have arisen, with them we will be able to improve current knowledge about human cognition. Here we focus on the multimodal integration of magnetoencephalography (MEG) and diffusion weighted magnetic resonance imaging. These techniques are sensitive to the magnetic fields emitted by the neuronal currents and to the white matter microstructure, respectively. The combination of such techniques could bring up evidences about structural-functional synergies in the brain information processing and which part of this synergy fails in specific neurological pathologies. In particular, we are interested in the relationship between functional and structural connectivity, and how two integrate this information. We quantify the functional connectivity by studying the phase synchronization or the amplitude correlation between time series obtained by MEG, and so we get an index indicating similarity between neuronal entities, i.e. brain regions. In addition we quantify structural connectivity by performing diffusion tensor estimation from the diffusion weighted images, thus obtaining an indicator of the integrity of the white matter or, if preferred, the strength of the structural connections between regions. These quantifications are then combined following three different approaches, from the lowest to the highest level of integration, in chapters 3, 4 and 5. We finally apply the fused information to the characterization or prediction of mild cognitive impairment, a clinical entity which is considered as an early step in the continuum pathological process of dementia. The dissertation is divided in six chapters. In chapter 1 I introduce connectomics within the fields of neuroimaging and neuroscience. Later in this chapter we describe the objectives of this thesis, and the specific objectives of each of the scientific publications that were produced as result of this work. In chapter 2 I describe the methods for each of the techniques that were employed, namely structural connectivity, resting state functional connectivity, complex brain networks and graph theory, and finally, I describe the clinical condition of mild cognitive impairment and the current state of the art in the search for early biomarkers. In chapters 3, 4 and 5 I have included the scientific publications that were generated along this work. They have been included in in their original format and they contain introduction, materials and methods, results and discussion. All methods that were employed in these papers have been described in chapter 2. Finally, in chapter 6 I summarize all the results from this thesis, both locally for each of the scientific publications and globally for the whole work.

Short DNA Fragments without Sequence Similarity Are Initiation Sites for Replication in the Chromosome of the Yeast Yarrowia lipolytica

We have previously shown that both a centromere (CEN) and a replication origin are necessary for plasmid maintenance in the yeast Yarrowia lipolytica (Vernis et al., 1997). Because of this requirement, only a small number of centromere-proximal replication origins have been isolated from Yarrowia. We used a CEN-based plasmid to obtain noncentromeric origins, and several new fragments, some unique and some repetitive sequences, were isolated. Some of them were analyzed by two-dimensional gel electrophoresis and correspond to actual sites of initiation (ORI) on the chromosome. We observed that a 125-bp fragment is sufficient for a functional ORI on plasmid, and that chromosomal origins moved to ectopic sites on the chromosome continue to act as initiation sites. These Yarrowia origins share an 8-bp motif, which is not essential for origin function on plasmids. The Yarrowia origins do not display any obvious common structural features, like bent DNA or DNA unwinding elements, generally present at or near eukaryotic replication origins. Y. lipolytica origins thus share features of those in the unicellular Saccharomyces cerevisiae and in multicellular eukaryotes: they are discrete and short genetic elements without sequence similarity.

Binding of bisubstrate analog promotes large structural changes in the unregulated catalytic trimer of aspartate transcarbamoylase: Implications for allosteric regulation

A central problem in understanding enzyme regulation is to define the conformational states that account for allosteric changes in catalytic activity. For Escherichia coli aspartate transcarbamoylase (ATCase; EC 2.1.3.2) the active, relaxed (R state) holoenzyme is generally assumed to be represented by the crystal structure of the complex of the holoenzyme with the bisubstrate analog N-phosphonacetyl-l-aspartate (PALA). It is unclear, however, which conformational differences between the unliganded, inactive, taut (T state) holoenzyme and the PALA complex are attributable to localized effects of inhibitor binding as contrasted to the allosteric transition. To define the conformational changes in the isolated, nonallosteric C trimer resulting from the binding of PALA, we determined the 1.95-Å resolution crystal structure of the C trimer–PALA complex. In contrast to the free C trimer, the PALA-bound trimer exhibits approximate threefold symmetry. Conformational changes in the C trimer upon PALA binding include ordering of two active site loops and closure of the hinge relating the N- and C-terminal domains. The C trimer–PALA structure closely resembles the liganded C subunits in the PALA-bound holoenzyme. This similarity suggests that the pronounced hinge closure and other changes promoted by PALA binding to the holoenzyme are stabilized by ligand binding. Consequently, the conformational changes attributable to the allosteric transition of the holoenzyme remain to be defined.

HyPaLib: a database of RNAs and RNA structural elements defined by hybrid patterns

The database, called HyPaLib (for Hybrid Pattern Library), contains annotated structural elements characteristic for certain classes of structural and/or functional RNAs. These elements are described in a language specifically designed for this purpose. The language allows convenient specification of hybrid patterns, i.e. motifs consisting of sequence features and structural elements together with sequence similarity and thermodynamic constraints. We are currently developing software tools that allow a user to search sequence databases for any pattern in HyPaLib, thus providing functionality which is similar to PROSITE, but dedicated to the more complex patterns in RNA sequences. HyPaLib is available at http://bibiserv.techfak.uni-bielefeld.de/HyPa/.

Treble clef finger—a functionally diverse zinc-binding structural motif

Detection of similarity is particularly difficult for small proteins and thus connections between many of them remain unnoticed. Structure and sequence analysis of several metal-binding proteins reveals unexpected similarities in structural domains classified as different protein folds in SCOP and suggests unification of seven folds that belong to two protein classes. The common motif, termed treble clef finger in this study, forms the protein structural core and is 25–45 residues long. The treble clef motif is assembled around the central zinc ion and consists of a zinc knuckle, loop, β-hairpin and an α-helix. The knuckle and the first turn of the helix each incorporate two zinc ligands. Treble clef domains constitute the core of many structures such as ribosomal proteins L24E and S14, RING fingers, protein kinase cysteine-rich domains, nuclear receptor-like fingers, LIM domains, phosphatidylinositol-3-phosphate-binding domains and His-Me finger endonucleases. The treble clef finger is a uniquely versatile motif adaptable for various functions. This small domain with a 25 residue structural core can accommodate eight different metal-binding sites and can have many types of functions from binding of nucleic acids, proteins and small molecules, to catalysis of phosphodiester bond hydrolysis. Treble clef motifs are frequently incorporated in larger structures or occur in doublets. Present analysis suggests that the treble clef motif defines a distinct structural fold found in proteins with diverse functional properties and forms one of the major zinc finger groups.

Mimotope/anti-mimotope probing of structural relationships in platelet glycoprotein Ib alpha.

A bacteriophage library displaying random decapeptides was used to characterize the binding preference of C-34, a monoclonal antibody originally raised against platelet-type von Willebrand disease platelets heterozygous for the mutation 23OWKQ (G --> V)233V234 in the alpha chain of glycoprotein Ib (GPIb alpha). Three rounds of biopanning C-34 against the library resulted in striking convergence upon the sequence WNWRYREYV. Since no portion of this sequence corresponds to a recognizable peptide sequence within human platelet GPIb alpha, it may be considered a "mimotope" of the naturally occurring C-34 epitope, presumably bearing similarity to it in three-dimensional structure. Synthetic AWNWRYREYV peptide preincubated with C-34 fully neutralized the ability of C-34 to inhibit platelet aggregation, with an IC50 of approximately 6 microg/ml. When biotinylated AWNWRYREYV was subsequently bioparmed against the original decapeptide library, the sole clone demonstrating inhibitory activity above background level in a functional platelet assay displayed the sequence RHVAWWRQGV, and chemically synthesized peptide fully inhibited ristocetin-induced aggregation, with an IC50 of 200-400 microg/ml. Synthesized RHVAWWKQGV peptide exerted only slight inhibition, whereas RHVAWWKQVV peptide showed potent inhibitory activity. Moreover, whereas synthesized wild-type 228YVWKQGVDVK237 GPIb alpha peptide was virtually without inhibitory activity, the 228YVWKQ(G -->V) 233VDVK237 peptide fully inhibited ristocetin-induced aggregation, with an IC50 of approximately 400 microg/ml. These studies raise the possibility of an intramolecular association of peptide regions within GPIb alpha that may play a role in the regulation of von Willebrand factor-dependent platelet aggregation.

Transcriptional activation by protein-induced DNA bending: evidence for a DNA structural transmission model.

Integration host factor (IHF) is a DNA-bending protein that binds to an upstream activating sequence (UAS1) and, on a negatively supercoiled DNA template, activates transcription from the ilvPG promoter of the ilvG-MEDA operon of Escherichia coli. The transcriptional initiation site of the ilvGMEDA operon is located 92 bp downstream of UAS1. Activation is still observed when the orientation of the upstream IHF binding site is reversed. This manipulation places the IHF binding site on the opposite face of the DNA helix, directs the IHF-induced DNA bend in the opposite direction, and presents the opposite face of the nonsymmetrical, heterodimeric, IHF molecule to the downstream RNA polymerase. Lymphoid enhancer-binding factor, LEF-1, is a DNA-bending, lymphoid-specific, mammalian transcription factor that shares no amino acid sequence similarity with IHF. When the IHF site in UAS1 is replaced with a LEF-1 site, LEF-1 activates transcription from the downstream ilvPG promoter in E. coli as well as it is activated by its natural activator, IHF. These results suggest that specific interactions between IHF and RNA polymerase are not required for activation. The results of DNA structural studies show that IHF forms a protein-DNA complex in the UAS1 region that, in the absence of RNA polymerase, alters the structure of the DNA helix in the -10 hexanucleotide region of the downstream ilvPG promoter. The results of in vitro abortive transcription assays show that IIIF also increases the apparent rate of RNA polymerase isomerization from a closed to an open complex. We suggest, therefore, that IHF activates transcription by forming a higher-order protein-DNA complex in the UAS1 region that structurally alters the DNA helix in a way that facilitates open complex formation at the downstream ilvPG promoter site.

Three-dimensional structure of recombinant human osteogenic protein 1: structural paradigm for the transforming growth factor beta superfamily.

We report the three-dimensional structure of osteogenic protein 1 (OP-1, also known as bone morphogenetic protein 7) to 2.8-A resolution. OP-1 is a member of the transforming growth factor beta (TGF-beta) superfamily of proteins and is able to induce new bone formation in vivo. Members of this superfamily share sequence similarity in their C-terminal regions and are implicated in embryonic development and adult tissue repair. Our crystal structure makes possible the structural comparison between two members of the TGF-beta superfamily. We find that although there is limited sequence identity between OP-1 and TGF-beta 2, they share a common polypeptide fold. These results establish a basis for proposing the OP-1/TGF-beta 2 fold as the primary structural motif for the TGF-beta superfamily as a whole. Detailed comparison of the OP-1 and TGF-beta 2 structures has revealed striking differences that provide insights into how these growth factors interact with their receptors.

Structural variation of the pseudoautosomal region between and within inbred mouse strains.

The pseudoautosomal region (PAR) is a segment of shared homology between the sex chromosomes. Here we report additional probes for this region of the mouse genome. Genetic and fluorescence in situ hybridization analyses indicate that one probe, PAR-4, hybridizes to the pseudoautosomal telomere and a minor locus at the telomere of chromosome 9 and that a PCR assay based on the PAR-4 sequence amplifies only the pseudoautosomal locus (DXYHgu1). The region detected by PAR-4 is structurally unstable; it shows polymorphism both between mouse strains and between animals of the same inbred strain, which implies an unusually high mutation rate. Variation occurs in the region adjacent to a (TTAGGG)n array. Two pseudoautosomal probes can also hybridize to the distal telomeres of chromosomes 9 and 13, and all three telomeres contain DXYMov15. The similarity between these telomeres may reflect ancestral telomere-telomere exchange.

Bacterial cell division protein FtsZ assembles into protofilament sheets and minirings, structural homologs of tubulin polymers.

The bacterial cell division protein FtsZ is a homolog of tubulin, but it has not been determined whether FtsZ polymers are structurally related to the microtubule lattice. In the present study, we have obtained high-resolution electron micrographs of two FtsZ polymers that show remarkable similarity to tubulin polymers. The first is a two-dimensional sheet of protofilaments with a lattice very similar to that of the microtubule wall. The second is a miniring, consisting of a single protofilament in a sharply curved, planar conformation. FtsZ minirings are very similar to tubulin rings that are formed upon disassembly of microtubules but are about half the diameter. This suggests that the curved conformation occurs at every FtsZ subunit, but in tubulin rings the conformation occurs at either beta- or alpha-tubulin subunits but not both. We conclude that the functional polymer of FtsZ in bacterial cell division is a long thin sheet of protofilaments. There is sufficient FtsZ in Escherichia coli to form a protofilament that encircles the cell 20 times. The similarity of polymers formed by FtsZ and tubulin implies that the protofilament sheet is an ancient cytoskeletal system, originally functioning in bacterial cell division and later modified to make microtubules.

Modes of expression and common structural features of the complete phenylalanine ammonia-lyase gene family in parsley.

We describe a complete gene family encoding phenylalanine ammonia-lyase (PAL; EC 4.3.1.5) in one particular plant species. In parsley (Petroselinum crispum), the PAL gene family comprises two closely related members, PAL1 and PAL2, whose TATA-proximal promoter and coding regions are almost identical, and two additional members, PAL3 and PAL4, with less similarity to one another and to the PAL1 and PAL2 genes. Using gene-specific probes derived from the 5' untranslated regions of PAL1/2, PAL3, and PAL4, we determined the respective mRNA levels in parsley leaves and cell cultures treated with UV light or fungal elicitor and in wounded leaves and roots. For comparison, the functionally closely related cinnamate 4-hydroxylase (C4H) and 4-coumarate:CoA ligase (4CL) mRNAs were measured in parallel. The results indicate various degrees of differential responsiveness of PAL4 relative to the other PAL gene family members, in contrast to a high degree of coordination in the overall expression of the PAL, C4H, and 4CL genes. The only significant sequence similarities shared by all four PAL gene promoters are a TATA-proximal set of three putative cis-acting elements (boxes P, A, and L). None of these elements alone, or the promoter region containing all of them together, conferred elicitor or light responsiveness on a reporter gene in transient expression assays. The elements appear to be necessary but not sufficient for elicitor- or light-mediated PAL gene activation, similar to the situation previously reported for 4CL.

Structural interpretation of the mutations in the beta-cardiac myosin that have been implicated in familial hypertrophic cardiomyopathy.

In 10-30% of hypertrophic cardiomyopathy kindreds, the disease is caused by > 29 missense mutations in the cardiac beta-myosin heavy chain (MYH7) gene. The amino acid sequence similarity between chicken skeletal muscle and human beta-cardiac myosin and the three-dimensional structure of the chicken skeletal muscle myosin head have provided the opportunity to examine the structural consequences of these naturally occurring mutations in human beta-cardiac myosin. This study demonstrates that the mutations are related to distinct structural and functional domains. Twenty-four are clustered around four specific locations in the myosin head that are (i) associated with the actin binding interface, (ii) around the nucleotide binding site, (iii) adjacent to the region that connects the two reactive cysteine residues, and (iv) in close proximity to the interface of the heavy chain with the essential light chain. The remaining five mutations are in the myosin rod. The locations of these mutations provide insight into the way they impair the functioning of this molecular motor and also into the mechanism of energy transduction.

Structural Analysis with ANSYS on Stone Constructions in the Historical Spanish Heritage

According to the importance of rehabilitation and recovery of Architectural Heritage in the live of people, this paper is aimed to strengthen the traditional methods of stone vaults calculation taking advantage of the technological characteristics of the powerful program ANSYS Workbench. As an example of this, it could find out the possible pathologies that could arise during the construction history of the building. To limit this research, the upper vault of the main chapel of the Santiago parish church in Orihuela -Alicante- is selected as a reference which is a Jeronimo Quijano´s important building work in the XVI century in the Renaissance. Moreover, it is an innovative stone masonry vault that consists of 8 double intercrossed arches with each other and braced by severies. During the seventeenth century there was a lantern in the central cap and it is unknown why it was removed. Its construction could justify the original constructive solution with intercrossed arches that freed the center to create a more enlightened and comfortable presbytery. By similarity with other Quijano’s works, it is considered a small lantern drilling the central spherical cap. It is proposed to carry out a comparative study of it with different architectural solutions from the same period and based on several common parameters such as: a vault of square plant with spherical surround, intercrossed arches, a possible lantern, the dimension of the permitted space, similar states of loads and compact limestone masonry. The three solutions are mainly differentiated by their size and the type of lantern and its comparison lets us know which one is the most resistant and stable. The other two building works maintain some connection with the Quijano's professional scope. It has selected the particular case of the Communion chapel of the Basilica in Elche (a large prismatic lantern with a large cylindrical drum that starts from the own arches and an upper hemispherical dome), for its conservation, its proximity to Orihuela and its implementation during the century XVIII. Finally, a significant Dome Spanish Renaissance complete the selection: a cross vault of the Benavides Chapel of the Saint Francisco Convent in Baeza - Jaén-, designed by Andres of Vandelvira in the sixteenth century (a large hemispherical dome that starts from the own arcs). To simplify the calculation and standardize the work that have to be contrasted, all of them were considered with some similar characteristics: 30 cm constant thickness, the intercrossed arches were specifically analyzed and had identical loads, Young's modulus and Poisson's ratio. Regarding the calculation solutions, in general terms, the compressive stresses predominate, influencing on it the joint collaboration of the filling material on the vault, the vault itself, the thick side walls, the buttresses and the top cover weight . In addition, the three solutions are suitable, being the Orihuela one the safest and the Baeza one the riskiest for its large dimensions. Thus, the idea of intercrossed arches with suitable thickness would allow carry out the heaviest lantern and this would confirm it as a Renaissance architectural typology built in stone.

Functional recycling of C2 domains throughout evolution: A comparative study of synaptotagmin, protein kinase C and phospholipase C by sequence, structural and modelling approaches

The C2 domain is one of the most frequent and widely distributed calcium-binding motifs. Its structure comprises an eight-stranded beta-sandwich with two structural types as if the result of a circular permutation. Combining sequence, structural and modelling information, we have explored, at different levels of granularity, the functional characteristics of several families of C2 domains. At the coarsest level,the similarity correlates with key structural determinants of the C2 domain fold and, at the finest level, with the domain architecture of the proteins containing them, highlighting the functional diversity between the various subfamilies. The functional diversity appears as different conserved surface patches throughout this common fold. In some cases, these patches are related to substrate-binding sites whereas in others they correspond to interfaces of presumably permanent interaction between other domains within the same polypeptide chain. For those related to substrate-binding sites, the predictions overlap with biochemical data in addition to providing some novel observations. For those acting as protein-protein interfaces' our modelling analysis suggests that slight variations between families are a result of not only complementary adaptations in the interfaces involved but also different domain architecture. In the light of the sequence and structural genomic projects, the work presented here shows that modelling approaches along with careful sub-typing of protein families will be a powerful combination for a broader coverage in proteomics. (C) 2003 Elsevier Ltd. All rights reserved.