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.

Analysis and application of an equilibrium model for in vitro bioassay systems with three components: Receptor, hormone and hormone-binding-protein

A simple theoretical framework is presented for bioassay studies using three component in vitro systems. An equilibrium model is used to derive equations useful for predicting changes in biological response after addition of hormone-binding-protein or as a consequence of increased hormone affinity. Sets of possible solutions for receptor occupancy and binding protein occupancy are found for typical values of receptor and binding protein affinity constants. Unique equilibrium solutions are dictated by the initial condition of total hormone concentration. According to the occupancy theory of drug action, increasing the affinity of a hormone for its receptor will result in a proportional increase in biological potency. However, the three component model predicts that the magnitude of increase in biological potency will be a small fraction of the proportional increase in affinity. With typical initial conditions a two-fold increase in hormone affinity for its receptor is predicted to result in only a 33% increase in biological response. Under the same conditions an Ii-fold increase in hormone affinity for receptor would be needed to produce a two-fold increase in biological potency. Some currently used bioassay systems may be unrecognized three component systems and gross errors in biopotency estimates will result if the effect of binding protein is not calculated. An algorithm derived from the three component model is used to predict changes in biological response after addition of binding protein to in vitro systems. The algorithm is tested by application to a published data set from an experimental study in an in vitro system (Lim et al., 1990, Endocrinology 127, 1287-1291). Predicted changes show good agreement (within 8%) with experimental observations. (C) 1998 Academic Press Limited.

Analysis of the substrate specificity of human sulfotransferases SULT1A1 and SULT1A3: Site-directed mutagenesis and kinetic studies

Sulfonation is an important metabolic process involved in the excretion and in some cases activation of various endogenous compounds and xenobiotics. This reaction is catalyzed by a family of enzymes named sulfotransferases. The cytosolic human sulfotransferases SULT1A1 and SULT1A3 have overlapping yet distinct substrate specificities. SULT1A1 favors simple phenolic substrates such as p-nitrophenol, whereas SULT1A3 prefers monoamine substrates such as dopamine. In this study we have used a variety of phenolic substrates to functionally characterize the role of the amino acid at position 146 in SULT1A1 and SULT1A3. First, the mutation A146E in SULT1A1 yielded a SULT1A3-like protein with respect to the Michaelis constant for simple phenols. The mutation E146A in SULT1A3 resulted in a SULT1A1-like protein with respect to the Michaelis constant for both simple phenols and monoamine compounds. When comparing the specificity of SULT1A3 toward tyramine with that for p-ethylphenol (which differs from tyramine in having no amine group on the carbon side chain), we saw a 200-fold preference for tyramine. The kinetic data obtained with the E146A mutant of SULT1A3 for these two substrates clearly showed that this protein preferred substrates without an amine group attached. Second, changing the glutamic acid at position 146 of SULT1A3 to a glutamine, thereby neutralizing the negative charge at this position, resulted in a 360-fold decrease in the specificity constant for dopamine. The results provide strong evidence that residue 146 is crucial in determining the substrate specificity of both SULT1A1 and SULT1A3 and suggest that there is a direct interaction between glutamic acid 146 in SULT1A3 and monoamine substrates.

Analysis of the genes flanking xabB: a methyltransferase gene is involved in albicidin biosynthesis in Xanthomonas albilineans

Transposon mutagenesis and complementation studies previously identified a gene (xabB) for a large (526 kDa) polyketide-peptide synthase required for biosynthesis of albicidin antibiotics and phytotoxins in the sugarcane leaf scald pathogen Xanthomonas albilineans. A cistron immediately downstream from xabB encodes a polypeptide of 343 aa containing three conserved motifs characteristic of a family of S-adenosyl-L-methionine (SAM)-dependent O-methyltransferases. Insertional mutagenesis and complementation indicate that the product of this cistron (designated xabC) is essential for albicidin production, and that there is no other required downstream cistron. The xab promoter region is bidirectional, and insertional mutagenesis of the first open reading frame (ORF) in the divergent gene also blocks albicidin biosynthesis. This divergent ORF (designated thp) encodes a protein of 239 aa displaying high similarity to several IS21-like transposition helper proteins. The thp cistron is not located in a recognizable transposon, and is probably a remnant from a past transposition event that may have contributed to the development of the albicidin biosynthetic gene cluster. Failure of 'in trans' complementation of rhp indicates that a downstream cistron transcribed with thp is required for albicidin biosynthesis. (C) 2000 Elsevier Science B.V. All rights reserved.

