Computational immunology: The coming of age

The explosive growth in biotechnology combined with major advancesin information technology has the potential to radically transformimmunology in the postgenomics era. Not only do we now have readyaccess to vast quantities of existing data, but new data with relevanceto immunology are being accumulated at an exponential rate. Resourcesfor computational immunology include biological databases and methodsfor data extraction, comparison, analysis and interpretation. Publiclyaccessible biological databases of relevance to immunologists numberin the hundreds and are growing daily. The ability to efficientlyextract and analyse information from these databases is vital forefficient immunology research. Most importantly, a new generationof computational immunology tools enables modelling of peptide transportby the transporter associated with antigen processing (TAP), modellingof antibody binding sites, identification of allergenic motifs andmodelling of T-cell receptor serial triggering.

Methods for prediction of peptide binding to MHC molecules: A comparative study

Background: A variety of methods for prediction of peptide binding to major histocompatibility complex (MHC) have been proposed. These methods are based on binding motifs, binding matrices, hidden Markov models (HMM), or artificial neural networks (ANN). There has been little prior work on the comparative analysis of these methods. Materials and Methods: We performed a comparison of the performance of six methods applied to the prediction of two human MHC class I molecules, including binding matrices and motifs, ANNs, and HMMs. Results: The selection of the optimal prediction method depends on the amount of available data (the number of peptides of known binding affinity to the MHC molecule of interest), the biases in the data set and the intended purpose of the prediction (screening of a single protein versus mass screening). When little or no peptide data are available, binding motifs are the most useful alternative to random guessing or use of a complete overlapping set of peptides for selection of candidate binders. As the number of known peptide binders increases, binding matrices and HMM become more useful predictors. ANN and HMM are the predictive methods of choice for MHC alleles with more than 100 known binding peptides. Conclusion: The ability of bioinformatic methods to reliably predict MHC binding peptides, and thereby potential T-cell epitopes, has major implications for clinical immunology, particularly in the area of vaccine design.

Bioinformatics for characterisation of allergens, allergenicity and allergic crossreactivity

Allergies are a major cause of chronic ill health in industrialised countries with the incidence of reported cases steadily increasing. This Research Focus details how bioinformatics is transforming the field of allergy through providing databases for management of allergen data, algorithms for characterisation of allergic crossreactivity, structural motifs and B- and T-cell epitopes, tools for prediction of allergenicity and techniques for genomic and proteomic analysis of allergens.

trans-Cyano(6-methyl-1,4,8,11-tetraazacyclotetradecan-6-amine)cobalt(III) bis(perchlorate) hydrate and trans-hydroxo(6-methyl-1,4,8,11-tetraazacyclotetradecan-6-amine)cobalt(III) bis(perchlorate)

The crystal structures of a pair of closely related macrocyclic cyano- and hydroxopentaaminecobalt(III) complexes, as their perchlorate salts, are reported. Although the two complexes, [Co(CN)(C11H27N5)](ClO4)2.H2O and [Co(OH)(C11H27N5)](ClO4)(2), exhibit similar conformations, significant differences in the Co-N bond lengths arise from the influence of the sixth ligand (cyano as opposed to hydroxo). The ensuing hydrogen-bonding patterns are also distinctly different. Disorder in the perchlorate anions was clearly resolved and this was rationalized on the basis of distinct hydrogen-bonding motifs involving the anion O atoms and the N-H and O-H donors.

Expression analysis of putative vitellogenin and lipophorin receptors in honey bee (Apis mellifera L.) queens and workers

