Structural and molecular basis of interaction of HCV non-structural protein 5A with human casein kinase 1 alpha and PKR

Background: Interaction of non-structural protein 5A (NS5A) of Hepatitis C virus (HCV) with human kinases namely, casein kinase 1 alpha (ck1 alpha) and protein kinase R (PKR) have different functional implications such as regulation of viral replication and evasion of interferon induced immune response respectively. Understanding the structural and molecular basis of interactions of the viral protein with two different human kinases can be useful in developing strategies for treatment against HCV. Results: Serine 232 of NS5A is known to be phosphorylated by human ck1 alpha. A structural model of NS5A peptide containing phosphoacceptor residue Serine 232 bound to ck1 alpha has been generated using the known 3-D structures of kinase-peptide complexes. The substrate interacting residues in ck1 alpha has been identified from the model and these are found to be conserved well in the ck1 family. ck1 alpha - substrate peptide complex has also been used to understand the structural basis of association between ck1 alpha and its other viral stress induced substrate, tumour suppressor p53 transactivation domain which has a crystal structure available. Interaction of NS5A with another human kinase PKR is primarily genotype specific. NS5A from genotype 1b has been shown to interact and inhibit PKR whereas NS5A from genotype 2a/3a are unable to bind and inhibit PKR efficiently. This is one of the main reasons for the varied response to interferon therapy in HCV patients across different genotypes. Using PKR crystal structure, sequence alignment and evolutionary trace analysis some of the critical residues responsible for the interaction of NS5A 1b with PKR have been identified. Conclusions: The substrate interacting residues in ck1 alpha have been identified using the structural model of kinase substrate peptide. The PKR interacting NS5A 1b residues have also been predicted using PKR crystal structure, NS5A sequence analysis along with known experimental results. Functional significance and nature of interaction of interferon sensitivity determining region and variable region 3 of NS5A in different genotypes with PKR which was experimentally shown are also supported by the findings of evolutionary trace analysis. Designing inhibitors to prevent this interaction could enable the HCV genotype 1 infected patients respond well to interferon therapy.

Site-specific N-Linked Glycosylation of Receptor Guanylyl Cyclase C Regulates Ligand Binding, Ligand-mediated Activation and Interaction with Vesicular Integral Membrane Protein 36, VIP36

Guanylyl cyclase C (GC-C) is a multidomain, membrane-associated receptor guanylyl cyclase. GC-C is primarily expressed in the gastrointestinal tract, where it mediates fluid-ion homeostasis, intestinal inflammation, and cell proliferation in a cGMP-dependent manner, following activation by its ligands guanylin, uroguanylin, or the heat-stable enterotoxin peptide (ST). GC-C is also expressed in neurons, where it plays a role in satiation and attention deficiency/hyperactive behavior. GC-C is glycosylated in the extracellular domain, and differentially glycosylated forms that are resident in the endoplasmic reticulum (130 kDa) and the plasma membrane (145 kDa) bind the ST peptide with equal affinity. When glycosylation of human GC-C was prevented, either by pharmacological intervention or by mutation of all of the 10 predicted glycosylation sites, ST binding and surface localization was abolished. Systematic mutagenesis of each of the 10 sites of glycosylation in GC-C, either singly or in combination, identified two sites that were critical for ligand binding and two that regulated ST-mediated activation. We also show that GC-C is the first identified receptor client of the lectin chaperone vesicular integral membrane protein, VIP36. Interaction with VIP36 is dependent on glycosylation at the same sites that allow GC-C to fold and bind ligand. Because glycosylation of proteins is altered in many diseases and in a tissue-dependent manner, the activity and/or glycan-mediated interactions of GC-C may have a crucial role to play in its functions in different cell types.

Computational studies on histidine kinase protein BaeS to target multidrug-resistant Salmonella

