mulPBA: an efficient multiple protein structure alignment method based on a structural alphabet

The increasing number of available protein structures requires efficient tools for multiple structure comparison. Indeed, multiple structural alignments are essential for the analysis of function, evolution and architecture of protein structures. For this purpose, we proposed a new web server called multiple Protein Block Alignment (mulPBA). This server implements a method based on a structural alphabet to describe the backbone conformation of a protein chain in terms of dihedral angles. This sequence-like' representation enables the use of powerful sequence alignment methods for primary structure comparison, followed by an iterative refinement of the structural superposition. This approach yields alignments superior to most of the rigid-body alignment methods and highly comparable with the flexible structure comparison approaches. We implement this method in a web server designed to do multiple structure superimpositions from a set of structures given by the user. Outputs are given as both sequence alignment and superposed 3D structures visualized directly by static images generated by PyMol or through a Jmol applet allowing dynamic interaction. Multiple global quality measures are given. Relatedness between structures is indicated by a distance dendogram. Superimposed structures in PDB format can be also downloaded, and the results are quickly obtained. mulPBA server can be accessed at www.dsimb.inserm.fr/dsimb_tools/mulpba/.

Metal nanoparticles as better protein-aggregation prevention agents than chaperones

Prevention or suppression of protein aggregation is of great importance in the context of protein storage, transportation and delivery. Traditionally chaperones or other chemically active agents are used to stop or diffuse native protein aggregation. We have used gold nanoparticles to prevent thermal aggregation of alcohol dehydrogenase (ADH), a protein that maintains the alcohol level in the liver and stomach. A light-scattering assay has been used to investigate the effect of gold nanoparticles on thermal aggregation of ADH and the result of our study has been summarized in Fig. 1. The scattered light intensity from the solution containing ADH decreases when 45 nm gold nanoparticles are added prior to heating (thermal denaturation) the solution, which indicates prevention of aggregation. The aggregation of the protein is suppressed to the extent of 96% with picomolar concentration of 45 nm gold nanoparticles while micromolar amounts of other proteins and biological substances are necessary to achieve the same effect. The extent varies with the size and the concentration of the gold NPs for the same protein concentration.

Molecular Dynamics Perspective on the Protein Thermal Stability: A Case Study Using SAICAR Synthetase

