Search

E1 and E2 (E1E2), the fusion proteins of Hepatitis C Virus (HCV), are unlike that of any other virus yet described, and the detailed molecular mechanisms of HCV entry/fusion remain unknown. Hypervariable region-1 (HVR-1) of E2 is a putative intrinsically disordered protein tail. Here, we demonstrate that HVR-1 has an autoinhibitory function that suppresses the activity of E1E2 on free virions; this is dependent on its conformational entropy. Thus, HVR-1 is akin to a safety catch that prevents premature triggering of E1E2 activity. Crucially, this mechanism is turned off by host receptor interactions at the cell surface to allow entry. Mutations that reduce conformational entropy in HVR-1, or genetic deletion of HVR-1, turn off the safety catch to generate hyper-reactive HCV that exhibits enhanced virus entry but is thermally unstable and acutely sensitive to neutralising antibodies. Therefore, the HVR-1 safety catch controls the efficiency of virus entry and maintains resistance to neutralising antibodies. This discovery provides an explanation for the ability of HCV to persist in the face of continual immune assault and represents a novel regulatory mechanism that is likely to be found in other viral fusion machinery.

We present flexiMAP (flexible Modeling of Alternative PolyAdenylation), a new beta-regression-based method implemented in R, for discovering differential alternative polyadenylation events in standard RNA-seq data. We show, using both simulated and real data, that flexiMAP exhibits a good balance between specificity and sensitivity and compares favourably to existing methods, especially at low fold changes. In addition, the tests on simulated data reveal some hitherto unrecognised caveats of existing methods. Importantly, flexiMAP allows modeling of multiple known covariates that often confound the results of RNA-seq data analysis. This repository contains scripts and all data required to reproduce data in the manuscript. _simulationData zip files contain extensive set of simulations which should serve as a useful resource for the development of similar methods in the future. Additionally, it contains 3’ sequencing data (PolyA-seq) and RNA-seq data from the Human Brain Reference and the Universal Human Reference MAQC samples (Bullard et al., 2010) that was used in the manuscript., {"references":["Krzysztof J. Szkop, David S. Moss and Irene Nobeli flexiMAP: A regression-based method for discovering differential alternative polyadenylation events in standard RNA-seq data"]}

Over 500 missense F8 mutations have been reported to cause non-severe haemophilia A. Some F8 genotypes appear to confer a higher risk of inhibitor formation than others and individuals with the same F8 genotype may have differing risks of inhibitor formation. We present an in silico strategy demonstrating the heterogeneity of factor VIII (FVIII)-derived antigen presentation whilst identifying patterns of human leucocyte antigen (HLA) peptide binding that might predict future inhibitor risk. A well-validated computational tool, NetMHCII, enabled large-scale comparison of predicted antigen presentation between endogenous, mutated FVIII-derived peptides and wild-type, therapeutic FVIII-derived peptides spanning all F8 missense mutation positions reported to The Haemophilia A Mutation, Structure and Resource Site (HADB). We identify 40 F8 genotypes to be 'low risk' at a 50% inhibitory concentration (IC50 )-binding threshold of 300 nmol/l (P = 0·00005), defined as absence of novel peptide-major histocompatibility complex (MHC) surfaces for all 14 common HLA-DR alleles assessed. Analysing each of the possible 7280 F8 genotype/HLA-DR permutations individually at an IC50 threshold of 300 nmol/l, 65% are predicted to not generate a novel peptide-MHC surface that would be necessary to engage T cell help for subsequent anti-FVIII antibody generation. This study demonstrates the future importance of interpreting F8 genotype in the context of an individual's HLA profile to personalize inhibitor risk prediction.

International audience; CPE (Clostridium perfringens enterotoxin) is the major virulence determinant for C. perfringens type-A food poisoning, the second most common bacterial food-borne illness in the UK and USA. After binding to its receptors, which include particular human claudins, the toxin forms pores in the cell membrane. The mature pore apparently contains a hexamer of CPE, claudin and, possibly, occludin. The combination of high binding specificity with cytotoxicity has resulted in CPE being investigated, with some success, as a targeted cytotoxic agent for oncotherapy. In this paper, we present the X-ray crystallographic structure of CPE in complex with a peptide derived from extracellular loop 2 of a modified, CPE-binding Claudin-2, together with high-resolution native and pore-formation mutant structures. Our structure provides the first atomic-resolution data on any part of a claudin molecule and reveals that claudin's CPE-binding fingerprint (NPLVP) is in a tight turn conformation and binds, as expected, in CPE's C-terminal claudin-binding groove. The leucine and valine residues insert into the binding groove while the first residue, asparagine, tethers the peptide via an interaction with CPE's aspartate 225 and the two prolines are required to maintain the tight turn conformation. Understanding the structural basis of the contribution these residues make to binding will aid in engineering CPE to target tumor cells.

