A tale of two symmetrical tails: structural and functional characteristics of palindromes in proteins

WAS IT A MATch I SAW? Approximate palindromes lead to overstated false match rates in benchmarks using reversed sequences

Background

Software for labeling biological sequences typically produces a theory-based statistic for each match (the E-value) that indicates the likelihood of seeing that match's score by chance. E-values accurately predict false match rate for comparisons of random (shuffled) sequences, and thus provide a reasoned mechanism for setting score thresholds that enable high sensitivity with low expected false match rate. This threshold-setting strategy is challenged by real biological sequences, which contain regions of local repetition and low sequence complexity that cause excess matches between non-homologous sequences. Knowing this, tool developers often develop benchmarks that use realistic-seeming decoy sequences to explore empirical tradeoffs between sensitivity and false match rate. A recent trend has been to employ reversed biological sequences as realistic decoys, because these preserve the distribution of letters and the existence of local repeats, while disrupting the original sequence's functional properties. However, we and others have observed that sequences appear to produce high scoring alignments to their reversals with surprising frequency, leading to overstatement of false match risk that may negatively affect downstream analysis.

Results

We demonstrate that an alignment between a sequence S and its (possibly mutated) reversal tends to produce higher scores than alignment between truly unrelated sequences, even when S is a shuffled string with no notable repetitive or low-complexity regions. This phenomenon is due to the unintuitive fact that (even randomly shuffled) sequences contain palindromes that are on average longer than the longest common substrings (LCS) shared between permuted variants of the same sequence. Though the expected palindrome length is only slightly larger than the expected LCS, the distribution of alignment scores involving reversed sequences is strongly right-shifted, leading to greatly increased frequency of high-scoring alignments to reversed sequences.

Impact

Overestimates of false match risk can motivate unnecessarily high score thresholds, leading to potentially reduced true match sensitivity. Also, when tool sensitivity is only reported up to the score of the first matched decoy sequence, a large decoy set consisting of reversed sequences can obscure sensitivity differences between tools. As a result of these observations, we advise that reversed biological sequences be used as decoys only when care is taken to remove positive matches in the original (un-reversed) sequences, or when overstatement of false labeling is not a concern. Though the primary focus of the analysis is on sequence annotation, we also demonstrate that the prevalence of internal palindromes may lead to an overstatement of the rate of false labels in protein identification with mass spectrometry.

An IL13Rα2 peptide exhibits therapeutic activity against metastatic colorectal cancer

Background

Interleukin 13 receptor α2 (IL13Rα2) is overexpressed in metastatic colorectal cancer. Here, we have developed novel strategies to block IL-13 binding to IL13Rα2 in order to reduce metastatic spread.

Methods

Synthetic IL13Rα2 D1 peptide (GSETWKTIITKN) was tested for the inhibition of IL-13 binding to IL13Rα2 using ELISA and different cellular assays. Peptide blocking effects on different cell signalling mediators were determined by western blot. An enantiomer version of the peptide (D-D1) was prepared to avoid proteolytic digestion. Nude mice were used for tumour growth and survival analysis after treatment with IL13Rα2 peptides.

Results

IL13Rα2 D1 peptide inhibited migration, invasion, and proliferation in metastatic colorectal and glioblastoma cancer cells treated with IL-13. Residues ⁸²K, ⁸³T, ⁸⁵I and ⁸⁶T were essential for blocking IL-13. IL13Rα2 peptide abolished ligand-mediated receptor internalisation and degradation, and substantially decreased IL-13 signalling capacity through IL13Rα2 to activate the FAK, PI3K/AKT and Src pathways as well as MT1-MMP expression. In addition, D1 significantly inhibited IL-13-mediated STAT6 activation through IL13Rα1. Nude mice treated with the enantiomer D-D1 peptide showed a remarkable survival increase.

Conclusions

We propose that the D-D1 peptide from IL13Rα2 represents a promising therapeutic agent to inhibit metastatic progression in colorectal cancer and, likely, other solid tumours.

KCNA4 deficiency leads to a syndrome of abnormal striatum, congenital cataract and intellectual disability

BACKGROUND

Voltage-gated potassium channels are highly diverse proteins representing the most complex class of voltage-gated ion channels from structural and functional perspectives. Deficiency of these channels usually results in various human disorders.

OBJECTIVES

To describe a novel autosomal recessive syndrome associated with KCNA4 deficiency leading to congenital cataract, abnormal striatum, intellectual disability and attention deficit hyperactivity disorder.

METHODS