The enzyme SAICAR synthetase ligates aspartate with CAIR (5'-phosphoribosyl-4-carboxy-5-aminoimidazole) forming SAICAR (5-amino-4-imidazole-N-succinocarboxamide ribonucleotide) in the presence of ATP. In continuation with our previous study on the thermostability of this enzyme in hyper-/thermophiles based on the structural aspects, here, we present the dynamic aspects that differentiate the mesophilic (E. coli, E. chaffeensis), thermophilic (G. kaustophilus), and hyperthermophilic (M. jannaschii, P. horikoshii) SAICAR synthetases by carrying out a total of 11 simulations. The five functional dimers from the above organisms were simulated using molecular dynamics for a period of 50 ns each at 300 K, 363 K, and an additional simulation at 333 K for the thermophilic protein. The basic features like root-mean-square deviations, root-mean-square fluctuations, surface accessibility, and radius of gyration revealed the instability of mesophiles at 363 K. Mean square displacements establish the reduced flexibility of hyper-/thermophiles at all temperatures. At the simulations time scale considered here, the long-distance networks are considerably affected in mesophilic structures at 363 K. In mesophiles, a comparatively higher number of short-lived (having less percent existence time) C alpha, hydrogen bonds, hydrophobic interactions are formed, and long-lived (with higher percentage existence time) contacts are lost. The number of time-averaged salt-bridges is at least 2-fold higher in hyperthermophiles at 363 K. The change in surface accessibility of salt-bridges at 363 K from 300 K is nearly doubled in mesophilic protein compared to proteins from other temperature classes.

Investigations of Ramachandran disallowed conformations in protein domain families

In peptide and protein structures, occurrence of (phi,psi.) angles in the disallowed region of the Ramachandran map almost always suggests local regions of error or poor accuracy. However, very rarely genuine disallowed conformations occur as noted in the current study in proteins of known structure available at ultra-high resolution (<= 1.2 (A) over circle). In the current work, extent of conservation of genuine disallowed conformations in homologous proteins of known structures has been analyzed. From a dataset of 124 protein domain families, with structure of at least one constituent member in each family available at a resolution of 1.2 (A) over circle or better, we have analyzed the conservation of 221 disallowed conformations. It is observed that the disallowed conformation is only moderately conservedin protein domain families. In the gross dataset no particular residue type adopting disallowed conformation elicit high conservation of residue type though there are alignment positions in the dataset with complete conservation of both the residue type and the disallowed conformation. Conserved disallowed conformation in protein domain families play biologically significant role in roughly 50% of the cases. The residues with the disallowed conformation or its flanking residues are often located within or around the functional site of the protein. (C) 2013 Elsevier B.V. All rights reserved.

Insulin like growth factor binding protein 4 promotes GBM progression and regulates key factors involved in EMT and invasion

Insulin like growth factor binding protein 4 (IGFBP4) regulates growth and development of tissues and organs by negatively regulating IGF signaling. Among most cancers, IGFBP4 has growth inhibitory role and reported as a down-regulated gene, except for renal cell carcinoma, wherein IGFBP4 promotes tumor progression. IGFBP4 expression has been shown to be higher in increasing grades of astrocytoma. However, the functional role of IGFBP4 in gliomas has not been explored. Surgical biopsies of 20 normal brain and 198 astrocytoma samples were analyzed for IGFBP4 expression by qRT-PCR. Highest expression of IGFBP4 mRNA was seen in GBM tumors compared to control brain tissues (median log2 of 2.035, p < 0.0001). Immunohistochemical analysis of 53 tissue samples revealed predominant nuclear staining of IGFBP4, seen maximally in GBMs when compared to DA and AA tumors (median LI = 29.12 +/- A 16.943, p < 0.001). Over expression of IGFBP4 in U343 glioma cells resulted in up-regulation of molecules involved in tumor growth, EMT and invasion such as pAkt, pErk, Vimentin, and N-cadherin and down-regulation of E-cadherin. Functionally, IGFBP4 over expression in these cells resulted in increased proliferation, migration and invasion as assessed by MTT, transwell migration, and Matrigel invasion assays. These findings were confirmed upon IGFBP4 knockdown in U251 glioma cells. Our data suggest a pro-tumorigenic role for IGFBP4 in glioma.

Transcriptional Repression of Tumor Suppressor CDC73, Encoding an RNA Polymerase II Interactor, by Wilms Tumor 1 Protein (WT1) Promotes Cell Proliferation IMPLICATION FOR CANCER THERAPEUTICS

The Wilms tumor 1 gene (WT1) can either repress or induce the expression of genes. Inconsistent with its tumor suppressor role, elevated WT1 levels have been observed in leukemia and solid tumors. WT1 has also been suggested to act as an oncogene by inducing the expression of MYC and BCL-2. However, these are only the correlational studies, and no functional study has been performed to date. Consistent with its tumor suppressor role, CDC73 binds to RNA polymerase II as part of a PAF1 transcriptional regulatory complex and causes transcriptional repression of oncogenes MYC and CCND1. It also represses beta-catenin-mediated transcription. Based on the reduced level of CDC73 in oral squamous cell carcinoma (OSCC) samples in the absence of loss-of-heterozygosity, promoter methylation, and mutations, we speculated that an inhibitory transcription factor is regulating its expression. The bioinformatics analysis predicted WT1 as an inhibitory transcription factor to regulate the CDC73 level. Our results showed that overexpression of WT1 decreased CDC73 levels and promoted proliferation of OSCC cells. ChIP and EMSA results demonstrated binding of WT1 to the CDC73 promoter. The 5-azacytidine treatment of OSCC cells led to an up-regulation of WT1 with a concomitant down-regulation of CDC73, further suggesting regulation of CDC73 by WT1. Exogenous CDC73 attenuated the protumorigenic activity of WT1 by apoptosis induction. An inverse correlation between expression levels of CDC73 and WT1 was observed in OSCC samples. These observations indicated that WT1 functions as an oncogene by repressing the expression of CDC73 in OSCC. We suggest that targeting WT1 could be a therapeutic strategy for cancer, including OSCC.

Filling-in Void and Sparse Regions in Protein Sequence Space by Protein-Like Artificial Sequences Enables Remarkable Enhancement in Remote Homology Detection Capability

Protein functional annotation relies on the identification of accurate relationships, sequence divergence being a key factor. This is especially evident when distant protein relationships are demonstrated only with three-dimensional structures. To address this challenge, we describe a computational approach to purposefully bridge gaps between related protein families through directed design of protein-like ``linker'' sequences. For this, we represented SCOP domain families, integrated with sequence homologues, as multiple profiles and performed HMM-HMM alignments between related domain families. Where convincing alignments were achieved, we applied a roulette wheel-based method to design 3,611,010 protein-like sequences corresponding to 374 SCOP folds. To analyze their ability to link proteins in homology searches, we used 3024 queries to search two databases, one containing only natural sequences and another one additionally containing designed sequences. Our results showed that augmented database searches showed up to 30% improvement in fold coverage for over 74% of the folds, with 52 folds achieving all theoretically possible connections. Although sequences could not be designed between some families, the availability of designed sequences between other families within the fold established the sequence continuum to demonstrate 373 difficult relationships. Ultimately, as a practical and realistic extension, we demonstrate that such protein-like sequences can be ``plugged-into'' routine and generic sequence database searches to empower not only remote homology detection but also fold recognition. Our richly statistically supported findings show that complementary searches in both databases will increase the effectiveness of sequence-based searches in recognizing all homologues sharing a common fold. (C) 2013 Elsevier Ltd. All rights reserved.

The N-terminal region containing the zinc finger domain of tobacco streak virus coat protein is essential for the formation of virus-like particles

Tobacco streak virus (TSV), a member of the genus Ilarvirus (family Bromoviridae), has a tripartite genome and forms quasi-isometric virions. All three viral capsids, encapsidating RNA 1, RNA 2 or RNA 3 and subgenomic RNA 4, are constituted of a single species of coat protein (CP). Formation of virus-like particles (VLPs) could be observed when the TSV CP gene was cloned and the recombinant CP (rCP) was expressed in E. coli. TSV VLPs were found to be stabilized by Zn2+ ions and could be disassembled in the presence of 500 mM CaCl2. Mutational analysis corroborated previous studies that showed that an N-terminal arginine-rich motif was crucial for RNA binding; however, the results presented here demonstrate that the presence of RNA is not a prerequisite for assembly of TSV VLPs. Instead, the N-terminal region containing the zinc finger domain preceding the arginine-rich motif is essential for assembly of these VLPs.

Evaluation of Brugia malayi sheath protein (Shp-1) as a diagnostic antigen for human lymphatic filariasis

Lymphatic filariasis is the second leading cause of permanent long-term disability globally and control of this disease needs efficient diagnostic methods. In this study, abundantly expressing microfilarial sheath protein (Shp-1) from Brugia malayi was characterized as a filarial diagnostic candidate using samples from different clinical population. Monoclonal antibodies were developed against E. coil expressed recombinant Shp-1 in order to assess its efficiency in filarial antigen detection assay system. Endemic Normal (EN, n = 170), asymptomatic microfilaeremics (MF, n = 65), symptomatic chronic pathology (CP, n = 45) and non endemic normal (NEN, n = 10) sera were analyzed by antigen capture enzyme-linked immunosorbent assay. Of the 290 individuals, all MF individuals (both brugian and bancroftian) were positive in this assay followed by CP and EN. When compared with SXP-1 and Og4C3 antigen assays, all assays detected Wb MF correctly, Bm MF was detected by Shp-1 and SXP-1 assays, and only Shp-1 was able to detect EN (12%) and CP (29%). Results showed that this assay may be useful for monitoring prior to mass drug administration. (c) 2014 Elsevier Inc. All rights reserved.

Specific Sequence of a Beta Turn in Human La Protein May Contribute to Species Specificity of Hepatitis C Virus

Human La protein is known to be an essential host factor for translation and replication of hepatitis C virus (HCV) RNA. Previously, we have demonstrated that residues responsible for interaction of human La protein with the HCV internal ribosomal entry site (IRES) around the initiator AUG within stem-loop IV form a beta-turn in the RNA recognition motif (RRM) structure. In this study, sequence alignment and mutagenesis suggest that the HCV RNA-interacting beta-turn is conserved only in humans and chimpanzees, the species primarily known to be infected by HCV. A 7-mer peptide corresponding to the HCV RNA-interacting region of human La inhibits HCV translation, whereas another peptide corresponding to the mouse La sequence was unable to do so. Furthermore, IRES-mediated translation was found to be significantly high in the presence of recombinant human La protein in vitro in rabbit reticulocyte lysate. We observed enhanced replication with HCV subgenomic and full-length replicons upon overexpression of either human La protein or a chimeric mouse La protein harboring a human La beta-turn sequence in mouse cells. Taken together, our results raise the possibility of creating an immunocompetent HCV mouse model using human-specific cell entry factors and a humanized form of La protein.

Chaperoning mitochondrial permeability transition: regulation of transition pore complex by a J-protein, DnaJC15

Mitochondria have a central role in the intrinsic pathway of apoptosis and involve activation of several transmembrane channels leading to release of death factors. Reduced expression of a mitochondrial J-protein DnaJC15 was associated with the development of chemoresistance in ovarian cancer cells. DnaJC15 was found to be a part of mitochondrial protein-transport machinery, though its connection with cell death mechanisms is still unclear. In the present study, we have provided evidence towards a novel function of DnaJC15 in regulation of mitochondrial permeability transition pore (MPTP) complex in normal and cancer cells. Overexpression of DnaJC15 resulted in MPTP opening and induction of apoptosis, whereas reduced amount of protein suppressed MPTP activation, upon cisplatin treatment. DnaJC15 was found to exert its proapoptotic function through the essential component of MPTP, cyclophilin D (CypD). Our results reveal a specific role of DnaJC15 in recruitment and coupling of CypD with mitochondrial permeability transition. In summary, our analysis provides first-time insights on the functional connection between mitochondrial inner membrane protein translocation machinery-associated J-protein DnaJC15 and regulation of cell death pathways.

A Novel C-Terminal Homologue of Aha1 Co-Chaperone Binds to Heat Shock Protein 90 and Stimulates Its ATPase Activity in Entamoeba histolytica

Cytosolic heat shock protein 90 (Hsp90) has been shown to be essential for many infectious pathogens and is considered a potential target for drug development. In this study, we have carried out biochemical characterization of Hsp90 from a poorly studied protozoan parasite of clinical importance, Entamoeba histolytica. We have shown that Entamoeba Hsp90 can bind to both ATP and its pharmacological inhibitor, 17-AAG (17-allylamino-17-demethoxygeldanamycin), with K-d values of 365.2 and 10.77 mu M, respectively, and it has a weak ATPase activity with a catalytic efficiency of 4.12 x 10(-4) min(-1) mu M-1. Using inhibitor 17-AAG, we have shown dependence of Entamoeba on Hsp90 for its growth and survival. Hsp90 function is regulated by various co-chaperones. Previous studies suggest a lack of several important co-chaperones in E. histolytica. In this study, we describe the presence of a novel homologue of co-chaperone Aha1 (activator of Hsp90 ATPase), EhAha1c, lacking a canonical Aha1 N-terminal domain. We also show that EhAha1c is capable of binding and stimulating ATPase activity of EhHsp90. In addition to highlighting the potential of Hsp90 inhibitors as drugs against amoebiasis, our study highlights the importance of E. histolytica in understanding the evolution of Hsp90 and its co-chaperone repertoire. (C) 2014 Elsevier Ltd. All rights reserved.

PLIC: protein-ligand interaction clusters

Most of the biological processes are governed through specific protein-ligand interactions. Discerning different components that contribute toward a favorable protein-ligand interaction could contribute significantly toward better understanding protein function, rationalizing drug design and obtaining design principles for protein engineering. The Protein Data Bank (PDB) currently hosts the structure of similar to 68 000 protein-ligand complexes. Although several databases exist that classify proteins according to sequence and structure, a mere handful of them annotate and classify protein-ligand interactions and provide information on different attributes of molecular recognition. In this study, an exhaustive comparison of all the biologically relevant ligand-binding sites (84 846 sites) has been conducted using PocketMatch: a rapid, parallel, in-house algorithm. PocketMatch quantifies the similarity between binding sites based on structural descriptors and residue attributes. A similarity network was constructed using binding sites whose PocketMatch scores exceeded a high similarity threshold (0.80). The binding site similarity network was clustered into discrete sets of similar sites using the Markov clustering (MCL) algorithm. Furthermore, various computational tools have been used to study different attributes of interactions within the individual clusters. The attributes can be roughly divided into (i) binding site characteristics including pocket shape, nature of residues and interaction profiles with different kinds of atomic probes, (ii) atomic contacts consisting of various types of polar, hydrophobic and aromatic contacts along with binding site water molecules that could play crucial roles in protein-ligand interactions and (iii) binding energetics involved in interactions derived from scoring functions developed for docking. For each ligand-binding site in each protein in the PDB, site similarity information, clusters they belong to and description of site attributes are provided as a relational database-protein-ligand interaction clusters (PLIC).

Predicting gene ontology annotations of orphan GWAS genes using protein-protein interactions

Background: The number of genome-wide association studies (GWAS) has increased rapidly in the past couple of years, resulting in the identification of genes associated with different diseases. The next step in translating these findings into biomedically useful information is to find out the mechanism of the action of these genes. However, GWAS studies often implicate genes whose functions are currently unknown; for example, MYEOV, ANKLE1, TMEM45B and ORAOV1 are found to be associated with breast cancer, but their molecular function is unknown. Results: We carried out Bayesian inference of Gene Ontology (GO) term annotations of genes by employing the directed acyclic graph structure of GO and the network of protein-protein interactions (PPIs). The approach is designed based on the fact that two proteins that interact biophysically would be in physical proximity of each other, would possess complementary molecular function, and play role in related biological processes. Predicted GO terms were ranked according to their relative association scores and the approach was evaluated quantitatively by plotting the precision versus recall values and F-scores (the harmonic mean of precision and recall) versus varying thresholds. Precisions of similar to 58% and similar to 40% for localization and functions respectively of proteins were determined at a threshold of similar to 30 (top 30 GO terms in the ranked list). Comparison with function prediction based on semantic similarity among nodes in an ontology and incorporation of those similarities in a k nearest neighbor classifier confirmed that our results compared favorably. Conclusions: This approach was applied to predict the cellular component and molecular function GO terms of all human proteins that have interacting partners possessing at least one known GO annotation. The list of predictions is available at http://severus.dbmi.pitt.edu/engo/GOPRED.html. We present the algorithm, evaluations and the results of the computational predictions, especially for genes identified in GWAS studies to be associated with diseases, which are of translational interest.