Multidrug-resistant Salmonella serovars have been a recent concern in curing infectious diseases like typhoid. Salmonella BaeS and BaeR are the two-component system (TCS) that signal transduction proteins found to play an important role in its multidrug resistance. A canonical TCS comprises a histidine kinase (HK) and its cognate partner response regulator (RR). The general approaches for therapeutic targeting are either the catalytic ATP-binding domain or the dimerization domain HisKA (DHp) of the HK, and in some cases, the receiver or the regulatory domain of the RR proteins. Earlier efforts of identifying novel drugs targeting the signal transduction protein have not been quite successful, as it shares similar ATP-binding domain with the key house keeping gene products of the mammalian GHL family. However, targeting the dimerization domain of HisKA through which the signals are received from the RR can be a better approach. In this article, we show stepwise procedure to specifically identify the key interacting residues involved in the dimerization with the RR along with effective targeting by ligands screened from the public database. We have found a few inhibitors which target effectively the important residues for the dimerization activity. Our results suggest a plausible de novo design of better DHp domain inhibitors.

Dramatic Structural Changes Resulting from the Loss of a Crucial Hydrogen Bond in the Hinge Region Involved in C-Terminal Helix Swapping in SurE: A Survival Protein from Salmonella typhimurium

Domain swapping is an interesting feature of some oligomeric proteins in which each protomer of the oligomer provides an identical surface for exclusive interaction with a segment or domain belonging to another protomer. Here we report results of mutagenesis experiments on the structure of C-terminal helix swapped dimer of a stationary phase survival protein from Salmonella typhimurium (StSurE). Wild type StSurE is a dimer in which a large helical segment at the C-terminus and a tetramerization loop comprising two beta strands are swapped between the protomers. Key residues in StSurE that might promote C-terminal helix swapping were identified by sequence and structural comparisons. Three mutants in which the helix swapping is likely to be avoided were constructed and expressed in E. coli. Three-dimensional X-ray crystal structures of the mutants H234A and D230A/H234A could be determined at 2.1 angstrom and 2.35 angstrom resolutions, respectively. Contrary to expectations, helix swapping was mostly retained in both the mutants. The loss of the crucial D230 OD2- H234 NE2 hydrogen bond (2.89 angstrom in the wild type structure) in the hinge region was compensated by new inter and intra-chain interactions. However, the two fold molecular symmetry was lost and there were large conformational changes throughout the polypeptide. In spite of these changes, the dimeric structure and an approximate tetrameric organization were retained, probably due to the interactions involving the tetramerization loop. Mutants were mostly functionally inactive, highlighting the importance of precise inter-subunit interactions for the symmetry and function of StSurE.

Tethering preferences of domain families co-occurring in multi-domain proteins

Genomic data of several organisms have revealed the presence of a vast repertoire of multi-domain proteins. The role played by individual domains in a multi-domain protein has a profound influence on the overall function of the protein. In the present analysis an attempt has been made to better understand the tethering preferences of domain families that occur in multi-domain proteins. The analysis has been carried out on an exhaustive dataset of 2 961 898 sequences of proteins from 930 organisms, where 741 274 proteins are comprised of at least two domain families. For every domain family, the number of other domain families with which it co-occurs within a protein in this dataset has been enumerated and is referred to as the tethering number of the domain family. It was found that, in the general dataset, the AAA ATPase family and the family of Ser/Thr kinases have the highest tethering numbers of 450 and 444 respectively. Further analysis reveals significant correlation between the number of members in a family and its tethering number. Positive correlation was also observed for the extent of a sequence and functional diversity within a family and the tethering numbers of domain families. Domain families that are present ubiquitously in diverse organisms tend to have large tethering numbers, while organism/kingdom-specific families have low tethering numbers. Thus, the analysis uncovers how domain families recombine and evolve to give rise to multi-domain proteins.

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 relationship between classification of multi-domain proteins using an alignment-free approach and their functions: a case study with immunoglobulins

