Refinement by shifting secondary structure elements improves sequence alignments

PC_ali: a tool for improved multiple alignments and evolutionary inference based on a hybrid protein sequence and structure similarity score

MOTIVATION

Evolutionary inference depends crucially on the quality of multiple sequence alignments (MSA), which is problematic for distantly related proteins. Since protein structure is more conserved than sequence, it seems natural to use structure alignments for distant homologs. However, structure alignments may not be suitable for inferring evolutionary relationships.

RESULTS

Here we examined four protein similarity measures that depend on sequence and structure (fraction of aligned residues, sequence identity, fraction of superimposed residues, and contact overlap), finding that they are intimately correlated but none of them provides a complete and unbiased picture of conservation in proteins. Therefore, we propose the new hybrid protein sequence and structure similarity score PC_sim based on their main principal component. The corresponding divergence measure PC_div shows the strongest correlation with divergences obtained from individual similarities, suggesting that it infers accurate evolutionary divergences. We developed the program PC_ali that constructs protein MSAs either de novo or modifying an input MSA, using a similarity matrix based on PC_sim. The program constructs a starting MSA based on the maximal cliques of the graph of these PAs and it refines it through progressive alignments along the tree reconstructed with PC_div. Compared with eight state-of-the-art multiple structure or sequence alignment tools, PC_ali achieves higher or equal aligned fraction and structural scores, sequence identity higher than structure aligners although lower than sequence aligners, highest score PC_sim, and highest similarity with the MSAs produced by other tools and with the reference MSA Balibase.

AVAILABILITY AND IMPLEMENTATION

https://github.com/ugobas/PC_ali.

New structural insights reveal an expanded reaction cycle for inositol pyrophosphate hydrolysis by human DIPP1

Nudix hydrolases attract considerable attention for their wide range of specialized activities in all domains of life. One particular group of Nudix phosphohydrolases (DIPPs), through their metabolism of diphosphoinositol polyphosphates (PP-InsPs), regulates the actions of these polyphosphates upon bioenergetic homeostasis. In the current study, we describe, at an atomic level, hitherto unknown properties of human DIPP1.We provide X-ray analysis of the catalytic core of DIPP1 in crystals complexed with either natural PP-InsPs, alternative PP-InsP stereoisomers, or non-hydrolysable methylene bisphosphonate analogs ("PCP-InsPs"). The conclusions that we draw from these data are interrogated by studying the impact upon catalytic activity upon mutagenesis of certain key residues. We present a picture of a V-shaped catalytic furrow with overhanging ridges constructed from flexible positively charged side chains; within this cavity, the labile phosphoanhydride bond is appropriately positioned at the catalytic site by an extensive series of interlocking polar contacts which we analogize as "suspension cables." We demonstrate functionality for a triglycine peptide within a β-strand which represents a non-canonical addition to the standard Nudix catalytic core structure. We describe pre-reaction enzyme/substrate states which we posit to reflect a role for electrostatic steering in substrate capture. Finally, through time-resolved analysis, we uncover a chronological sequence of DIPP1/product post-reaction states, one of which may rationalize a role for InsP₆ as an inhibitor of catalytic activity.

Asp1 from Schizosaccharomyces pombe binds a [2Fe-2S](2+) cluster which inhibits inositol pyrophosphate 1-phosphatase activity

Iron-sulfur (Fe-S) clusters are widely distributed protein cofactors that are vital to cellular biochemistry and the maintenance of bioenergetic homeostasis, but to our knowledge, they have never been identified in any phosphatase. Here, we describe an iron-sulfur cluster in Asp1, a dual-function kinase/phosphatase that regulates cell morphogenesis in Schizosaccharomyces pombe. Full-length Asp1, and its phosphatase domain (Asp1(371-920)), were each heterologously expressed in Escherichia coli. The phosphatase activity is exquisitely specific: it hydrolyzes the 1-diphosphate from just two members of the inositol pyrophosphate (PP-InsP) signaling family, namely, 1-InsP7 and 1,5-InsP8. We demonstrate that Asp1 does not hydrolyze either InsP6, 2-InsP7, 3-InsP7, 4-InsP7, 5-InsP7, 6-InsP7, or 3,5-InsP8. We also recorded 1-phosphatase activity in a human homologue of Asp1, hPPIP5K1, which was heterologously expressed in Drosophila S3 cells with a biotinylated N-terminal tag, and then isolated from cell lysates with avidin beads. Purified, recombinant Asp1(371-920) contained iron and acid-labile sulfide, but the stoichiometry (0.8 atoms of each per protein molecule) indicates incomplete iron-sulfur cluster assembly. We reconstituted the Fe-S cluster in vitro under anaerobic conditions, which increased the stoichiometry to approximately 2 atoms of iron and acid-labile sulfide per Asp1 molecule. The presence of a [2Fe-2S](2+) cluster in Asp1(371-920) was demonstrated by UV-visible absorption, resonance Raman spectroscopy, and electron paramagnetic resonance spectroscopy. We determined that this [2Fe-2S](2+) cluster is unlikely to participate in redox chemistry, since it rapidly degraded upon reduction by dithionite. Biochemical and mutagenic studies demonstrated that the [2Fe-2S](2+) cluster substantially inhibits the phosphatase activity of Asp1, thereby increasing its net kinase activity.

SFESA: a web server for pairwise alignment refinement by secondary structure shifts

Background

Protein sequence alignment is essential for a variety of tasks such as homology modeling and active site prediction. Alignment errors remain the main cause of low-quality structure models. A bioinformatics tool to refine alignments is needed to make protein alignments more accurate.

Results

We developed the SFESA web server to refine pairwise protein sequence alignments. Compared to the previous version of SFESA, which required a set of 3D coordinates for a protein, the new server will search a sequence database for the closest homolog with an available 3D structure to be used as a template. For each alignment block defined by secondary structure elements in the template, SFESA evaluates alignment variants generated by local shifts and selects the best-scoring alignment variant. A scoring function that combines the sequence score of profile-profile comparison and the structure score of template-derived contact energy is used for evaluation of alignments. PROMALS pairwise alignments refined by SFESA are more accurate than those produced by current advanced alignment methods such as HHpred and CNFpred. In addition, SFESA also improves alignments generated by other software.

Conclusions

SFESA is a web-based tool for alignment refinement, designed for researchers to compute, refine, and evaluate pairwise alignments with a combined sequence and structure scoring of alignment blocks. To our knowledge, the SFESA web server is the only tool that refines alignments by evaluating local shifts of secondary structure elements. The SFESA web server is available at http://prodata.swmed.edu/sfesa.