We used SNP arrays, linkage analyses, autozygosity mapping, whole-exome sequencing, RT-PCR and two-electrode voltage-clamp recording.

RESULTS

We identified a missense variant (p.Arg89Gln) in KCNA4 in four patients from a consanguineous family manifesting a novel syndrome of congenital cataract, abnormal striatum, intellectual disability and attention deficit hyperactivity disorder. The variant was fully segregated with the disease and absent in 747 ethnically matched exomes. Xenopus oocytes were injected with human Kv1.4 wild-type mRNA, R89Q and WT/R89Q channels. The wild type had mean current amplitude that was significantly greater than those recorded from the cells expressing the same amount of mutant mRNA. Co-expression of the wild type and mutant mRNAs resulted in mean current amplitude that was significantly different from that of the wild type. RT-PCR indicated that KCNA4 is present in mouse brain, lens and retina. KCNA4 interacts with several molecules including synaptotagmin I, DLG1 and DLG2. The channel co-localises with cholinergic amacrine and rod bipolar cells in rats and is widely distributed in the central nervous system. Based on previous studies, the channel is highly expressed in outer retina, rod inner segments, hippocampus and concentrated in axonal membranes.

CONCLUSION

KCNA4 (Kv1.4) is implicated in a novel syndrome characterised by striatal thinning, congenital cataract and attention deficit hyperactivity disorder. Our study highlights potassium channels' role in ocular and neuronal genetics.

Analyses of the Sequence and Structural Properties Corresponding to Pentapeptide and Large Palindromes in Proteins

The analyses of 3967 representative proteins selected from the Protein Data Bank revealed the presence of 2803 pentapeptide and large palindrome sequences with known secondary structure conformation. These represent 2014 unique palindrome sequences. 60% palindromes are not associated with any regular secondary structure and 28% are in helix conformation, 11% in strand conformation and 1% in the coil conformation. The average solvent accessibility values are in the range between 0–155.28 Ų suggesting that the palindromes in proteins can be either buried, exposed to the solvent or share an intermittent property. The number of residue neighborhood contacts defined by interactions ≤ 3.2 Ǻ is in the range between 0–29 residues. Palindromes of the same length in helix, strand and coil conformation are associated with different amino acid residue preferences at the individual positions. Nearly, 20% palindromes interact with catalytic/active site residues, ligand or metal ions in proteins and may therefore be important for function in the corresponding protein. The average hydrophobicity values for the pentapeptide and large palindromes range between -4.3 to +4.32 and the number of palindromes is almost equally distributed between the negative and positive hydrophobicity values. The palindromes represent 107 different protein families and the hydrolases, transferases, oxidoreductases and lyases contain relatively large number of palindromes.

A conservative mutant of a proteolytic fragment produced during fibril formation enhances fibrillogenesis

The fibrillogenesis of a peptide corresponding to residues 35-51 of human α-lactalbumin (¹GYDTQAIVENNESTEYG¹⁷) can be dramatically enhanced by the addition of a tetrapeptide TDYG homologous to its C-terminus (TEYG). Generation of spontaneous hydrolytic products similar to this peptide was demonstrated by mass-spectrometry analysis of GYDTQAIVENNESTEYG peptide solution components during fibrillogenesis. Possible mechanisms and roles of short peptides in protein metabolism are discussed.

Unravelling the biodiversity of nanoscale signatures of spider silk fibres

Living organisms are masters at designing outstanding self-assembled nanostructures through a hierarchical organization of modular proteins. Protein-based biopolymers improved and selected by the driving forces of molecular evolution are among the most impressive archetypes of nanomaterials. One of these biomacromolecules is the myriad of compound fibroins of spider silks, which combine surprisingly high tensile strength with great elasticity. However, no consensus on the nano-organization of spider silk fibres has been reached. Here we explore the biodiversity of spider silk fibres, focusing on nanoscale characterization with high-resolution atomic force microscopy. Our results reveal an evolution of the nanoroughness, nanostiffness, nanoviscoelastic, nanotribological and nanoelectric organization of microfibres, even when they share similar sizes and shapes. These features are related to unique aspects of their molecular structures. The results show that combined nanoscale analyses of spider silks may enable the screening of appropriate motifs for bioengineering synthetic fibres from recombinant proteins.

PPS: A computing engine to find Palindromes in all Protein sequences