Establishing functional relationships between multi-domain protein sequences is a non-trivial task. Traditionally, delineating functional assignment and relationships of proteins requires domain assignments as a prerequisite. This process is sensitive to alignment quality and domain definitions. In multi-domain proteins due to multiple reasons, the quality of alignments is poor. We report the correspondence between the classification of proteins represented as full-length gene products and their functions. Our approach differs fundamentally from traditional methods in not performing the classification at the level of domains. Our method is based on an alignment free local matching scores (LMS) computation at the amino-acid sequence level followed by hierarchical clustering. As there are no gold standards for full-length protein sequence classification, we resorted to Gene Ontology and domain-architecture based similarity measures to assess our classification. The final clusters obtained using LMS show high functional and domain architectural similarities. Comparison of the current method with alignment based approaches at both domain and full-length protein showed superiority of the LMS scores. Using this method we have recreated objective relationships among different protein kinase sub-families and also classified immunoglobulin containing proteins where sub-family definitions do not exist currently. This method can be applied to any set of protein sequences and hence will be instrumental in analysis of large numbers of full-length protein sequences.

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.

THE KNICKKOPF DOMON DOMAIN IS ESSENTIAL FOR CUTICLE DIFFERENTIATION IN Drosophila melanogaster

The dopamine monoxygenase N-terminal (DOMON) domain is found in extracellular proteins across several eukaryotic and prokaryotic taxa. It has been proposed that this domain binds to heme or sugar moieties. Here, we have analyzed the role of four highly conserved amino acids in the DOMON domain of the Drosophila melanogaster Knickkopf protein that is inserted into the apical plasma membrane and assists extracellular chitin organization. In principal, we generated Knickkopf versions with exchanged residues tryptophan(299,) methionine(333), arginine(401), or histidine(437), and scored for the ability of the respective engineered protein to normalize the knickkopf mutant phenotype. Our results confirm the absolute necessity of tryptophan(299,) methionine(333), and histidine(437) for Knickkopf function and stability, the latter two being predicted to be critical for heme binding. In contrast, arginine(401) is required for full efficiency of Knickkopf activity. Taken together, our genetic data support the prediction of these residues to mediate the function of Knickkopf during cuticle differentiation in insects. Hence, the DOMON domain is apparently an essential factor contributing to the construction of polysaccharide-based extracellular matrices.

Accurate prediction of interfacial residues in two-domain proteins using evolutionary information: Implications for three-dimensional modeling

With the preponderance of multidomain proteins in eukaryotic genomes, it is essential to recognize the constituent domains and their functions. Often function involves communications across the domain interfaces, and the knowledge of the interacting sites is essential to our understanding of the structure-function relationship. Using evolutionary information extracted from homologous domains in at least two diverse domain architectures (single and multidomain), we predict the interface residues corresponding to domains from the two-domain proteins. We also use information from the three-dimensional structures of individual domains of two-domain proteins to train naive Bayes classifier model to predict the interfacial residues. Our predictions are highly accurate (approximate to 85%) and specific (approximate to 95%) to the domain-domain interfaces. This method is specific to multidomain proteins which contain domains in at least more than one protein architectural context. Using predicted residues to constrain domain-domain interaction, rigid-body docking was able to provide us with accurate full-length protein structures with correct orientation of domains. We believe that these results can be of considerable interest toward rational protein and interaction design, apart from providing us with valuable information on the nature of interactions. Proteins 2014; 82:1219-1234. (c) 2013 Wiley Periodicals, Inc.

The HORMA domain: an evolutionarily conserved domain discovered in chromatin-associated proteins, has unanticipated diverse functions

