Computational tools for the study of allergens

Allergy is a major cause of morbidity worldwide. The number of characterized allergens and related information is increasing rapidly creating demands for advanced information storage, retrieval and analysis. Bioinformatics provides useful tools for analysing allergens and these are complementary to traditional laboratory techniques for the study of allergens. Specific applications include structural analysis of allergens, identification of B- and T-cell epitopes, assessment of allergenicity and cross-reactivity, and genome analysis. In this paper, the most important bioinformatic tools and methods with relevance to the study of allergy have been reviewed.

Determination of the solution structures of conantokin-G and conantokin-T by CD and NMR spectroscopy

Conantokin-G and conantokin-T are two paralytic polypeptide toxins originally isolated from the venom of the fish-hunting cone snails of the genus Conus. Conantokin-G and conantokin-T are the only naturally occurring peptidic compounds which possess N-methyl-D-aspartate receptor antagonist activity, produced by a selective non-competitive antagonism of polyamine responses, They are also structurally unusual in that they contain a disproportionately large number of acid labile post-translational gamma-carboxyglutamic acid (Gla) residues, Although no precise structural information has previously been published for these peptides, early spectroscopic measurements have indicated that both conantokin-G and conantokin-T form alpha-helical structures, although there is some debate whether the presence of calcium ions is required for these peptides to adopt this fold, We now report a detailed structural study of synthetic conantokin-G and conantokin-T in a range of solution conditions using CD and H-1 NMR spec troscopy. The three-dimensional structures of conantokin-T and conantokin-G were calculated from H-1 NMR-derived distance and dihedral restraints. Both conantokins were found to contain a mixture of alpha- and 3(10) helix, that give rise to curved and straight helical conformers. Conantokin-G requires the presence of divalent cations (Zn2+, Ca2+, Cu2+, Or Mg2+) to form a stable iv-helix, while conantokin-T adopts a stable alpha-helical structure in aqueous conditions, in the presence or absence of divalent cations (Zn2+, Ca2+, Cu2+, Or Mg2+).

The Relationships between West Nile and Kunjin Viruses

Until recently, West Nile (WN) and Kunjin (KUN) viruses were classified as distinct types in the Flavivirus genus. However, genetic and antigenic studies on isolates of these two viruses indicate that the relationship between them is more complex. To better define this relationship, we performed sequence analyses on 32 isolates of KUN virus and 28 isolates of WN virus from different geographic areas, including a WN isolate from the recent outbreak in New York. Sequence comparisons showed that the KUN virus isolates from Australia were tightly grouped but that the WN virus isolates exhibited substantial divergence and could be differentiated into four district groups. KUN virus isolates from Australia were antigenically homologous and distinct from the WN isolates and a Malaysian KUN virus. Our results suggest that KUN and WN viruses comprise a group of closely related viruses that can be differentiated into subgroups on the basis of genetic and antigenic analyses.

Recurrent gains and losses of large (84-109 bp) repeats in the rDNA internal transcribed spacer 2 (ITS2) of rhipicephaline ticks

We studied the internal transcribed spacer 2 (ITS2) in twenty-two spp. of ticks from the subfamily Rhipicephalinae. A 104-109 base pair (bp) region was Imperfectly repeated In most ticks studied. Mapping the number of repeat copies on to a phylogeny from the ITS2 showed that there have been many Independent gains and losses of repeats. Comparison of the sequences of the repeat copies Indicated that in most taxa concerted evolution had played little if any role in the evolution of these regions, as the copies clustered by sequence position rather than species, In our putative secondary structure, each repeat copy can fold into a distinct and almost identical stem-loop complex; a gain or loss of a repeat copy apparently does not impair the function of the ITS2 in these ticks.

Alterations in Plasmodium falciparum genotypes during sequential infections suggest the presence of strain specific immunity

Many of the asexual stage Plasmodium falciparum proteins that are the targets of host protective responses are markedly polymorphic. The full repertoire of diversity is not defined for any antigen. Most studies have focused on the genes encoding merozoite surface proteins 1 and 2 (MSP1, MSP2). We explored the extent of diversity of some of the less studied merozoite surface antigens and analyzed the degree of complexity of malaria field isolates by deriving nucleotide sequences of several antigens. We have determined the genotype of apical membrane antigen 1 (AMA1) in a group of 30 field samples, collected over 29 months, from individuals living in an area of intense malaria transmission in Irian Jaya, identifying 14 different alleles. AMA1 genotyping was combined with previously determined MSP2 typing. AMA1 had the greatest power in distinguishing between isolates but methodological problems, especially when mixed infections are present, suggest it is not an ideal typing target. MSP1, MSP3, and glutamate-rich protein genotypes were also determined from a smaller group of samples, and all results were combined to derive an extended antigenic haplotype. Within this subset of 10 patients, nine different genotypes could be discerned; however, five patients were all infected with the same strain. This strain was present in individuals from two separate villages and was still present 12 months later. This strain was predominant at the first time point but had disappeared at the fourth time point. This significant change in malaria genotypes could be due to strain-specific immunity developing in this population.

