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.

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.

Thermal Characteristics of Polybia scutellaris Nests (Hymenoptera: Vespidae) Using Computational Fluid Dynamics: A Possible Adaptation to Tropical Climates

Polybia scutellaris constructs huge nests characterized by numerous spinal projections on the surface. We investigated the thermal characteristics of P scutellaris nests in order to determine whether their nest temperature is homeothermically maintained and whether the spines play a role in the thermoregulation of the nests. In order to examine these hypotheses, we measured the nest temperature in a active nest and in an abandoned nest. The temperature in the active nest was almost stable at 27 degrees C, whereas that of the abandoned nest varied with changes in the ambient temperature, suggesting that nest temperature was maintained by the thermogenesis of colony individuals. In order to predict the thermal properties of the spines, a numerical simulation was employed. To construct a 3D-model of a P scutellaris nest, the nest architecture was simplified into an outer envelope and the surface spines, for both of which the initial temperature was set at 27 degrees C. The physical properties of the simulated nest were regarded to be those of wood since the nest of this species is constructed from plant materials. When the model was exposed to cool air (12 degrees C), the temperature was lower in the models with more spines. On the other hand, when the nest was heated (42 degrees C), the temperature increase was smaller in models with more spines. It is suggested that the spines act as a heat radiator, not as an insulator, against the changes in ambient temperature.

Generation of Naphthoquinone Radical Anions by Electrospray Ionization: Solution, Gas-Phase, and Computational Chemistry Studies

Radical anions are present in several chemical processes, and understanding the reactivity of these species may be described by their thermodynamic properties. Over the last years, the formation of radical ions in the gas phase has been an important issue concerning electrospray ionization mass spectrometry studies. In this work, we report on the generation of radical anions of quinonoid compounds (Q) by electrospray ionization mass spectrometry. The balance between radical anion formation and the deprotonated molecule is also analyzed by influence of the experimental parameters (gas-phase acidity, electron affinity, and reduction potential) and solvent system employed. The gas-phase parameters for formation of radical species and deprotonated species were achieved on the basis of computational thermochemistry. The solution effects on the formation of radical anion (Q(center dot-)) and dianion (Q(2-)) were evaluated on the basis of cyclic voltammetry analysis and the reduction potentials compared with calculated electron affinities. The occurrence of unexpected ions [Q + 15](-) was described as being a reaction between the solvent system and the radical anion, Q(center dot-).The gas-phase chemistry of the electrosprayed radical anions was obtained by collisional-induced dissociation and compared to the relative energy calculations. These results are important for understanding the formation and reactivity of radical anions and to establish their correlation with the reducing properties by electrospray ionization analyses.

A Computational Study of Tetrafluoro-[2.2]Cyclophanes

A computational study of the isomers of tetrafluorinated [2.2]cyclophanes persubstituted in one ring, namely F-4-[2.2]paracyclophane (4), F-4-anti-[2.2]metacyclophane (5a), F-4-syn-[2.2]metacyclophane (5b), and F-4-[2.2]metaparacyclophane (6a and 6b), was carried out. The effects of fluorination on the geometries, relative energies, local and global aromaticity, and strain energies of the bridges and rings were investigated. An analysis of the electron density by B3PW91/6-31+G(d,p), B3LYP/6-31+G(d,p), and MP2/6-31+G(d,p) was carried out using the natural bond orbitals (NBO), natural steric analysis (NSA), and atoms in molecules (AIM) methods. The analysis of frontier molecular orbitals (MOs) was also employed. The results indicated that the molecular structure of [2.2]paracyclophane is the most affected by the fluorination. Isodesmic reactions showed that the fluorinated rings are more strained than the nonfluorinated ones. The NICS, HOMA, and PDI criteria evidenced that the fluorination affects the aromaticity of both the fluorinated and the nonfluorinated rings. The NBO and NSA analyses gave an indication that the fluorination increases not only the number of through-space interactions but also their magnitude. The AIM analysis suggested that the through-space interactions are restricted to the F-4-[2.2]metacyclophanes. In addition, the atomic properties, computed over the atomic basins, shave evidence that not only the substitution, but also the position of the bridges could affect the atomic charges. the first atomic moments, and the atomic volumes.

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.

Solution structure of the sodium channel antagonist conotoxin GS: A new molecular caliper for probing sodium channel geometry

