Prediction of subcellular localization using sequence-biased recurrent networks

Motivation: Targeting peptides direct nascent proteins to their specific subcellular compartment. Knowledge of targeting signals enables informed drug design and reliable annotation of gene products. However, due to the low similarity of such sequences and the dynamical nature of the sorting process, the computational prediction of subcellular localization of proteins is challenging. Results: We contrast the use of feed forward models as employed by the popular TargetP/SignalP predictors with a sequence-biased recurrent network model. The models are evaluated in terms of performance at the residue level and at the sequence level, and demonstrate that recurrent networks improve the overall prediction performance. Compared to the original results reported for TargetP, an ensemble of the tested models increases the accuracy by 6 and 5% on non-plant and plant data, respectively.

The applicability of recurrent neural networks for biological sequence analysis

Selection of machine learning techniques requires a certain sensitivity to the requirements of the problem. In particular, the problem can be made more tractable by deliberately using algorithms that are biased toward solutions of the requisite kind. In this paper, we argue that recurrent neural networks have a natural bias toward a problem domain of which biological sequence analysis tasks are a subset. We use experiments with synthetic data to illustrate this bias. We then demonstrate that this bias can be exploitable using a data set of protein sequences containing several classes of subcellular localization targeting peptides. The results show that, compared with feed forward, recurrent neural networks will generally perform better on sequence analysis tasks. Furthermore, as the patterns within the sequence become more ambiguous, the choice of specific recurrent architecture becomes more critical.

Detecting and sorting targeting peptides wtih neural networks and support vector machines

This paper presents a composite multi-layer classifier system for predicting the subcellular localization of proteins based on their amino acid sequence. The work is an extension of our previous predictor PProwler v1.1 which is itself built upon the series of predictors SignalP and TargetP. In this study we outline experiments conducted to improve the classifier design. The major improvement came from using Support Vector machines as a "smart gate" sorting the outputs of several different targeting peptide detection networks. Our final model (PProwler v1.2) gives MCC values of 0.873 for non-plant and 0.849 for plant proteins. The model improves upon the accuracy of our previous subcellular localization predictor (PProwler v1.1) by 2% for plant data (which represents 7.5% improvement upon TargetP).

Exploiting sequence dependencies in the prediction of peroxisomal proteins

Prediction of peroxisomal matrix proteins generally depends on the presence of one of two distinct motifs at the end of the amino acid sequence. PTS1 peroxisomal proteins have a well conserved tripeptide at the C-terminal end. However, the preceding residues in the sequence arguably play a crucial role in targeting the protein to the peroxisome. Previous work in applying machine learning to the prediction of peroxisomal matrix proteins has failed W capitalize on the full extent of these dependencies. We benchmark a range of machine learning algorithms, and show that a classifier - based on the Support Vector Machine - produces more accurate results when dependencies between the conserved motif and the preceding section are exploited. We publish an updated and rigorously curated data set that results in increased prediction accuracy of most tested models.

DNA chip designed antisense oligodeoxynucleotides targeting EGFR MRNA for brain tumour therapy

Glioblastoma multiforme (GBM) is a malignant brain tumour for which there is currently no effective treatment regime. It is thought to develop due to the overexpression of a number of genes, including the epidermal growth factor receptor (EGFR), which is found in over 40% of GBM. Novel forms of treatment such as antisense therapy may allow for the specific inhibition of aberrant genes and thus they are optimistic therapies for future treatment of GBM. Oligodeoxynucleotides (ODNs) are small pieces of DNA that are often modified to increase their stability to nucleases and can be targeted to the aberrant gene in order to inhibit it and thus prevent its transcription into protein. By specifically binding to mRNA in an antisense manner, they can bring about its degradation by a variety of mechanisms including the activation of RNase H and thus have great potential as therapeutic agents. One of the main drawbacks to the utilisation of this therapy so far is the lack of techniques that can successfully predict accessible regions on the target mRNA that the ODNs can bind to. DNA chip technology has been utilised here to predict target sequences on the EGFR mRNA and these ODNs (AS 1 and AS2) have been tested in vitro for their stability, uptake into cells and their efficacy on cellular growth, EGFR protein and mRNA. Studies showed that phosphorothioate and 2'O-methyl ODNs were significantly more stable than phosphodiester ODNs both in serum and serum-free conditions and that the mechanism of uptake into A431 cells was temperature dependent and more efficient with the use of optimised lipofectin. Efficacy results show that AS 1 and AS2 phosphorothioate antisense ODNs were capable of inhibiting cell proliferation by 69% ±4% and 65% ±4.5% respectively at 500nM in conjunction with a non-toxic dose of lipofectinTM used to enhance cellular delivery. Furthermore, control ODN sequences, 2' O-methyl derivatives and a third ODN sequence, that was found not to be capable of binding efficiently to the EGFR mRNA by DNA chip technology, showed no significant effect on cell proliferation. AS 1 almost completely inhibited EGFR protein levels within 48 hours with two doses of 500nM AS 1 with no effect on other EGFR family member proteins or by control sequences. RNA analysis showed a decrease in mRNA levels of 32.4% ±0.8% but techniques require further optimisation to confirm this. As there are variations found between human glioblastoma in situ and those developed as xenografts, analysis of effect of AS 1 and AS2 was performed on primary tumour cell lines derived from glioma patients. ODN treatment showed a specific knockdown of cell growth compared to any of the controls used. Furthermore, combination therapies were tested on A431 cell growth to determine the advantage of combining different antisense approaches and that of conventional drugs. Results varied between the combination treatments but indicated that with optimisation of treatment regimes and delivery techniques that combination therapies utilising antisense therapies would be plausible.

