TOPOFIT-DB, a database of protein structural alignments based on the TOPOFIT method

TOPOFIT-DB (T-DB) is a public web-based database of protein structural alignments based on the TOPOFIT method, providing a comprehensive resource for comparative analysis of protein structure families. The TOPOFIT method is based on the discovery of a saturation point on the alignment curve (topomax point) which presents an ability to objectively identify a border between common and variable parts in a protein structural family, providing additional insight into protein comparison and functional annotation. TOPOFIT also effectively detects non-sequential relations between protein structures. T-DB provides users with the convenient ability to retrieve and analyze structural neighbors for a protein; do one-to-all calculation of a user provided structure against the entire current PDB release with T-Server, and pair-wise comparison using the TOPOFIT method through the T-Pair web page. All outputs are reported in various web-based tables and graphics, with automated viewing of the structure-sequence alignments in the Friend software package for complete, detailed analysis. T-DB presents researchers with the opportunity for comprehensive studies of the variability in proteins and is publicly available at

Molecular Model of the Microvillar Cytoskeleton and Organization of the Brush Border

BACKGROUND. Brush border microvilli are ~1-µm long finger-like projections emanating from the apical surfaces of certain, specialized absorptive epithelial cells. A highly symmetric hexagonal array of thousands of these uniformly sized structures form the brush border, which in addition to aiding in nutrient absorption also defends the large surface area against pathogens. Here, we present a molecular model of the protein cytoskeleton responsible for this dramatic cellular morphology. METHODOLOGY/PRINCIPAL FINDINGS. The model is constructed from published crystallographic and microscopic structures reported by several groups over the last 30+ years. Our efforts resulted in a single, unique, self-consistent arrangement of actin, fimbrin, villin, brush border myosin (Myo1A), calmodulin, and brush border spectrin. The central actin core bundle that supports the microvillus is nearly saturated with fimbrin and villin cross-linkers and has a density similar to that found in protein crystals. The proposed model accounts for all major proteinaceous components, reproduces the experimentally determined stoichiometry, and is consistent with the size and morphology of the biological brush border membrane. CONCLUSIONS/SIGNIFICANCE. The model presented here will serve as a structural framework to explain many of the dynamic cellular processes occurring over several time scales, such as protein diffusion, association, and turnover, lipid raft sorting, membrane deformation, cytoskeletal-membrane interactions, and even effacement of the brush border by invading pathogens. In addition, this model provides a structural basis for evaluating the equilibrium processes that result in the uniform size and structure of the highly dynamic microvilli.Boston University (Graduate Student Research Fellowship); National Institutes of Health (GM62886

Structure SNP (StSNP): a web server for mapping and modeling nsSNPs on protein structures with linkage to metabolic pathways

SNPs located within the open reading frame of a gene that result in an alteration in the amino acid sequence of the encoded protein [nonsynonymous SNPs (nsSNPs)] might directly or indirectly affect functionality of the protein, alone or in the interactions in a multi-protein complex, by increasing/decreasing the activity of the metabolic pathway. Understanding the functional consequences of such changes and drawing conclusions about the molecular basis of diseases, involves integrating information from multiple heterogeneous sources including sequence, structure data and pathway relations between proteins. The data from NCBI's SNP database (dbSNP), gene and protein databases from Entrez, protein structures from the PDB and pathway information from KEGG have all been cross referenced into the StSNP web server, in an effort to provide combined integrated, reports about nsSNPs. StSNP provides ‘on the fly’ comparative modeling of nsSNPs with links to metabolic pathway information, along with real-time visual comparative analysis of the modeled structures using the Friend software application. The use of metabolic pathways in StSNP allows a researcher to examine possible disease-related pathways associated with a particular nsSNP(s), and link the diseases with the current available molecular structure data. The server is publicly available at http://glinka.bio.neu.edu/StSNP/

Tools for integrated sequence-structure analysis with UCSF Chimera

BACKGROUND: Comparing related structures and viewing the structures in the context of sequence alignments are important tasks in protein structure-function research. While many programs exist for individual aspects of such work, there is a need for interactive visualization tools that: (a) provide a deep integration of sequence and structure, far beyond mapping where a sequence region falls in the structure and vice versa; (b) facilitate changing data of one type based on the other (for example, using only sequence-conserved residues to match structures, or adjusting a sequence alignment based on spatial fit); (c) can be used with a researcher's own data, including arbitrary sequence alignments and annotations, closely or distantly related sets of proteins, etc.; and (d) interoperate with each other and with a full complement of molecular graphics features. We describe enhancements to UCSF Chimera to achieve these goals. RESULTS: The molecular graphics program UCSF Chimera includes a suite of tools for interactive analyses of sequences and structures. Structures automatically associate with sequences in imported alignments, allowing many kinds of crosstalk. A novel method is provided to superimpose structures in the absence of a pre-existing sequence alignment. The method uses both sequence and secondary structure, and can match even structures with very low sequence identity. Another tool constructs structure-based sequence alignments from superpositions of two or more proteins. Chimera is designed to be extensible, and mechanisms for incorporating user-specific data without Chimera code development are also provided. CONCLUSION: The tools described here apply to many problems involving comparison and analysis of protein structures and their sequences. Chimera includes complete documentation and is intended for use by a wide range of scientists, not just those in the computational disciplines. UCSF Chimera is free for non-commercial use and is available for Microsoft Windows, Apple Mac OS X, Linux, and other platforms from

An AP Endonuclease 1–DNA Polymerase β Complex: Theoretical Prediction of Interacting Surfaces

Abasic (AP) sites in DNA arise through both endogenous and exogenous mechanisms. Since AP sites can prevent replication and transcription, the cell contains systems for their identification and repair. AP endonuclease (APEX1) cleaves the phosphodiester backbone 5′ to the AP site. The cleavage, a key step in the base excision repair pathway, is followed by nucleotide insertion and removal of the downstream deoxyribose moiety, performed most often by DNA polymerase beta (pol-β). While yeast two-hybrid studies and electrophoretic mobility shift assays provide evidence for interaction of APEX1 and pol-β, the specifics remain obscure. We describe a theoretical study designed to predict detailed interacting surfaces between APEX1 and pol-β based on published co-crystal structures of each enzyme bound to DNA. Several potentially interacting complexes were identified by sliding the protein molecules along DNA: two with pol-β located downstream of APEX1 (3′ to the damaged site) and three with pol-β located upstream of APEX1 (5′ to the damaged site). Molecular dynamics (MD) simulations, ensuring geometrical complementarity of interfaces, enabled us to predict interacting residues and calculate binding energies, which in two cases were sufficient (∼−10.0 kcal/mol) to form a stable complex and in one case a weakly interacting complex. Analysis of interface behavior during MD simulation and visual inspection of interfaces allowed us to conclude that complexes with pol-β at the 3′-side of APEX1 are those most likely to occur in vivo. Additional multiple sequence analyses of APEX1 and pol-β in related organisms identified a set of correlated mutations of specific residues at the predicted interfaces. Based on these results, we propose that pol-β in the open or closed conformation interacts and makes a stable interface with APEX1 bound to a cleaved abasic site on the 3′ side. The method described here can be used for analysis in any DNA-metabolizing pathway where weak interactions are the principal mode of cross-talk among participants and co-crystal structures of the individual components are available

MultiSeq: unifying sequence and structure data for evolutionary analysis

BACKGROUND: Since the publication of the first draft of the human genome in 2000, bioinformatic data have been accumulating at an overwhelming pace. Currently, more than 3 million sequences and 35 thousand structures of proteins and nucleic acids are available in public databases. Finding correlations in and between these data to answer critical research questions is extremely challenging. This problem needs to be approached from several directions: information science to organize and search the data; information visualization to assist in recognizing correlations; mathematics to formulate statistical inferences; and biology to analyze chemical and physical properties in terms of sequence and structure changes. RESULTS: Here we present MultiSeq, a unified bioinformatics analysis environment that allows one to organize, display, align and analyze both sequence and structure data for proteins and nucleic acids. While special emphasis is placed on analyzing the data within the framework of evolutionary biology, the environment is also flexible enough to accommodate other usage patterns. The evolutionary approach is supported by the use of predefined metadata, adherence to standard ontological mappings, and the ability for the user to adjust these classifications using an electronic notebook. MultiSeq contains a new algorithm to generate complete evolutionary profiles that represent the topology of the molecular phylogenetic tree of a homologous group of distantly related proteins. The method, based on the multidimensional QR factorization of multiple sequence and structure alignments, removes redundancy from the alignments and orders the protein sequences by increasing linear dependence, resulting in the identification of a minimal basis set of sequences that spans the evolutionary space of the homologous group of proteins. CONCLUSION: MultiSeq is a major extension of the Multiple Alignment tool that is provided as part of VMD, a structural visualization program for analyzing molecular dynamics simulations. Both are freely distributed by the NIH Resource for Macromolecular Modeling and Bioinformatics and MultiSeq is included with VMD starting with version 1.8.5. The MultiSeq website has details on how to download and use the software

A comprehensive analysis of non-sequential alignments between all protein structures

<p>Abstract</p> <p>Background</p> <p>The majority of relations between proteins can be represented as a conventional sequential alignment. Nevertheless, unusual non-sequential alignments with different connectivity of the aligned fragments in compared proteins have been reported by many researchers. It is interesting to understand those non-sequential alignments; are they unique, sporadic cases or they occur frequently; do they belong to a few specific folds or spread among many different folds, as a common feature of protein structure. We present here a comprehensive large-scale study of non-sequential alignments between available protein structures in Protein Data Bank.</p> <p>Results</p> <p>The study has been conducted on a non-redundant set of 8,865 protein structures aligned with the aid of the TOPOFIT method. It has been estimated that between 17.4% and 35.2% of all alignments are non-sequential depending on variations in the parameters. Analysis of the data revealed that non-sequential relations between proteins do occur systematically and in large quantities. Various sizes and numbers of non-sequential fragments have been observed with all possible complexities of fragment rearrangements found for alignments consisting of up to 12 fragments. It has been found that non-sequential alignments are not limited to proteins of any particular fold and are present in more than two hundred of them. Moreover, many of them are found between proteins with different fold assignments. It has been shown that protein structure symmetry does not explain non-sequential alignments. Therefore, compelling evidences have been provided that non-sequential alignments between proteins are systematic and widespread across the protein universe.</p> <p>Conclusion</p> <p>The phenomenon of the widespread occurrence of non-sequential alignments between proteins might represent a missing rule of protein structure organization. More detailed study of this phenomenon will enhance our understanding of protein stability, folding, and evolution.</p

Skipping of Exons by Premature Termination of Transcription and Alternative Splicing within Intron-5 of the Sheep SCF Gene: A Novel Splice Variant

Stem cell factor (SCF) is a growth factor, essential for haemopoiesis, mast cell development and melanogenesis. In the hematopoietic microenvironment (HM), SCF is produced either as a membrane-bound (−) or soluble (+) forms. Skin expression of SCF stimulates melanocyte migration, proliferation, differentiation, and survival. We report for the first time, a novel mRNA splice variant of SCF from the skin of white merino sheep via cloning and sequencing. Reverse transcriptase (RT)-PCR and molecular prediction revealed two different cDNA products of SCF. Full-length cDNA libraries were enriched by the method of rapid amplification of cDNA ends (RACE-PCR). Nucleotide sequencing and molecular prediction revealed that the primary 1519 base pair (bp) cDNA encodes a precursor protein of 274 amino acids (aa), commonly known as ‘soluble’ isoform. In contrast, the shorter (835 and/or 725 bp) cDNA was found to be a ‘novel’ mRNA splice variant. It contains an open reading frame (ORF) corresponding to a truncated protein of 181 aa (vs 245 aa) with an unique C-terminus lacking the primary proteolytic segment (28 aa) right after the D175G site which is necessary to produce ‘soluble’ form of SCF. This alternative splice (AS) variant was explained by the complete nucleotide sequencing of splice junction covering exon 5-intron (5)-exon 6 (948 bp) with a premature termination codon (PTC) whereby exons 6 to 9/10 are skipped (Cassette Exon, CE 6–9/10). We also demonstrated that the Northern blot analysis at transcript level is mediated via an intron-5 splicing event. Our data refine the structure of SCF gene; clarify the presence (+) and/or absence (−) of primary proteolytic-cleavage site specific SCF splice variants. This work provides a basis for understanding the functional role and regulation of SCF in hair follicle melanogenesis in sheep beyond what was known in mice, humans and other mammals

Importancia biológica de la Enzima Arginina Deiminasa en el proceso de crecimiento y diferenciación del protozoario intestinal Giardia Lamblia.

Tesis (Grado Doctor en Ciencias Biológicas)--Universidad Nacional de Córdoba. Facultad de Ciencias Exactas, Físicas y Naturales. Lugar de Trabajo: Laboratorio de Microbiología e Inmunología. Instituto de Investigación Médica Mercedes y Martín Ferreyra. INIMEC- CONICET- Universidad Nacional de Córdoba. 2014 - 155 h. + CD. tabls.; figus.; grafs. Contiene Referencia Bibliográfica Y Publicaciones Derivadas de la Tesis. Abstract en español e inglés.En este Trabajo de Tesis se estudiaron SUMOilación y citrulinación de proteínas en el parásito Giardia lamblia. Se demostró que tanto la proteína SUMO como las enzimas que participarían en la vía de SUMOilación se encuentran presentes en Giardia, identificándose a la enzima Arginina Deiminasa (ADI) como una proteína SUMOilada. ADI participa en la producción de energía y en el proceso de variación antigénica en el parásito, y durante el enquistamiento es translocada desde el citoplasma a los núcleos regulando negativamente la expresión de genes específicos del enquistamiento. Mutaciones puntuales en el sitio de SUMOilación de ADI mostraron que la translocación de la enzima al núcleo disminuye en las etapas tempranas de la diferenciación a quiste, aunque la enzima localiza en los núcleos en las últimas etapas del enquistamiento, participando en este proceso las señales de localización nuclear presentes en la secuencia de la proteína. Dentro del núcleo, ADI posiblemente citrulina histonas, encontrándose involucrada esta modificación en la menor expresión de las proteínas de enquistamiento y en la formación del quiste. Estos resultados indican que ADI desempeña un rol regulatorio en el proceso de enquistamiento, siendo las modificaciones post-traduccionales, mecanismos claves en este proceso

Effects of metal ions on the structural and biochemical properties of trypanosomatid phosphoglycerate mutases

Flagellate protozoa from the order Trypanosomatida have developed a range of strategies to survive in their mammalian hosts. A consequence is that the glycolytic pathway has assumed an important role, especially in bloodstream-form Trypanosoma brucei, where it is essential as the sole producer of ATP. The seventh enzyme in the pathway, 2,3-bisphosphoglycerate-independent phosphoglycerate mutase (iPGAM) is particularly attractive as a drug target because it shares no common properties with the corresponding enzyme in humans. This enzyme catalyses the conversion of 3PGA to 2PGA, with the requirement for metal ions to assist the catalytic function. In this study, two important biochemical and structural aspects of the enzyme were investigated: i) The in vitro and in vivo requirements for biologically relevant metal ions to support the activity of iPGAM, and ii) The ability of trypanosomatid iPGAM to exist in multiple conformations and oligomeric states in solution. The maximum activity of iPGAM in vitro requires Co2+, but this cannot be the case in vivo where ICP-OES analyses confirmed that Co2+ was essentially undetectable in T. brucei cytosolic fractions. The activity of iPGAM in vivo is therefore one of the lowest among the glycolytic enzymes. By contrast, Mg2+ and Zn2+ were found to be the most abundant metals in both cytosolic fractions and in purified bacterially expressed iPGAM. Our newly-developed multimode-plate reader discontinuous assay further revealed that of the biologically relevant metals, only Mg2+ can support iPGAM activity, but at less than 50% of the level of Co2+. By contrast, Zn2+ strongly inhibits iPGAM. This assay which was developed with minimal metal interference on the coupling enzymes, also showed that in solution, the ratio of the concentrations of 3PGA:2PGA (substrate:product) at equilibrium is not 1:1 as observed in the crystal structure, but is in fact 12:1, which may be due to the tighter binding of 2PGA to the enzyme. A series of biophysical analyses, notably by SEC-MALS showed that iPGAM from Leishmania mexicana, another trypanosomatid protozoan parasite exists in different forms and oligomeric states in solution, either as the closed-form monomer, openiii form monomer, or closed/open-form dimer which can be successfully separated by ion-exchange chromatography. The open-form LmiPGAM is particularly relevant for drug development, as the catalytic site in the closed-form structure is poorly inaccessible. Both virtual and high-throughput screening approaches were used to identify novel potential inhibitors. Out of a collection of 11 compounds tested at 1 mM, two showed substantial inhibition with 49% and 14% remaining activity. Taken together, the findings from this study demonstrated the potential of iPGAM to be a key modulator in controlling glycolytic flux in trypanosomes, and thus further validated it as an important drug target