Sampling phylogenetic tree space with the generalized Gibbs sampler

The generalized Gibbs sampler (GGS) is a recently developed Markov chain Monte Carlo (MCMC) technique that enables Gibbs-like sampling of state spaces that lack a convenient representation in terms of a fixed coordinate system. This paper describes a new sampler, called the tree sampler, which uses the GGS to sample from a state space consisting of phylogenetic trees. The tree sampler is useful for a wide range of phylogenetic applications, including Bayesian, maximum likelihood, and maximum parsimony methods. A fast new algorithm to search for a maximum parsimony phylogeny is presented, using the tree sampler in the context of simulated annealing. The mathematics underlying the algorithm is explained and its time complexity is analyzed. The method is tested on two large data sets consisting of 123 sequences and 500 sequences, respectively. The new algorithm is shown to compare very favorably in terms of speed and accuracy to the program DNAPARS from the PHYLIP package.

Solution structure of robustoxin, the lethal neurotoxin from the funnel-web spider Atrax robustus

The solution structure of robustoxin, the lethal neurotoxin from the Sydney funnel-web spider Atrax robustus, has been determined from 2D H-1 NMR data, Robustoxin is a polypeptide of 42 residues cross-linked by four disulphide bonds, the connectivities of which were determined from NMR data and trial structure calculations to be 1-15, 8-20, 14-31 and 16-42 (a 1-4/2-6/3-7/5-8 pattern), The structure consists of a small three-stranded, anti-parallel beta-sheet and a series of interlocking gamma-turns at the C-terminus. It also contains a cystine knot, thus placing it in the inhibitor cystine knot motif family of structures, which includes the omega-conotoxins and a number of plant and animal toxins and protease inhibitors. Robustoxin contains three distinct charged patches on its surface, and an extended loop that includes several aromatic and non-polar residues, Both of these structural features may play a role in its binding to the voltage-gated sodium channel. (C) 1997 Federation of European Biochemical Societies.

Sequence and tissue distribution of chicken cellular nucleic acid binding protein cDNA

We sequenced cDNAs coding for chicken cellular nucleic acid binding protein (CNBP). Two slightly different variations of the open reading frame were found, each of which translates into a protein with seven zinc finger domains. The longest transcript contains an in-frame insert of 3 bp. The sequence conservation between chick CNBP cDNAs with human, rat and mouse CNBP cDNAs is extreme, especially in the coding region, where the deduced amino acid sequence identity with human, rat and mouse CNBP is 99%. CNBP-like transcripts were also found in various tissues from insect, shrimp, fish and lizard. Regions with remarkable nucleotide conservation were also found in the 3' untranslated region, indicating important functions for these regions. Quantitative reverse transcription polymerase chain reaction (RT-PCR) indicated that in the chick, CNBP is present in all tissues examined in approximately equal ratios to total RNA. RT-PCR of total RNA isolated from different phyla indicate CNBP-like proteins art widespread throughout the animal kingdom. The extraordinary level of conservation suggests an important physiological role for CNBP. (C) 1997 Elsevier Science Inc.

Solution structure and proposed binding mechanism of a novel potassium channel toxin kappa-conotoxin PVIIA

Background: kappa-PVIIA is a 27-residue polypeptide isolated from the venom of Conus purpurascens and is the first member of a new class of conotoxins that block potassium channels. By comparison to other ion channels of eukaryotic cell membranes, voltage-sensitive potassium channels are relatively simple and methodology has been developed for mapping their interactions with small-peptide toxins, PVIIA, therefore, is a valuable new probe of potassium channel structure. This study of the solution structure and mode of channel binding of PVIIA forms the basis for mapping the interacting residues at the conotoxin-ion channel interface. Results: The three-dimensional structure of PVIIA resembles the triple-stranded beta sheet/cystine-knot motif formed by a number of toxic and inhibitory peptides. Subtle structural differences, predominantly in loops 2 and 4, are observed between PVIIA and other conotoxins with similar structural frameworks, however. Electrophysiological binding data suggest that PVIIA blocks channel currents by binding in a voltage-sensitive manner to the external vestibule and occluding the pore, Comparison of the electrostatic surface of PVIIA with that of the well-characterised potassium channel blocker charybdotoxin suggests a likely binding orientation for PVIIA, Conclusions: Although the structure of PVIIA is considerably different to that of the alpha K scorpion toxins, it has a similar mechanism of channel blockade. On the basis of a comparison of the structures of PVIIA and charybdotoxin, we suggest that Lys19 of PVIIA is the residue which is responsible for physically occluding the pore of the potassium channel.

