Molecular determinants of improved cathepsin B inhibition by new cystatins obtained by DNA shuffling

Background: Cystatins are inhibitors of cysteine proteases. The majority are only weak inhibitors of human cathepsin B, which has been associated with cancer, Alzheimer's disease and arthritis. Results: Starting from the sequences of oryzacystatin-1 and canecystatin-1, a shuffling library was designed and a hybrid clone obtained, which presented higher inhibitory activity towards cathepsin B. This clone presented two unanticipated point mutations as well as an N-terminal deletion. Reversing each point mutation independently or both simultaneously abolishes the inhibitory activity towards cathepsin B. Homology modeling together with experimental studies of the reverse mutants revealed the likely molecular determinants of the improved inhibitory activity to be related to decreased protein stability. Conclusion: A combination of experimental approaches including gene shuffling, enzyme assays and reverse mutation allied to molecular modeling has shed light upon the unexpected inhibitory properties of certain cystatin mutants against Cathepsin B. We conclude that mutations disrupting the hydrophobic core of phytocystatins increase the flexibility of the N-terminus, leading to an increase in inhibitory activity. Such mutations need not affect the inhibitory site directly but may be observed distant from it and manifest their effects via an uncoupling of its three components as a result of increased protein flexibility.

Lossless filter for multiple repeats with bounded edit distance

Background: Identifying local similarity between two or more sequences, or identifying repeats occurring at least twice in a sequence, is an essential part in the analysis of biological sequences and of their phylogenetic relationship. Finding such fragments while allowing for a certain number of insertions, deletions, and substitutions, is however known to be a computationally expensive task, and consequently exact methods can usually not be applied in practice. Results: The filter TUIUIU that we introduce in this paper provides a possible solution to this problem. It can be used as a preprocessing step to any multiple alignment or repeats inference method, eliminating a possibly large fraction of the input that is guaranteed not to contain any approximate repeat. It consists in the verification of several strong necessary conditions that can be checked in a fast way. We implemented three versions of the filter. The first is simply a straightforward extension to the case of multiple sequences of an application of conditions already existing in the literature. The second uses a stronger condition which, as our results show, enable to filter sensibly more with negligible (if any) additional time. The third version uses an additional condition and pushes the sensibility of the filter even further with a non negligible additional time in many circumstances; our experiments show that it is particularly useful with large error rates. The latter version was applied as a preprocessing of a multiple alignment tool, obtaining an overall time (filter plus alignment) on average 63 and at best 530 times smaller than before (direct alignment), with in most cases a better quality alignment. Conclusion: To the best of our knowledge, TUIUIU is the first filter designed for multiple repeats and for dealing with error rates greater than 10% of the repeats length.

TESTING STATISTICAL HYPOTHESIS ON RANDOM TREES AND APPLICATIONS TO THE PROTEIN CLASSIFICATION PROBLEM

Efficient automatic protein classification is of central importance in genomic annotation. As an independent way to check the reliability of the classification, we propose a statistical approach to test if two sets of protein domain sequences coming from two families of the Pfam database are significantly different. We model protein sequences as realizations of Variable Length Markov Chains (VLMC) and we use the context trees as a signature of each protein family. Our approach is based on a Kolmogorov-Smirnov-type goodness-of-fit test proposed by Balding et at. [Limit theorems for sequences of random trees (2008), DOI: 10.1007/s11749-008-0092-z]. The test statistic is a supremum over the space of trees of a function of the two samples; its computation grows, in principle, exponentially fast with the maximal number of nodes of the potential trees. We show how to transform this problem into a max-flow over a related graph which can be solved using a Ford-Fulkerson algorithm in polynomial time on that number. We apply the test to 10 randomly chosen protein domain families from the seed of Pfam-A database (high quality, manually curated families). The test shows that the distributions of context trees coming from different families are significantly different. We emphasize that this is a novel mathematical approach to validate the automatic clustering of sequences in any context. We also study the performance of the test via simulations on Galton-Watson related processes.

THE LOWER CENTRAL AND DERIVED SERIES OF THE BRAID GROUPS OF THE SPHERE

In this paper, we determine the lower central and derived series for the braid groups of the sphere. We are motivated in part by the study of Fadell-Neuwirth short exact sequences, but the problem is important in its own right. The braid groups of the 2-sphere S(2) were studied by Fadell, Van Buskirk and Gillette during the 1960s, and are of particular interest due to the fact that they have torsion elements (which were characterised by Murasugi). We first prove that for all n epsilon N, the lower central series of the n-string braid group B(n)(S(2)) is constant from the commutator subgroup onwards. We obtain a presentation of Gamma(2)(Bn(S(2))), from which we observe that Gamma(2)(B(4)(S(2))) is a semi-direct product of the quaternion group Q(8) of order 8 by a free group F(2) of rank 2. As for the derived series of Bn(S(2)), we show that for all n >= 5, it is constant from the derived subgroup onwards. The group Bn(S(2)) being finite and soluble for n <= 3, the critical case is n = 4 for which the derived subgroup is the above semi-direct product Q(8) (sic) F(2). By proving a general result concerning the structure of the derived subgroup of a semi-direct product, we are able to determine completely the derived series of B(4)(S(2)) which from (B(4)(S(2)))(4) onwards coincides with that of the free group of rank 2, as well as its successive derived series quotients.

Environmental stresses inhibit splicing in the aquatic fungus Blastocladiella emersonii

