A simulated annealing algorithm for finding consensus sequences

Motivation: A consensus sequence for a family of related sequences is, as the name suggests, a sequence that captures the features common to most members of the family. Consensus sequences are important in various DNA sequencing applications and are a convenient way to characterize a family of molecules. Results: This paper describes a new algorithm for finding a consensus sequence, using the popular optimization method known as simulated annealing. Unlike the conventional approach of finding a consensus sequence by first forming a multiple sequence alignment, this algorithm searches for a sequence that minimises the sum of pairwise distances to each of the input sequences. The resulting consensus sequence can then be used to induce a multiple sequence alignment. The time required by the algorithm scales linearly with the number of input sequences and quadratically with the length of the consensus sequence. We present results demonstrating the high quality of the consensus sequences and alignments produced by the new algorithm. For comparison, we also present similar results obtained using ClustalW. The new algorithm outperforms ClustalW in many cases.

Heterogeneous nuclear ribonucleoprotein A3, a novel RNA trafficking response element-binding protein

The cis-acting response element, A2RE, which is sufficient for cytoplasmic mRNA trafficking in oligodendrocytes, binds a small group of rat brain proteins. Predominant among these is heterogeneous nuclear ribonucleoprotein (hnRNP) A2, a trans-acting factor for cytoplasmic trafficking of RNAs bearing A2RE-like sequences. We have now identified the other A2RE-binding proteins as hnRNP A1/A1(B), hnRNP B1, and four isoforms of hnRNP A3. The rat and human hnRNP A3 cDNAs have been sequenced, revealing the existence of alternatively spliced mRNAs. In Western blotting, 38-, 39-, 41 -, and 41.5-kDa components were all recognized by antibodies against a peptide in the glycine-rich region of hnRNP A3, but only the 41- and 41.5-kDa bands bound antibodies to a 15-residue N-terminal peptide encoded by an alternatively spliced part of exon 1. The identities of these four proteins were verified by Edman sequencing and mass spectral analysis of tryptic fragments generated from electrophoretically separated bands. Sequence-specific binding of bacterially expressed hnRNP A3 to A2RE has been demonstrated by biosensor and UV cross-linking electrophoretic mobility shift assays. Mutational analysis and confocal microscopy data support the hypothesis that the hnRNP A3 isoforms have a role in cytoplasmic trafficking of RNA.

Characterization of mouse profilaggrin: Evidence for nuclear engulfment and translocation of the profilaggrin B-domain during epidermal differentiation

Filaggrin is a keratin filament associated protein that is expressed in granular layer keratinocytes and derived by sequential proteolysis from a polyprotein precursor termed profilaggrin. Depending on the species, each profilaggrin molecule contains between 10 and 20 filaggrin subunits organized as tandem repeats with a calcium-binding domain at the N-terminal end. We now report the characterization of the complete mouse gene. The structural organization of the mouse gene is identical to the human profilaggrin gene and consists of three exons with a 4 kb intron within the 5' noncoding region and a 1.7 kb intron separating the sequences encoding the calcium-binding EF-hand motifs. A processed pseudogene was found embedded within the second intron. The third and largest exon encodes the second EF-hand, a basic domain (designated the B-domain) followed by 12 filaggrin repeats and a unique C-terminal tail domain. A polyclonal anti-body raised against the conceptually translated sequence of the B-domain specifically stained keratohyalin granules and colocalized with a filaggrin antibody in granular layer cells. In upper granular layer cells, B-domain containing keratohyalin granules were in close apposition to the nucleus and, in some cells, appeared to be completely engulfed by the nucleus. In transition layer cells, B-domain staining was evident in the nucleus whereas filaggrin staining remained cytoplasmic. Nuclear staining of the B-domain was also observed in primary mouse keratinocytes induced to differentiate. This study has also revealed significant sequence homology between the mouse and human promoter sequences and in the calcium-binding domain but the remainder of the protein-coding region shows substantial divergence.

Cytochrome P450(cin) (CYP176A), isolation, expression, and characterization

