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.

Comparative Study of Protein Unfolding in Aqueous Urea and Dimethyl Sulfoxide Solutions: Surface Polarity, Solvent Specificity, and Sequence of Secondary Structure Melting

Elucidation of possible pathways between folded (native) and unfolded states of a protein is a challenging task, as the intermediates are often hard to detect. Here, we alter the solvent environment in a controlled manner by choosing two different cosolvents of water, urea, and dimethyl sulfoxide (DMSO) and study unfolding of four different proteins to understand the respective sequence of melting by computer simulation methods. We indeed find interesting differences in the sequence of melting of alpha helices and beta sheets in these two solvents. For example, in 8 M urea solution, beta-sheet parts of a protein are found to unfold preferentially, followed by the unfolding of alpha helices. In contrast, 8 M DMSO solution unfolds alpha helices first, followed by the separation of beta sheets for the majority of proteins. Sequence of unfolding events in four different alpha/beta proteins and also in chicken villin head piece (HP-36) both in urea and DMSO solutions demonstrate that the unfolding pathways are determined jointly by relative exposure of polar and nonpolar residues of a protein and the mode of molecular action of a solvent on that protein.

Bi-Domain Protein Tyrosine Phosphatases Reveal an Evolutionary Adaptation to Optimize Signal Transduction

Significance: The bi-domain protein tyrosine phosphatases (PTPs) exemplify functional evolution in signaling proteins for optimal spatiotemporal signal transduction. Bi-domain PTPs are products of gene duplication. The catalytic activity, however, is often localized to one PTP domain. The inactive PTP domain adopts multiple functional roles. These include modulation of catalytic activity, substrate specificity, and stability of the bi-domain enzyme. In some cases, the inactive PTP domain is a receptor for redox stimuli. Since multiple bi-domain PTPs are concurrently active in related cellular pathways, a stringent regulatory mechanism and selective cross-talk is essential to ensure fidelity in signal transduction. Recent Advances: The inactive PTP domain is an activator for the catalytic PTP domain in some cases, whereas it reduces catalytic activity in other bi-domain PTPs. The relative orientation of the two domains provides a conformational rationale for this regulatory mechanism. Recent structural and biochemical data reveal that these PTP domains participate in substrate recruitment. The inactive PTP domain has also been demonstrated to undergo substantial conformational rearrangement and oligomerization under oxidative stress. Critical Issues and Future Directions: The role of the inactive PTP domain in coupling environmental stimuli with catalytic activity needs to be further examined. Another aspect that merits attention is the role of this domain in substrate recruitment. These aspects have been poorly characterized in vivo. These lacunae currently restrict our understanding of neo-functionalization of the inactive PTP domain in the bi-domain enzyme. It appears likely that more data from these research themes could form the basis for understanding the fidelity in intracellular signal transduction.

A Genetic Analysis of the Functional Interactions within Mycobacterium tuberculosis Single-Stranded DNA Binding Protein

Single-stranded DNA binding proteins (SSBs) are vital in all organisms. SSBs of Escherichia coli (EcoSSB) and Mycobacterium tuberculosis (MtuSSB) are homotetrameric. The N-terminal domains (NTD) of these SSBs (responsible for their tetramerization and DNA binding) are structurally well defined. However, their C-terminal domains (CTD) possess undefined structures. EcoSSB NTD consists of beta 1-beta 1'-beta 2-beta 3-alpha-beta 4-beta 45(1)-beta 45(2)-beta 5 secondary structure elements. MtuSSB NTD includes an additional beta-strand (beta 6) forming a novel hook-like structure. Recently, we observed that MtuSSB complemented an E. coli Delta ssb strain. However, a chimeric SSB (m beta 4-beta 5), wherein only the terminal part of NTD (beta 4-beta 5 region possessing L-45 loop) of EcoSSB was substituted with that from MtuSSB, failed to function in E. coli in spite of its normal DNA binding and oligomerization properties. Here, we designed new chimeras by transplanting selected regions of MtuSSB into EcoSSB to understand the functional significance of the various secondary structure elements within SSB. All chimeric SSBs formed homotetramers and showed normal DNA binding. The m beta 4-beta 6 construct obtained by substitution of the region downstream of beta 5 in m beta 4-beta 5 SSB with the corresponding region (beta 6) of MtuSSB complemented the E. coli strain indicating a functional interaction between the L-45 loop and the beta 6 strand of MtuSSB.

PROTEIN DISULFIDE ANALYSIS AND DESIGN