Highlights • Etx mutant Y30A-Y196A showed markedly reduced cytotoxicity towards MDCK.2 cells. • Y30A-Y196A is inactive in mice after intraperitoneal administration. • Y30A-Y196A is able to induce a specific antibody response in rabbits. • Y30A-Y196A polyclonal antibody is able to induce protective immunity in vitro. • Y30A-Y196A could form the basis of a recombinant vaccine against enterotoxemia., Epsilon toxin (Etx) is a β-pore-forming toxin produced by Clostridium perfringens toxinotypes B and D and plays a key role in the pathogenesis of enterotoxemia, a severe, often fatal disease of ruminants that causes significant economic losses to the farming industry worldwide. This study aimed to determine the potential of a site-directed mutant of Etx (Y30A-Y196A) to be exploited as a recombinant vaccine against enterotoxemia. Replacement of Y30 and Y196 with alanine generated a stable variant of Etx with significantly reduced cell binding and cytotoxic activities in MDCK.2 cells relative to wild type toxin (>430-fold increase in CT50) and Y30A-Y196A was inactive in mice after intraperitoneal administration of trypsin activated toxin at 1000× the expected LD50 dose of trypsin activated wild type toxin. Moreover, polyclonal antibody raised in rabbits against Y30A-Y196A provided protection against wild type toxin in an in vitro neutralisation assay. These data suggest that Y30A-Y196A mutant could form the basis of an improved recombinant vaccine against enterotoxemia.

Recently, a number of broad-spectrum human antibodies binding to the stalk region of influenza A haemagglutinin (HA) have been isolated. As this region tends to develop substitutions at a slower rate than other regions of HA, a vaccine eliciting such antibodies could have a longer effective life. But this begs a question: is the stalk resistant to change even in the face of evolutionary pressure? In this paper, we analysed the known epitopes in the H3 stalk and, utilizing a collection of 3440 sequences, present a novel approach for detecting putative B-cell epitopes in regions such as this, in which mutations occur infrequently. We concluded that there have been periods of activity in the stalk that are consistent with the evolution of antigenic escape. This work casts light on the presence of stalk-binding antibodies in the population as a whole and, through the analysis of antigenically active regions in the stalk, may contribute to the identification of epitopes that are refractive to change and hence useful for vaccine development.

Clostridium perfringens epsilon toxin (Etx) is a pore-forming toxin responsible for a severe and rapidly fatal enterotoxemia of ruminants. The toxin is classified as a category B bioterrorism agent by the U.S. Government Centres for Disease Control and Prevention (CDC), making work with recombinant toxin difficult. To reduce the hazard posed by work with recombinant Etx, we have used a variant of Etx that contains a H149A mutation (Etx-H149A), previously reported to have reduced, but not abolished, toxicity. The three-dimensional structure of H149A prototoxin shows that the H149A mutation in domain III does not affect organisation of the putative receptor binding loops in domain I of the toxin. Surface exposed tyrosine residues in domain I of Etx-H149A (Y16, Y20, Y29, Y30, Y36 and Y196) were mutated to alanine and mutants Y30A and Y196A showed significantly reduced binding to MDCK.2 cells relative to Etx-H149A that correlated with their reduced cytotoxic activity. Thus, our study confirms the role of surface exposed tyrosine residues in domain I of Etx in binding to MDCK cells and the suitability of Etx-H149A for further receptor binding studies. In contrast, binding of all of the tyrosine mutants to ACHN cells was similar to that of Etx-H149A, suggesting that Etx can recognise different cell surface receptors. In support of this, the crystal structure of Etx-H149A identified a glycan (β-octyl-glucoside) binding site in domain III of Etx-H149A, which may be a second receptor binding site. These findings have important implications for developing strategies designed to neutralise toxin activity.