Cytochromes P450 are members of a superfamily of hemoproteins involved in the oxidative metabolism of various physiologic and xenobiotic compounds in eukaryotes and prokaryotes. Studies on bacterial P450s, particularly those involved in monoterpene oxidation, have provided an integral contribution to our understanding of these proteins, away from the problems encountered with eukaryotic forms. We report here a novel cytochrome P450 (P450(cin), CYP176A1) purified from a strain of Citrobacter braakii that is capable of using cineole 1 as its sole source of carbon and energy. This enzyme has been purified to homogeneity and the amino acid sequences of three tryptic peptides determined. By using this information, a PCR-based cloning strategy was developed that allowed the isolation of a 4-kb DNA fragment containing the cytochrome P450(cin) gene (cinA). Sequencing revealed three open reading frames that were identified on the basis of sequence homology as a cytochrome P450, an NADPH-dependent flavodoxin/ferrodoxin reductase, and a flavodoxin. This arrangement suggests that P450(cin) may be the first isolated P450 to use a flavodoxin as its natural redox partner. Sequencing also identified the unprecedented substitution of a highly conserved, catalytically, important active site threonine with an asparagine residue. The P450 gene was subcloned and heterologously expressed in Escherichia coli at similar to2000 nmol/liter of original culture, and purification was achieved by standard protocols. Postulating the native E. coli flavodoxin/flavodoxin reductase system might mimic the natural redox partners of P450,in, it was expressed in E. coli in the presence of cineole 1. A product was formed in vivo that was tentatively identified by gas chromatography-mass spectrometry as 2-hydroxycineole 2. Examination of P450(cin) by UV-visible spectroscopy revealed typical spectra characteristic of P450s, a high affinity for cineole 1 (K-D = 0.7 mum), and a large spin state change of the heme iron associated with binding of cineole 1. These facts support the hypothesis that cineole 1 is the natural substrate for this enzyme and that P450(cin) catalyzes the initial monooxygenation of cineole 1 biodegradation. This constitutes the first characterization of an enzyme involved in this pathway.

Molecular phylogenetics of Lentibulariaceae inferred from plastid rps16 Intron and trnL-F DNA sequences: implications for character evolution and biogeography

Phylogenetic relationships among 75 species of Lentibulariaceae, representing the three recognized genera, were assessed by cladistic analysis of DNA sequences from the plastid rps16 intron and the trnL-F region. Sequence data from the two loci were analyzed both separately and in combination. Consensus trees from all analyses are congruent, and parsimony jackknife results demonstrate strong support for relationships both between and within each of the three demonstrably monophyletic genera. The genus Pinguicula is sister to a Genlisea-Utricularia clade, the phylogenetic structure within this clade closely follows Taylor's recent sectional delimitations based on morphology. Three principal clades are shown within Utricularia, with the basal sections Polypoinpholyx and Pleiochasia together forming the sister lineage of the remaining Utricularia species. Of the fundamental morphological specializations, the stoloniferous growth form apparently arose independently within Genlisea and Utricularia three times, and within Utricularia itself, perhaps more than once. The epiphytic habit has evolved independently at least three times, in Pinguicula, in Utricularia section Phyllaria, and within the two sections Orchidioides and Iperua (in the latter as bromeliad tank-epiphytes). The suspended aquatic habit may have evolved independently within sections Utricularia and Vesiculina. Biogeographic optimization on the phylogeny demonstrates patterns commonly associated with the boreotropics hypothesis and limits the spatial origin of Lentibulariaceae to temperate Eurasia or tropical America.

Catalytically active Dengue virus NS3 protease forms aggregates that are separable by size exclusion chromatography

An active form of the Dengue virus protease NS3 (CF40.Gly.NS3pro) was expressed in Escherichia coli. This construct consists of a critical 40 amino acid cofactor domain from NS2B fused to the N-terminal 184 amino acid protease domain of NS3 via a flexible, covalent linker (Gly(4)SerGly(4)). The recombinantly produced protein is soluble and has a hexa-histidine tag engineered at the N-terminus for ease of purification using metal affinity chromatography. However, the presence of lower molecular weight impurities after affinity chromatography indicated the need for additional purification steps. The consistent appearance of these impurities suggested that they may be the products of proteolysis and/or auto-proteolysis. The latter possibility was subsequently excluded by the observation of the same impurities in a purified, catalytically inactive form of the recombinant protease (CF40.Gly.NS3pro.SA). Further analysis indicated that these impurities may represent premature translation termination products. Regardless of their origin, they were shown to form various sized aggregates with full-length CF40.Gly.NS3pro that can be separated by size exclusion chromatography, yielding fractions of active protease of sufficient purity for crystallisation trials. The ultimate goal of these studies is to obtain a crystal structure of a catalytically active form of the Dengue virus NS3 protease for structure-based drug design. (C) 2002 Elsevier Science (USA). All rights reserved.

Inferring genome trees by using a filter to eliminate phylogenetically discordant sequences and a distance matrix based on mean normalized BLASTP scores