Protein fold recognition without Boltzmann statistics or explicit physical basis

We present a fast method for finding optimal parameters for a low-resolution (threading) force field intended to distinguish correct from incorrect folds for a given protein sequence. In contrast to other methods, the parameterization uses information from >10(7) misfolded structures as well as a set of native sequence-structure pairs. In addition to testing the resulting force field's performance on the protein sequence threading problem, results are shown that characterize the number of parameters necessary for effective structure recognition.

Characterization of calcium-activated chloride channels in patches excised from the dendritic knob of mammalian olfactory receptor neurons

We investigated the properties of calcium-activated chloride channels in inside-out membrane patches from the dendritic knobs of acutely dissociated rat olfactory receptor neurons. Patches typically contained large calcium-activated currents, with total conductances in the range 30-75 nS. The dose response curve for calcium exhibited an EC50 of about 26 mu M. In symmetrical NaCl solutions, the current-voltage relationship reversed at 0 mV and was linear between -80 and +70 mV. When the intracellular NaCl concentration was progressively reduced from 150 to 25 mM, the reversal potential changed in a manner consistent with a chloride-selective conductance. Indeed, modeling these data with the Goldman-Hodgkin-Katz equation revealed a P-Na/P-Cl of 0.034. The halide permeability sequence was P-Cl > P-F > P-I > P-Br indicating that permeation through the channel was dominated by ion binding sites with a high field strength. The channels were also permeable to the large organic anions, SCN-, acetate(-), and gluconate(-), with the permeability sequence P-Cl > P-SCN > gluconaie. Significant permeation to gluconate ions suggested that the channel pore had a minimum diameter of at least 5.8 Angstrom.

MHCPEP, a database of MHC-binding peptides: update 1997