The ultimate aim of the EU-funded ImmunoGrid project is to develop a natural-scale model of the human immune system—that is, one that reflects both the diversity and the relative proportions of the molecules and cells that comprise it—together with the grid infrastructure necessary to apply this model to specific applications in the field of immunology. These objectives present the ImmunoGrid Consortium with formidable challenges in terms of complexity of the immune system, our partial understanding about how the immune system works, the lack of reliable data and the scale of computational resources required. In this paper, we explain the key challenges and the approaches adopted to overcome them. We also consider wider implications for the present ambitious plans to develop natural-scale, integrated models of the human body that can make contributions to personalized health care, such as the European Virtual Physiological Human initiative. Finally, we ask a key question: How long will it take us to resolve these challenges and when can we expect to have fully functional models that will deliver health-care benefits in the form of personalized care solutions and improved disease prevention?

In this paper we undertake an analysis of the antigenicity of influenza A virus hemagglutinin. We developed a novel computational approach to the identification of antigenically active regions and showed that the amino acid substitutions between successive predominant seasonal strains form clusters that are consistent, in terms of both their location and their size, with the properties of B-cell epitopes in general and with those epitopes that have been identified experimentally in influenza A virus hemagglutinin to date. Such an interpretation provides a biologically plausible framework for an understanding of the location of antigenically important substitutions that is more specific than the canonical “antigenic site” model and provides an effective basis for deriving models that predict antigenic escape in the H3N2 subtype. Our results support recent indications that antibodies binding to the “stalk” region of hemagglutinin are found in the human population and exert evolutionary pressure on the virus. Our computational approach provides a possible method for identifying antigenic escape through evolution in this region, which in some cases will not be identified by the hemagglutinin inhibition assay.

We are used to thinking of proteins as beneficial, so it is surprising to realize that the most toxic substances known to man are also protein molecules. These are the bacterial exotoxins, proteins secreted by pathogenic bacteria. Toxicity is measured by the median lethal dose (LD50). An LD50 value is defined as the mass of toxin per kg of body weight required to wipe out half of an animal population. Whereas classic poisons, such as potassium cyanide or arsenic trioxide, have LD50 values in the range 5–15 mg/ kg, the causative agent of botulism, botulinum toxin, has an LD50 in the range 1–3 ng/kg, a million times more toxic!

T cell activation is the final step in a complex pathway through which pathogen-derived peptide fragments can elicit an immune response. For it to occur, peptides must form stable complexes with Major Histocompatibility Complex (MHC) molecules and be presented on the cell surface. Computational predictors of MHC binding are often used within in silico vaccine design pathways. We have previously shown that, paradoxically, most bacterial proteins known experimentally to elicit an immune response in disease models are depleted in peptides predicted to bind to human MHC alleles. The results presented here, derived using software proven through benchmarking to be the most accurate currently available, show that vaccine antigens contain fewer predicted MHC-binding peptides than control bacterial proteins from almost all subcellular locations with the exception of cell wall and some cytoplasmic proteins. This effect is too large to be explained from the undoubted lack of precision of the software or from the amino acid composition of the antigens. Instead, we propose that pathogens have evolved under the influence of the host immune system so that surface proteins are depleted in potential MHC-binding peptides, and suggest that identification of a protein likely to contain a single immuno-dominant epitope is likely to be a productive strategy for vaccine design.

We have developed a computational Grid that enables us to exploit through a single interface a range of local, national and international resources. It insulates the user as far as possible from issues concerning administrative boundaries, passwords and different operating system features. This work has been undertaken as part of the European Union ImmunoGrid project whose aim is to develop simulations of the immune system at the molecular, cellular and organ levels. The ImmunoGrid consortium has members with computational resources on both sides of the Atlantic. By making extensive use of existing Grid middleware, our Grid has enabled us to exploit consortium and publicly available computers in a unified way, notwithstanding the diverse local software and administrative environments. We took 40 000 polypeptide sequences from 4000 avian and mammalian influenza strains and used a neural network for class I T-cell epitope prediction tools for 120 class I alleles and haplotypes to generate over 14 million high-quality protein–peptide binding predictions that we are mapping onto the three-dimensional structures of the proteins. By contrast, the Grid is also being used for developing new methods for class T-cell epitope predictions, where we have running batches of 120 molecular dynamics free-energy calculations.

