3D protein structure matching by patch signatures

For determining functionality dependencies between two proteins, both represented as 3D structures, it is an essential condition that they have one or more matching structural regions called patches. As 3D structures for proteins are large, complex and constantly evolving, it is computationally expensive and very time-consuming to identify possible locations and sizes of patches for a given protein against a large protein database. In this paper, we address a vector space based representation for protein structures, where a patch is formed by the vectors within the region. Based on our previews work, a compact representation of the patch named patch signature is applied here. A similarity measure of two patches is then derived based on their signatures. To achieve fast patch matching in large protein databases, a match-and-expand strategy is proposed. Given a query patch, a set of small k-sized matching patches, called candidate patches, is generated in match stage. The candidate patches are further filtered by enlarging k in expand stage. Our extensive experimental results demonstrate encouraging performances with respect to this biologically critical but previously computationally prohibitive problem.

Overcoming barriers to membrane protein structure determination

After decades of slow progress, the pace of research on membrane protein structures is beginning to quicken thanks to various improvements in technology, including protein engineering and microfocus X-ray diffraction. Here we review these developments and, where possible, highlight generic new approaches to solving membrane protein structures based on recent technological advances. Rational approaches to overcoming the bottlenecks in the field are urgently required as membrane proteins, which typically comprise ~30% of the proteomes of organisms, are dramatically under-represented in the structural database of the Protein Data Bank.

Solving Maximum Clique Problem for Protein Structure Similarity

Computing the similarity between two protein structures is a crucial task in molecular biology, and has been extensively investigated. Many protein structure comparison methods can be modeled as maximum weighted clique problems in specific k-partite graphs, referred here as alignment graphs. In this paper we present both a new integer programming formulation for solving such clique problems and a dedicated branch and bound algorithm for solving the maximum cardinality clique problem. Both approaches have been integrated in VAST, a software for aligning protein 3D structures largely used in the National Center for Biotechnology Information, an original clique solver which uses the well known Bron and Kerbosch algorithm (BK). Our computational results on real protein alignment instances show that our branch and bound algorithm is up to 116 times faster than BK.

The accuracy of protein structure alignment servers

Background: Protein structural alignment is one of the most fundamental and crucial areas of research in the domain of computational structural biology. Comparison of a protein structure with known structures helps to classify it as a new or belonging to a known group of proteins. This, in turn, is useful to determine the function of protein, its evolutionary relationship with other protein molecules and grasping principles underlying protein architecture and folding. Results: A large number of protein structure alignment methods are available. Each protein structure alignment tool has its own strengths andweaknesses that need to be highlighted.We compared and presented results of six most popular and publically available servers for protein structure comparison. These web-based servers were compared with the respect to functionality (features provided by these servers) and accuracy (how well the structural comparison is performed). The CATH was used as a reference. The results showed that overall CE was top performer. DALI and PhyreStorm showed similar results whereas PDBeFold showed the lowest performance. In case of few secondary structural elements, CE, DALI and PhyreStorm gave 100% success rate. Conclusion: Overall none of the structural alignment servers showed 100% success rate. Studies of overall performance, effect of mainly alpha and effect of mainly beta showed consistent performance. CE, DALI, FatCat and PhyreStorm showed more than 90% success rate.

Oligomerization, Membrane Association, And In Vivo Phosphorylation Of Sugarcane Udp-glucose Pyrophosphorylase.

Sugarcane is a monocot plant that accumulates sucrose to levels of up to 50% of dry weight in the stalk. The mechanisms that are involved in sucrose accumulation in sugarcane are not well understood, and little is known with regard to factors that control the extent of sucrose storage in the stalks. UDP-glucose pyrophosphorylase (UGPase; EC 2.7.7.9) is an enzyme that produces UDP-glucose, a key precursor for sucrose metabolism and cell wall biosynthesis. The objective of this work was to gain insights into the ScUGPase-1 expression pattern and regulatory mechanisms that control protein activity. ScUGPase-1 expression was negatively correlated with the sucrose content in the internodes during development, and only slight differences in the expression patterns were observed between two cultivars that differ in sucrose content. The intracellular localization of ScUGPase-1 indicated partial membrane association of this soluble protein in both the leaves and internodes. Using a phospho-specific antibody, we observed that ScUGPase-1 was phosphorylated in vivo at the Ser-419 site in the soluble and membrane fractions from the leaves but not from the internodes. The purified recombinant enzyme was kinetically characterized in the direction of UDP-glucose formation, and the enzyme activity was affected by redox modification. Preincubation with H2O2 strongly inhibited this activity, which could be reversed by DTT. Small angle x-ray scattering analysis indicated that the dimer interface is located at the C terminus and provided the first structural model of the dimer of sugarcane UGPase in solution.

