Maximizing the number of polychronous groups in spiking networks

In this paper we investigate the effect of biasing the axonal connection delay values in the number of polychronous groups produced for a spiking neuron network model. We use an estimation of distribution algorithm (EDA) that learns tree models to search for optimal delay configurations. Our results indicate that the introduced approach can be used to considerably increase the number of such groups.

Gravitational actions upon a tether in a non-uniform gravity field with arbitrary number of zonal harmonics

We develop general closed-form expressions for the mutual gravitational potential, resultant and torque acting upon a rigid tethered system moving in a non-uniform gravity field produced by an attracting body with revolution symmetry, such that an arbitrary number of zonal harmonics is considered. The final expressions are series expansion in two small parameters related to the reference radius of the primary and the length of the tether, respectively, each of which are scaled by the mutual distance between their centers of mass. A few numerical experiments are performed to study the convergence behavior of the final expressions, and conclude that for high precision applications it might be necessary to take into account additional perturbation terms, which come from the mutual Two-Body interaction.

A simple mathematical model that describes the growth of the area and the number of total and viable cells in yeast colonies

We propose a model, based on the Gompertz equation, to describe the growth of yeasts colonies on agar medium. This model presents several advantages: (i) one equation describes the colony growth, which previously needed two separate ones (linear increase of radius and of the squared radius); (ii) a similar equation can be applied to total and viable cells, colony area or colony radius, because the number of total cells in mature colonies is proportional to their area; and (iii) its parameters estimate the cell yield, the cell concentration that triggers growth limitation and the effect of this limitation on the specific growth rate. To elaborate the model, area, total and viable cells of 600 colonies of Saccharomyces cerevisiae, Debaryomyces fabryi, Zygosaccharomyces rouxii and Rhodotorula glutinis have been measured. With low inocula, viable cells showed an initial short exponential phase when colonies were not visible. This phase was shortened with higher inocula. In visible or mature colonies, cell growth displayed Gompertz-type kinetics. It was concluded that the cells growth in colonies is similar to liquid cultures only during the first hours, the rest of the time they grow, with near-zero specific growth rates, at least for 3 weeks.

Mechanical and chemical unfolding of a single protein: A comparison

Is the mechanical unraveling of protein domains by atomic force microscopy (AFM) just a technological feat or a true measurement of their unfolding? By engineering a protein made of tandem repeats of identical Ig modules, we were able to get explicit AFM data on the unfolding rate of a single protein domain that can be accurately extrapolated to zero force. We compare this with chemical unfolding rates for untethered modules extrapolated to 0 M denaturant. The unfolding rates obtained by the two methods are the same. Furthermore, the transition state for unfolding appears at the same position on the folding pathway when assessed by either method. These results indicate that mechanical unfolding of a single protein by AFM does indeed reflect the same event that is observed in traditional unfolding experiments. The way is now open for the extensive use of AFM to measure folding reactions at the single-molecule level. Single-molecule AFM recordings have the added advantage that they define the reaction coordinate and expose rare unfolding events that cannot be observed in the absence of chemical denaturants.

The optimization principle in phylogenetic analysis tends to give incorrect topologies when the number of nucleotides or amino acids used is small

In the maximum parsimony (MP) and minimum evolution (ME) methods of phylogenetic inference, evolutionary trees are constructed by searching for the topology that shows the minimum number of mutational changes required (M) and the smallest sum of branch lengths (S), respectively, whereas in the maximum likelihood (ML) method the topology showing the highest maximum likelihood (A) of observing a given data set is chosen. However, the theoretical basis of the optimization principle remains unclear. We therefore examined the relationships of M, S, and A for the MP, ME, and ML trees with those for the true tree by using computer simulation. The results show that M and S are generally greater for the true tree than for the MP and ME trees when the number of nucleotides examined (n) is relatively small, whereas A is generally lower for the true tree than for the ML tree. This finding indicates that the optimization principle tends to give incorrect topologies when n is small. To deal with this disturbing property of the optimization principle, we suggest that more attention should be given to testing the statistical reliability of an estimated tree rather than to finding the optimal tree with excessive efforts. When a reliability test is conducted, simplified MP, ME, and ML algorithms such as the neighbor-joining method generally give conclusions about phylogenetic inference very similar to those obtained by the more extensive tree search algorithms.

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.

The number of unidentified amacrine cells in the mammalian retina

