A general additive distance with time-reversibility and rate variation among nucleotide sites.

As additivity is a very useful property for a distance measure, a general additive distance is proposed under the stationary time-reversible (SR) model of nucleotide substitution or, more generally, under the stationary, time-reversible, and rate variable (SRV) model, which allows rate variation among nucleotide sites. A method for estimating the mean distance and the sampling variance is developed. In addition, a method is developed for estimating the variance-covariance matrix of distances, which is useful for the statistical test of phylogenies and molecular clocks. Computer simulation shows (i) if the sequences are longer than, say, 1000 bp, the SR method is preferable to simpler methods; (ii) the SR method is robust against deviations from time-reversibility; (iii) when the rate varies among sites, the SRV method is much better than the SR method because the distance is seriously underestimated by the SR method; and (iv) our method for estimating the sampling variance is accurate for sequences longer than 500 bp. Finally, a test is constructed for testing whether DNA evolution follows a general Markovian model.

Variable efficacy of repeated annual influenza vaccination

Conclusions have differed in studies that have compared vaccine efficacy in groups receiving influenza vaccine for the first time to efficacy in groups vaccinated more than once. For example, the Hoskins study [Hoskins, T. W., Davis, J. R., Smith, A. J., Miller, C. L. & Allchin, A. (1979) Lancet i, 33–35] concluded that repeat vaccination was not protective in the long term, whereas the Keitel study [Keitel, W. A., Cate, T. R., Couch, R. B., Huggins, L. L. & Hess, K. R. (1997) Vaccine 15, 1114–1122] concluded that repeat vaccination provided continual protection. We propose an explanation, the antigenic distance hypothesis, and test it by analyzing seven influenza outbreaks that occurred during the Hoskins and Keitel studies. The hypothesis is that variation in repeat vaccine efficacy is due to differences in antigenic distances among vaccine strains and between the vaccine strains and the epidemic strain in each outbreak. To test the hypothesis, antigenic distances were calculated from historical hemagglutination inhibition assay tables, and a computer model of the immune response was used to predict the vaccine efficacy of individuals given different vaccinations. The model accurately predicted the observed vaccine efficacies in repeat vaccinees relative to the efficacy in first-time vaccinees (correlation 0.87). Thus, the antigenic distance hypothesis offers a parsimonious explanation of the differences between and within the Hoskins and Keitel studies. These results have implications for the selection of influenza vaccine strains, and also for vaccination strategies for other antigenically variable pathogens that might require repeated vaccination.

A neuronal model of a global workspace in effortful cognitive tasks

A minimal hypothesis is proposed concerning the brain processes underlying effortful tasks. It distinguishes two main computational spaces: a unique global workspace composed of distributed and heavily interconnected neurons with long-range axons, and a set of specialized and modular perceptual, motor, memory, evaluative, and attentional processors. Workspace neurons are mobilized in effortful tasks for which the specialized processors do not suffice. They selectively mobilize or suppress, through descending connections, the contribution of specific processor neurons. In the course of task performance, workspace neurons become spontaneously coactivated, forming discrete though variable spatio-temporal patterns subject to modulation by vigilance signals and to selection by reward signals. A computer simulation of the Stroop task shows workspace activation to increase during acquisition of a novel task, effortful execution, and after errors. We outline predictions for spatio-temporal activation patterns during brain imaging, particularly about the contribution of dorsolateral prefrontal cortex and anterior cingulate to the workspace.

Mouse strain differences determine severity of iron accumulation in Hfe knockout model of hereditary hemochromatosis

Hereditary hemochromatosis (HH) is a common disorder of iron metabolism caused by mutation in HFE, a gene encoding an MHC class I-like protein. Clinical studies demonstrate that the severity of iron loading is highly variable among individuals with identical HFE genotypes. To determine whether genetic factors other than Hfe genotype influence the severity of iron loading in the murine model of HH, we bred the disrupted murine Hfe allele onto three different genetically defined mouse strains (AKR, C57BL/6, and C3H), which differ in basal iron status and sensitivity to dietary iron loading. Serum transferrin saturations (percent saturation of serum transferrin with iron), hepatic and splenic iron concentrations, and hepatocellular iron distribution patterns were compared for wild-type (Hfe +/+), heterozygote (Hfe +/−), and knockout (Hfe −/−) mice from each strain. Although the Hfe −/− mice from all three strains demonstrated increased transferrin saturations and liver iron concentrations compared with Hfe +/+ mice, strain differences in severity of iron accumulation were striking. Targeted disruption of the Hfe gene led to hepatic iron levels in Hfe −/− AKR mice that were 2.5 or 3.6 times higher than those of Hfe −/− C3H or Hfe −/− C57BL/6 mice, respectively. The Hfe −/− mice also demonstrated strain-dependent differences in transferrin saturation, with the highest values in AKR mice and the lowest values in C3H mice. These observations demonstrate that heritable factors markedly influence iron homeostasis in response to Hfe disruption. Analysis of mice from crosses between C57BL/6 and AKR mice should allow the mapping and subsequent identification of genes modifying the severity of iron loading in this murine model of HH.

