Predicting sub-cellular localization of tRNA synthetases from their primary structures

Since endo-symbiotic events occur, all genes of mitochondrial aminoacyl tRNA synthetase (AARS) were lost or transferred from ancestral mitochondrial genome into the nucleus. The canonical pattern is that both cytosolic and mitochondrial AARSs coexist in the nuclear genome. In the present scenario all mitochondrial AARSs are nucleus-encoded, synthesized on cytosolic ribosomes and post-translationally imported from the cytosol into the mitochondria in eukaryotic cell. The site-based discrimination between similar types of enzymes is very challenging because they have almost same physico-chemical properties. It is very important to predict the sub-cellular location of AARSs, to understand the mitochondrial protein synthesis. We have analyzed and optimized the distinguishable patterns between cytosolic and mitochondrial AARSs. Firstly, support vector machines (SVM)-based modules have been developed using amino acid and dipeptide compositions and achieved Mathews correlation coefficient (MCC) of 0.82 and 0.73, respectively. Secondly, we have developed SVM modules using position-specific scoring matrix and achieved the maximum MCC of 0.78. Thirdly, we developed SVM modules using N-terminal, intermediate residues, C-terminal and split amino acid composition (SAAC) and achieved MCC of 0.82, 0.70, 0.39 and 0.86, respectively. Finally, a SVM module was developed using selected attributes of split amino acid composition (SA-SAAC) approach and achieved MCC of 0.92 with an accuracy of 96.00%. All modules were trained and tested on a non-redundant data set and evaluated using fivefold cross-validation technique. On the independent data sets, SA-SAAC based prediction model achieved MCC of 0.95 with an accuracy of 97.77%. The web-server 'MARSpred' based on above study is available at http://www.imtech.res.in/raghava/marspred/.     

Voro3D: 3D Voronoi tessellations applied to protein structures

SUMMARY: Voro3D is an original easy-to-use tool, which provides a brand new point of view on protein structures through the three-dimensional (3D) Voronoi tessellations. To construct the Voronoi cells associated with each amino acid by a number of different tessellation methods, Voro3D uses a protein structure file in the PDB format as an input. After calculation, different structural properties of interest like secondary structures assignment, environment accessibility and exact contact matrices can be derived without any geometrical cut-off. Voro3D provides also a visualization of these tessellations superimposed on the associated protein structure, from which it is possible to model a polygonal protein surface using a model solvent or to quantify, for instance, the contact areas between a protein and a ligand. AVAILABILITY: The software executable file for PC using Windows 98, 2000, NT, XP can be freely downloaded at http://www.lmcp.jussieu.fr/~mornon/voronoi.html CONTACT: franck.dupuis@sanofi-aventis.com; jean-paul-mornon@imcp.jussieu.fr.     

Protein Peeling 3D: new tools for analyzing protein structures

We present an improved version of our Protein Peeling web server dedicated to the analysis of protein structure architecture through the identification of protein units produced by an iterative splitting algorithm. New features include identification of structural domains, detection of unstructured terminal elements and evaluation of the stability of protein unit structures. AVAILABILITY: The website is free and open to all users with no login requirements at http://www.dsimb.inserm.fr/dsimb-tools/peeling3.     

ProSAT: functional annotation of protein 3D structures

ProSAT (for Protein Structure Annotation Tool) is a tool to facilitate interactive visualization of non-structure-based functional annotations in protein 3D structures. It performs automated mapping of the functional annotations onto the protein structure and allows functional sites to be readily identified upon visualization. The current version of ProSAT can be applied to large datasets of protein structures for fast visual identification of active and other functional sites derived from the SwissProt and Prosite databases.     

Pre-docking filter for protein and ligand 3D structures

Virtual drug screening using protein-ligand docking techniques is a time-consuming process, which requires high computational power for binding affinity calculation. There are millions of chemical compounds available for docking. Eliminating compounds that are unlikely to exhibit high binding affinity from the screening set should speed-up the virtual drug screening procedure. We performed docking of 6353 ligands against twenty-one protein X-ray crystal structures. The docked ligands were ranked according to their calculated binding affinities, from which the top five hundred and the bottom five hundred were selected. We found that the volume and number of rotatable bonds of the top five hundred docked ligands are similar to those found in the crystal structures and corresponded with the volume of the binding sites. In contrast, the bottom five hundred set contains ligands that are either too large to enter the binding site, or too small to bind with high specificity and affinity to the binding site. A pre-docking filter that takes into account shapes and volumes of the binding sites as well as ligand volumes and flexibilities can filter out low binding affinity ligands from the screening sets. Thus, the virtual drug screening procedure speed is increased.     