The fourth extracellular loop of the alpha subunit of Na,K-ATPase. Functional evidence for close proximity with the second extracellular loop.

Na,K-ATPase is the main active transport system that maintains the large gradients of Na(+) and K(+) across the plasma membrane of animal cells. The crystal structure of a K(+)-occluding conformation of this protein has been recently published, but the movements of its different domains allowing for the cation pumping mechanism are not yet known. The structure of many more conformations is known for the related calcium ATPase SERCA, but the reliability of homology modeling is poor for several domains with low sequence identity, in particular the extracellular loops. To better define the structure of the large fourth extracellular loop between the seventh and eighth transmembrane segments of the alpha subunit, we have studied the formation of a disulfide bond between pairs of cysteine residues introduced by site-directed mutagenesis in the second and the fourth extracellular loop. We found a specific pair of cysteine positions (Y308C and D884C) for which extracellular treatment with an oxidizing agent inhibited the Na,K pump function, which could be rapidly restored by a reducing agent. The formation of the disulfide bond occurred preferentially under the E2-P conformation of Na,K-ATPase, in the absence of extracellular cations. Using recently published crystal structure and a distance constraint reproducing the existence of disulfide bond, we performed an extensive conformational space search using simulated annealing and showed that the Tyr(308) and Asp(884) residues can be in close proximity, and simultaneously, the SYGQ motif of the fourth extracellular loop, known to interact with the extracellular domain of the beta subunit, can be exposed to the exterior of the protein and can easily interact with the beta subunit.

Regulation of Nedd4-2 self-ubiquitination and stability by a PY motif located within its HECT-domain.

Ubiquitin ligases play a pivotal role in substrate recognition and ubiquitin transfer, yet little is known about the regulation of their catalytic activity. Nedd4 (neural-precursor-cell-expressed, developmentally down-regulated 4)-2 is an E3 ubiquitin ligase composed of a C2 domain, four WW domains (protein-protein interaction domains containing two conserved tryptophan residues) that bind PY motifs (L/PPXY) and a ubiquitin ligase HECT (homologous with E6-associated protein C-terminus) domain. In the present paper we show that the WW domains of Nedd4-2 bind (weakly) to a PY motif (LPXY) located within its own HECT domain and inhibit auto-ubiquitination. Pulse-chase experiments demonstrated that mutation of the HECT PY-motif decreases the stability of Nedd4-2, suggesting that it is involved in stabilization of this E3 ligase. Interestingly, the HECT PY-motif mutation does not affect ubiquitination or down-regulation of a known Nedd4-2 substrate, ENaC (epithelial sodium channel). ENaC ubiquitination, in turn, appears to promote Nedd4-2 self-ubiquitination. These results support a model in which the inter- or intra-molecular WW-domain-HECT PY-motif interaction stabilizes Nedd4-2 by preventing self-ubiquitination. Substrate binding disrupts this interaction, allowing self-ubiquitination of Nedd4-2 and subsequent degradation, resulting in down-regulation of Nedd4-2 once it has ubiquitinated its target. These findings also point to a novel mechanism employed by a ubiquitin ligase to regulate itself differentially compared with substrate ubiquitination and stability.

HCF-1 is cleaved in the active site of O-GlcNAc transferase.

Host cell factor-1 (HCF-1), a transcriptional co-regulator of human cell-cycle progression, undergoes proteolytic maturation in which any of six repeated sequences is cleaved by the nutrient-responsive glycosyltransferase, O-linked N-acetylglucosamine (O-GlcNAc) transferase (OGT). We report that the tetratricopeptide-repeat domain of O-GlcNAc transferase binds the carboxyl-terminal portion of an HCF-1 proteolytic repeat such that the cleavage region lies in the glycosyltransferase active site above uridine diphosphate-GlcNAc. The conformation is similar to that of a glycosylation-competent peptide substrate. Cleavage occurs between cysteine and glutamate residues and results in a pyroglutamate product. Conversion of the cleavage site glutamate into serine converts an HCF-1 proteolytic repeat into a glycosylation substrate. Thus, protein glycosylation and HCF-1 cleavage occur in the same active site.

Essential role of CD8 palmitoylation in CD8 coreceptor function.