The primary structure of a protein molecule comprises a linear chain of amino acid residues. Certain parts of this linear chain are unique in nature and function. They can be classified under different categories and their roles studied in detail. Two such unique categories are the palindromic sequences and the Single Amino Acid Repeats (SAARs), which plays a major role in the structure, function and evolution of the protein molecule. In spite of their presence in various protein sequences, palindromes have not yet been investigated in detail. Thus, to enable a comprehensive understanding of these sequences, a computing engine, PPS, has been developed. The users can search the occurrences of palindromes and SAARs in all the protein sequences available in various databases and can view the three-dimensional structures (in case it is available in the known three-dimensional protein structures deposited to the Protein Data Bank) using the graphics plug-in Jmol. The proposed server is the first of its kind and can be freely accessed through the World Wide Web.

AVAILABILITY

URL http://pranag.physics.iisc.ernet.in/pps/

RepEx: repeat extractor for biological sequences

Genomic sequences are far from being random but are made up of systematically ordered and information rich patterns. These repeated sequence patterns have been vastly utilized for their fundamental importance in understanding the genome function and organization. To this end, a comprehensive toolkit, RepEx, has been developed which extracts repeat (inverted, everted and mirror) patterns from the given genome sequence(s) without any constraints. The toolkit can also be used to fetch the inverted repeats present in the protein sequence(s). Further, it is capable of extracting exact and degenerate repeats with a user defined spacer intervals. It is remarkably more precise and sensitive when compared to the existing tools. An example with comprehensive case studies and a performance evaluation of the proposed toolkit has been presented to authenticate its efficiency and accuracy.

An algorithm to find all palindromic sequences in proteins

A palindrome is a set of characters that reads the same forwards and backwards. Since the discovery of palindromic peptide sequences two decades ago, little effort has been made to understand its structural, functional and evolutionary significance. Therefore, in view of this, an algorithm has been developed to identify all perfect palindromes (excluding the palindromic subset and tandem repeats) in a single protein sequence. The proposed algorithm does not impose any restriction on the number of residues to be given in the input sequence. This avant-garde algorithm will aid in the identification of palindromic peptide sequences of varying lengths in a single protein sequence.

Flexible and efficient genome tiling design with penalized uniqueness score

Background

As a powerful tool in whole genome analysis, tiling array has been widely used in the answering of many genomic questions. Now it could also serve as a capture device for the library preparation in the popular high throughput sequencing experiments. Thus, a flexible and efficient tiling array design approach is still needed and could assist in various types and scales of transcriptomic experiment.

Results

In this paper, we address issues and challenges in designing probes suitable for tiling array applications and targeted sequencing. In particular, we define the penalized uniqueness score, which serves as a controlling criterion to eliminate potential cross-hybridization, and a flexible tiling array design pipeline. Unlike BLAST or simple suffix array based methods, computing and using our uniqueness measurement can be more efficient for large scale design and require less memory. The parameters provided could assist in various types of genomic tiling task. In addition, using both commercial array data and experiment data we show, unlike previously claimed, that palindromic sequence exhibiting relatively lower uniqueness.

Conclusions

Our proposed penalized uniqueness score could serve as a better indicator for cross hybridization with higher sensitivity and specificity, giving more control of expected array quality. The flexible tiling design algorithm incorporating the penalized uniqueness score was shown to give higher coverage and resolution. The package to calculate the penalized uniqueness score and the described probe selection algorithm are implemented as a Perl program, which is freely available at http://www1.fbn-dummerstorf.de/en/forschung/fbs/fb3/paper/2012-yang-1/OTAD.v1.1.tar.gz.

Analysis of compensatory substitution and gene evolution on the MAGEA/CSAG-palindrome of the primate X chromosomes

The human X chromosome contains a large number of inverted repeat DNA palindromes. Although arbitrary substitutions destroyed the inverted repeat structure of MAGEA/CSAG-palindrome during the evolutionary process of the primates, most of the substitutions are compensatory. Using maximum parsimony, it is demonstrated that the compensatory substitutions are prone to occur between bases with similar structures on the human, chimpanzee and orangutan MAGEA/CSAG-palindromes. Furthermore, it is found that MAGEA/CSAG genes also exist in orangutan and rhesus monkey palindromes by homologous searching. This suggests that the MAGEA/CSAG-palindrome might predate the divergence of human and other primate lineages. Comparative sequence analysis of the arms and genes on the primate MAGEA/CSAG-palindromes provides possible evidence of subsequently arm to arm gene conversion. These compensatory substitutions on the MAGEA/CSAG-palindrome of the primate X chromosomes play an important role in maintaining their structural symmetry during the process of formation.