The significant contribution of naturally occurring disulfide bonds to protein stability has encouraged development of methods to engineer non-native disulfides in proteins. These have yielded mixed results. We summarize applications of the program MODIP for disulfide engineering. The program predicts sites in proteins where disulfides can be stably introduced. The program has also been used as an aid in conformational analysis of naturally occurring disulfides in a-helices, antiparallel and parallel beta-strands. Disulfides in a-helices occur only at N-termini, where the first cysteine residue is the N-cap residue of the helix. The disulfide occurs as a CXXC motif and can possess redox activity. In antiparallel beta-strands, disulfides occur exclusively at non-hydrogen bonded (NHB) registered pairs of antiparallel beta-sheets with only 1 known natural example occurring at a hydrogen bonded (HB) registered pair. Conformational analysis suggests that disulfides between HB residue pairs are under torsional strain. A similar analysis to characterize disulfides in parallel beta-strands was carried out. We observed that only 9 instances of cross-strand disulfides exist in a non-redundant dataset. Stereochemical analysis shows that while tbe chi(square) angles are similar to those of other disulfides, the chi(1) and chi(2) angles show more variation and that one of tbe strands is generally an edge strand.

Sparsely populated residue conformations in protein structures: Revisiting ``experimental'' Ramachandran maps

The Ramachandran map clearly delineates the regions of accessible conformational (phi-) space for amino acid residues in proteins. Experimental distributions of phi, values in high-resolution protein structures, reveal sparsely populated zones within fully allowed regions and distinct clusters in apparently disallowed regions. Conformational space has been divided into 14 distinct bins. Residues adopting these relatively rare conformations are presented and amino acid propensities for these regions are estimated. Inspection of specific examples in a completely arid, fully allowed region in the top left quadrant establishes that side-chain and backbone interactions may provide the energetic compensation necessary for populating this region of phi- space. Asn, Asp, and His residues showed the highest propensities in this region. The two distinct clusters in the bottom right quadrant which are formally disallowed on strict steric considerations correspond to the gamma turn (C7 axial) conformation (Bin 12) and the i + 1 position of Type II turns (Bin 13). Of the 516 non-Gly residues in Bin 13, 384 occupied the i + 1 position of Type II turns. Further examination of these turn segments revealed a high propensity to occur at the N-terminus of helices and as a tight turn in hairpins. The strand-helix motif with the Type II turn as a connecting element was also found in as many as 57 examples. Proteins 2014; 82:1101-1112. (c) 2013 Wiley Periodicals, Inc.

Trans-spliced Heat Shock Protein 90 Modulates Encystation in Giardia lamblia

Background: Hsp90 from Giardia lamblia is expressed by splicing of two independently transcribed RNA molecules, coded by genes named HspN and HspC located 777 kb apart. The reasons underlying such unique trans-splicing based generation of GlHsp90 remain unclear. Principle Finding: In this study using mass-spectrometry we identify the sequence of the unique, junctional peptide contributed by the 5' UTR of HspC ORF. This peptide is critical for the catalytic function of Hsp90 as it harbours an essential ``Arg'' in its sequence. We also show that full length GlHsp90 possesses all the functional hall marks of a canonical Hsp90 including its ability to bind and hydrolyze ATP. Using qRT-PCR as well as western blotting approach we find the reconstructed Hsp90 to be induced in response to heat shock. On the contrary we find GlHsp90 to be down regulated during transition from proliferative trophozoites to environmentally resistant cysts. This down regulation of GlHsp90 appears to be mechanistically linked to the encystation process as we find pharmacological inhibition of GlHsp90 function to specifically induce encystation. Significance: Our results implicate the trans-spliced GlHsp90 from Giardia lamblia to regulate an essential stage transition in the life cycle of this important human parasite.

Suramin is a potent and selective inhibitor of Mycobacterium tuberculosis RecA protein and the SOS response: RecA as a potential target for antibacterial drug discovery

In eubacteria, RecA is essential for recombinational DNA repair and for stalled replication forks to resume DNA synthesis. Recent work has implicated a role for RecA in the development of antibiotic resistance in pathogenic bacteria. Consequently, our goal is to identify and characterize small-molecule inhibitors that target RecA both in vitro and in vivo. We employed ATPase, DNA strand exchange and LexA cleavage assays to elucidate the inhibitory effects of suramin on Mycobacterium tuberculosis RecA. To gain insights into the mechanism of suramin action, we directly visualized the structure of RecA nucleoprotein filaments by atomic force microscopy. To determine the specificity of suramin action in vivo, we investigated its effect on the SOS response by pull-down and western blot assays as well as for its antibacterial activity. We show that suramin is a potent inhibitor of DNA strand exchange and ATPase activities of bacterial RecA proteins with IC50 values in the low micromolar range. Additional evidence shows that suramin inhibits RecA-catalysed proteolytic cleavage of the LexA repressor. The mechanism underlying such inhibitory actions of suramin involves its ability to disassemble RecA-single-stranded DNA filaments. Notably, suramin abolished ciprofloxacin-induced recA gene expression and the SOS response and augmented the bactericidal action of ciprofloxacin. Our findings suggest a strategy to chemically disrupt the vital processes controlled by RecA and hence the promise of small molecules for use against drug-susceptible as well as drug-resistant strains of M. tuberculosis for better infection control and the development of new therapies.

The Glycosylated Rv1860 Protein of Mycobacterium tuberculosis Inhibits Dendritic Cell Mediated TH1 and TH17 Polarization of T Cells and Abrogates Protective Immunity Conferred by BCG