To investigate the molecular basis that makes heterodimeric CD8alphabeta a more efficient coreceptor than homodimeric CD8alphaalpha, we used various CD8 transfectants of T1.4 T cell hybridomas, which are specific for H-2Kd, and a photoreactive derivative of the Plasmodium berghei circumsporozoite peptide PbCS 252-260 (SYIPSAEKI). We demonstrate that CD8 is palmitoylated at the cytoplasmic tail of CD8beta and that this allows partitioning of CD8alphabeta, but not of CD8alphaalpha, in lipid rafts. Localization of CD8 in rafts is crucial for its coreceptor function. First, association of CD8 with the src kinase p56lck takes place nearly exclusively in rafts, mainly due to increased concentration of both components in this compartment. Deletion of the cytoplasmic domain of CD8beta abrogated localization of CD8 in rafts and association with p56lck. Second, CD8-mediated cross-linking of p56lck by multimeric Kd-peptide complexes or by anti-CD8 Ab results in p56lck activation in rafts, from which the abundant phosphatase CD45 is excluded. Third, CD8-associated activated p56lck phosphorylates CD3zeta in rafts and hence induces TCR signaling and T cell activation. This study shows that palmitoylation of CD8beta is required for efficient CD8 coreceptor function, mainly because it dramatically increases CD8 association with p56lck and CD8-mediated activation of p56lck in lipid rafts.

Sequence conservation in Plasmodium falciparum alpha-helical coiled coil domains proposed for vaccine development.

BACKGROUND: The availability of the P. falciparum genome has led to novel ways to identify potential vaccine candidates. A new approach for antigen discovery based on the bioinformatic selection of heptad repeat motifs corresponding to alpha-helical coiled coil structures yielded promising results. To elucidate the question about the relationship between the coiled coil motifs and their sequence conservation, we have assessed the extent of polymorphism in putative alpha-helical coiled coil domains in culture strains, in natural populations and in the single nucleotide polymorphism data available at PlasmoDB. METHODOLOGY/PRINCIPAL FINDINGS: 14 alpha-helical coiled coil domains were selected based on preclinical experimental evaluation. They were tested by PCR amplification and sequencing of different P. falciparum culture strains and field isolates. We found that only 3 out of 14 alpha-helical coiled coils showed point mutations and/or length polymorphisms. Based on promising immunological results 5 of these peptides were selected for further analysis. Direct sequencing of field samples from Papua New Guinea and Tanzania showed that 3 out of these 5 peptides were completely conserved. An in silico analysis of polymorphism was performed for all 166 putative alpha-helical coiled coil domains originally identified in the P. falciparum genome. We found that 82% (137/166) of these peptides were conserved, and for one peptide only the detected SNPs decreased substantially the probability score for alpha-helical coiled coil formation. More SNPs were found in arrays of almost perfect tandem repeats. In summary, the coiled coil structure prediction was rarely modified by SNPs. The analysis revealed a number of peptides with strictly conserved alpha-helical coiled coil motifs. CONCLUSION/SIGNIFICANCE: We conclude that the selection of alpha-helical coiled coil structural motifs is a valuable approach to identify potential vaccine targets showing a high degree of conservation.

The gamma-crystallins and human cataracts: a puzzle made clearer.

Despite the fact that cataracts constitute the leading cause of blindness worldwide, the mechanisms of lens opacification remain unclear. We recently mapped the aculeiform cataract to the gamma-crystallin locus (CRYG) on chromosome 2q33-35, and mutational analysis of the CRYG-genes cluster identified the aculeiform-cataract mutation in exon 2 of gamma-crystallin D (CRYGD). This mutation occurred in a highly conserved amino acid and could be associated with an impaired folding of CRYGD. During our study, we observed that the previously reported Coppock-like-cataract mutation, the first human cataract mutation, in the pseudogene CRYGE represented a polymorphism seen in 23% of our control population. Further analysis of the original Coppock-like-cataract family identified a missense mutation in a highly conserved segment of exon 2 of CRYGC. These mutations were not seen in a large control population. There is no direct evidence, to date, that up-regulation of a pseudogene causes cataracts. To our knowledge, these findings are the first evidence of an involvement of CRYGC and support the role of CRYGD in human cataract formation.

TWEAK-FN14 signaling induces lysosomal degradation of a cIAP1-TRAF2 complex to sensitize tumor cells to TNFalpha.