Glycosylation controls sodium-calcium exchanger 3 sub-cellular localization during cell cycle

The Na⁺/Ca²⁺ exchanger (NCX) is a membrane antiporter that has been identified in the plasma membrane, the inner membrane of the nuclear envelope and in the membrane of the endoplasmic reticulum (ER). In humans, three genes have been identified, encoding unique NCX proteins. Although extensively studied, the NCX’s sub-cellular localization and mechanisms regulating the activity of different subtypes are still ambiguous. Here we investigated the subcellular localization of the NCX subtype 3 (NCX3) and its impact on the cell cycle. Two phenotypes, switching from one to the other during the cell cycle, were detected. One phenotype was NCX3 in the plasma membrane during S and M phase, and the other was NCX3 in the ER membrane during resting and interphase. Glycosylation of NCX3 at the N45 site was required for targeting the protein to the plasma membrane, and the N45 site functioned as an on-off switch for the translocation of NCX3 to either the plasma membrane or the membrane of the ER. Introduction of an N-glycosylation deficient NCX3 mutant led to an arrest of cells in the G0/G1 phase of the cell cycle. This was accompanied by accumulation of de-glycosylated NCX3 in the cytosol (that is in the ER), where it transported calcium ions (Ca²⁺) from the cytosol to the ER. These results, obtained in transfected HEK293T and HeLa and confirmed endogenously in SH-SY5Y cells, suggest that cells can use a dynamic Ca²⁺ signaling toolkit in which the NCX3 sub-cellular localization changes in synchrony with the cell cycle.     

LGA: A method for finding 3D similarities in protein structures

We present the LGA (Local-Global Alignment) method, designed to facilitate the comparison of protein structures or fragments of protein structures in sequence dependent and sequence independent modes. The LGA structure alignment program is available as an online service at http://PredictionCenter.llnl.gov/local/lga. Data generated by LGA can be successfully used in a scoring function to rank the level of similarity between two structures and to allow structure classification when many proteins are being analyzed. LGA also allows the clustering of similar fragments of protein structures.     

Efficient shape descriptors for feature extraction in 3D protein structures

Structural Genomics initiatives are generating an increasing number of protein structures with very limited biochemical characterization. Characterization of a protein's function and understanding the specific nature of a protein's binding is a critical part of both protein engineering and structure-based drug discovery. The accurate detection of binding site in these protein structures can be valuable in determining its function. As shape plays a crucial role in bimolecular recognition and function, the development of shape analysis techniques is important for understanding protein structure-function relationships. This paper describes the use of the continuous wavelet transforms (CWT) for characterizing shape features of 3D protein structures. The goal is to explore the CWT as a multiscale tool to generate rotation- and translation-invariant shape features.     

Reconstruction of 3D structures from protein contact maps

The prediction of the protein tertiary structure from solely its residue sequence (the so called Protein Folding Problem) is one of the most challenging problems in Structural Bioinformatics. We focus on the protein residue contact map. When this map is assigned it is possible to reconstruct the 3D structure of the protein backbone. The general problem of recovering a set of 3D coordinates consistent with some given contact map is known as a unit-disk-graph realization problem and it has been recently proven to be NP-Hard. In this paper we describe a heuristic method (COMAR) that is able to reconstruct with an unprecedented rate (3-15 seconds) a 3D model that exactly matches the target contact map of a protein. Working with a non-redundant set of 1760 proteins, we find that the scoring efficiency of finding a 3D model very close to the protein native structure depends on the threshold value adopted to compute the protein residue contact map. Contact maps whose threshold values range from 10 to 18 Ångstroms allow reconstructing 3D models that are very similar to the proteins native structure.     

Library of disordered patterns in 3D protein structures

Intrinsically disordered regions serve as molecular recognition elements, which play an important role in the control of many cellular processes and signaling pathways. It is useful to be able to predict positions of disordered regions in protein chains. The statistical analysis of disordered residues was done considering 34,464 unique protein chains taken from the PDB database. In this database, 4.95% of residues are disordered (i.e. invisible in X-ray structures). The statistics were obtained separately for the N- and C-termini as well as for the central part of the protein chain. It has been shown that frequencies of occurrence of disordered residues of 20 types at the termini of protein chains differ from the ones in the middle part of the protein chain. Our systematic analysis of disordered regions in PDB revealed 109 disordered patterns of different lengths. Each of them has disordered occurrences in at least five protein chains with identity less than 20%. The vast majority of all occurrences of each disordered pattern are disordered. This allows one to use the library of disordered patterns for predicting the status of a residue of a given protein to be ordered or disordered. We analyzed the occurrence of the selected patterns in three eukaryotic and three bacterial proteomes.     