The mitochondrial 12S gene is a suitable marker of populations of Sarcoptes scabiei from wombats, dogs and humans in Australia

We sequenced part of the mitochondrial 12S ribosomal RNA gene of 23 specimens of Sarcoptes scabiei from eight wombats, one dog and three humans. Twelve of the 326 nucleotide positions varied among these mites and there were nine haplotypes (sequences) that differed by 1-8 nucleotides. Phylogenetic analyses indicated that these mites were from two lineages: (1) mites from wombats from Victoria, Australia, and mites from the humans and dog from the Northern Territory, Australia (haplotypes 1-4, 9); and (2) mites from the humans and dog from the Northern Territory (haplotypes 5-8). Mites from the three different hosts (wombats, a dog and humans) had not diverged phylogenetically; rather, these mites had similar 12S sequences. Thus, we conclude that these mites from wombats, humans and a dog are closely related, and that they diverged from a common ancestor relatively recently. This conclusion is consistent with the argument that people and/or their dogs introduced to Australia the S. scabiei mites that infect wombats Australia. So, S. scabiei, which has been blamed for the extinction of populations of wombats in Australia, may be a parasitic mite that was introduced to Australia with people and/or their dogs. These data show that the mitochondrial 12S rRNA gene may be a suitable population marker of S. scabiei from wombats, dogs and humans in Australia.

Conformationally constrained macrocycles that mimic tripeptide beta-strands in water and aprotic solvents

The beta-strand conformation is unknown for short peptides in aqueous solution, yet it is a fundamental building block in proteins and the crucial recognition motif for proteolytic enzymes that enable formation and turnover of all proteins. To create a generalized scaffold as a peptidomimetic that is preorganized in a beta-strand, we individually synthesized a series of 15-22-membered macrocyclic analogues of tripeptides and analyzed their structures. Each cycle is highly constrained by two trans amide bonds and a planar aromatic ring with a short nonpeptidic linker between them. A measure of this ring strain is the restricted rotation of the component tyrosinyl aromatic ring (DeltaG(rot) 76.7 kJ mol(-1) (16-membered ring), 46.1 kJ mol(-1) (17-membered ring)) evidenced by variable temperature proton NMR spectra (DMF-d(7), 200-400 K). Unusually large amide coupling constants ((3)J(NH-CHalpha) 9-10 Hz) corresponding to large dihedral angles were detected in both protic and aprotic solvents for these macrocycles, consistent with a high degree of structure in solution. The temperature dependence of all amide NH chemical shifts (Deltadelta/T7-12 ppb/deg) precluded the presence of transannular hydrogen bonds that define alternative turn structures. Whereas similar sized conventional cyclic peptides usually exist in solution as an equilibrium mixture of multiple conformers, these macrocycles adopt a well-defined beta-strand structure even in water as revealed by 2-D NMR spectral data and by a structure calculation for the smallest (15-membered) and most constrained macrocycle. Macrocycles that are sufficiently constrained to exclusively adopt a beta-strand-mimicking structure in water may be useful pre-organized and generic templates for the design of compounds that interfere with beta-strand recognition in biology.

The value of idiosyncratic markers and changes to conserved tRNA sequences from the mitochondrial genome of hard ticks (Acari: Ixodida: Ixodidae) for phylogenetic inference