We previously reported interferon gamma secretion by human CD4(+) and CD8(+) T cells in response to recombinant E. coli-expressed Rv1860 protein of Mycobacterium tuberculosis (MTB) as well as protection of guinea pigs against a challenge with virulent MTB following prime-boost immunization with DNA vaccine and poxvirus expressing Rv1860. In contrast, a Statens Serum Institute Mycobacterium bovis BCG (BCG-SSI) recombinant expressing MTB Rv1860 (BCG-TB1860) showed loss of protective ability compared to the parent BCG strain expressing the control GFP protein (BCG-GFP). Since Rv1860 is a secreted mannosylated protein of MTB and BCG, we investigated the effect of BCG-TB1860 on innate immunity. Relative to BCG-GFP, BCG-TB1860 effected a significant near total reduction both in secretion of cytokines IL-2, IL-12p40, IL-12p70, TNF-alpha, IL-6 and IL-10, and up regulation of co-stimulatory molecules MHC-II, CD40, CD54, CD80 and CD86 by infected bone marrow derived dendritic cells (BMDC), while leaving secreted levels of TGF-beta unchanged. These effects were mimicked by BCG-TB1860His which carried a 6-Histidine tag at the C-terminus of Rv1860, killed sonicated preparations of BCG-TB1860 and purified H37Rv-derived Rv1860 glycoprotein added to BCG-GFP, but not by E. coli-expressed recombinant Rv1860. Most importantly, BMDC exposed to BCG-TB1860 failed to polarize allogeneic as well as syngeneic T cells to secrete IFN-gamma and IL-17 relative to BCG-GFP. Splenocytes from mice infected with BCG-SSI showed significantly less proliferation and secretion of IL-2, IFN-gamma and IL-17, but secreted higher levels of IL-10 in response to in vitro restimulation with BCG-TB1860 compared to BCG-GFP. Spleens from mice infected with BCG-TB1860 also harboured significantly fewer DC expressing MHC-II, IL-12, IL-2 and TNF-alpha compared to mice infected with BCG-GFP. Glycoproteins of MTB, through their deleterious effects on DC may thus contribute to suppress the generation of a TH1- and TH17-dominated adaptive immune response that is vital for protection against tuberculosis.

Direct control of type IIA topoisomerase activity by a chromosomally encoded regulatory protein

Precise control of supercoiling homeostasis is critical to DNA-dependent processes such as gene expression, replication, and damage response. Topoisomerases are central regulators of DNA supercoiling commonly thought to act independently in the recognition and modulation of chromosome superstructure; however, recent evidence has indicated that cells tightly regulate topoisomerase activity to support chromosome dynamics, transcriptional response, and replicative events. How topoisomerase control is executed and linked to the internal status of a cell is poorly understood. To investigate these connections, we determined the structure of Escherichia coil gyrase, a type HA topoisomerase bound to YacG, a recently identified chromosomally encoded inhibitor protein. Phylogenetic analyses indicate that YacG is frequently associated with coenzyme A (CoA) production enzymes, linking the protein to metabolism and stress. The structure, along with supporting solution studies, shows that YacG represses gyrase by sterically occluding the principal DNA-binding site of the enzyme. Unexpectedly, YacG acts by both engaging two spatially segregated regions associated with small-molecule inhibitor interactions (fluoroquinolone antibiotics and the newly reported antagonist GSK299423) and remodeling the gyrase holo enzyme into an inactive, ATP-trapped configuration. This study establishes a new mechanism for the protein-based control of topoisomerases, an approach that may be used to alter supercoiling levels for responding to changes in cellular state.

HupB, a Nucleoid-Associated Protein of Mycobacterium tuberculosis, Is Modified by Serine/Threonine Protein Kinases In Vivo

HU, a widely conserved bacterial histone-like protein, regulates many genes, including those involved in stress response and virulence. Whereas ample data are available on HU-DNA communication, the knowledge on how HU perceives a signal and transmit it to DNA remains limited. In this study, we identify HupB, the HU homolog of the human pathogen Mycobacterium tuberculosis, as a component of serine/threonine protein kinase (STPK) signaling. HupB is extracted in its native state from the exponentially growing cells of M. tuberculosis H37Ra and is shown to be phosphorylated on both serine and threonine residues. The STPKs capable of modifying HupB are determined in vitro and the residues modified by the STPKs are identified for both in vivo and the in vitro proteins through mass spectrometry. Of the identified phosphosites, Thr(65) and Thr(74) in the DNA-embracing beta-strand of the N-terminal domain of HupB (N-HupB) are shown to be crucial for its interaction with DNA. In addition, Arg(55) is also identified as an important residue for N-HupB-DNA interaction. N-HupB is shown to have a diminished interaction with DNA after phosphorylation. Furthermore, hupB is shown to be maximally expressed during the stationary phase in M. tuberculosis H37Ra, while HupB kinases were found to be constitutively expressed (PknE and PknF) or most abundant during the exponential phase (PknB). In conclusion, HupB, a DNA-binding protein, with an ability to modulate chromatin structure is proposed to work in a growth-phase-dependent manner through its phosphorylation carried out by the mycobacterial STPKs.