Information services of the European Bioinformatics Institute

The Bioinformatics Program at Washington University's Bernard Becker Medical Library: making it happen

The need for a Bioinformatics Program became apparent after repeated requests from researchers for bioinformatics resources and the University's announcement of BioMed 21, a 300 million dollars biomedical science initiative to bring "new knowledge of the human genetic blueprint to the patient's bedside and change how illnesses are understood, diagnosed, and successfully treated." A Bioinformatics Specialist with a strong background in the fields of Molecular Biology and Biostatistics was hired. Program development started in 2003. Initially, three core courses were developed: Sequence Similarity Search, Genetic Variation, and Human Genome Resources. Each of these courses was offered twice per semester. Additionally, partnerships were established within the University Medical School Community (e.g., the Genome Sequencing Center) to license and teach Spotfire's DecisionSite for Functional Genomics, a software package used to analyze microarray data. From March to May 2005, seven Spotfire classes were taught. Each Spotfire class consisted of seven hours of classroom work. Also, in-depth consultations were scheduled with faculty and researchers to address their specific needs. These consultations led to requests for other software packages to purchase and manage, including Lasergene and CSD (Cambridge Crystal Structure Database). Efforts to reach outside of the University Community were made through the development of a Bioinformatics Web site.

Allosteric determinants in guanine nucleotide-binding proteins

Members of the G protein superfamily contain nucleotide-dependent switches that dictate the specificity of their interactions with binding partners. Using a sequence-based method termed statistical coupling analysis (SCA), we have attempted to identify the allosteric core of these proteins, the network of amino acid residues that couples the domains responsible for nucleotide binding and protein-protein interactions. One-third of the 38 residues identified by SCA were mutated in the G protein Gs alpha, and the interactions of guanosine 5'-3-O-(thio)triphosphate- and GDP-bound mutant proteins were tested with both adenylyl cyclase (preferential binding to GTP-Gs alpha) and the G protein beta gamma subunit complex (preferential binding to GDP-Gs alpha). A two-state allosteric model predicts that mutation of residues that control the equilibrium between GDP- and GTP-bound conformations of the protein will cause the ratio of affinities of these species for adenylyl cyclase and G beta gamma to vary in a reciprocal fashion. Observed results were consistent with this prediction. The network of residues identified by the SCA appears to comprise a core allosteric mechanism conferring nucleotide-dependent switching; the specific features of different G protein family members are built on this core.

Molecular modeling of the amyloid-beta-peptide using the homology to a fragment of triosephosphate isomerase that forms amyloid in vitro

The main component of the amyloid senile plaques found in Alzheimer's brain is the amyloid-beta-peptide (A beta), a proteolytic product of a membrane precursor protein. Previous structural studies have found different conformations for the A beta peptide depending on the solvent and pH used. In general, they have suggested an alpha-helix conformation at the N-terminal domain and a beta-sheet conformation for the C-terminal domain. The structure of the complete A beta peptide (residues 1-40) solved by NMR has revealed that only helical structure is present in A beta. However, this result cannot explain the large beta-sheet A beta aggregates known to form amyloid under physiological conditions. Therefore, we investigated the structure of A beta by molecular modeling based on extensive homology using the Smith and Waterman algorithm implemented in the MPsrch program (Blitz server). The results showed a mean value of 23% identity with selected sequences. Since these values do not allow a clear homology to be established with a reference structure in order to perform molecular modeling studies, we searched for detailed homology. A 28% identity with an alpha/beta segment of a triosephosphate isomerase (TIM) from Culex tarralis with an unsolved three-dimensional structure was obtained. Then, multiple sequence alignment was performed considering A beta, TIM from C.tarralis and another five TIM sequences with known three-dimensional structures. We found a TIM segment with secondary structure elements in agreement with previous experimental data for A beta. Moreover, when a synthetic peptide from this TIM segment was studied in vitro, it was able to aggregate and to form amyloid fibrils, as established by Congo red binding and electron microscopy. The A beta model obtained was optimized by molecular dynamics considering ionizable side chains in order to simulate A beta in a neutral pH environment. We report here the structural implications of this study.

MIPS: a database for protein sequences, homology data and yeast genome information

The MIPS group (Martinsried Institute for Protein Sequences) at the Max-Planck-Institute for Biochemistry, Martinsried near Munich, Germany, collects, processes and distributes protein sequence data within the framework of the tripartite association of the PIR-International Protein Sequence Database (,). MIPS contributes nearly 50% of the data input to the PIR-International Protein Sequence Database. The database is distributed on CD-ROM together with PATCHX, an exhaustive supplement of unique, unverified protein sequences from external sources compiled by MIPS. Through its WWW server (http://www.mips.biochem.mpg.de/ ) MIPS permits internet access to sequence databases, homology data and to yeast genome information. (i) Sequence similarity results from the FASTA program () are stored in the FASTA database for all proteins from PIR-International and PATCHX. The database is dynamically maintained and permits instant access to FASTA results. (ii) Starting with FASTA database queries, proteins have been classified into families and superfamilies (PROT-FAM). (iii) The HPT (hashed position tree) data structure () developed at MIPS is a new approach for rapid sequence and pattern searching. (iv) MIPS provides access to the sequence and annotation of the complete yeast genome (), the functional classification of yeast genes (FunCat) and its graphical display, the 'Genome Browser' (). A CD-ROM based on the JAVA programming language providing dynamic interactive access to the yeast genome and the related protein sequences has been compiled and is available on request.