Analysis of Mouse Keratin 6a Regulatory Sequences in Transgenic Mice Reveals Constitutive, Tissue-Specific Expression by a Keratin 6a Minigene

The analysis of keratin 6 expression is complicated by the presence of multiple isoforms that are expressed constitutively in a number of internal stratified epithelia, in palmoplantar epidermis, and in the companion cell layer of the hair follicle. In addition, keratin 6 expression is inducible in interfollicular epidermis and the outer root sheath of the follicle, in response to wounding stimuli, phorbol esters, or retinoic acid. In order to establish the critical regions involved in the regulation of keratin 6a (the dominant isoform in mice), we generated transgenic mice with two different-sized mouse keratin 6a constructs containing either 1.3 kb or 0.12 kb of 5' flanking sequence linked to the lacZ reporter gene. Both constructs also contained the first intron and the 3' flanking sequence of mouse keratin 6a. Ectopic expression of either transgene was not observed. Double-label immunofluorescence analyses demonstrated expression of the reporter gene in keratin 6 expressing tissues, including the hair follicle, tongue, footpad, and nail bed, showing that both transgenes retained keratinocyte-specific expression. Quantitative analysis of beta -galactosidase activity verified that both the 1.3 and 0.12 kb keratin 6a promoter constructs produced similar levels of the reporter. Notably, both constructs were constitutively expressed in the outer root sheath and interfollicular epidermis in the absence of any activating stimulus, suggesting that they lack the regulatory elements that normally silence transcription in these cells. This study has revealed that a keratin 6a minigene contains critical cis elements that mediate tissue-specific expression and that the elements regulating keratin 6 induction lie distal to the 1.3 kb promoter region.

Analysis of thermodynamic non-ideality in terms of protein solvation

The effects of thermodynamic non-ideality on the forms of sedimentation equilibrium distributions for several isoelectric proteins have been analysed on the statistical-mechanical basis of excluded volume to obtain an estimate of the extent of protein solvation. Values of the effective solvation. parameter delta are reported for ellipsoidal as well as spherical models of the proteins, taken to be rigid, impenetrable macromolecular structures. The dependence of the effective solvated radius upon protein molecular mass exhibits reasonable agreement with the relationship calculated for a model in which the unsolvated protein molecule is surrounded by a 0.52-nm solvation shell. Although the observation that this shell thickness corresponds to a double layer of water molecules may be of questionable relevance to mechanistic interpretation of protein hydration, it augurs well for the assignment of magnitudes to the second virial coefficients of putative complexes in the quantitative characterization of protein-protein interactions under conditions where effects of thermodynamic non-ideality cannot justifiably be neglected. (C) 2001 Elsevier Science B.V. All rights reserved.

Analysis of the mouse transcriptome based on functional annotation of 60,770 full-length cDNAs

Only a small proportion of the mouse genome is transcribed into mature messenger RNA transcripts. There is an international collaborative effort to identify all full-length mRNA transcripts from the mouse, and to ensure that each is represented in a physical collection of clones. Here we report the manual annotation of 60,770 full-length mouse complementary DNA sequences. These are clustered into 33,409 'transcriptional units', contributing 90.1% of a newly established mouse transcriptome database. Of these transcriptional units, 4,258 are new protein-coding and 11,665 are new non-coding messages, indicating that non-coding RNA is a major component of the transcriptome. 41% of all transcriptional units showed evidence of alternative splicing. In protein-coding transcripts, 79% of splice variations altered the protein product. Whole-transcriptome analyses resulted in the identification of 2,431 sense-antisense pairs. The present work, completely supported by physical clones, provides the most comprehensive survey of a mammalian transcriptome so far, and is a valuable resource for functional genomics.