Structural Analysis and Molecular Model of a Self-Incompatibility RNase from Wild Tomato1

Self-incompatibility RNases (S-RNases) are an allelic series of style glycoproteins associated with rejection of self-pollen in solanaceous plants. The nucleotide sequences of S-RNase alleles from several genera have been determined, but the structure of the gene products has only been described for those from Nicotiana alata. We report on the N-glycan structures and the disulfide bonding of the S3-RNase from wild tomato (Lycopersicon peruvianum) and use this and other information to construct a model of this molecule. The S3-RNase has a single N-glycosylation site (Asn-28) to which one of three N-glycans is attached. S3-RNase has seven Cys residues; six are involved in disulfide linkages (Cys-16-Cys-21, Cys-46-Cys-91, and Cys-166-Cys-177), and one has a free thiol group (Cys-150). The disulfide-bonding pattern is consistent with that observed in RNase Rh, a related RNase for which radiographic-crystallographic information is available. A molecular model of the S3-RNase shows that four of the most variable regions of the S-RNases are clustered on one surface of the molecule. This is discussed in the context of recent experiments that set out to determine the regions of the S-RNase important for recognition during the self-incompatibility response.

Evolutionarily conserved Galphabetagamma binding surfaces support a model of the G protein-receptor complex.

The pivotal role of G proteins in sensory, hormonal, inflammatory, and proliferative responses has provoked intense interest in understanding how they interact with their receptors and effectors. Nonetheless, the locations of the receptors and effector binding sites remain poorly characterized, although nearly complete structures of the alphabetagamma heterotrimeric complex are available. Here we apply evolutionary trace (ET) analysis [Lichtarge, O., Bourne, H. R. & Cohen, F. E. (1996) J. Mol. Biol. 257, 342-358] to propose plausible locations for these sites. On each subunit, ET identifies evolutionarily selected surfaces composed of residues that do not vary within functional subgroups and that form spatial clusters. Four clusters correctly identify subunit interfaces, and additional clusters on Galpha point to likely receptor or effector binding sites. Our results implicate the conformationally variable region of Galpha in an effector binding role. Furthermore the range of predicted interactions between the receptor and Galphabetagamma, is sufficiently limited that we can build a low resolution and testable model of the receptor-G protein complex.

Novel unconventional binding site in the variable region of immunoglobulins.

The variable immunoglobulin (Ig) domains contain hypervariable regions that are involved in the formation of the antigen binding site. Besides the canonical antigen binding site, so-called unconventional sites also reside in the variable region that bind bacterial and viral proteins. Docking to these unconventional sites does not typically interfere with antigen binding, which suggests that these sites may be a part of the biological functions of Igs. Herein, a novel unconventional binding site is described. The site is detected with 8-azidopurine nucleotide photoaffinity probes that label antibodies efficiently and under mild conditions. Tryptic peptides were isolated from photolabeled monoclonal antibodies and aligned with the variable antibody domains of heavy and light chains. The structure of a variable Ig fragment was used to model the binding of the purine nucleotide to invariant residues in a hydrophobic pocket of the Ig molecule at a location distant from the antigen binding site. Monoclonal and polyclonal antibodies were biotinylated with the photoaffinity linker and used in fluorescence-activated cell sorter and ELISA analyses. The data support the utility of this site for tethering diagnostic and therapeutic agents to the variable Ig fragment region without impairing the structural and functional integrity of antibodies.

Tertiary structure of an amyloid immunoglobulin light chain protein: a proposed model for amyloid fibril formation.

An immunoglobulin light chain protein was isolated from the urine of an individual (BRE) with systemic amyloidosis. Complete amino acid sequence of the variable region of the light chain (VL) protein established it as a kappa I, which when compared with other kappa I amyloid associated proteins had unique residues, including Ile-34, Leu-40, and Tyr-71. To study the tertiary structure, BRE VL was expressed in Escherichia coli by using a PCR product amplified from the patient BRE's bone marrow DNA. The PCR product was ligated into pCZ11, a thermal-inducible replication vector. Recombinant BRE VL was isolated, purified to homogeneity, and crystallized by using ammonium sulfate as the precipitant. Two crystal forms were obtained. In crystal form I the BRE VL kappa domain crystallizes as a dimer with unit cell constants isomorphous to previously published kappa protein structures. Comparison with a nonamyloid VL kappa domain from patient REI, identified significant differences in position of residues in the hypervariable segments plus variations in framework region (FR) segments 40-46 (FR2) and 66-67 (FR3). In addition, positional differences can be seen along the two types of local diads, corresponding to the monomer-monomer and dimer-dimer interfaces. From the packing diagram, a model for the amyloid light chain (AL) fibril is proposed based on a pseudohexagonal spiral structure with a rise of approximately the width of two dimers per 360 degree turn. This spiral structure could be consistent with the dimensions of amyloid fibrils as determined by electron microscopy.