Two members of the low density lipoprotein receptor (LDLR) family were identified as putative orthologs for a vitellogenin receptor (Amvgr) and a lipophorin receptor (Amlpr) in the Apis mellifera genome. Both receptor sequences have the structural motifs characteristic of LDLR family members and show a high degree of similarity with sequences of other insects. RT-PCR analysis of Amvgr and Amlpr expression detected the presence of both transcripts in different tissues of adult female (ovary, fat body, midgut, head and specifically hypopharyngeal gland), as well as in embryos. In the head RNA samples we found two variant forms of AmLpR: a full length one and a shorter one lacking 29 amino acids in the O-linked sugar domain. In ovaries the expression levels of the two honey bee LDLR members showed opposing trends: whereas Amvgr expression was upregulated as the ovaries became activated, Amlpr transcript levels gradually declined. In situ hybridization analysis performed on ovaries detected Amvgr mRNA exclusively in germ line cells and corroborated the qPCR results showing an increase in Amvgr gene expression concomitant with follicle growth. (C) 2008 Elsevier Ltd. All rights reserved.

Identification of a juvenile hormone esterase-like gene in the honey bee, Apis mellifera L. - Expression analysis and functional assays

Tight control over circulating juvenile hormone (JH) levels is of prime importance in an insect`s life cycle. Consequently, enzymes involved in JH metabolism, especially juvenile hormone esterases (JHEs), play major roles during metamorphosis and reproduction. In the highly eusocial Hymenoptera, JH has been co-opted into additional functions, primarily in the development of the queen and worker castes and in age-related behavioral development of workers. Within a set of 21 carboxylesterases predicted in the honey bee genome we identified one gene (Amjhe-like) that contained the main functional motifs of insect JHEs. Its transcript levels during larval development showed a maximum at the switch from feeding to spinning behavior, coinciding with a JH titer minimum. In adult workers, the highest levels were observed in nurse bees, where a low JH titer is required to prevent the switch to foraging. Functional assays showed that Amjhe-like expression is induced by JH-III and suppressed by 20-hydroxyecdysone. RNAi-mediated silencing of Amjhe-like gene function resulted in a six-fold increase in the JH titer in adult worker bees. The temporal profile of Amjhe-like expression in larval and adult workers, the pattern of hormonal regulation and the knockdown phenotype are consistent with the function of this gene as an authentic JHE. (C) 2008 Elsevier Inc. All rights reserved.

Computational methods for prediction of T-cell epitopes - a framework for modelling, testing, and applications

Computational models complement laboratory experimentation for efficient identification of MHC-binding peptides and T-cell epitopes. Methods for prediction of MHC-binding peptides include binding motifs, quantitative matrices, artificial neural networks, hidden Markov models, and molecular modelling. Models derived by these methods have been successfully used for prediction of T-cell epitopes in cancer, autoimmunity, infectious disease, and allergy. For maximum benefit, the use of computer models must be treated as experiments analogous to standard laboratory procedures and performed according to strict standards. This requires careful selection of data for model building, and adequate testing and validation. A range of web-based databases and MHC-binding prediction programs are available. Although some available prediction programs for particular MHC alleles have reasonable accuracy, there is no guarantee that all models produce good quality predictions. In this article, we present and discuss a framework for modelling, testing, and applications of computational methods used in predictions of T-cell epitopes. (C) 2004 Elsevier Inc. All rights reserved.

The use of electrospray ionization tandem mass spectrometry on the structural characterization of novel asymmetric metallo-organic supermolecules, based on pentafluorophenylporphyrins and ruthenium complexes

The novel asymmetric metallo-organic triads cis- and trans-[B(4-py)BPFPH(2){Ru(3)O(Ac)(6)(py)(2)}(Ru(bpy)(2)Cl}](PF(6))(2) (5a,b) for which cis- and trans-B(4-py)BPFPH(2)=5,10-bis(pentafluorophenyl)-15,20-bis(4-pyridyl)porphyrin and 5,15-bis(pentafluorophenyl)-10,20-bis(4-pyridyl)porphyrin, respectively; Ac = acetate; py = pyridine and bpy = 2,2`-bipyridine, as well as their corresponding monosubstituted dyads cis- and trans-[B(4-py)BPFPH(2){Ru(3)O(Ac)(6)(py)(2)}]PF(6) (4a,b) have been structurally characterized via electrospray ionization mass spectrometry (ESI-MS and ESI-MS/MS). The ESI-MS of dyads 4a,b display two characteristic Ru-multicomponent clusters of isotopologue ions corresponding to singly charged ions 4a,b(+) of m/z 1629 and doubly charged ions [4a,b+H](2+) of m/z 815 and the triads 5a,b are detected by ESI-MS as the intact doubly charged cluster of isotopologue ions of m/z 1039 [5a,b](2+). The ESI-MS/MS of 4a,b(+), [4a,b+H](2+) and [5a,b](2+) reveal characteristic dissociation pathways, which confirm the structural assignments providing additional information on the intrinsic binding strengths of the gaseous ions. Although the gas-phase behavior of each pair of isomers was rather similar, the less symmetric dyads 4a,b are distinguished via the (1)H NMR spectral profile of the pyrrolic signals. Exploratory photophysical assays have shown that both modifying motifs alter the porphyrinic core emission profile, opening the possibility to use these asymmetric systems as photophysical devices. (C) 2008 Elsevier Ltd. All rights reserved.