Synthetic inhibitor of apoptosis (IAP) antagonists induce degradation of IAP proteins such as cellular IAP1 (cIAP1), activate nuclear factor kappaB (NF-kappaB) signaling, and sensitize cells to tumor necrosis factor alpha (TNFalpha). The physiological relevance of these discoveries to cIAP1 function remains undetermined. We show that upon ligand binding, the TNF superfamily receptor FN14 recruits a cIAP1-Tnf receptor-associated factor 2 (TRAF2) complex. Unlike IAP antagonists that cause rapid proteasomal degradation of cIAP1, signaling by FN14 promotes the lysosomal degradation of cIAP1-TRAF2 in a cIAP1-dependent manner. TNF-like weak inducer of apoptosis (TWEAK)/FN14 signaling nevertheless promotes the same noncanonical NF-kappaB signaling elicited by IAP antagonists and, in sensitive cells, the same autocrine TNFalpha-induced death occurs. TWEAK-induced loss of the cIAP1-TRAF2 complex sensitizes immortalized and minimally passaged tumor cells to TNFalpha-induced death, whereas primary cells remain resistant. Conversely, cIAP1-TRAF2 complex overexpression limits FN14 signaling and protects tumor cells from TWEAK-induced TNFalpha sensitization. Lysosomal degradation of cIAP1-TRAF2 by TWEAK/FN14 therefore critically alters the balance of life/death signals emanating from TNF-R1 in immortalized cells.

Progesterone down-regulates the open probability of the amiloride-sensitive epithelial sodium channel via a Nedd4-2-dependent mechanism.

Activation of the mitogen-activated protein (MAP) kinase cascade by progesterone in Xenopus oocytes leads to a marked down-regulation of activity of the amiloride-sensitive epithelial sodium channel (ENaC). Here we have studied the signaling pathways involved in progesterone effect on ENaC activity. We demonstrate that: (i) the truncation of the C termini of the alphabetagammaENaC subunits results in the loss of the progesterone effect on ENaC; (ii) the effect of progesterone was also suppressed by mutating conserved tyrosine residues in the Pro-X-X-Tyr (PY) motif of the C termini of the beta and gamma ENaC subunits (beta(Y618A) and gamma(Y628A)); (iii) the down-regulation of ENaC activity by progesterone was also suppressed by co-expression ENaC subunits with a catalytically inactive mutant of Nedd4-2, a ubiquitin ligase that has been previously demonstrated to decrease ENaC cell-surface expression via a ubiquitin-dependent internalization/degradation mechanism; (iv) the effect of progesterone was significantly reduced by suppression of consensus sites (beta(T613A) and gamma(T623A)) for ENaC phosphorylation by the extracellular-regulated kinase (ERK), a MAP kinase previously shown to facilitate the binding of Nedd4 ubiquitin ligases to ENaC; (v) the quantification of cell-surface-expressed ENaC subunits revealed that progesterone decreases ENaC open probability (whole cell P(o), wcP(o)) and not its cell-surface expression. Collectively, these results demonstrate that the binding of active Nedd4-2 to ENaC is a crucial step in the mechanism of ENaC inhibition by progesterone. Upon activation of ERK, the effect of Nedd4-2 on ENaC open probability can become more important than its effect on ENaC cell-surface expression.

Structural prediction of peptides bound to MHC class I.

An ab initio structure prediction approach adapted to the peptide-major histocompatibility complex (MHC) class I system is presented. Based on structure comparisons of a large set of peptide-MHC class I complexes, a molecular dynamics protocol is proposed using simulated annealing (SA) cycles to sample the conformational space of the peptide in its fixed MHC environment. A set of 14 peptide-human leukocyte antigen (HLA) A0201 and 27 peptide-non-HLA A0201 complexes for which X-ray structures are available is used to test the accuracy of the prediction method. For each complex, 1000 peptide conformers are obtained from the SA sampling. A graph theory clustering algorithm based on heavy atom root-mean-square deviation (RMSD) values is applied to the sampled conformers. The clusters are ranked using cluster size, mean effective or conformational free energies, with solvation free energies computed using Generalized Born MV 2 (GB-MV2) and Poisson-Boltzmann (PB) continuum models. The final conformation is chosen as the center of the best-ranked cluster. With conformational free energies, the overall prediction success is 83% using a 1.00 Angstroms crystal RMSD criterion for main-chain atoms, and 76% using a 1.50 Angstroms RMSD criterion for heavy atoms. The prediction success is even higher for the set of 14 peptide-HLA A0201 complexes: 100% of the peptides have main-chain RMSD values < or =1.00 Angstroms and 93% of the peptides have heavy atom RMSD values < or =1.50 Angstroms. This structure prediction method can be applied to complexes of natural or modified antigenic peptides in their MHC environment with the aim to perform rational structure-based optimizations of tumor vaccines.