Understanding and predicting protein assemblies with 3D structures

Protein interactions are central to most biological processes, and are currently the subject of great interest. Yet despite the many recently developed methods for interaction discovery, little attention has been paid to one of the best sources of data: complexes of known three-dimensional (3D) structure. Here we discuss how such complexes can be used to study and predict protein interactions and complexes, and to interrogate interaction networks proposed by methods such as two-hybrid screens or affinity purifications.     

SOSUI-GramN: high performance prediction for sub-cellular localization of proteins in gram-negative bacteria

A predictive software system, SOSUI-GramN, was developed for assessing the subcellular localization of proteins in Gram-negative bacteria. The system does not require the sequence homology data of any known sequences; instead, it uses only physicochemical parameters of the N- and C-terminal signal sequences, and the total sequence. The precision of the prediction system for subcellular localization to extracellular, outer membrane, periplasm, inner membrane and cytoplasmic medium was 92.3%, 89.4%, 86.4%, 97.5% and 93.5%, respectively, with corresponding recall rates of 70.3%, 87.5%, 76.0%, 97.5% and 88.4%, respectively. The overall performance for precision and recall obtained using this method was 92.9% and 86.7%, respectively. The comparison of performance of SOSUI-GramN with that of other methods showed the performance of prediction for extracellular proteins, as well as inner and outer membrane proteins, was either superior or equivalent to that obtained with other systems. SOSUI-GramN particularly improved the accuracy for predictions of extracellular proteins which is an area of weakness common to the other methods.     

Monoclonal antibodies against human translation termination factor eRF3 and their utilization for sub-cellular localization of eRF3

Eukaryotic translation termination is triggered by peptide release factors eRF1 and eRF3. eRF1 recognizes the stop codon and promotes nascent peptide chain release, while eRF3 facilitates this peptide release in a GTP-dependent manner. In addition to its role in termination, eRF3 is involved in normal and nonsense-mediated mRNA decay. Despite extensive investigation, the complete understanding of eRF3 function have been hampered by the lack of specific anti-eRF3 monoclonal antibodies (Mabs). The purpose of the study was production of recombinant eRF3a/GSPT1, development of anti-eRF3a/GSPT1 Mabs and their utilization for eRF3a/GSPT1 sub-cellular localization. Plasmid encoding C-terminal part of human GSPT1/eRF3a was constructed. Purified protein, which was predominantly present in the inclusion bodies, was used for the development of Mabs. Characterization of the regions recognized by Mabs using GSPT1/eRF3a mutants and its visualization in the 3D space suggested that Mabs recognize different epitopes. Consistent with its function in translational termination, immunostaining of the cells with developed Mabs revealed that the endogenous GSPT1/eRF3a localized in endoplasmic reticulum. Taking into account the important role of eRF3 for the fundamental research one can suggests that developed Mabs have great prospective to be used as a research reagent in a wide range of applications.     

ESyPred3D: Prediction of proteins 3D structures

MOTIVATION: Homology or comparative modeling is currently the most accurate method to predict the three-dimensional structure of proteins. It generally consists in four steps: (1) databanks searching to identify the structural homolog, (2) target-template alignment, (3) model building and optimization, and (4) model evaluation. The target-template alignment step is generally accepted as the most critical step in homology modeling. RESULTS: We present here ESyPred3D, a new automated homology modeling program. The method gets benefit of the increased alignment performances of a new alignment strategy. Alignments are obtained by combining, weighting and screening the results of several multiple alignment programs. The final three-dimensional structure is build using the modeling package MODELLER. ESyPred3D was tested on 13 targets in the CASP4 experiment (Critical Assessment of Techniques for Proteins Structural Prediction). Our alignment strategy obtains better results compared to PSI-BLAST alignments and ESyPred3D alignments are among the most accurate compared to those of participants having used the same template. AVAILABILITY: ESyPred3D is available through its web site at http://www.fundp.ac.be/urbm/bioinfo/esypred/ CONTACT: christophe.lambert@fundp.ac.be; http://www.fundp.ac.be/~lambertc     