The three largest known populations of amacrine cells in the rabbit retina were stained with fluorescent probes in whole mounts and counted at a series of retinal eccentricities. The retinas were counterstained using a fluorescent DNA-binding molecule and the total number of nuclei in the inner nuclear layer were counted in confocal sections. From the total number of inner nuclear layer cells and the known fraction of them occupied by amacrine cells, the fraction of amacrine cells made up by the stained populations could be calculated. Starburst cells made up 3%, indoleamine-accumulating cells made up 4%, and AII cells made up 11% of all amacrine cells. By referring four smaller populations of amacrine cells to the number of indoleamine-accumulating cells, they were estimated to make up 4% of all amacrine cells. Thus, 78% of all amacrine cells in the rabbit’s retina are known only from isolated examples, if at all. This proportion is similar in the retinas of the mouse, cat, and monkey. It is likely that a substantial fraction of the local circuit neurons present in other regions of the central nervous system are also invisible as populations to current techniques.

Telomeric repeat amplification, without shortening or lengthening of the telomerase products: a method to analyze the processivity of telomerase enzyme

The Telomeric Repeat Amplification Protocol (TRAP) and its modified versions (including ours, TP-TRAP) change the size and/or the ratio of the telomerase products in the amplification stage of the assay. Based on our recently published method we developed a new TRAP. This method ensures that the number of telomeric repeats present in the original telomerase products does not change on PCR amplification. The usefulness of the method was proved with amplification of chemically synthesized telomerase products and a newly designed telomerase substrate oligonucleotide. This is the first report in which the PCR products directly reflect the size distribution of telomerase products generated by the enzyme.

Heat effects on DNA repair after ionising radiation: hyperthermia commonly increases the number of non-repaired double-strand breaks and structural rearrangements

After ionising radiation double-strand breaks (dsb) are lethal if not repaired or misrepaired. Cell killing is greatly enhanced by hyperthermia and it is questioned here whether heat not only affects dsb repair capacity but also fidelity in a chromosomal context. dsb repair experiments were designed so as to mainly score non-homologous end joining, while homologous recombination was largely precluded. Human male G0 fibroblasts were either preheated (45°C, 20 min) or not before X-irradiation. dsb induction and repair were measured by conventional gel electrophoresis and an assay combining restriction digestion using a rare cutting enzyme (NotI) and Southern hybridisation, which detects large chromosomal rearrangements (>100 kb). dsb induction rate in an X-chromosomal NotI fragment was 4.8 × 10–3 dsb/Gy/Mb. Similar values were found for the genome overall and also when cells were preheated. After 50 Gy, fibroblasts were competent to largely restore the original restriction fragment size. Five per cent of dsb remained non-rejoined and 14% were misrejoined. Correct restitution of restriction fragments occurred preferably during the first hour but continued at a slow rate for 12–16 h. In addition, dsb appeared to misrejoin throughout the entire repair period. After hyperthermia the fractions of non-rejoined and misrejoined dsb were similarly increased to 13 and 51%, respectively. It is suggested that heat increases the probability of dsb being incorrectly rejoined but it is not likely to interfere with one dsb repair pathway in particular.

Telomeric repeats (TTAGGC)n are sufficient for chromosome capping function in Caenorhabditis elegans.

Telomeres are specialized structures located at the ends of linear eukaryotic chromosomes that ensure their complete replication and protect them from fusion and degradation. We report here the characterization of the telomeres of the nematode Caenorhabditis elegans. We show that the chromosomes terminate in 4-9 kb of tandem repeats of the sequence TTAGGC. Furthermore, we have isolated clones corresponding to 11 of the 12 C. elegans telomeres. Their subtelomeric sequences are all different from each other, demonstrating that the terminal TTAGGC repeats are sufficient for general chromosomal capping functions. Finally, we demonstrate that the me8 meiotic mutant, which is defective in X chromosome crossing over and segregation, bears a terminal deficiency, that was healed by the addition of telomeric repeats, presumably by the activity of a telomerase enzyme. The 11 cloned telomeres represent an important advance for the completion of the physical map and for the determination of the entire sequence of the C. elegans genome.

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.

Multiplex genotype determination at a large number of gene loci.

To facilitate large-scale genotype analysis, an efficient PCR-based multiplex approach has been developed. For simultaneously amplifying the target sequences at a large number of genetic loci, locus-specific primers containing 5' universal tails are used. Attaching the universal tails to the target sequences in the initial PCR steps allows replacement of all specific primers with a pair of primers identical to the universal tails and converts the multiplex amplification into "uniplex." Simultaneous amplification of 26 genetic loci with this approach is described. The multiplex amplification can be coupled with genotype determination. By incorporating a single-base mismatch between a primer and the template into the target sequences, a polymorphic site can be converted into a desirable restriction fragment length polymorphism when it is necessary. In this way, the allelic PCR products for the polymorphic loci can be discriminated by gel electrophoresis after restriction enzyme digestion. In this study, 32 loci were typed in such a multiplex way.