Biological determinants of extinction risk: why are smaller species less vulnerable?

It is becoming increasingly clear that species of smaller body size tend to be less vulnerable to contemporary extinction threats than larger species, but few studies have examined the mechanisms underlying this pattern. In this paper, data for the Australian terrestrial mammal fauna are used to ask whether higher reproductive output or smaller home ranges can explain the reduced extinction risk of smaller species. Extinct and endangered species do indeed have smaller litters and larger home ranges for their body size than expected under a null model. In multiple regressions, however, only litter size is a significant predictor of extinction risk once body size and phylogeny are controlled for. Larger litters contribute to fast population growth, and are probably part of the reason that smaller species are less extinction-prone. The effect of litter size varies between the mesic coastal regions and the and interior of Australia, indicating that the environment a species inhabits mediates the effect of biology on extinction risk. These results suggest that predicting extinction risk from biological traits is likely to be a complex task which must consider explicitly interactions between biology and environment.

Analysis of RNA-protein interactions by flow cytometry

Flow cytometry, in combination with advances in bead coding technologies, is maturing as a powerful high-throughput approach for analyzing molecular interactions. Applications of this technology include antibody assays and single nucleotide polymorphism mapping. This review describes the recent development of a microbead flow cytometric approach to analyze RNA-protein interactions and discusses emerging bead coding strategies that together will allow genome-wide identification of RNA-protein complexes. The microbead flow cytometric approach is flexible and provides new opportunities for functional genomic studies and small-molecule screening.

Analysis of the role of pglI in pilin glycosylation of Neisseria meningitidis

Pilin is the major subunit of the essential virulence factor pili and is glycosylated at Ser63. In this study we investigated the gene pglI to determine whether it is involved in the biosynthesis of the pilin-linked glycan of Neisseria meningitidis strain C311#3. A N. meningitidis C311#3pglI mutant resulted in a change of apparent molecular weight in SDS-PAGE and altered binding of antisera, consistent with a role in the biosynthesis of the pilin-linked glycan. These data, in conjunction with homology with well-characterised acyltransferases suggests a specific role for pglI in the biosynthesis of the basal 2,4-diacetamido-2,4,6-trideoxyhexose residue of the pilin-linked glycan. (C) 2004 Published by Elsevier B.V. on behalf of the Federation of European Microbiological Societies.

The use and implications of ribosomal DNA sequencing for the discrimination of digenean species

In just over a decade, the use of molecular approaches for the recognition of parasites has become commonplace. For trematodes, the internal transcribed spacer region of ribosomal DNA (ITS rDNA) has become the default region of choice. Here, we review the findings of 63 studies that report ITS rDNA sequence data for about 155 digenean species from 19 families, and then review the levels of variation that have been reported and how the variation has been interpreted. Overall, complete ITS sequences (or ITS1 or ITS2 regions alone) usually distinguish trematode species clearly, including combinations for which morphology gives ambiguous results. Closely related species may have few base differences and in at least one convincing case the ITS2 sequences of two good species are identical. In some cases, the ITS1 region gives greater resolution than the ITS2 because of the presence of variable repeat units that are generally lacking in the ITS2. Intraspecific variation is usually low and frequently apparently absent. Information on geographical variation of digeneans is limited but at least some of the reported variation probably reflects the presence of multiple species. Despite the accepted dogma that concerted evolution makes the individual representative of the entire species, a significant number of studies have reported at least some intraspecific variation. The significance of such variation is difficult to assess a posteriori, but it seems likely that identification and sequencing errors account for some of it and failure to recognise separate species may also be significant. Some reported variation clearly requires further analysis. The use of a yardstick to determine when separate species should be recognised is flawed. Instead, we argue that consistent genetic differences that are associated with consistent morphological or biological traits should be considered the marker for separate species. We propose a generalised approach to the use of rDNA to distinguish trematode species.