Background: Exposure of cells to environmental stress conditions can lead to the interruption of several intracellular processes, in particular those performed by macromolecular complexes such as the spliceosome. Results: During nucleotide sequencing of cDNA libraries constructed using RNA isolated from B. emersonii cells submitted to heat shock and cadmium stress, a large number of ESTs with retained introns was observed. Among the 6,350 ESTs obtained through sequencing of stress cDNA libraries, 181 ESTs presented putative introns (2.9%), while sequencing of cDNA libraries from unstressed B. emersonii cells revealed only 0.2% of ESTs containing introns. These data indicate an enrichment of ESTs with introns in B. emersonii stress cDNA libraries. Among the 85 genes corresponding to the ESTs that retained introns, 19 showed more than one intron and three showed three introns, with intron length ranging from 55 to 333 nucleotides. Canonical splicing junctions were observed in most of these introns, junction sequences being very similar to those found in introns from genes previously characterized in B. emersonii, suggesting that inhibition of splicing during stress is apparently a random process. Confirming our observations, analyses of gpx3 and hsp70 mRNAs by Northern blot and S1 protection assays revealed a strong inhibition of intron splicing in cells submitted to cadmium stress. Conclusion: In conclusion, data indicate that environmental stresses, particularly cadmium treatment, inhibit intron processing in B. emersonii, revealing a new adaptive response to cellular exposure to this heavy metal.

Cross-Talk Between Mitochondria and NADPH Oxidase: Effects of Mild Mitochondrial Dysfunction on Angiotensin II-Mediated Increase in Nox Isoform Expression and Activity in Vascular Smooth Muscle Cells

Mitochondria and NADPH oxidase activation are concomitantly involved in pathogenesis of many vascular diseases. However, possible cross-talk between those ROS-generating systems is unclear. We induced mild mitochondrial dysfunction due to mitochondrial DNA damage after 24 h incubation of rabbit aortic smooth muscle (VSMC) with 250 ng/mL ethidium bromide (EtBr). VSMC remained viable and had 29% less oxygen consumption, 16% greater baseline hydrogen peroxide, and unchanged glutathione levels. Serum-stimulated proliferation was unaltered at 24 h. Although PCR amplification of several mtDNA sequences was preserved, D-Loop mtDNA region showed distinct amplification of shorter products after EtBr. Such evidence for DNA damage was further enhanced after angiotensin-II (AngII) incubation. Remarkably, the normally observed increase in VSMC membrane fraction NADPH oxidase activity after AngII was completely abrogated after EtBr, together with failure to upregulate Nox1 mRNA expression. Conversely, basal Nox4 mRNA expression increased 1.6-fold, while being unresponsive to AngII. Similar loss in AngII redox response occurred after 24 h antimycin-A incubation. Enhanced Nox4 expression was unassociated with endoplasmic reticulum stress markers. Protein disulfide isomerase, an NADPH oxidase regulator, exhibited increased expression and inverted pattern of migration to membrane fraction after EtBr. These results unravel functionally relevant cross-talk between mitochondria and NADPH oxidase, which markedly affects redox responses to AngII. Antioxid Redox Signal 11, 1265-1278.

Haplotype distribution of five nuclear genes based on network genealogies and Bayesian inference indicates that Trypanosoma cruzi hybrid strains are polyphyletic

Chagas disease is still a major public health problem in Latin America. Its causative agent, Trypanosoma cruzi, can be typed into three major groups, T. cruzi I, T. cruzi II and hybrids. These groups each have specific genetic characteristics and epidemiological distributions. Several highly virulent strains are found in the hybrid group; their origin is still a matter of debate. The null hypothesis is that the hybrids are of polyphyletic origin, evolving independently from various hybridization events. The alternative hypothesis is that all extant hybrid strains originated from a single hybridization event. We sequenced both alleles of genes encoding EF-1 alpha, actin and SSU rDNA of 26 T. cruzi strains and DHFR-TS and TR of 12 strains. This information was used for network genealogy analysis and Bayesian phylogenies. We found T. cruzi I and T. cruzi II to be monophyletic and that all hybrids had different combinations of T. cruzi I and T. cruzi II haplotypes plus hybrid-specific haplotypes. Bootstrap values (networks) and posterior probabilities (Bayesian phylogenies) of clades supporting the monophyly of hybrids were far below the 95% confidence interval, indicating that the hybrid group is polyphyletic. We hypothesize that T. cruzi I and T. cruzi II are two different species and that the hybrids are extant representatives of independent events of genome hybridization, which sporadically have sufficient fitness to impact on the epidemiology of Chagas disease.

GPIomics: global analysis of glycosylphosphatidylinositol-anchored molecules of Trypanosoma cruzi

Glycosylphosphatidylinositol (GPI) anchoring is a common, relevant posttranslational modification of eukaryotic surface proteins. Here, we developed a fast, simple, and highly sensitive (high attomole-low femtomole range) method that uses liquid chromatography-tandem mass spectrometry (LC-MS(n)) for the first large-scale analysis of GPI-anchored molecules (i.e., the GPIome) of a eukaryote, Trypanosoma cruzi, the etiologic agent of Chagas disease. Our genome-wise prediction analysis revealed that approximately 12% of T. cruzi genes possibly encode GPI-anchored proteins. By analyzing the GPIome of T. cruzi insect-dwelling epimastigote stage using LC-MS(n), we identified 90 GPI species, of which 79 were novel. Moreover, we determined that mucins coded by the T. cruzi small mucin-like gene (TcSMUG S) family are the major GPI-anchored proteins expressed on the epimastigote cell surface. TcSMUG S mucin mature sequences are short (56-85 amino acids) and highly O-glycosylated, and contain few proteolytic sites, therefore, less likely susceptible to proteases of the midgut of the insect vector. We propose that our approach could be used for the high throughput GPIomic analysis of other lower and higher eukaryotes. Molecular Systems Biology 7 April 2009; doi:10.1038/msb.2009.13