Background: The venoms of Conus snails contain small, disulfide-rich inhibitors of voltage-dependent sodium channels. Conotoxin GS is a 34-residue polypeptide isolated from Conus geographus that interacts with the extracellular entrance of skeletal muscle sodium channels to prevent sodium ion conduction. Although conotoxin GS binds competitively with mu conotoxin GIIIA to the sodium channel surface, the two toxin types have little sequence identity with one another, and conotoxin GS has a four-loop structural framework rather than the characteristic three-loop mu-conotoxin framework. The structural study of conotoxin GS will form the basis for establishing a structure-activity relationship and understanding its interaction with the pore region of sodium channels. Results: The three-dimensional structure of conotoxin GS was determined using two-dimensional NMR spectroscopy. The protein exhibits a compact fold incorporating a beta hairpin and several turns. An unusual feature of conotoxin GS is the exceptionally high proportion (100%) of cis-imide bond geometry for the three proline or hydroxyproline residues. The structure of conotoxin GS bears little resemblance to the three-loop mu conotoxins, consistent with the low sequence identity between the two toxin types and their different structural framework. However, the tertiary structure and cystine-knot motif formed by the three disulfide bonds is similar to that present in several other polypeptide ion channel inhibitors. Conclusions: This is the first three-dimensional structure of a 'four-loop' sodium channel inhibitor, and it represents a valuable new structural probe for the pore region of voltage-dependent sodium channels. The distribution of amino acid sidechains in the structure creates several polar and charged patches, and comparison with the mu conotoxins provides a basis for determining the binding surface of the conotoxin GS polypeptide.

Implications of prevalent and incident diabetes mellitus on left ventricular geometry and function in the ageing heart: The MONICA/KORA Augsburg cohort study

Background and Aim: It is unclear to what extent diabetes modulates the ageing-related adaptations of cardiac geometry and function. Methods and Results: We examined 1005 adults, aged 25-74 years, from a population-based survey at baseline in 1994/5 and at follow-up in 2004/5. We compared persistently non-diabetic individuals (ND; no diabetes at baseline and at follow-up, n = 833) with incident (ID; non-diabetic at baseline and diabetic at follow-up, n = 36) and with prevalent diabetics (PD; diabetes at baseline and follow-up examination, n = 21). Left ventricular (LV) geometry and function were evaluated by echocardiography. Statistical analyses were performed with multivariate linear regression models. Over ten years the PD group displayed a significantly stronger relative increase of LV mass (+9.34% vs. +23.7%) that was mediated by a more pronounced increase of LV end-diastolic diameter (+0% vs. +6.95%) compared to the ND group. In parallel, LA diameter increased (+4.50% vs. +12.7%), whereas ejection fraction decreased (+3.02% vs. -4.92%) more significantly in the PD group. Moreover, at the follow-up examination the PD and ID groups showed a significantly worse diastolic function, indicated by a higher E/EM ratio compared with the ND group (11.6 and 11.8 vs. 9.79, respectively). Conclusions: Long-standing diabetes was associated with an acceleration of age-related changes of left ventricular geometry accumulating in an eccentric remodelling of the left ventricle. Likewise, echocardiographic measures of systolic and diastolic ventricular function deteriorated more rapidly in individuals with diabetes. (C) 2009 Elsevier B.V. All rights reserved.

Evaluation of computational methods for the reconstruction of HLA haplotypes

Human leukocyte antigen (HLA) haplotypes are frequently evaluated for population history inferences and association studies. However, the available typing techniques for the main HLA loci usually do not allow the determination of the allele phase and the constitution of a haplotype, which may be obtained by a very time-consuming and expensive family-based segregation study. Without the family-based study, computational inference by probabilistic models is necessary to obtain haplotypes. Several authors have used the expectation-maximization (EM) algorithm to determine HLA haplotypes, but high levels of erroneous inferences are expected because of the genetic distance among the main HLA loci and the presence of several recombination hotspots. In order to evaluate the efficiency of computational inference methods, 763 unrelated individuals stratified into three different datasets had their haplotypes manually defined in a family-based study of HLA-A, -B, -DRB1 and -DQB1 segregation, and these haplotypes were compared with the data obtained by the following three methods: the Expectation-Maximization (EM) and Excoffier-Laval-Balding (ELB) algorithms using the arlequin 3.11 software, and the PHASE method. When comparing the methods, we observed that all algorithms showed a poor performance for haplotype reconstruction with distant loci, estimating incorrect haplotypes for 38%-57% of the samples considering all algorithms and datasets. We suggest that computational haplotype inferences involving low-resolution HLA-A, HLA-B, HLA-DRB1 and HLA-DQB1 haplotypes should be considered with caution.