PREDBALB/c: a system for the prediction of peptide binding to H2(d) molecules, a haplotype of the BALB/c mouse

PREDBALB/c is a computational system that predicts peptides binding to the major histocompatibility complex-2 (H2(d)) of the BALB/c mouse, an important laboratory model organism. The predictions include the complete set of H2(d) class I ( H2-K-d, H2-L-d and H2-D-d) and class II (I-E-d and I-A(d)) molecules. The prediction system utilizes quantitative matrices, which were rigorously validated using experimentally determined binders and non-binders and also by in vivo studies using viral proteins. The prediction performance of PREDBALB/c is of very high accuracy. To our knowledge, this is the first online server for the prediction of peptides binding to a complete set of major histocompatibility complex molecules in a model organism (H2(d) haplotype). PREDBALB/c is available at http://antigen.i2r.a-star.edu.sg/predBalbc/.

The myosin-I-binding protein Acan125 binds the SH3 domain and belongs to the superfamily of leucine-rich repeat proteins

The SH3 domains of src and other nonreceptor tyrosine kinases have been shown to associate with the motif PXXP, where P and X stand for proline and an unspecified amino acid, but a motif that binds to the SH3 domain of myosin has thus far not been characterized. We previously showed that the SH3 domain of Acanthamoeba myosin-IC interacts with the protein Acan125. We now report that the Acan125 protein sequence contains two tandem consensus PXXP motifs near the C terminus. To test for binding, we expressed a polypeptide, AD3p, which includes 344 residues of native C-terminal sequence and a mutant polypeptide, AD3 Delta 977-994p, which lacks the sequence RPKPVPPPRGAKPAPPPR containing both PXXP motifs. The SH3 domain of Acanthamoeba myosin-IC bound AD3p and not AD3 Delta 977-994p, showing that the PXXP motifs are required for SH3 binding. The sequence of Acan125 is related overall to a protein of unknown function coded by Caenorhabditis elegans gene K07G5.1. The K07G5.1 gene product contains a proline-rich segment similar to the SH3 binding motif found in Acan125. The aligned sequences show considerable conservation of leucines and other hydrophobic residues, including the spacing of these residues, which matches a motif for leucine-rich repeats (LRRs). LRR domains have been demonstrated to be sites for ligand binding. Having an LRR domain and an SH3-binding domain, Acan125 and the C. elegans homologue define a novel family of bifunctional binding proteins.

Peginterferon alfa-2b and Ribavirin: Effective in Patients With Hepatitis C Who Failed Interferon alfa/Ribavirin Therapy