The HORMA domain (for Hop1p, Rev7p and MAD2) was discovered in three chromatin-associated proteins in the budding yeast Saccharomyces cerevisiae. This domain has also been found in proteins with similar functions in organisms including plants, animals and nematodes. The HORMA domain containing proteins are thought to function as adaptors for meiotic checkpoint protein signaling and in the regulation of meiotic recombination. Surprisingly, new work has disclosed completely unanticipated and diverse functions for the HORMA domain containing proteins. A. M. Villeneuve and colleagues (Schvarzstein et al., 2013) show that meiosis-specific HORMA domain containing proteins plays a vital role in preventing centriole disengagement during Caenorhabditis elegans spermatocyte meiosis. Another recent study reveals that S. cerevisiae Atg13 HORMA domain acts as a phosphorylation-dependent conformational switch in the cellular autophagic process. (C) 2014 Elsevier B.V. All rights reserved.

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.

Groundnut Bud Necrosis Virus Encoded NSm Associates with Membranes via Its C-Terminal Domain

Groundnut Bud Necrosis Virus (GBNV) is a tripartite ambisense RNA plant virus that belongs to serogroup IV of Tospovirus genus. Non-Structural protein-m (NSm), which functions as movement protein in tospoviruses, is encoded by the M RNA. In this communication, we demonstrate that despite the absence of any putative transmembrane domain, GBNV NSm associates with membranes when expressed in E. coli as well as in N. benthamiana. Incubation of refolded NSm with liposomes ranging in size from 200-250 nm resulted in changes in the secondary and tertiary structure of NSm. A similar behaviour was observed in the presence of anionic and zwitterionic detergents. Furthermore, the morphology of the liposomes was found to be modified in the presence of NSm. Deletion of coiled coil domain resulted in the inability of in planta expressed NSm to interact with membranes. Further, when the C-terminal coiled coil domain alone was expressed, it was found to be associated with membrane. These results demonstrate that NSm associates with membranes via the C-terminal coiled coil domain and such an association may be important for movement of viral RNA from cell to cell.

Targeting Mycobacterium tuberculosis nucleoid-associated protein HU with structure-based inhibitors

The nucleoid-associated protein HU plays an important role in maintenance of chromosomal architecture and in global regulation of DNA transactions in bacteria. Although HU is essential for growth in Mycobacterium tuberculosis (Mtb), there have been no reported attempts to perturb HU function with small molecules. Here we report the crystal structure of the N-terminal domain of HU from Mtb. We identify a core region within the HU-DNA interface that can be targeted using stilbene derivatives. These small molecules specifically inhibit HU-DNA binding, disrupt nucleoid architecture and reduce Mtb growth. The stilbene inhibitors induce gene expression changes in Mtb that resemble those induced by HU deficiency. Our results indicate that HU is a potential target for the development of therapies against tuberculosis.

The Non-catalytic ``Cap Domain'' of a Mycobacterial Metallophosphoesterase Regulates Its Expression and Localization in the Cell

Despite highly conserved core catalytic domains, members of the metallophosphoesterase (MPE) superfamily perform diverse and crucial functions ranging from nucleotide and nucleic acid metabolism to phospholipid hydrolysis. Unique structural elements outside of the catalytic core called ``cap domains'' are thought to provide specialization to these enzymes; however, no directed study has been performed to substantiate this. The cap domain of Rv0805, an MPE from Mycobacterium tuberculosis, is located C-terminal to its catalytic domain and is dispensable for the catalytic activity of this enzyme in vitro. We show here that this C-terminal extension (CTE) mediates in vivo localization of the protein to the cell membrane and cell wall as well as modulates expression levels of Rv0805 in mycobacteria. We also demonstrate that Rv0805 interacts with the cell wall of mycobacteria, possibly with the mycolyl-arabinogalactan-peptidoglycan complex, by virtue of its C terminus, a hitherto unknown property of this MPE. Using a panel of mutant proteins, we identify interactions between active site residues of Rv0805 and the CTE that determine its association with the cell wall. Finally, we show that Rv0805 and a truncated mutant devoid of the CTE produce different phenotypic effects when expressed in mycobacteria. Our study thus provides a detailed dissection of the functions of the cap domain of an MPE and suggests that the repertoire of cellular functions of MPEs cannot be understood without exploring the modulatory effects of these subdomains.