Identification of protein functions from a molecular surface database,eF-site

A bioinformatics method was developed to identify the protein surface around the functional site and to estimate the biochemical function, using a newly constructed molecular surface database named the eF-site (electrostatic surface of Functional site. Molecular surfaces of protein molecules were computed based on the atom coordinates, and the eF-site database was prepared by adding the physical properties on the constructed molecular surfaces. The electrostatic potential on each molecular surface was individually calculated solving the Poisson–Boltzmann equation numerically for the precise continuum model, and the hydrophobicity information of each residue was also included. The eF-site database is accessed by the internet (http://pi.protein.osaka-u.ac.jp/eF-site/). We have prepared four different databases, eF-site/antibody, eF-site/prosite, eF-site/P-site, and eF-site/ActiveSite, corresponding to the antigen binding sites of antibodies with the same orientations, the molecular surfaces for the individual motifs in PROSITE database, the phosphate binding sites, and the active site surfaces for the representatives of the individual protein family, respectively. An algorithm using the clique detection method as an applied graph theory was developed to search of the eF-site database, so as to recognize and discriminate the characteristic molecular surfaces of the proteins. The method identifies the active site having the similar function to those of the known proteins.     

Expanding the nitrogen regulatory protein superfamily: Homology detection at below random sequence identity

Nitrogen regulatory (PII) proteins are signal transduction molecules involved in controlling nitrogen metabolism in prokaryots. PII proteins integrate the signals of intracellular nitrogen and carbon status into the control of enzymes involved in nitrogen assimilation. Using elaborate sequence similarity detection schemes, we show that five clusters of orthologs (COGs) and several small divergent protein groups belong to the PII superfamily and predict their structure to be a (betaalphabeta)(2) ferredoxin-like fold. Proteins from the newly emerged PII superfamily are present in all major phylogenetic lineages. The PII homologs are quite diverse, with below random (as low as 1%) pairwise sequence identities between some members of distant groups. Despite this sequence diversity, evidence suggests that the different subfamilies retain the PII trimeric structure important for ligand-binding site formation and maintain a conservation of conservations at residue positions important for PII function. Because most of the orthologous groups within the PII superfamily are composed entirely of hypothetical proteins, our remote homology-based structure prediction provides the only information about them. Analogous to structural genomics efforts, such prediction gives clues to the biological roles of these proteins and allows us to hypothesize about locations of functional sites on model structures or rationalize about available experimental information. For instance, conserved residues in one of the families map in close proximity to each other on PII structure, allowing for a possible metal-binding site in the proteins coded by the locus known to affect sensitivity to divalent metal ions. Presented analysis pushes the limits of sequence similarity searches and exemplifies one of the extreme cases of reliable sequence-based structure prediction. In conjunction with structural genomics efforts to shed light on protein function, our strategies make it possible to detect homology between highly diverse sequences and are aimed at understanding the most remote evolutionary connections in the protein world.     

Pharmacophore-based molecular docking to account for ligand flexibility

Rapid computational mining of large 3D molecular databases is central to generating new drug leads. Accurate virtual screening of large 3D molecular databases requires consideration of the conformational flexibility of the ligand molecules. Ligand flexibility can be included without prohibitively increasing the search time by docking ensembles of precomputed conformers from a conformationally expanded database. A pharmacophore-based docking method whereby conformers of the same or different molecules are overlaid by their largest 3D pharmacophore and simultaneously docked by partial matches to that pharmacophore is presented. The method is implemented in DOCK 4.0.     

Mining of assembled expressed sequence tag (EST) data for protein families: application to the G protein-coupled receptor superfamily

The availability of large expressed sequence tag (EST) databases has led to a revolution in the way new genes are identified. Mining of these databases using known protein sequences as queries is a powerful technique for discovering orthologous and paralogous genes. The scientist is often confronted, however, by an enormous amount of search output owing to the inherent redundancy of EST data. In addition, high search sensitivity often cannot be achieved using only a single member of a protein superfamily as a query. In this paper a technique for addressing both of these issues is described. Assembled EST databases are queried with every member of a protein superfamily, the results are integrated and false positives are pruned from the set. The result is a set of assemblies enriched in members of the protein superfamily under consideration. The technique is applied to the G protein-coupled receptor (GPCR) superfamily in the construction of a GPCR Resource. A novel full-length human GPCR identified from the GPCR Resource is presented, illustrating the utility of the method.     

HOMSTRAD: a database of protein structure alignments for homologous families.

We describe a database of protein structure alignments for homologous families. The database HOMSTRAD presently contains 130 protein families and 590 aligned structures, which have been selected on the basis of quality of the X-ray analysis and accuracy of the structure. For each family, the database provides a structure-based alignment derived using COMPARER and annotated with JOY in a special format that represents the local structural environment of each amino acid residue. HOMSTRAD also provides a set of superposed atomic coordinates obtained using MNYFIT, which can be viewed with a graphical user interface or used for comparative modeling studies. The database is freely available on the World Wide Web at: http://www-cryst.bioc.cam. ac.uk/-homstrad/, with search facilities and links to other databases.     

Structure, dynamics and mapping of membrane-binding residues of micelle-bound antimicrobial peptides by natural abundance C NMR spectroscopy

Worldwide bacterial resistance to traditional antibiotics has drawn much research attention to naturally occurring antimicrobial peptides (AMPs) owing to their potential as alternative antimicrobials. Structural studies of AMPs are essential for an in-depth understanding of their activity, mechanism of action, and in guiding peptide design. Two-dimensional solution proton NMR spectroscopy has been the major tool. In this article, we describe the applications of natural abundance ¹³C NMR spectroscopy that provides complementary information to 2D ¹H NMR. The correlation of ¹³Cα secondary shifts with both 3D structure and heteronuclear ¹⁵N NOE values indicates that natural abundance carbon chemical shifts are useful probes for backbone structure and dynamics of membrane peptides. Using human LL-37-derived peptides (GF-17, KR-12, and RI-10), as well as amphibian antimicrobial and anticancer peptide aurein 1.2 and its analog LLAA, as models, we show that the cross peak intensity plots of 2D ¹H-¹³Cα HSQC spectra versus residue number present a wave-like pattern (HSQC wave) where key hydrophobic residues of micelle-bound peptides are located in the troughs with weaker intensities, probably due to fast exchange between the free and bound forms. In all the cases, the identification of aromatic phenylalanines as a key membrane-binding residue is consistent with previous intermolecular Phe-lipid NOE observations. Furthermore, mutation of one of the key hydrophobic residues of KR-12 to Ala significantly reduced the antibacterial activity of the peptide mutants. These results illustrate that natural abundance heteronuclear-correlated NMR spectroscopy can be utilized to probe backbone structure and dynamics, and perhaps to map key membrane-binding residues of peptides in complex with micelles. ¹H-¹³Cα HSQC wave, along with other NMR waves such as dipolar wave and chemical shift wave, offers novel insights into peptide-membrane interactions from different angles.     

YTPdb: A wiki database of yeast membrane transporters

Membrane transporters constitute one of the largest functional categories of proteins in all organisms. In the yeast Saccharomyces cerevisiae, this represents about 300 proteins (∼ 5% of the proteome). We here present the Yeast Transport Protein database (YTPdb), a user-friendly collaborative resource dedicated to the precise classification and annotation of yeast transporters. YTPdb exploits an evolution of the MediaWiki web engine used for popular collaborative databases like Wikipedia, allowing every registered user to edit the data in a user-friendly manner. Proteins in YTPdb are classified on the basis of functional criteria such as subcellular location or their substrate compounds. These classifications are hierarchical, allowing queries to be performed at various levels, from highly specific (e.g. ammonium as a substrate or the vacuole as a location) to broader (e.g. cation as a substrate or inner membranes as location). Other resources accessible for each transporter via YTPdb include post-translational modifications, K_m values, a permanently updated bibliography, and a hierarchical classification into families. The YTPdb concept can be extrapolated to other organisms and could even be applied for other functional categories of proteins. YTPdb is accessible at http://homes.esat.kuleuven.be/ytpdb/.