Darwin's paradigm holds that the diversity of present-day organisms has arisen via a process of genetic descent with modification, as on a bifurcating tree. Evidence is accumulating that genes are sometimes transferred not along lineages but rather across lineages. To the extent that this is so, Darwin's paradigm can apply only imperfectly to genomes, potentially complicating or perhaps undermining attempts to reconstruct historical relationships among genomes (i.e., a genome tree). Whether most genes in a genome have arisen via treelike (vertical) descent or by lateral transfer across lineages can be tested if enough complete genome sequences are used. We define a phylogenetically discordant sequence (PDS) as an open reading frame (ORF) that exhibits patterns of similarity relationships statistically distinguishable from those of most other ORFs in the same genome. PDSs represent between 6.0 and 16.8% (mean, 10.8%) of the analyzable ORFs in the genomes of 28 bacteria, eight archaea, and one eukaryote (Saccharomyces cerevisiae). In this study we developed and assessed a distance-based approach, based on mean pairwise sequence similarity, for generating genome trees. Exclusion of PDSs improved bootstrap support for basal nodes but altered few topological features, indicating that there is little systematic bias among PDSs. Many but not all features of the genome tree from which PDSs were excluded are consistent with the 16S rRNA tree.

Behaviour of the black flying fox Pteropus alecto: 2. Territoriality and courtship

Flying foxes are commonly thought of as highly social mammals, yet little is known about the dynamics of their social interactions at a day roost. The aim of the present study was to examine the nature of the seasonal activities of territoriality and courtship amongst wild flying foxes in Australia. Focal observations were conducted at two permanent roosts of black flying foxes Pteropus alecto during periods of peak social interaction in the summers of 1999 and 2000 in urban Brisbane, Queensland. Observations of male territoriality were conducted at dawn and began eight weeks prior to the commencement of mating. The majority of defense bouts (87%) consisted of ritualised pursuit, while 13% of bouts involved physical contact expressed as either wrestling or hooking. One male with an unusually large territory took significantly longer to defend it than other males with less territory to defend. Observations of courtship revealed repetitive courtship sequences, including pre-copulatory approaches by the males, copulation attempts and grooming/resting periods. Thirty-four complete courtship sequences incorporating 135 copulation attempts were recorded over two seasons. Females actively resisted courtship approaches by males, forcing males to display a continuous determination to mate over time where determination can be considered an indicator of 'fitness'. The courtship bout length of females with suckling young was significantly longer ((x) over bar +/- SE; 230.9 +/- 22.16 s) than that of females unencumbered by large pups (158.5 +/- 9.69 s), although the length of copulations within those courtships was not (45.6 +/- 5.19 versus 36.2 +/- 3.43 s).

Production of the baculovirus-expressed dengue virus glycoprotein NS1 can be improved dramatically with optimised regimes for fed-batch cultures and the addition of the insect moulting hormone, 20-Hydroxyecdysone

A perennial problem in recombinant protein expression is low yield of the product of interest. A strategy which has been shown to increase the production of baculovirus-expressed proteins is to utilise fed-batch cultures. One disadvantage of this approach is the time-consuming task of optimising the feeding strategy. Previously, a statistical optimisation routine was applied to develop a feeding strategy that increased the yield of beta-Galactosidase (beta-Gal) by 2.4-fold (Biotechnol. Bioeng, 59 (1998) 178). This involves the single addition of nutrient concentrates (amino acids, lipids. glucose and yeastolate ultrafiltrate) into Sf9 cell cultures grown in SF900II medium. In this study, it is demonstrated that this optimised fed-batch strategy developed for a high-yielding intracellular product beta-Gal could be applied successfully to a relatively low-yielding glycosylated and secreted product such as the dengue virus glycoprotein NS1. Optimised batch infections yielded 4 mug/ml of NS1 at a peak cell density of 4.2 x 10(6) cells/ml. In contrast. optimised fed-batch infections exhibited a 3-fold improvement in yield, with 12 mug ml of NS1 produced at a peak cell density of 11.3 x 10(6) cells/ml. No further improvements in yield were recorded when the feed volumes were doubled and the peak cell density was increased to 23 x 10(6) cells/ml, unless the cultures were stimulated by the addition of 4 mug/ml of 20-Hydroxyecdysone (an insect moulting hormone). In this case, the NS1 yield was increased to 20 mug/ml. which was nearly 5-fold higher than optimised batch cultures. (C) 2002 Elsevier Science B.V. All rights reserved.

In vivo DNA electrotransfer