Idiosyncratic markers are features of genes and genomes that are so unusual that it is unlikely that they evolved more than once in a lineage of organisms. Here we explore further the potential of idiosyncratic markers and changes to typically conserved tRNA sequences for phylogenetic inference. Hard ticks were chosen as the model group because their phylogeny has been studied extensively. Fifty-eight candidate markers from hard ticks ( family Ixodidae) and 22 markers from the subfamily Rhipicephalinae sensu lato were mapped onto phylogenies of these groups. Two of the most interesting markers, features of the secondary structure of two different tRNAs, gave strong support to the hypothesis that species of the Prostriata ( Ixodes spp.) are monophyletic. Previous analyses of genes and morphology did not strongly support this relationship, instead suggesting that the Prostriata is paraphyletic with respect to the Metastriata ( the rest of the hard ticks). Parallel or convergent evolution was not found in the arrangements of mitochondrial genes in ticks nor were there any reversals to the ancestral arthropod character state. Many of the markers identified were phylogenetically informative, whereas others should be informative with study of additional taxa. Idiosyncratic markers and changes to typically conserved nucleotides in tRNAs that are phylogenetically informative were common in this data set, and thus these types of markers might be found in other organisms.

Insights into the structural basis for zinc inhibition of the glycine receptor

Histidines 107 and 109 in the glycine receptor ( GlyR) alpha(1) subunit have previously been identified as determinants of the inhibitory zinc-binding site. Based on modeling of the GlyR alpha(1) subunit extracellular domain by homology to the acetylcholine-binding protein crystal structure, we hypothesized that inhibitory zinc is bound within the vestibule lumen at subunit interfaces, where it is ligated by His(107) from one subunit and His(109) from an adjacent subunit. This was tested by co-expressing alpha(1) subunits containing the H107A mutation with alpha(1) subunits containing the H109A mutation. Although sensitivity to zinc inhibition is markedly reduced when either mutation is individually incorporated into all five subunits, the GlyRs formed by the co-expression of H107A mutant subunits with H109A mutant subunits exhibited an inhibitory zinc sensitivity similar to that of the wild type alpha(1) homomeric GlyR. This constitutes strong evidence that inhibitory zinc is coordinated at the interface between adjacent alpha(1) subunits. No evidence was found for beta subunit involvement in the coordination of inhibitory zinc, indicating that a maximum of two zinc-binding sites per alpha(1)beta receptor is sufficient for maximal zinc inhibition. Our data also show that two zinc-binding sites are sufficient for significant inhibition of alpha(1) homomers. The binding of zinc at the interface between adjacent alpha(1) subunits could restrict intersubunit movements, providing a feasible mechanism for the inhibition of channel activation by zinc.

The production, purification and crystallization of a pocilloporin pigment from a reef-forming coral

Reef-building corals contain fluorescent pigments termed pocilloporins that function by regulating the light environment of coral and acting as a photoprotectant in excessive sunlight. These pocilloporins are related to the monomeric green fluorescent protein and the tetrameric DsRed fluorescent proteins, which have widespread use as biotechnological tools. An intensely blue-coloured pocilloporin, termed Rtms5, was expressed in Escherichia coli, purified and crystallized. Rtms5 was shown to be tetrameric, with deep blue crystals that diffract to 2.2 Angstrom resolution and belong to space group I4(1)22. The colour of this pocilloporin was observed to be sensitive to pH and a yellow (pH 3.5) and a red form (pH 4.5) of Rtms5 were also crystallized. These crystals belong to space group P4(2)22 and diffract to 2.4 Angstrom resolution or better.

Searching oligo sets of human chromosome 12 using evolutionary strategies

DNA Microarray is a powerful tool to measure the level of a mixed population of nucleic acids at one time, which has great impact in many aspects of life sciences research. In order to distinguish nucleic acids with very similar composition by hybridization, it is necessary to design microarray probes with high specificities and sensitivities. Highly specific probes correspond to probes having unique DNA sequences; whereas highly sensitive probes correspond to those with melting temperature within a desired range and having no secondary structure. The selection of these probes from a set of functional DNA sequences (exons) constitutes a computationally expensive discrete non-linear search problem. We delegate the search task to a simple yet effective Evolution Strategy algorithm. The computational efficiency is also greatly improved by making use of an available bioinformatics tool.

Correlating gene expression with larval competence and the effect of age and parentage on metamorphosis in the tropical abalone Haliotis asinina