High-level expression, polyclonal antibody preparation and sub-cellular localization analysis of mouse Rhox5 protein

Mouse reproductive homeobox on the X chromosome (Rhox) is a novel homeobox gene cluster. Rhox5, also called Pem, belongs to the beta subcluster of Rhox. Codon analysis indicated that the cDNA contains 16% of codons rarely used in Escherichia coli. To achieve high-level expression of Rhox5, the coding sequence of Rhox5 was amplified and subcloned into the prokaryotic expression vector pET22b (+) in order to produce 6His-tagged fusion protein in the modified BL21 (DE3) cells, namely Rosetta2 (DE3) cells. The 6His-tagged Rhox5 was expressed efficiently in Rosetta2 (DE3), compared with marginal expression in BL21 (DE3). The fusion protein amounted to 16% of the total bacterial proteins after induction with 0.4mM IPTG for 1.5h at 37 degrees C. After purification, Rhox5-6His was used to immunize New Zealand white rabbits following standard protocol. The homemade antiserum could detect both endogenous Rhox5 protein expressed in eukaryotic cells (Cos-7) and exogenous GFP-Rhox5 protein. Furthermore, the antiserum was used to determine the localization of Rhox5 in NIH3T3 cells using an immunofluorescence technique. The results demonstrated that Rhox5 was localized predominantly in the nucleus. Preparation of the anti-Rhox5 polyclonal antibody will facilitate further functional study of Rhox5.     

Exploiting 3D structural templates for detection of metal-binding sites in protein structures

High throughput structural genomics efforts have been making the structures of proteins available even before their function has been fully characterized. Therefore, methods that exploit the structural knowledge to provide evidence about the functions of proteins would be useful. Such methods would be needed to complement the sequence-based function annotation approaches. The current study describes generation of 3D-structural motifs for metal-binding sites from the known metalloproteins. It then scans all the available protein structures in the PDB database for putative metal-binding sites. Our analysis predicted more than 1000 novel metal-binding sites in proteins using three-residue templates, and more than 150 novel metal-binding sites using four-residue templates. Prediction of metal-binding site in a yeast protein YDR533c led to the hypothesis that it might function as metal-dependent amidopeptidase. The structural motifs identified by our method present novel metal-binding sites that reveal newer mechanisms for a few well-known proteins.     

Assessing a novel approach for predicting local 3D protein structures from sequence

We developed a novel approach for predicting local protein structure from sequence. It relies on the Hybrid Protein Model (HPM), an unsupervised clustering method we previously developed. This model learns three-dimensional protein fragments encoded into a structural alphabet of 16 protein blocks (PBs). Here, we focused on 11-residue fragments encoded as a series of seven PBs and used HPM to cluster them according to their local similarities. We thus built a library of 120 overlapping prototypes (mean fragments from each cluster), with good three-dimensional local approximation, i.e., a mean accuracy of 1.61 A Calpha root-mean-square distance. Our prediction method is intended to optimize the exploitation of the sequence-structure relations deduced from this library of long protein fragments. This was achieved by setting up a system of 120 experts, each defined by logistic regression to optimize the discrimination from sequence of a given prototype relative to the others. For a target sequence window, the experts computed probabilities of sequence-structure compatibility for the prototypes and ranked them, proposing the top scorers as structural candidates. Predictions were defined as successful when a prototype <2.5 A from the true local structure was found among those proposed. Our strategy yielded a prediction rate of 51.2% for an average of 4.2 candidates per sequence window. We also proposed a confidence index to estimate prediction quality. Our approach predicts from sequence alone and will thus provide valuable information for proteins without structural homologs. Candidates will also contribute to global structure prediction by fragment assembly.     

1,25-Dihydroxyvitamin D3-induced calcium-binding protein: localization in organ-cultured embryonic chick duodenum

1,25-Dihydroxyvitamin D3 (1,25(OH)2D3) induces de novo biosynthesis of a specific calcium-binding protein (CaBP) in embryonic chick duodenum in organ culture. Using a highly sensitive and specific, peroxidase-antiperoxidase immunocytochemical procedure, 1,25(OH)2D3-induced CaBP in the organ-cultured duodenum was found only in the cytoplasm of absorptive cells, corresponding to its localization in rachitic chick duodenal cells after a single injection of 1,25(OH)2D3 in vivo. This observation, along with evidence correlating CaBP with calcium transport, strongly supports the use of the embryonic chick duodenal organ culture system as a physiologically relevant model of the vitamin D-dependent calcium absorptive mechanism.