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.

Identification of heat shock factor binding protein in Plasmodium falciparum

Background: Heat shock factor binding protein (HSBP) was originally discovered in a yeast two-hybrid screen as an interacting partner of heat shock factor (HSF). It appears to be conserved in all eukaryotes studied so far, with yeast being the only exception. Cell biological analysis of HSBP in mammals suggests its role as a negative regulator of heat shock response as it appears to interact with HSF only during the recovery phase following exposure to heat stress. While the identification of HSF in the malaria parasite is still eluding biologists, this study for the first time, reports the presence of a homologue of HSBP in Plasmodium falciparum. Methods: PfHSBP was cloned and purified as his-tag fusion protein. CD (Circular dichroism) spectroscopy was performed to predict the secondary structure. Immunoblots and immunofluorescence approaches were used to study expression and localization of HSBP in P. falciparum. Cellular fractionation was performed to examine subcellular distribution of PfHSBP. Immunoprecipitation was carried out to identify HSBP interacting partner in P. falciparum. Results: PfHSBP is a conserved protein with a high helical content and has a propensity to form homo-oligomers. PfHSBP was cloned, expressed and purified. The in vivo protein expression profile shows maximal expression in trophozoites. The protein was found to exist in oligomeric form as trimer and hexamer. PfHSBP is predominantly localized in the parasite cytosol, however, upon heat shock, it translocates to the nucleus. This study also reports the interaction of PfHSBP with PfHSP70-1 in the cytoplasm of the parasite. Conclusions: This study emphasizes the structural and biochemical conservation of PfHSBP with its mammalian counterpart and highlights its potential role in regulation of heat shock response in the malaria parasite. Analysis of HSBP may be an important step towards identification of the transcription factor regulating the heat shock response in P. falciparum.

PocketMatch (version 2.0): A parallel algorithm for the detection of structural similarities between protein ligand binding-sites

Knowledge of protein-ligand interactions is essential to understand several biological processes and important for applications ranging from understanding protein function to drug discovery and protein engineering. Here, we describe an algorithm for the comparison of three-dimensional ligand-binding sites in protein structures. A previously described algorithm, PocketMatch (version 1.0) is optimised, expanded, and MPI-enabled for parallel execution. PocketMatch (version 2.0) rapidly quantifies binding-site similarity based on structural descriptors such as residue nature and interatomic distances. Atomic-scale alignments may also be obtained from amino acid residue pairings generated. It allows an end-user to compute database-wide, all-to-all comparisons in a matter of hours. The use of our algorithm on a sample dataset, performance-analysis, and annotated source code is also included.

The Presence of Multiple Cellular Defects Associated with a Novel G50E Iron-Sulfur Cluster Scaffold Protein ( ISCU) Mutation Leads to Development of Mitochondrial Myopathy

Background: Muscle-specific deficiency of iron-sulfur (Fe-S) cluster scaffold protein (ISCU) leads to myopathy. Results: Cells carrying the myopathy-associated G50E ISCU mutation demonstrate impaired Fe-S cluster biogenesis and mitochondrial dysfunction. Conclusion: Reduced mitochondrial respiration as a result of diminished Fe-S cluster synthesis results in muscle weakness in myopathy patients. Significance: The molecular mechanism behind disease progression should provide invaluable information to combat ISCU myopathy. Iron-sulfur (Fe-S) clusters are versatile cofactors involved in regulating multiple physiological activities, including energy generation through cellular respiration. Initially, the Fe-S clusters are assembled on a conserved scaffold protein, iron-sulfur cluster scaffold protein (ISCU), in coordination with iron and sulfur donor proteins in human mitochondria. Loss of ISCU function leads to myopathy, characterized by muscle wasting and cardiac hypertrophy. In addition to the homozygous ISCU mutation (g.7044GC), compound heterozygous patients with severe myopathy have been identified to carry the c.149GA missense mutation converting the glycine 50 residue to glutamate. However, the physiological defects and molecular mechanism associated with G50E mutation have not been elucidated. In this report, we uncover mechanistic insights concerning how the G50E ISCU mutation in humans leads to the development of severe ISCU myopathy, using a human cell line and yeast as the model systems. The biochemical results highlight that the G50E mutation results in compromised interaction with the sulfur donor NFS1 and the J-protein HSCB, thus impairing the rate of Fe-S cluster synthesis. As a result, electron transport chain complexes show significant reduction in their redox properties, leading to loss of cellular respiration. Furthermore, the G50E mutant mitochondria display enhancement in iron level and reactive oxygen species, thereby causing oxidative stress leading to impairment in the mitochondrial functions. Thus, our findings provide compelling evidence that the respiration defect due to impaired biogenesis of Fe-S clusters in myopathy patients leads to manifestation of complex clinical symptoms.