Background & Aims: Treatment with peginterferon alfa and ribavirin produces a sustained virologic response (SVR) in approximately 60% of hepatitis C virus (HCV)-infected patients. Alternate options are needed for patients who relapse or do not respond to therapy. Methods: This prospective, international, multicenter, open-label study evaluated efficacy and safety of peginterferon alfa-2b (1.5 mu g/kg/wk) plus weight-based ribavirin (800-1400 mg/day) in 2333 chronic HCV-infected patients with significant fibrosis/cirrhosis whose previous interferon alfa/ribavirin therapy failed. Patients with undetectable HCV-RNA at treatment week (TW) 12 received 48 weeks of therapy; patients with detectable HCV-RNA at TW12 could enter maintenance studies at TW18; 188 patients with low/detectable HCV-RNA at TW12 continued therapy at the investigator`s request. Results: Overall, 22% of the patients attained SVR (56% with undetectable HCV-RNA and 12% with low/detectable HCV-RNA at TW12). SVR was better in relapsers (38%) than nonresponders (14%), regardless of previous treatment, and in patients previously treated with interferon-alfa/ribavirin (25%) than peginterferon alfa-ribavirin (17%). Predictors of response in patients with undetectable HCV-RNA at TW12 were genotype (2/3 vs 1, respectively; odds ratio [OR] 2.4; P < .0001), fibrosis score (F2 vs F4; OR, 2.2; F3 vs F4; OR, 1.7; P < .0001), and baseline viral load (<= 600,000 vs >600,000 IU/mL; OR, 1.4; P = .0223). These factors plus previous treatment and response were overall predictors of SVR. Safety was similar among fibrosis groups. Conclusions: Peginterferon alfa-2b plus weight-based ribavirin is effective and safe in patients who failed interferon alfa/ribavirin therapy. Genotype, baseline viral load, and fibrosis stage were predictors of response.

Association between a 15q25 gene variant, smoking quantity and tobacco-related cancers among 17 000 individuals

Methods We performed a detailed analysis of one 15q single nucleotide polymorphism (SNP) (rs16969968) with smoking behaviour and cancer risk in a total of 17 300 subjects from five LC studies and four upper aerodigestive tract (UADT) cancer studies. Results Subjects with one minor allele smoked on average 0.3 cigarettes per day (CPD) more, whereas subjects with the homozygous minor AA genotype smoked on average 1.2 CPD more than subjects with a GG genotype (P < 0.001). The variant was associated with heavy smoking (> 20 CPD) [odds ratio (OR) = 1.13, 95% confidence interval (CI) 0.96-1.34, P = 0.13 for heterozygotes and 1.81, 95% CI 1.39-2.35 for homozygotes, P < 0.0001]. The strong association between the variant and LC risk (OR = 1.30, 95% CI 1.23-1.38, P = 1 x 10(-18)), was virtually unchanged after adjusting for this smoking association (smoking adjusted OR = 1.27, 95% CI 1.19-1.35, P = 5 x 10(-13)). Furthermore, we found an association between the variant allele and an earlier age of LC onset (P = 0.02). The association was also noted in UADT cancers (OR = 1.08, 95% CI 1.01-1.15, P = 0.02). Genome wide association (GWA) analysis of over 300 000 SNPs on 11 219 subjects did not identify any additional variants related to smoking behaviour. Conclusions This study confirms the strong association between 15q gene variants and LC and shows an independent association with smoking quantity, as well as an association with UADT cancers.

Mucosal and systemic anti-GAG immunity induced by neonatal immunization with HIV LAMP/gag DNA vaccine in mice

Vaccines capable of inducing mucosal immunity in early postnatal life until adulthood, protecting early sexual initiation, should be considered as strategies to vaccination against HIV. The HIV-1 GAG protein as a chimera with the lysosome-associated membrane protein (LAMP/gag), encoded by a DNA vaccine, is targeted to the endosomal/lysosomal compartment that contains class II MHC molecules and has been shown to be immunogenic in adult mice. Assuming that one such strategy could help to overcome the immunological immaturity in the early postnatal period, we have evaluated the systemic and mucosal immunogenicity of LAMP/gag immunization in neonatal mice. Intranasal immunization with LAMP/gag vaccine induced higher levels of sIgA and IgG anti-GAG antibodies in intestinal washes than did the gag vaccine. The combination of ID injections and the IN protocol with the chimeric vaccine promoted the increase of Ab levels in sera. Both vaccines induced splenic IFN-gamma- secreting cells against GAG peptide pools, as well as in vivo cytotoxic T lymphocyte (CTL) function, and increased the percentage of CD8+ T cells to the immunodominant class I peptide in gut and spleen. However, only the chimeric vaccine was able to enhance Th1/Th2 cytokine secretion in response to class II GAG peptide and to enhance IL-4-secreting cells against GAG peptides and p24 protein stimuli. Long-lasting humoral and cellular responses were detected until adult age, following neonatal immunization with the chimeric vaccine. The LAMP/gag vaccination was able to induce potent GAG-specific T and B cell immune responses in early life which are essential to elicit sustained and long-lasting mucosal and systemic humoral response. (C) 2010 Elsevier GmbH. All rights reserved.

Global Analysis of Protein Palmitoylation in African Trypanosomes

Many eukaryotic proteins are posttranslationally modified by the esterification of cysteine thiols to long-chain fatty acids. This modification, protein palmitoylation, is catalyzed by a large family of palmitoyl acyltransferases that share an Asp-His-His-Cys Cys-rich domain but differ in their subcellular localizations and substrate specificities. In Trypanosoma brucei, the flagellated protozoan parasite that causes African sleeping sickness, protein palmitoylation has been observed for a few proteins, but the extent and consequences of this modification are largely unknown. We undertook the present study to investigate T. brucei protein palmitoylation at both the enzyme and substrate levels. Treatment of parasites with an inhibitor of total protein palmitoylation caused potent growth inhibition, yet there was no effect on growth by the separate, selective inhibition of each of the 12 individual T. brucei palmitoyl acyltransferases. This suggested either that T. brucei evolved functional redundancy for the palmitoylation of essential palmitoyl proteins or that palmitoylation of some proteins is catalyzed by a noncanonical transferase. To identify the palmitoylated proteins in T. brucei, we performed acyl biotin exchange chemistry on parasite lysates, followed by streptavidin chromatography, two-dimensional liquid chromatography-tandem mass spectrometry protein identification, and QSpec statistical analysis. A total of 124 palmitoylated proteins were identified, with an estimated false discovery rate of 1.0%. This palmitoyl proteome includes all of the known palmitoyl proteins in procyclic-stage T. brucei as well as several proteins whose homologues are palmitoylated in other organisms. Their sequences demonstrate the variety of substrate motifs that support palmitoylation, and their identities illustrate the range of cellular processes affected by palmitoylation in these important pathogens.

Molecular cloning and characterization of a ripening-induced polygalacturonase related to papaya fruit softening

Pulp softening is one of the most remarkable changes during ripening of papaya (Carica papaya) fruit and it is a major cause for post-harvest losses. Although cell wall catabolism has a major influence on papaya fruit, quality information on the gene products involved in this process is limited. A full-length polygalacturonase cDNA (cpPG) was isolated from papaya pulp and used to study gene expression and enzyme activity during normal and ethylene-induced ripening and after exposure of the fruit to 1-MCP. Northern-blot analysis demonstrated that cpPG transcription was strongly induced during ripening and was highly ethylene-dependent. The accumulation of cpPG transcript was paralleled by enzyme activity, and inversely correlated to the pulp firmness. Preliminary in silica analysis of the cpPG genomic sequence revealed the occurrence of putative regulatory motifs in the promoter region that may help to explain the effects of plant hormones and non-abiotic stresses on papaya fruit firmness. This newly isolated cpPG is an important candidate for functional characterization and manipulation to control the process of pulp softening during papaya ripening. (C) 2009 Elsevier Masson SAS. All rights reserved.