Theory of a mode-locked atom laser with toroidal geometry

We consider a possible technique for mode locking an atom laser, based on the generation of a dark soliton in a ring-shaped Bose-Einstein condensate, with repulsive atomic interactions. The soliton is a kink, with angular momentum per particle equal to (h) over bar /2. It emerges naturally when the condensate is stirred at the soliton velocity and cleansed with a periodic out coupler. The result is a replicating coherent field inside the atom laser, stabilized by topology. We give a numerical demonstration of the generation and stabilization of the soliton.

The effect of channel geometry and wall boundary conditions on the formation of extrusion surface instabilities for LLDPE

It is believed that surface instabilities can occur during the extrusion of linear low density polyethylene due to high extensional stresses at the exit of the die. Local crack development can occur at a critical stress level when melt rupture is reached. This high extensional stress results from the rearrangement of the flow at the boundary transition between the wall exit and the free surface. The stress is highest at the extrudate surface and decreases into the bulk of the material. The location of the region where the critical level is reached can determine the amplitude of the extrudate surface distortion, This paper studies the effect of wall slip on the numerically simulated extensional stress level at the die exit and correlates this to the experimentally determined amplitude of the surface instability. The effect of die exit radius and die wall roughness on extrusion surface instabilities is also correlated to the exit stress level in the same way. Whereas full slip may completely suppress the surface instability, a reduction in the exit stress level and instability amplitude is also shown for a rounded die exit and a slight increase in instability is shown to result from a rough die wall. A surface instability map demonstrates how the shear rate for onset of extrusion surface instabilities can be predicted on the basis of melt strength measurements and simulated stress peaks at the exit of the die. (C) 2001 Elsevier Science B.V. All rights reserved.

Linear ketenimines. Variable structures of C,C-dicyanoketenimines and C,C-bis-sulfonylketenimines

C,C-Dicyanoketenimines 10a-c were generated by flash vacuum thermolysis of ketene NS-acetals 9a-c or by thermal or photochemical decomposition of alpha-azido-,beta-cyanocinnamonitrile 11. In the latter reaction, 3,3-dicyano-2-phenyl-1-azirine 12 is also formed. IR spectroscopy of the keteniminines isolated in Ar matrixes or as neat films, NMR spectroscopy of 10c, and theoretical calculations (B3LYP/6-31G*) demonstrate that these ketenimines have variable geometry, being essentially linear along the CCN-R framework in polar media (neat films and solution), but in the gas phase or Ar matrix they are bent, as is usual for ketenimines. Experiments and calculations agree that a single CN substituent as in 13 is not enough to enforce linearity, and sulfonyl groups are less effective that cyano groups in causing linearity. C,C-Bis(methylsulfonyl)ketenimines 4-5 and a C-cyano-C-(methylsulfonyl)ketenimine 15 are not linear. The compound p-O2NC6H4N=C= C(COOMe)2 previously reported in the literature is probably somewhat linearized along the CCNR moiety. A computational survey (B3LYP/6-31G*) of the inversion barrier at nitrogen indicates that electronegative C-substituents dramatically lower the barrier; this is also true of N-acyl substituents. Increasing polarity causes lower barriers. Although N-alkylbis(methylsulfonyl)ketenimines are not calculated to be linear, the barriers are so low that crystal lattice forces can induce planarity in N-methylbis(methylsulfonyl)ketenimine 3.

Influence of wall stress and left ventricular geometry on the accuracy of Dobutamine stress Echocardiography

OBJECTIVES The goal of this study was to determine whether wall stress at rest and during stress could explain the influence of left ventricular (LV) morphology on the accuracy of dobutamine stress echocardiography (DSE). BACKGROUND The sensitivity of DSE appears to be reduced in patients with concentric remodeling, but the cause of this finding is unclear. METHODS We studied 161 patients without resting wall motion abnormalities who underwent DSE and coronary angiography. Patients were classified into four groups according to relative wan thickness (normal