Sequence and Structural Determinants of Specificity Differences between Paralogous Transcription Factors

Transcription factors (TFs) control the temporal and spatial expression of target genes by interacting with DNA in a sequence-specific manner. Recent advances in high throughput experiments that measure TF-DNA interactions in vitro and in vivo have facilitated the identification of DNA binding sites for thousands of TFs. However, it remains unclear how each individual TF achieves its specificity, especially in the case of paralogous TFs that recognize distinct target genomic sites despite sharing very similar DNA binding motifs. In my work, I used a combination of high throughput in vitro protein-DNA binding assays and machine-learning algorithms to characterize and model the binding specificity of 11 paralogous TFs from 4 distinct structural families. My work proves that even very closely related paralogous TFs, with indistinguishable DNA binding motifs, oftentimes exhibit differential binding specificity for their genomic target sites, especially for sites with moderate binding affinity. Importantly, the differences I identify in vitro and through computational modeling help explain, at least in part, the differential in vivo genomic targeting by paralogous TFs. Future work will focus on in vivo factors that might also be important for specificity differences between paralogous TFs, such as DNA methylation, interactions with protein cofactors, or the chromatin environment. In this larger context, my work emphasizes the importance of intrinsic DNA binding specificity in targeting of paralogous TFs to the genome.

Development and Characterization of Monovalent and Bivalent RNA Aptamers Targeting the Common Pathway of Coagulation

Anticoagulant agents are commonly used drugs to reduce blood coagulation in acute and chronic clinical settings. Many of these drugs target the common pathway of coagulation because it is critical for thrombin generation and disruption of this portion of the pathway has profound effects on the hemostatic process. Currently available drugs for these indications struggle with balancing desired activity with immunogenicity and poor reversibility or irreversibility in the event of hemorrhage. While improvements are being made with the current drugs, new drugs with better therapeutic indices are needed for surgical intervention and chronic indications to prevent thrombosis from occurring.

A class of therapeutics known as aptamers may be able to meet the need for safer anticoagulant agents. Aptamer are short single-stranded RNA oligonucleotides that adopt specific secondary and tertiary structures based upon their sequence. They can be generated to both enzymes and cofactors because they derive their inhibitory activity by blocking protein-protein interactions, rather than active site inhibition. They inhibit their target proteins with a high level of specificity and bind with high affinity to their target. Additionally, they can be reversed using two different antidote approaches, specific oligonucleotide antidotes, or with cationic, “universal” antidotes. The reversal of their activity is both rapid and durable.

The ability of aptamers to be generated to cofactors has been conclusively proven by generating an aptamer targeting the common pathway coagulation cofactor, Factor V (FV). We developed two aptamers with anticoagulant ability that bind to both FV and FVa, the active cofactor. Both aptamers were truncated to smaller functional sizes and had specific point mutant aptamers developed for use as controls. The anticoagulant activity of both aptamer-mutant pairs was characterized using plasma-based clotting assays and whole blood assays. The mechanism of action resulting in anticoagulant activity was assessed for one aptamer. The aptamer was found to block FVa docking to membrane surfaces, a mechanism not previously observed in any of our other anticoagulant aptamers.