MHCPEP (http://wehih.wehi.edu.au/mhcpep/) is a curated database comprising over 13 000 peptide sequences known to bind MHC molecules, Entries are compiled from published reports as well as from direct submissions of experimental data, Each entry contains the peptide sequence, its MHC specificity and where available, experimental method, observed activity, binding affinity, source protein and anchor positions, as well as publication references, The present format of the database allows text string matching searches but can easily be converted for use in conjunction with sequence analysis packages. The database can be accessed via Internet using WWW or FTP.

Calcium is essential for the structural integrity of the cysteine-rich, ligand-binding repeat of the low-density lipoprotein receptor

Seven cysteine-rich repeats form the ligand-binding region of the low-density lipoprotein (LDL) receptor. Each of these repeats is assumed to bind a calcium ion, which is needed for association of the receptor with its ligands, LDL and beta-VLDL. The effects of metal ions on the folding of the reduced N-terminal cysteine-rich repeat have been examined by using reverse-phase high-performance liquid chromatography to follow the formation of fully oxidized isomers with different disulfide connectivities. in the absence of calcium many of the 15 possible isomers formed on oxidation, whereas in its presence the predominant product at equilibrium had the native disulfide bond connectivities. Other metals were far less effective at directing disulfide bond formation: Mn2+ partly mimicked the action of Ca2+, but Ba2+, Sr2+, and Mg2+ had little effect. This metal-ion specificity was also observed in two-dimensional H-1 NMR spectral studies: only Ca2+ induced the native three-dimensional fold. The two paramagnetic ions, Gd3+ and Mn2+, and Cd2+ did not promote adoption of a well-defined structure, and the two paramagnetic ions did not displace calcium ions. The location of calcium ion binding sites in the repeat was also explored by NMR spectroscopy. The absence of chemical shift changes for the side chain proton resonances of Asp26, Asp36, and Glu37 from pH 3.9 to 6.8 in the presence of calcium ions and their proximal location in the NMR structures implicated these side chains as calcium ligands. Deuterium exchange NMR experiments also revealed a network of hydrogen bonds that stabilizes the putative calcium-binding loop.

Structural and functional characterisation of human sulfotransferases

The human aryl sulfotransferases HAST4 and HAST4v vary by only two amino acids but exhibit markedly different affinity towards the sulfonate acceptor p-nitrophenol and the sulfonate donor 3'-phosphoadenosine-5'-phosphosulfate (PAPS). To determine the importance of each of these amino acid differences, chimeric constructs were made of HAST4 and HAST4v. By attaching the last 120 amino acids of HAST-4v to HAST4 (changing Thr235 to Asn235) we have been able to produce a protein that has a K-m for PAPS similar to HAST4v. The reverse construct, HAST4v/4 produces a protein with a K-m for PAPS similar to HAST4. These data suggests that the COOH-terminal of sulfotransferases is involved in co-factor binding. (C) 1998 Elsevier Science Ireland Ltd. All rights reserved.

Bacterial expression of two human aryl sulfotransferases

The effect of replacing a single codon in the N-terminal of human aryl sulfotransferase (HAST) 1 and 3 with one that is more commonly found in E. coli genes was assessed. The pKK233-2 E. coli expression vector was employed and the polymerase chain reaction (PCR) was used to introduce the 5' nucleotide substitution, at the same time maintaining the fidelity of the amino acid sequence. The data indicates that this change had a minimal effect on protein production, subcellular localization or, in the case of HAST3, catalytic activity. In general, the pKK233-2 E. coli vector has been less than optimal for expressing human sulfotransferase cDNAs. (C) 1998 Elsevier Science Ireland Ltd. All rights reserved.

Localisation of aryl sulfotransferase expression in human tissues using hybridisation histochemistry and immunohistochemistry

To date, the laboratory has cloned seven unique human sulfotransferases; five aryl sulfotransferases (HAST1, HAST2, HAST3, HAST4 and HAST4v), an estrogen sulfotransferase and a dehydroepiandrosterone sulfotransferase. The cellular distribution of human aryl sulfotransferases in human hepatic and extrahepatic tissues has been determined using the techniques of hybridization histochemistry and immunohistochemistry. Human aryl sulfotransferase expression was detected in liver, epithelial cells of the gastrointestinal mucosal layer, epithelial cells lining bronchioles and in mammary duct epithelial cells. (C) 1998 Elsevier Science Ireland Ltd. All rights reserved.

Carrageenans with complex substitution patterns from red algae of the genus Erythroclonium

Extracellular polysaccharides from three Erythroclonium spp. were shown, by a combination of compositional, linkage analyses, and Fourier transform infrared and C-13-nuclear magnetic resonance spectroscopy, to be highly substituted carrageenans with at least five types of repeating disaccharide units. These are the carrabiose 2,4'-disulfate of iota-carrageenan, carrabiose 2-sulfate of alpha-carrageenan, the 6'-O-methylated counterparts of each of these repeating units, and 4',6'-O-(1-carboxyethylidene)carrabiose 2-sulfate. The polysaccharides also contain significant amounts of unsubstituted, 4-linked galactopyranose and small amounts of 4-linked 3-O-methylgalactopyranose and terminal glycosyl residues. The carrageenan preparations of the three species are similar, differing only in the proportions of some components. (C) 1998 Elsevier Science Ltd.

Molecules and morphology in phylogenetic studies of the Hemiuroidea (Digenea : Trematoda : Platyhelminthes)

Phylogenies of trematodes based on characters derived from morphology and life cycles have been controversial. Here, we add molecular data to the phylogenetic study of a group of trematodes, members of the superfamily Hemiuroidea Looss, 1899. DNA sequences from the V4 domain of the nuclear small subunit (18S) rRNA gene and a matrix of morphological characters modified from a previous study were used. There was no significant incongruence between the molecular and the morphological data. However, this was probably due largely to the limited resolving power of the morphological data. Analyses support a monophyletic Hemiuroidea containing at least the families Accacoeliidae, Derogenidae, Didymozoidae, Hirudinellidae, Sclerodistomidae, Syncoeliidae, Isoparorchiidae, Lecithasteridae, and Hemiuridae. These families fall into two principal clades. One contains the first six families and the other the Hemiuridae and lecithasterine lecithasterids. The positions of the hysterolecithine lecithasterids and the Isoparorchiidae were poorly resolved. The Ptychogonimidae may be the sister group of the remaining Hemiuroidea, but there was no support from the molecular data for the placement of the Azygiidae within the superfamily. (C) 1998 Academic Press.