Integration of multiple repeats of geminiviral DNA into the nuclear genome of tobacco during evolution.

Integration of viral DNA into the host nuclear genome, although not unusual in bacterial and animal systems, has surprisingly not been reported for plants. We have discovered geminvirus-related DNA (GRD) sequences, in the form of distinct sets of multiple direct repeats comprising three related repeat classes, situated in a unique locus in the Nicotiana tabacum (tobacco) nuclear genome. The organization of these sequences is similar or identical in eight different tobacco cultivars we have examined. DNA sequence analysis reveals that each repeat has sequences most resembling those of the New World geminiviral DNA replication origin plus the adjacent AL1 gene, encoding the viral replication protein. We believe these GRD sequences originated quite recently in Nicotiana evolution through integration of geminiviral DNA by some combination of the processes of illegitimate recombination, amplification, deletions, and rearrangements. These events must have occurred in plant tissue that was subsequently able to contribute to meristematic tissue yielding gametes. GRD may have been retained in tobacco by selection or by random fixation in a small evolving population. Although we cannot detect transcription of these sequences, this does not exclude the possibility that they may originally have been expressed.

Incorporation of glutamine repeats makes protein oligomerize: implications for neurodegenerative diseases.

Many transcription factors and some other proteins contain glutamine repeats; their abnormal expansion has been linked to several dominantly inherited neuro-degenerative diseases. Having found that poly(L-glutamine) alone forms beta-strands held together by hydrogen bonds between their amide groups, we surmised that glutamine repeats may form polar zippers, an unusual motif for protein-protein interactions. To test this hypothesis, we have engineered a Gly-Gln10-Gly peptide into the inhibitory loop of truncated chymotrypsin inhibitor 2 (CI2), a small protein from barley seeds, by both insertion and replacement. Gel filtration resolved both mutant inhibitors into at least three fractions, which analytical ultracentrifugation identified as monomers, dimers, and trimers of the recombinant protein; the truncated wild-type CI2 formed only monomers. CD difference spectra of the dimers and trimers versus wild type indicated that their glutamine repeats formed beta-pleated sheets, while those of the monomers versus wild type were more suggestive of type I beta-turns. The CD spectra of all three fractions remained unchanged even after incubation at 70 degrees C; neither the dimers nor the trimers dissociated at this temperature. We argue that the stability of all three fractions is due to the multiplicity of hydrogen bonds between extended strands of glutamine repeats in the oligomers or within a beta-hairpin formed by the single glutamine repeat of each monomer. Pathological effects may arise when expanded glutamine repeats cause proteins to acquire excessively high affinities for each other or for other proteins with glutamine repeats.

Mutational analysis of the conserved region 2 site of adenovirus E1A and its effect on binding to the retinoblastoma gene product: use of the "double-tagging" assay.

We have explored the feasibility of using a "double-tagging" assay for assessing which amino acids of a protein are responsible for its binding to another protein. We have chosen the adenovirus E1A-retinoblastoma gene product (pRB) proteins for a model system, and we focused on the high-affinity conserved region 2 of adenovirus E1A (CR2). We used site-specific mutagenesis to generate a mutant E1A gene with a lysine instead of an aspartic acid at position 121 within the CR2 site. We demonstrated that this mutant exhibited little binding to pRB by the double-tagging assay. We also have shown that this lack of binding is not due to any significant decrease in the level of expression of the beta-galactosidase-E1A fusion protein. We then created a "library" of phage expressing beta-galactosidase-E1A fusion proteins with a variety of different mutations within CR2. This library of E1A mutations was used in a double-tagging screening to identify mutant clones that bound to pRB. Three classes of phage were identified: the vast majority of clones were negative and exhibited no binding to pRB. Approximately 1 in 10,000 bound to pRB but not to E1A ("true positives"). A variable number of clones appeared to bind equally well to both pRB and E1A ("false positives"). The DNA sequence of 10 true positive clones yielded the following consensus sequence: DLTCXEX, where X = any amino acid. The recovery of positive clones with only one of several allowed amino acids at each position suggests that most, if not all, of the conserved residues play an important role in binding to pRB. On the other hand, the DNA sequence of the negative clones appeared random. These results are consistent with those obtained from other sources. These data suggest that a double-tagging assay can be employed for determining which amino acids of a protein are important for specifying its interaction with another protein if the complex forms within bacteria. This assay is rapid and up to 1 x 10(6) mutations can be screened at one time.