To explore development of aptamers as anticoagulant agents targeting the common pathway for surgical interventions, we fused two anticoagulant aptamers targeting Factor X and prothrombin into a single molecule. The bivalent aptamer was truncated to a minimal size while maintaining robust anticoagulant activity. Characterization of the bivalent aptamer in plasma-based clotting assays indicated we had generated a very robust anticoagulant therapeutic. Furthermore, we were able to simultaneously reverse the activity of both aptamers with a single oligonucleotide antidote. This rapid and complete reversal of anticoagulant activity is not available in the antithrombotic agents currently used in surgery.

Genome Sequence Of Streptomyces Wadayamensis Strain A23, An Endophytic Actinobacterium From Citrus Reticulata.

The actinobacterium Streptomyces wadayamensis A23 is an endophyte of Citrus reticulata that produces the antimycin and mannopeptimycin antibiotics, among others. The strain has the capability to inhibit Xylella fastidiosa growth. The draft genome of S. wadayamensis A23 has ~7.0 Mb and 6,006 protein-coding sequences, with a 73.5% G+C content.

Genome Sequence Of Bacillus Safensis Cfa06, Isolated From Biodegraded Petroleum In Brazil.

Bacillus safensis is a microorganism recognized for its biotechnological and industrial potential due to its interesting enzymatic portfolio. Here, as a means of gathering information about the importance of this species in oil biodegradation, we report a draft genome sequence of a strain isolated from petroleum.

Overlapped Sequence Types (sts) And Serogroups Of Avian Pathogenic (apec) And Human Extra-intestinal Pathogenic (expec) Escherichia Coli Isolated In Brazil.

Avian pathogenic Escherichia coli (APEC) strains belong to a category that is associated with colibacillosis, a serious illness in the poultry industry worldwide. Additionally, some APEC groups have recently been described as potential zoonotic agents. In this work, we compared APEC strains with extraintestinal pathogenic E. coli (ExPEC) strains isolated from clinical cases of humans with extra-intestinal diseases such as urinary tract infections (UTI) and bacteremia. PCR results showed that genes usually found in the ColV plasmid (tsh, iucA, iss, and hlyF) were associated with APEC strains while fyuA, irp-2, fepC sitDchrom, fimH, crl, csgA, afa, iha, sat, hlyA, hra, cnf1, kpsMTII, clpVSakai and malX were associated with human ExPEC. Both categories shared nine serogroups (O2, O6, O7, O8, O11, O19, O25, O73 and O153) and seven sequence types (ST10, ST88, ST93, ST117, ST131, ST155, ST359, ST648 and ST1011). Interestingly, ST95, which is associated with the zoonotic potential of APEC and is spread in avian E. coli of North America and Europe, was not detected among 76 APEC strains. When the strains were clustered based on the presence of virulence genes, most ExPEC strains (71.7%) were contained in one cluster while most APEC strains (63.2%) segregated to another. In general, the strains showed distinct genetic and fingerprint patterns, but avian and human strains of ST359, or ST23 clonal complex (CC), presented more than 70% of similarity by PFGE. The results demonstrate that some zoonotic-related STs (ST117, ST131, ST10CC, ST23CC) are present in Brazil. Also, the presence of moderate fingerprint similarities between ST359 E. coli of avian and human origin indicates that strains of this ST are candidates for having zoonotic potential.

Draft Genome Sequence Of Ft9, A Novel Bacillus Cereus Strain Isolated From A Brazilian Thermal Spring.

A Bacillus cereus strain, FT9, isolated from a hot spring in the midwest region of Brazil, had its entire genome sequenced.

Phtx-ii A Basic Myotoxic Phospholipase A2 From Porthidium Hyoprora Snake Venom, Pharmacological Characterization And Amino Acid Sequence By Mass Spectrometry.