The use of electrotransfer for DNA delivery to prokaryotic cells, and eukaryotic cells in vitro, has been well known and widely used for many years. However, it is only recently that electric fields have been used to enhance DNA transfer to animal cells in vivo, and this is known as DNA electrotransfer or in vivo DNA electroporation. Some of the advantages of this method of somatic cell gene transfer are that it is a simple method that can be used to transfer almost any DNA construct to animal cells and tissues in vivo; multiple constructs can be co-transfected; it is equally applicable to dividing and nondividing cells; the DNA of interest does not need to be subeloned into a specific viral transfer vector and there is no need for the production of high titre viral stocks; and, as no viral genes are expressed there is less chance of an adverse immunologic reaction to vector sequences. The ease with which efficient in vivo gene transfer can be achieved with in vivo DNA electrotransfer is now allowing genetic analysis to be applied to a number of classic animal model systems where transgenic and embryonic stem cell techniques are not well developed, but for which a wealth of detailed descriptive embryological information is available, or surgical manipulation is much more feasible. As well as exciting applications in developmental biology, in vivo DNA electrotransfer is also being used to transfer genes to skeletal muscle and drive expression of therapeutically active proteins, and to examine exogenous gene and protein function in normal adult cells situated within the complex environment of a tissue and organ system in vivo. Thus, in effect providing the in vivo equivalent of the in vitro transient transfection assay. As the widespread use of in vivo electroporation has really only just begun, it is likely that the future will hold many more applications for this technology in basic research, biotechnology and clinical research areas.

A family of textilinin genes, two of which encode proteins with antihaemorrhagic properties

Two peptides, textilinins 1 and 2, isolated from the venom of the Australian common brown snake, Pseudonaja textilis textilis, are effective in preventing blood loss. To further investigate the potential of textilinins as anti-haemorrhagic agents, we cloned cDNAs encoding these proteins. The isolated full-length cDNA (430 bp in size) was shown to code for a 59 amino acid protein, corresponding in size to the native peptide, plus an additional 24 amino acid propeptide. Six such cDNAs were identified, differing in nucleotide sequence in the coding region but with an identical propeptide. All six sequences predicted peptides containing six conserved cysteines common to Kunitz-type serine protease inhibitors. When expressed as glutathione S-transferase (GST) fusion proteins and released by cleavage with thrombin, only those peptides corresponding to textilinin 1 and 2 were active in inhibiting plasmin with K-i values similar to those of their native counterparts and in binding to plasmin less tightly than aprotinin by two orders of magnitude. Similarly, in the mouse tail vein blood loss model only recombinant textilinin 1 and 2 were effective in reducing blood loss. These recombinant textilinins have potential as therapeutic agents for reducing blood loss in humans, obviating the need for reliance on aprotinin, a bovine product with possible risk of transmissible disease, and compromising the fibrinolytic system in a less irreversible manner.

MUC1 epithelial mucin (CD227) is expressed by activated dendritic cells

The MUC1 mucin (CD227) is a cell surface mucin originally thought to be restricted to epithelial tissues. We report that CD227 is expressed on human blood dendritic cells (DC) and monocyte-derived DC following in vitro activation. Freshly isolated murine splenic DC had very low levels of CD227; however, all DC expressed CD227 following in vitro culture. In the mouse spleen, CD227 was seen on clusters within the red pulp and surrounding the marginal zone in the white pulp. Additionally, we confirm CD227 expression by activated human T cells and show for the first time that the CD227 cytoplasmic domain is tyrosine-phosphorylated in activated T cells and DC and is associated with other phosphoproteins, indicating a role in signaling. The function of CD227 on DC and T cells requires further elucidation.

Prolactin-releasing peptide in ewe: cDNA cloning, mRNA distribution and effects on prolactin secretion in vitro and in vivo

RT-PCR followed by 5'- and 3'- rapid amplification of cDNA ends was used to clone and sequence ovine prolactin-releasing peptide (PrRP). The cDNA was characterised by short 5'- and 3'-untranslated regions and a GC-rich (71%) coding region. The nucleotide and deduced amino acid sequences for the coding region showed 95.6 and 94.9% identity with bovine PrRP but the amino acid sequence of PrRP31 was conserved between these species. Northern blot analysis and RT-PCR showed that, as in the rat, the peptide was more abundantly expressed in the brainstem than the hypothalamus. However, in the ovine hypothalamus, PrRP mRNA expression was more widespread than in the rat, with expression detected in both rostral and caudal parts of the mediobasal hypothalamus. The effects of synthetic ovine PrRP on prolactin secretion both in vitro and in vivo were also examined. In primary cultures of sheep pituitary cells, PrRP significantly (P