The non-geniculate crustose coralline alga (CCA) Mastophora pacifica can induce the metamorphosis of competent Haliotis asinina (Vetigastropoda) larvae. The ability to respond to this natural cue varies considerably with larval age, with a higher proportion of older larvae (e.g. 90 h) able to metamorphose in response to M. pacifica than younger larvae (e.g. 66 h). Here we document the variation in time to acquisition of competence within a larval age class. For example, after 18 h of exposure to M. pacifica, approximately 15 and 36% of 84 and 90-h-old H. asinina larvae had initiated metamorphosis, respectively. This age-dependent response to M. pacifica is also observed when different aged larvae are exposed to CCA for varying periods. A higher proportion of older larvae require shorter periods of exposure to CCA than younger larvae in order to initiate metamorphosis. In this experiment, as in the previous, a small proportion of young larvae were able to respond to brief periods of CCA exposure, suggesting that they had developed the same state of competency as the majority of their older counterparts. Comparisons of the proportions of larvae undergoing metamorphosis between families reveals that parentage also has a significant (P < 0.05) affect on whether an individual will initiate metamorphosis at a given age. These familial differences are more pronounced when younger, largely pre-competent larvae (i.e. 66 h old) are exposed to M. pacifica, with proportions of larvae undergoing metamorphosis differing by as much as 10 fold between families. As these data suggest that variation in the rate of development of the competent state has a genetic basis, and as a first step towards identifying the molecular basis to this variation, we have identified numerous genes that are differentially expressed later in larval development using a differential display approach. Spatial expression analysis of these genes suggests that they may be directly involved in the acquisition of competence, or may play a functional role in the postlarva following metamorphosis.

Accurate prediction of scorpion toxin functional properties from primary structures

Scorpion toxins are common experimental tools for studies of biochemical and pharmacological properties of ion channels. The number of functionally annotated scorpion toxins is steadily growing, but the number of identified toxin sequences is increasing at much faster pace. With an estimated 100,000 different variants, bioinformatic analysis of scorpion toxins is becoming a necessary tool for their systematic functional analysis. Here, we report a bioinformatics-driven system involving scorpion toxin structural classification, functional annotation, database technology, sequence comparison, nearest neighbour analysis, and decision rules which produces highly accurate predictions of scorpion toxin functional properties. (c) 2005 Elsevier Inc. All rights reserved.

Novel mitochondrial gene content and gene arrangement indicate illegitimate inter-mtDNA recombination in the chigger mite, Leptotrombidium pallidum

To better understand the evolution of mitochondrial (mt) genomes in the Acari (mites and ticks), we sequenced the mt genome of the chigger mite, Leptotrombidium pallidum (Arthropoda: Acari: Acariformes). This genome is highly rearranged relative to that of the hypothetical ancestor of the arthropods and the other species of Acari studied. The mt genome of L. pallidum has two genes for large subunit rRNA, a pseudogene for small subunit rRNA, and four nearly identical large noncoding regions. Nineteen of the 22 tRNAs encoded by this genome apparently lack either a T-arm or a D-arm. Further, the mt genome of L. pallidum has two distantly separated sections with identical sequences but opposite orientations of transcription. This arrangement cannot be accounted for by homologous recombination or by previously known mechanisms of mt gene rearrangement. The most plausible explanation for the origin of this arrangement is illegitimate inter-mtDNA recombination, which has not been reported previously in animals. In light of the evidence from previous experiments on recombination in nuclear and mt genomes of animals, we propose a model of illegitimate inter-mtDNA recombination to account for the novel gene content and gene arrangement in the mt genome of L. pallidum.

hnRNP A2, a potential ssDNA/RNA molecular adapter at the telomere

The heterogeneous nuclear ribonucleoprotein (hnRNP) A2 is a multi-tasking protein that acts in the cytoplasm and nucleus. We have explored the possibility that this protein is associated with telomeres and participates in their maintenance. Rat brain hnRNP A2 was shown to have two nucleic acid binding sites. In the presence of heparin one site binds single-stranded oligodeoxyribonucleotides irrespective of sequence but not the corresponding oligoribonucleotides. Both the hnRNP A2-binding cis-acting element for the cytoplasmic RNA trafficking element, A2RE, and the ssDNA telomere repeat match a consensus sequence for binding to a second sequence-specific site identified by mutational analysis. hnRNP A2 protected the telomeric repeat sequence, but not the complementary sequence, against DNase digestion: the glycine-rich domain was found to be necessary, but not sufficient, for protection. The N-terminal RRM (RNA recognition motif) and tandem RRMs of hnRNP A2 also bind the single-stranded, template-containing segment of telomerase RNA. hnRNP A2 colocalizes with telomeric chromatin in the subset of PML bodies that are a hallmark of ALT cells, reinforcing the evidence for hnRNPs having a role in telomere maintenance. Our results support a model in which hnRNP A2 acts as a molecular adapter between single-stranded telomeric repeats, or telomerase RNA, and another segment of ssDNA.