iPBA: a tool for protein structure comparison using sequence alignment strategies

With the immense growth in the number of available protein structures, fast and accurate structure comparison has been essential. We propose an efficient method for structure comparison, based on a structural alphabet. Protein Blocks (PBs) is a widely used structural alphabet with 16 pentapeptide conformations that can fairly approximate a complete protein chain. Thus a 3D structure can be translated into a 1D sequence of PBs. With a simple Needleman-Wunsch approach and a raw PB substitution matrix, PB-based structural alignments were better than many popular methods. iPBA web server presents an improved alignment approach using (i) specialized PB Substitution Matrices (SM) and (ii) anchor-based alignment methodology. With these developments, the quality of similar to 88% of alignments was improved. iPBA alignments were also better than DALI, MUSTANG and GANGSTA(+) in > 80% of the cases. The webserver is designed to for both pairwise comparisons and database searches. Outputs are given as sequence alignment and superposed 3D structures displayed using PyMol and Jmol. A local alignment option for detecting subs-structural similarity is also embedded. As a fast and efficient `sequence-based' structure comparison tool, we believe that it will be quite useful to the scientific community. iPBA can be accessed at http://www.dsimb.inserm.fr/dsimb_tools/ipba/.

PRUNE and PROBE-two modular web services for protein-protein docking

The protein-protein docking programs typically perform four major tasks: (i) generation of docking poses, (ii) selecting a subset of poses, (iii) their structural refinement and (iv) scoring, ranking for the final assessment of the true quaternary structure. Although the tasks can be integrated or performed in a serial order, they are by nature modular, allowing an opportunity to substitute one algorithm with another. We have implemented two modular web services, (i) PRUNE: to select a subset of docking poses generated during sampling search (http://pallab.serc.iisc.ernet.in/prune) and (ii) PROBE: to refine, score and rank them (http://pallab.serc.iisc.ernet.in/probe). The former uses a new interface area based edge-scoring function to eliminate > 95% of the poses generated during docking search. In contrast to other multi-parameter-based screening functions, this single parameter based elimination reduces the computational time significantly, in addition to increasing the chances of selecting native-like models in the top rank list. The PROBE server performs ranking of pruned poses, after structure refinement and scoring using a regression model for geometric compatibility, and normalized interaction energy. While web-service similar to PROBE is infrequent, no web-service akin to PRUNE has been described before. Both the servers are publicly accessible and free for use.

Helical structures: The geometry of protein helices and nanotubes

In nature, helical structures arise when identical structural subunits combine sequentially, the orientational and translational relation between each unit and its predecessor remaining constant. A helical structure is thus generated by the repeated action of a screw transformation acting on a subunit. A plane hexagonal lattice wrapped round a cylinder provides a useful starting point for describing the helical conformations of protein molecules, for investigating the geometrical properties of carbon nanotubes, and for certain types of dense packings of equal spheres.

Mapping of B-cell epitopes of hemagglutinin protein of rinderpest virus

Monoclonal antibodies (mAbs) against secreted hemagglutinin (H) protein of rinderpest virus (RPV) expressed by a recombinant baculovirus were generated to characterize the antigenic sites on H protein and regions of functional significance. Three of the mAbs displayed hemagglutination inhibition activity and these mAbs were unable to neutralize virus infectivity. Western immunoblot analysis of overlapping deletion mutants indicated that three mAbs recognize antigenic regions at the extreme carboxy terminus (between amino acids 569 and 609) and the fourth mAb between amino acids 512 and 568. Using synthetic peptides, aa 569-577 and 575-583 were identified as the epitopes for E2G4 and D2F4, respectively. The epitopic domains of A12A9 and E2B6 mAbs were mapped to regions encompassing aa 527-554 and 588-609. Two epitopes spanning the extreme carboxy terminal region of aa 573 to 587 and 588 to 609 were shown to be immunodominant employing a competitive ELISA with polyclonal sera form vaccinated cattle. The D2F4 mAb which recognizes a unique epitope on RPV-H is not present on the closely related peste des petits ruminant virus FIN protein and this mAb could serve as a tool in the seromonitoring program after rinderpest vaccination. (C) 2002 Elsevier Science (USA).

PDB Goodies - a web-based GUI to manipulate the Protein Data Bank file

PDB Goodies is a web-based graphical user interface (GUI) to manipulate the Protein Data Bank file containing the three-dimensional atomic coordinates of protein structures. The program also allows users to save the manipulated three-dimensional atomic coordinate file on their local client system. These fragments are used in various stages of structure elucidation and analysis. This software is incorporated with all the three-dimensional protein structures available in the Protein Data Bank, which presently holds approximately 18 000 structures. In addition, this program works on a three-dimensional atomic coordinate file (Protein Data Bank format) uploaded from the client machine. The program is written using CGI/PERL scripts and is platform independent. The program PDB Goodies can be accessed over the World Wide Web at http:// 144.16.71.11/pdbgoodies/.

The repertoire of protein kinases encoded in the draft version of the human genome: atypical variations and uncommon domain combinations

Background: Phosphorylation by protein kinases is central to cellular signal transduction. Abnormal functioning of kinases has been implicated in developmental disorders and malignancies. Their activity is regulated by second messengers and by the binding of associated domains, which are also influential in translocating the catalytic component to their substrate sites, in mediating interaction with other proteins and carrying out their biological roles. Results: Using sensitive profile-search methods and manual analysis, the human genome has been surveyed for protein kinases. A set of 448 sequences, which show significant similarity to protein kinases and contain the critical residues essential for kinase function, have been selected for an analysis of domain combinations after classifying the kinase domains into subfamilies. The unusual domain combinations in particular kinases suggest their involvement in ubiquitination pathways and alternative modes of regulation for mitogen-activated protein kinase kinases (MAPKKs) and cyclin-dependent kinase (CDK)-like kinases. Previously unexplored kinases have been implicated in osteoblast differentiation and embryonic development on the basis of homology with kinases of known functions from other organisms. Kinases potentially unique to vertebrates are involved in highly evolved processes such as apoptosis, protein translation and tyrosine kinase signaling. In addition to coevolution with the kinase domain, duplication and recruitment of non-catalytic domains is apparent in signaling domains such as the PH, DAG-PE, SH2 and SH3 domains. Conclusions: Expansion of the functional repertoire and possible existence of alternative modes of regulation of certain kinases is suggested by their uncommon domain combinations. Experimental verification of the predicted implications of these kinases could enhance our understanding of their biological roles.

The requirements for Ca2+, protein phosphorylation and concurrent protein synthesis for zeatin signaling of acidic chitinase transcript accumulation in Cucumis sativus L

We have previously reported that both Ca2+ and staurosporine-sensitive protein kinase(s) are involved in the cytokinin zeatin induction of cucumber chitinase activity and its protein content (Barwe et al. 2001). To further characterize signal transduction events involved in this cytokinin induction of chitinase gene expression, Northern hybridizations of total RNAs prepared from excised, dark-grown cucumber cotyledons treated with cytokinins and/or various agonists and antagonists of signal transduction components, were carried out using a cucumber acidic chitinase (CACHT) cDNA probe (Metraux et al. 1989). CACHT mRNA increased by approximately 5- to 6-fold in response to exogenous zeatin (Z), zeatin riboside (ZR), and benzyladenine (BA) treatment, but failed to accumulate in response to kinetin (K). Among the cytokinins tested, Z was most effective. The Z-induced accumulation of CACHT mRNA was inhibited by a plasma membrane Ca2+ channel blocker verapamil. Treatment of cotyledons with exogenous CaCl2 and calcium ionophore A23187 in the presence and absence of cytokinin enhanced CACHT mRNA accumulation. These two observations suggest the participation of extracellular calcium in signaling Z-induction. Furthermore, the presence of staurosporine (an inhibitor of protein kinase) in Z treatment reduced CACHT mRNA, suggesting the involvement of phosphorylation of one or more cellular proteins. In addition, we provide evidence that the Z-induction of CACHT mRNA is blocked by protein synthesis inhibitor cycloheximide treatment. Taken together, these results suggest that Ca2+ influx from extracellular space, protein phosphorylation, and concurrent protein synthesis events participate in cytokinin signaling during Z-induced CACHT transcript accumulation.

Development of recombinant coat protein antibody based IC-RT-PCR for detection and discrimination of sugarcane streak mosaic virus isolates from Southern India

Sugarcane streak mosaic virus (SCSMV), causes mosaic disease of sugarcane and is thought to belong to a new undescribed genus in the family Potyviridae. The coat protein (CP) gene from the Andhra Pradesh (AP) isolate of SCSMV (SCSMV AP) was cloned and expressed in Escherichia coli. The recombinant coat protein was used to raise high quality antiserum. The CP antiserum was used to develop an immunocapture reverse transcription-polymerase chain reaction (IC-RT-PCR) based assay for the detection and discrimination of SCSMV isolates in South India. The sequence of the cloned PCR products encoding 3'untranslated region (UTR) and CP regions of the virus isolates from three different locations in South India viz. Tanuku (Coastal Andhra Pradesh), Coimbatore (Tamil Nadu) and Hospet (Karnataka) was compared with that of SCSMV AP The analysis showed that they share 89.4, 89.5 and 90% identity respectively at the nucleotide level. This suggests that the isolates causing mosaic disease of sugarcane in South India are indeed strains of SCSMV In addition, the sensitivity of the IC-RT-PCR was compared with direct antigen coating-enzyme linked immunosorbent assay (DAC-ELISA) and dot-blot immunobinding assays and was found to be more sensitive and hence could be used to detect the presence of virus in sugarcane breeding, germplasm centres and in quarantine programs.

Denaturant mediated unfolding of both native and molten globule states of maltose binding protein are accompanied by large [Delta]Cp's

Maltose binding protein (MBP) is a large, monomeric two domain protein containing 370 amino acids. In the absence of denaturant at neutral pH, the protein is in the native state, while at pH 3.0 it forms a molten globule. The molten globule lacks a tertiary circular dichroism signal but has secondary structure similar to that of the native state. The molten globule binds 8-anilino-1-naphthalene sulfonate (ANS). The unfolding thermodynamics of MBP at both pHs were measured by carrying out a series of isothermal urea melts at temperatures ranging from 274–329 K. At 298 K, values of [Delta]G°, [Delta]Cp, and Cm were 3.1 ± 0.2 kcal mol−1, 5.9 ± 0.8 kcal mol−1 K−1 (15.9 cal (mol-residue)−1 K−1), and 0.8 M, respectively, at pH 3.0 and 14.5 ± 0.4 kcal mol−1, 8.3 ± 0.7 kcal mol−1 K−1 (22.4 kcal (mol-residue)−1 K−1), and 3.3 M, respectively, at pH 7.1. Guanidine hydrochloride denaturation at pH 7.1 gave values of [Delta]G° and [Delta]Cp similar to those obtained with urea. The m values for denaturation are strongly temperature dependent, in contrast to what has been previously observed for small globular proteins. The value of [Delta]Cp per mol-residue for the molten globule is comparable to corresponding values of [Delta]Cp for the unfolding of typical globular proteins and suggests that it is a highly ordered structure, unlike molten globules of many small proteins. The value of [Delta]Cp per mol-residue for the unfolding of the native state is among the highest currently known for any protein.

PocketAlign A Novel Algorithm for Aligning Binding Sites in Protein Structures

A fundamental task in bioinformatics involves a transfer of knowledge from one protein molecule onto another by way of recognizing similarities. Such similarities are obtained at different levels, that of sequence, whole fold, or important substructures. Comparison of binding sites is important to understand functional similarities among the proteins and also to understand drug cross-reactivities. Current methods in literature have their own merits and demerits, warranting exploration of newer concepts and algorithms, especially for large-scale comparisons and for obtaining accurate residue-wise mappings. Here, we report the development of a new algorithm, PocketAlign, for obtaining structural superpositions of binding sites. The software is available as a web-service at http://proline.physicslisc.emetin/pocketalign/. The algorithm encodes shape descriptors in the form of geometric perspectives, supplemented by chemical group classification. The shape descriptor considers several perspectives with each residue as the focus and captures relative distribution of residues around it in a given site. Residue-wise pairings are computed by comparing the set of perspectives of the first site with that of the second, followed by a greedy approach that incrementally combines residue pairings into a mapping. The mappings in different frames are then evaluated by different metrics encoding the extent of alignment of individual geometric perspectives. Different initial seed alignments are computed, each subsequently extended by detecting consequential atomic alignments in a three-dimensional grid, and the best 500 stored in a database. Alignments are then ranked, and the top scoring alignments reported, which are then streamed into Pymol for visualization and analyses. The method is validated for accuracy and sensitivity and benchmarked against existing methods. An advantage of PocketAlign, as compared to some of the existing tools available for binding site comparison in literature, is that it explores different schemes for identifying an alignment thus has a better potential to capture similarities in ligand recognition abilities. PocketAlign, by finding a detailed alignment of a pair of sites, provides insights as to why two sites are similar and which set of residues and atoms contribute to the similarity.

Dynamics of water at the interface of a small protein, enterotoxin

Fully atomistic molecular dynamics simulations have been carried out to investigate the correlation of biological activity with dynamics of water molecules in an aqueous protein solution of the toxic domain of enterotoxin (PDB ID: 1ETN). This is a small protein of 13 amino acid residues. Our study of this water soluble protein clearly reveals that water dynamics slows down in the hydration layer. Despite this general slowing down, water molecules in the vicinity of the second beta turn of this protein exhibit faster dynamics than those near other regions of the protein. Since this beta turn is believed to play a critical role in the receptor binding of this protein, the faster dynamics of water near the beta turn m ay have biological significance. The collective orientational dynamics of the water molecules in the protein solution exhibits a characteristic long time component of 27 ps, which agrees well with dielectric relaxation experiments.

Insulin Growth Factor-2 Binding Protein 3 (IGF2BP3) Is a Glioblastoma-specific Marker That Activates Phosphatidylinositol 3-Kinase/Mitogen-activated Protein Kinase (PI3K/MAPK) Pathways by Modulating IGF-2

Glioblastoma is the most common and malignant form of primary astrocytoma. Upon investigation of the insulin-like growth factor (IGF) pathway, we found the IGF2BP3/IMP3 transcript and protein to be up-regulated in GBMs but not in lower grade astrocytomas (p<0.0001). IMP3 is an RNA binding protein known to bind to the 5'-untranslated region of IGF-2 mRNA, thereby activating its translation. Overexpression-and knockdown-based studies establish a role for IMP3 in promoting proliferation, anchorage-independent growth, invasion, and chemoresistance. IMP3 overexpressing B16F10 cells also showed increased tumor growth, angiogenesis, and metastasis, resulting in poor survival in a mouse model. Additionally, the infiltrating front, perivascular, and subpial regions in a majority of the GBMs stained positive for IMP3. Furthermore, two different murine glioma models were used to substantiate the above findings. In agreement with the translation activation functions of IMP3, we also found increased IGF-2 protein in the GBM tumor samples without a corresponding increase in its transcript levels. Also, in vitro IMP3 overexpression/knockdown modulated the IGF-2 protein levels without altering its transcript levels. Additionally, IGF-2 neutralization and supplementation studies established that the proproliferative effects of IMP3 were indeed mediated through IGF-2. Concordantly, PI3K and MAPK, the downstream effectors of IGF-2, are activated by IMP3 and are found to be essential for IMP3-induced cell proliferation. Thus, we have identified IMP3 as a GBM-specific proproliferative and proinvasive marker acting through IGF-2 resulting in the activation of oncogenic PI3K and MAPK pathways.

Heat shock protein 70 as a biomarker of heat stress in a simulated hot cockpit

Background: Fighter pilots are frequently exposed to high temperatures during high-speed low-level flight. Heat strain can result in temporary impairment of cognitive functions and when severe, loss of consciousness and consequent loss of life and equipment. Induction of stress proteins is a highly conserved stress response mechanism from bacteria to humans. induced stress protein levels are known to be cytoprotective and have been correlated with stress tolerance. Although many studies on the heat shock response mechanisms have been performed in cell culture and animal model systems, there is very limited information on stress protein induction in human subjects. Hypothesis: Heat shock proteins (Hsp), especially Hsp70, may be induced in human subjects exposed to high temperatures in a hot cockpit designed to simulate heat stress experienced in low flying sorties. Methods: Six healthy volunteers were subjected to heat stress at 55degreesC in a high temperature cockpit simulator for a period of 1 h at 30% humidity. Physiological parameters such as oral and skin temperatures, heart rate, and sweat rate were monitored regularly during this time. The level of Hsp70 in leukocytes was examined before and after the heat exposure in each subject. Conclusions: Hsp70 was found to be significantly induced in all the six subjects exposed to heat stress. The level of induced Hsp70 appears to correlate with other strain indicators such as accumulative circulatory strain and Craig's modified index. The usefulness of Hsp70 as a molecular marker of heat stress in humans is discussed.