For many infectious diseases, protective immunity can be elicited by vaccination with pathogen-derived proteins. Peptides derived from these proteins are bound to major histocompatibility complex (MHC) molecules and presented to T-cell receptors to stimulate an immune response. We show here that, paradoxically, bacterial proteins known experimentally to elicit a protective immune response are relatively depleted in peptides predicted to bind to human MHC alleles. We propose three nonconflicting reasons for this: the lack of precision of current predictive software, the low incidence of hydrophobic residues in vaccine antigens or evolutionary pressure exerted on bacteria by the immune system. We suggest that there is little value in predicting candidate vaccines based on high MHC-binding epitope density.

Epitopes mediated by T cells lie at the heart of the adaptive immune response and form the essential nucleus of anti-tumour peptide or epitope-based vaccines. Antigenic T cell epitopes are mediated by major histocompatibility complex (MHC) molecules, which present them to T cell receptors. Calculating the affinity between a given MHC molecule and an antigenic peptide using experimental approaches is both difficult and time consuming, thus various computational methods have been developed for this purpose. A server has been developed to allow a structural approach to the problem by generating specific MHC:peptide complex structures and providing configuration files to run molecular modelling simulations upon them. A system has been produced which allows the automated construction of MHC:peptide structure files and the corresponding configuration files required to execute a molecular dynamics simulation using NAMD. The system has been made available through a web-based front end and stand-alone scripts. Previous attempts at structural prediction of MHC:peptide affinity have been limited due to the paucity of structures and the computational expense in running large scale molecular dynamics simulations. The MHCsim server (http://igrid-ext.cryst.bbk.ac.uk/MHCsim) allows the user to rapidly generate any desired MHC:peptide complex and will facilitate molecular modelling simulation of MHC complexes on an unprecedented scale.

Human leukocyte antigen (HLA)-DM is a critical participant in antigen presentation that catalyzes the dissociation of the Class II-associated Invariant chain-derived Peptide (CLIP) from the major histocompatibility complex (MHC) Class II molecules. There is competition amongst peptides for access to an MHC Class II groove and it has been hypothesised that DM functions as a 'peptide editor' that catalyzes the replacement of one peptide for another within the groove. It is established that the DM catalyst interacts directly with the MHC Class II but the precise location of the interface is unknown. Here, we combine previously described mutational data with molecular docking and energy minimisation simulations to identify a putative interaction site of >4000A2 which agrees with known point mutational data for both the DR and DM molecule. The docked structure is validated by comparison with experimental data and previously determined properties of protein-protein interfaces. A possible dissociation mechanism is suggested by the presence of an acidic cluster near the N terminus of the bound peptide.

In this review, we focus upon X-ray studies of diffuse scattering from molecular crystals and proteins, but reference is also made to other crystals of relevance to the discussion.

Antigenic peptide is presented to a T-cell receptor through the formation of a stable complex with a Major Histocompatibility Complex (MHC) molecule. Various predictive algorithms have been developed to estimate a peptide’s capacity to form a stable complex with a given MHC Class II allele, a technique integral to the strategy of vaccine design. These have previously incorporated such computational techniques as quantitative matrices and neural networks. We have developed a novel predictive technique that uses molecular modeling of predetermined crystal structures to estimate the stability of an MHC Class II peptide complex. This is the 1st structure-based technique, as previous methods have been based on binding data. ROC curves are used to quantify the accuracy of the molecular modeling technique. The novel predictive technique is found to be comparable with the best predictive software currently available.

The extent of the role of N-linked glycans (N-glycans) in shielding influenza A hemagglutinin (HA) against host antibodies has proved controversial, with different authors making widely different assumptions. One common assumption is that N-glycans physically shield surface residues that are near to glycosylation sites, thereby preventing antibodies from binding to them. However, it is unclear, from existing experimental evidence, whether antibodies that bind close to N-glycans are a rare or commonplace feature of human herd immune responses to influenza AHA. The aim of this paper is to present a computational analysis of mutations in the vicinity of N-glycans that will facilitate a better understanding of their protective role. We identify, from an analysis of over 6000 influenza A H3N2 sequences, a set of residues adjacent to N-glycosylation sites that are highly likely to be involved in antigenic escape from host antibodies. Fifteen of these residues occur within 10 A of an N-glycosylation site. Hence, we conclude that it is relatively common for antibodies to bind in close proximity to N-glycans on the surface ofHA, with any shielding effect largely attributable to the inability of host antibodies to bind across an N-glycan attachment site, rather than to the physical masking of neighboring residues.

Differences in the biological properties of the Clostridium perfringens phospholipase C (alpha-toxin) and the C. bifermentans phospholipase C (Cbp) have been attributed to differences in their carboxy-terminal domains. Three residues in the carboxy-terminal domain of alpha-toxin, which have been proposed to play a role in membrane recognition (D269, Y331 and F334), are not conserved in Cbp (Y, L and I respectively). We have characterised D269Y, Y331L and F334I variant forms of alpha-toxin. Variant D269Y had reduced phospholipase C activity towards aggregated egg yolk phospholipid but increased haemolytic and cytotoxic activity. Variants Y331L and F334I showed a reduction in phospholipase C, haemolytic and cytotoxic activities indicating that these substitutions contribute to the reduced haemolytic and cytotoxic activity of Cbp.

A panel of random mutants within the DNA encoding the carboxy-terminal domain of Clostridium perfringens alpha-toxin was constructed. Three mutants were identified which encoded alpha-toxin variants (Lys330Glu, Asp305Gly, and Asp293Ser) with reduced hemolytic activity. These variants also had diminished phospholipase C activity toward aggregated egg yolk phospholipid and reduced cytotoxic and myotoxic activities. Asp305Gly showed a significantly increased enzymatic activity toward the monodisperse substrate rhoNPPC, whereas Asp293Ser displayed a reduced activity toward this phospholipid analogue. In addition, Asp293Ser showed an increased dependence on calcium for enzymatic activity toward aggregated phospholipid and appeared calcium-depleted in PAGE band-shift assays. In contrast, neither Lys330Glu nor Asp305Gly showed altered dependence on calcium for enzymatic activity toward aggregated phospholipid. Asp305 is located in the interface between the amino- and carboxy-terminal domains, whereas Asp293 and Lys330 are surface exposed residues which may play a role in the recognition of membrane phospholipids.

The last five years have seen a large increase in the use of cross validation in the refinement of macromolecular structures using X-ray data. In this technique a test set of reflections is set aside from the working set and the progress of the refinement is monitored by the calculation of a free R factor which is based only on the excluded reflections. This paper gives estimates for the ratio of the free R factor to the R factor calculated from the working set for both unrestrained and restrained refinement. It is assumed that both the X-­ray and restraint observations have been weighted correctly and that there is no correlation of errors between the test and working sets. It is also shown that the least-squares weights that minimize the variances of the refined parameters, also approximately minimize the free R factor. The estimated free R-factor ratios are compared with those reported for structures in the Protein Data Bank.

Faster workstations with larger memories are making error estimation from full-matrix least-squares refinement a more practicable technique in protein crystallography. Using minimum variance weighting, estimated standard deviations of atomic positions have been calculated for two eye lens proteins from the inverse of a least-squares normal matrix which was full with respect to the coordinate parameters. gammaB-crystallin, refined at 1.49 A yielded average errors in atomic positions which ranged from 0.05 A for main-chain atoms to 0.27 A for unrestrained water molecules. The second structure used in this work was that of betaB2-crystallin refined at 2.1 A resolution where the corresponding average errors were 0.08 and 0.35 A, respectively. The relative errors in atomic positions are dependent on the number and kinds of restraints used in the refinements. It is also shown that minimum variance weighting leads to mean-square deviations from target geometry in the refined structures which are smaller than the variances used in the distance weighting.

The 3-dimensional organisation of crystallin polypeptides into globular proteins and their interactions into higher order structures are important factors governing optical functions related to refraction, accommodation and transparency. Single crystal X-ray diffraction studies have revealed the tertiary and quaternary structural organisation of β-, γ- and δ-crystallins. Regions of the lens with high refractive index contain high levels of monomeric y-crystallins while the accommodating, hydrated avian lens has largely replaced γ-crystallins with δ-crystallin. The βγ-crystallins form a superfamily of proteins of high symmetry and great diversity in which the basic building block is a 10 kD pseudo-symmetrical 2-Greek key domain. A γ-crystallin comprises two of these β-sheet domains, joined by a linker, and a short C-terminal extension. In β-crystallins the linker has an extended conformation resulting in dimer formation by a mechanism known as domain swapping. Crystallographic analysis of engineered single domains of γ-crystallins, analogous to the ancestral domain, has indicated the importance of the short C-terminal extension in directing domain pairing. γ-crystallins have numerous cysteine residues, some are conserved in the core of the protein molecule and some are variable on the protein surface. The structure of γB-crystallin has been determined at very high resolution using cryo-crystallography allowing the visualisation of the complete protein-protein and protein-water structure at the surface. β-crystallins are seen as tetramers in the crystal structures but their long sequence extensions are harder to visualise in the electron density of the hydrated crystal lattice structure. In one tight packing lattice of βB2 crystallin the N-terminal extension is seen to mediate protein-protein interactions between tetramers to form a 42 helix. The X-ray structure of the taxon-restricted avian δ-crystallin shows that the 50 kD subunit contains 22 helices that form three α-helical domains which dimerise followed by a dimer-dimer interaction to form a tetramer with a 20-helix bundle at the centre. Analysis of the spatial disposition of the sequence conserved regions showed the location of the active site cleft of the superfamily of enzymes related to δ-crystallin and argininosuccinate lyase. A different crystal structure of δ-crystallin solved under more physiological conditions revealed that tetramers assembled as higher order supramolecular helices and that the N-terminal extension may be involved. Combining the observations of higher order helical structures in both the oligomeric β-crystallin and δ-crystallin crystal lattices, we have proposed a highly speculative model for crystallin assembly in the lens fibre cells.

Background: Clostridium perfringens toxin NetB is a key factor in avian necrotic enteritis. Results: NetB forms heptameric pores structurally similar to Staphylococcus aureus toxins but lacks a phosphocholine binding pocket. NetB activity is enhanced by cholesterol. Conclusion: NetB has distinct binding specificity, and cholesterol may act as a receptor. Significance: The structure of NetB will facilitate development of control measures against necrotic enteritis., NetB is a pore-forming toxin produced by Clostridium perfringens and has been reported to play a major role in the pathogenesis of avian necrotic enteritis, a disease that has emerged due to the removal of antibiotics in animal feedstuffs. Here we present the crystal structure of the pore form of NetB solved to 3.9 Å. The heptameric assembly shares structural homology to the staphylococcal α-hemolysin. However, the rim domain, a region that is thought to interact with the target cell membrane, shows sequence and structural divergence leading to the alteration of a phosphocholine binding pocket found in the staphylococcal toxins. Consistent with the structure we show that NetB does not bind phosphocholine efficiently but instead interacts directly with cholesterol leading to enhanced oligomerization and pore formation. Finally we have identified conserved and non-conserved amino acid positions within the rim loops that significantly affect binding and toxicity of NetB. These findings present new insights into the mode of action of these pore-forming toxins, enabling the design of more effective control measures against necrotic enteritis and providing potential new tools to the field of bionanotechnology.

Clostridium perfringens epsilon toxin (Etx) is a pore-forming toxin responsible for a severe and rapidly fatal enterotoxemia of ruminants. The toxin is classified as a category B bioterrorism agent by the U.S. Government Centres for Disease Control and Prevention (CDC), making work with recombinant toxin difficult. To reduce the hazard posed by work with recombinant Etx, we have used a variant of Etx that contains a H149A mutation (Etx-H149A), previously reported to have reduced, but not abolished, toxicity. The three-dimensional structure of H149A prototoxin shows that the H149A mutation in domain III does not affect organisation of the putative receptor binding loops in domain I of the toxin. Surface exposed tyrosine residues in domain I of Etx-H149A (Y16, Y20, Y29, Y30, Y36 and Y196) were mutated to alanine and mutants Y30A and Y196A showed significantly reduced binding to MDCK.2 cells relative to Etx-H149A that correlated with their reduced cytotoxic activity. Thus, our study confirms the role of surface exposed tyrosine residues in domain I of Etx in binding to MDCK cells and the suitability of Etx-H149A for further receptor binding studies. In contrast, binding of all of the tyrosine mutants to ACHN cells was similar to that of Etx-H149A, suggesting that Etx can recognise different cell surface receptors. In support of this, the crystal structure of Etx-H149A identified a glycan (β-octyl-glucoside) binding site in domain III of Etx-H149A, which may be a second receptor binding site. These findings have important implications for developing strategies designed to neutralise toxin activity.

The magnitude of the static and dynamic disorder which is usually associated with macromolecular crystals often severely limits the resolution of their Bragg diffraction patterns. The balance of the diffracted energy from these crystals occurs as diffuse scattering which may be an order of magnitude weaker than the Bragg reflections. The measurement of this diffuse scattering has become much easier with the advent of intense well-collimated X-ray sources provided by synchrotron radiation. Diffuse intensity measurements of macromolecular diffraction on films, and latterly image plates, have made possible the acquisition of information on the correlation of atomic displacements which has not been possible from Bragg reflections. Models of correlated displacements have been proposed and verified by comparison of calculated and observed diffracted patterns. Future work will focus on the systematic methods for measuring diffuse scattering from macromolecules and the interpretation of this data in terms of the correlated motions of domains or secondary structural elements of proteins and nucleic acids. Methods of displaying this information and relating it to biological function will also need to be developed. The weakness of much diffuse scattering underlines the important role that synchrotrons will have to play in these tasks.

beta-Crystallins are oligomeric eye lens proteins that are related to monomeric gamma-crystallins. The main sequence difference between the two families is the presence of sequence extensions in the beta-crystallins. A major question concerns the role that these extensions play in mediating interactions at the high protein concentrations found in the lens. The predominant beta-crystallin polypeptide, beta B2, can be crystallized in two different space groups, I222 and C222. The I222 crystal structure revealed that the protein packed as a tetramer with perfect 222 symmetry but that the extensions were disordered. The X-ray structure of the C222 lattice of beta B2 has now been refined at 3.3 A, the structure analysed and compared with the I222 lattice. The protein is also a tetramer with 222 symmetry in the C222 lattice but differs in that parts of the N-terminal extensions have been visualized. In the asymmetric unit of the C222 lattice there are four subunits, each comprising a single polypeptide chain, in which certain flexible loops in the N-terminal domains and the N-terminal extensions have various conformations. The tetramers in the C222 lattice are more tightly packed than in the I222 form. Analysis of the tetramer contacts shows that the sites of interaction break the 222 symmetry of the tetramers. The N-terminal extensions play a major role in directing interactions between tetramers. One of the N-terminal extensions interacts with a hydrophobic patch on the N-terminal domain of another tetramer. These crystallographic observations obtained over a physiological concentration range indicate how, in beta-crystallin oligomers, the N-terminal extensions of beta B2 can switch from interacting with water to interacting with protein depending on their relative concentrations. This could be useful in maintaining a gradient of refractive index.

An algorithm is described for automatically generating protein topology cartoons. This algorithm optimally places circles and triangles depicting alpha-helices and beta-strands respectively giving a pictorial topological summary of any protein structure. beta-Sheets, sandwiches and barrels are automatically identified and represented using special templates. The output from this algorithm may be controlled by adjustment of variable weights during the optimization step giving a preferred result. The rules for generating protein toplogy cartoons, including consideration of the handedness of local structure motifs, are discussed. The design of this algorithm is completely general and is easily adapted to include further rules that dictate the generation of the cartoons.

Clostridium perfringens enterotoxin (CPE) is a major cause of food poisoning and antibiotic-associated diarrhea. Upon its release from C. perfringens spores, CPE binds to its receptor, claudin, at the tight junctions between the epithelial cells of the gut wall and subsequently forms pores in the cell membranes. A number of different complexes between CPE and claudin have been observed, and the process of pore formation has not been fully elucidated. We have determined the three-dimensional structure of the soluble form of CPE in two crystal forms by X-ray crystallography, to a resolution of 2.7 and 4.0 A, respectively, and found that the N-terminal domain shows structural homology with the aerolysin-like β-pore-forming family of proteins. We show that CPE forms a trimer in both crystal forms and that this trimer is likely to be biologically relevant but is not the active pore form. We use these data to discuss models of pore formation.

The PROCHECK suite of programs provides a detailed check on the stereochemistry of a protein structure. Its outputs comprise a number of plots in PostScript format and a comprehensive residue-by-residue listing. These give an assessment of the overall quality of the structure as compared with well refined structures of the same resolution and also highlight regions that may need further investigation. The PROCHECK programs are useful for assessing the quality not only of protein structures in the process of being solved but also of existing structures and of those being modelled on known structures.

The molecular structure of calf γB-crystallin (previously called γII), a lens-specific protein, has been refined to a crystallographic R factor of 18.1% for all reflection data, between 8.0 and 1.47 A, 25 959 hkl measured at 293 (1) K. 230 water molecules have been defined by difference Fourier techniques and included in a restrained least-squares refinement. Difference Fourier maps clearly indicated the presence of multiple sites for the sulfur atoms of Cys 18 and Cys 22 which were therefore given coupled second-site occupancies during the refinement. The sulfur atom in the major position of Cys 22 is in the reduced state. Either of the Cys 18 sites can form a high-energy disulfide bridge with the minor position of Cys 22. The position of the carboxy terminus and many other surface side chains have been further defined including the RGD signal peptide. The hydration of the backbone and the interdomain region has been analysed. 27 water molecules make extensive contacts to a single protein molecule and thus contribute to its stability.

Motivation: Modelling antigenic shift in influenza A H3N2 can help to predict the efficiency of vaccines. The virus is known to exhibit sudden jumps in antigenic distance, and prediction of such novel strains from amino acid sequence differences remains a challenge. Results: From analysis of 6624 amino acid sequences of wild-type H3, we propose updates to the frequently referenced list of 131 amino acids located at or near the five identified antibody binding regions in haemagglutinin (HA). We introduce a class of predictive models based on the analysis of amino acid changes in these binding regions, and extend the principle to changes in HA1 as a whole by dividing the molecule into regional bands. Our results show that a range of simple models based on banded changes give better predictive performance than models based on the established five canonical regions and can identify a higher proportion of vaccine escape candidates among novel strains than a current state-of-the-art model. Contact: wlees01@mail.cryst.bbk.ac.uk Supplementary information: Supplementary Data is available at Bioinformatics online.

International audience; IMGT, the international ImMunoGeneTics information system (http://imgt.cines.fr), is the reference in immunogenetics and immunoinformatics. IMGT standardizes and manages the complex immunogenetic data that include the immunoglobulins (IG) or antibodies, the T cell receptors (TR), the major histocompatibility complex (MHC) and the related proteins of the immune system (RPI), which belong to the immunoglobulin superfamily (IgSF) and the MHC superfamily (MhcSF). The accuracy and consistency of IMGT data and the coherence between the different IMGT components (databases, tools and Web resources) are based on IMGT-ONTOLOGY, the first ontology for immunogenetics and immunoinformatics. IMGT-ONTOLOGY manages the immunogenetics knowledge through diverse facets relying on seven axioms, 'IDENTIFICATION', 'DESCRIPTION', 'CLASSIFICATION', 'NUMEROTATION', 'LOCALIZATION', 'ORIENTATION' and 'OBTENTION', that postulate that objects, processes and relations have to be identified, described, classified, numerotated, localized, orientated, and that the way they are obtained has to be determined. These axioms constitute the Formal IMGT-ONTOLOGY, also designated as IMGT-Kaleidoscope. These axioms have been essential for the conceptualization of the molecular immunogenetics knowledge and for the creation of IMGT. Indeed all the components of the IMGT integrated system have been developed, based on standardized concepts and relations, thus allowing IMGT to bridge biological and computational spheres in bioinformatics. The same axioms can be used to generate concepts for multi-scale level approaches at the molecule, cell, tissue, organ, organism or population level, emphasizing the generalization of the application domain. In that way the Formal IMGT-ONTOLOGY represents a paradigm for the elaboration of ontologies in system biology.

A range of protein single-crystal diffraction patterns, recorded with intense collimated synchrotron Xradiation are presented. These patterns illustrate the rich and varied nature of the diffuse scattering from this kind of crystal and suggest that the continuous background diffraction will be a source of specific information on the molecular dynamics and flexibility of particular proteins. The relative contributions of the background X-ray diffuse scattering from the solvent, the glass capillary and the sample have been quantified for two of the samples; the so-called solvent ring is shown to be due principally to protein disorder in the crystal both because of the intensity and the marked anisotropy of the ring in specific cases. The form of the acoustic scattering associated with the Bragg peaks was studied in ribonuclease, and is shown, at low resolution at least, to be explained in terms of single-phonon interactions.