A monomeric basic PLA2 (PhTX-II) of 14149.08 Da molecular weight was purified to homogeneity from Porthidium hyoprora venom. Amino acid sequence by in tandem mass spectrometry revealed that PhTX-II belongs to Asp49 PLA2 enzyme class and displays conserved domains as the catalytic network, Ca2+-binding loop and the hydrophobic channel of access to the catalytic site, reflected in the high catalytic activity displayed by the enzyme. Moreover, PhTX-II PLA2 showed an allosteric behavior and its enzymatic activity was dependent on Ca2+. Examination of PhTX-II PLA2 by CD spectroscopy indicated a high content of alpha-helical structures, similar to the known structure of secreted phospholipase IIA group suggesting a similar folding. PhTX-II PLA2 causes neuromuscular blockade in avian neuromuscular preparations with a significant direct action on skeletal muscle function, as well as, induced local edema and myotoxicity, in mice. The treatment of PhTX-II by BPB resulted in complete loss of their catalytic activity that was accompanied by loss of their edematogenic effect. On the other hand, enzymatic activity of PhTX-II contributes to this neuromuscular blockade and local myotoxicity is dependent not only on enzymatic activity. These results show that PhTX-II is a myotoxic Asp49 PLA2 that contributes with toxic actions caused by P. hyoprora venom.

The Putative Leishmania Telomerase Rna (leishter) Undergoes Trans-splicing And Contains A Conserved Template Sequence.

Telomerase RNAs (TERs) are highly divergent between species, varying in size and sequence composition. Here, we identify a candidate for the telomerase RNA component of Leishmania genus, which includes species that cause leishmaniasis, a neglected tropical disease. Merging a thorough computational screening combined with RNA-seq evidence, we mapped a non-coding RNA gene localized in a syntenic locus on chromosome 25 of five Leishmania species that shares partial synteny with both Trypanosoma brucei TER locus and a putative TER candidate-containing locus of Crithidia fasciculata. Using target-driven molecular biology approaches, we detected a ∼2,100 nt transcript (LeishTER) that contains a 5' spliced leader (SL) cap, a putative 3' polyA tail and a predicted C/D box snoRNA domain. LeishTER is expressed at similar levels in the logarithmic and stationary growth phases of promastigote forms. A 5'SL capped LeishTER co-immunoprecipitated and co-localized with the telomerase protein component (TERT) in a cell cycle-dependent manner. Prediction of its secondary structure strongly suggests the existence of a bona fide single-stranded template sequence and a conserved C[U/C]GUCA motif-containing helix II, representing the template boundary element. This study paves the way for further investigations on the biogenesis of parasite TERT ribonucleoproteins (RNPs) and its role in parasite telomere biology.

Molecular Genetic Variability Of Commercial And Wild Accessions Of Passion Fruit (passiflora Spp.) Targeting Ex Situ Conservation And Breeding.

Passiflora species are distributed throughout Latin America, and Brazil and Colombia serve as the centers of diversity for this genus. We performed cross-species amplification to evaluate 109 microsatellite loci in 14 Passiflora species and estimated the diversity and genetic structure of Passiflora cincinnata, Passiflora setaceae and Passiflora edulis. A total of 127 accessions, including 85 accessions of P. edulis, a commercial species, and 42 accessions of 13 wild species, were examined. The cross-species amplification was effective for obtaining microsatellite loci (average cross-amplification of 70%). The average number of alleles per locus (five) was relatively low, and the average diversity ranged from 0.52 in P. cincinnata to 0.32 in P. setacea. The Bayesian analyses indicated that the P. cincinnata and P. setacea accessions were distributed into two groups, and the P. edulis accessions were distributed into five groups. Private alleles were identified, and suggestions for core collections are presented. Further collections are necessary, and the information generated may be useful for breeding and conservation.

Timing and sequence of primary tooth eruption in children with cleft lip and palate

OBJECTIVE: To determine the timing and sequence of eruption of primary teeth in children with complete bilateral cleft lip and palate. MATERIAL AND METHODS: This cross-sectional study was conducted at the Hospital for Rehabilitation of Craniofacial Anomalies of the University of São Paulo, Bauru, SP, Brazil, with a sample of 395 children (128 girls and 267 boys) aged 0 to 48 months, with complete bilateral cleft lip and palate. RESULTS: Children with complete bilateral clefts presented a higher mean age of eruption of all primary teeth for both arches and both genders, compared to children without clefts. This difference was statistically signifcant for all teeth, except for the maxillary first molar. Mean age of eruption of most teeth was lower for girls compared to boys. The greatest delay was found for the maxillary lateral incisor, which was the eighth tooth of children with clefts of both genders. Analyzing by gender, the maxillary lateral incisor was the eighth tooth to erupt in girls and the last in boys. CONCLUSION: The results suggest an interference of the cleft on the timing and sequence of eruption of primary teeth.