Equilibrium distributions of microsatellite repeat length resulting from a balance between slippage events and point mutations

We describe and test a Markov chain model of microsatellite evolution that can explain the different distributions of microsatellite lengths across different organisms and repeat motifs. Two key features of this model are the dependence of mutation rates on microsatellite length and a mutation process that includes both strand slippage and point mutation events. We compute the stationary distribution of allele lengths under this model and use it to fit DNA data for di-, tri-, and tetranucleotide repeats in humans, mice, fruit flies, and yeast. The best fit results lead to slippage rate estimates that are highest in mice, followed by humans, then yeast, and then fruit flies. Within each organism, the estimates are highest in di-, then tri-, and then tetranucleotide repeats. Our estimates are consistent with experimentally determined mutation rates from other studies. The results suggest that the different length distributions among organisms and repeat motifs can be explained by a simple difference in slippage rates and that selective constraints on length need not be imposed.

Defining the minimal length of sequence homology required for selective gene isolation by TAR cloning

The transformation-associated recombination (TAR) cloning technique allows selective and accurate isolation of chromosomal regions and genes from complex genomes. The technique is based on in vivo recombination between genomic DNA and a linearized vector containing homologous sequences, or hooks, to the gene of interest. The recombination occurs during transformation of yeast spheroplasts that results in the generation of a yeast artificial chromosome (YAC) containing the gene of interest. To further enhance and refine the TAR cloning technology, we determined the minimal size of a specific hook required for gene isolation utilizing the Tg.AC mouse transgene as a targeted region. For this purpose a set of vectors containing a B1 repeat hook and a Tg.AC-specific hook of variable sizes (from 20 to 800 bp) was constructed and checked for efficiency of transgene isolation by a radial TAR cloning. When vectors with a specific hook that was ≥60 bp were utilized, ∼2% of transformants contained circular YACs with the Tg.AC transgene sequences. Efficiency of cloning dramatically decreased when the TAR vector contained a hook of 40 bp or less. Thus, the minimal length of a unique sequence required for gene isolation by TAR is ∼60 bp. No transgene-positive YAC clones were detected when an ARS element was incorporated into a vector, demonstrating that the absence of a yeast origin of replication in a vector is a prerequisite for efficient gene isolation by TAR cloning.

Enhancing the catalytic repertoire of nucleic acids: a systematic study of linker length and rigidity

The incorporation of potentially catalytic groups in DNA is of interest for the in vitro selection of novel deoxyribozymes. A series of 10 C5-modified analogues of 2′-deoxyuridine triphosphate have been synthesised that possess side chains of differing flexibility and bearing a primary amino or imidazole functionality. For each series of nucleotide analogues differing degrees of flexibility of the C5 side chain was achieved through the use of alkynyl, alkenyl and alkyl moieties. The imidazole function was conjugated to these C5-amino-modified nucleotides using either imidazole 4-acetic acid or imidazole 4-acrylic acid (urocanic acid). The substrate properties of the nucleotides (fully replacing dTTP) with Taq polymerase during PCR have been investigated in order to evaluate their potential applications for in vitro selection experiments. 5-(3-Aminopropynyl)dUTP and 5-(E-3-aminopropenyl)dUTP and their imidazole 4-acetic acid- and urocanic acid-modified conjugates were found to be substrates. In contrast, C5-amino-modified dUTPs with alkane or Z-alkene linkers and their corresponding conjugates were not substrates. The incorporation of these analogues during PCR has been confirmed by inhibition of restriction enzyme digestion using XbaI and by mass spectrometry of the PCR products.

A modified gambler's ruin model of polyethylene chains in the amorphous region.

Polyethylene chains in the amorphous region between two crystalline lamellae M unit apart are modeled as random walks with one-step memory on a cubic lattice between two absorbing boundaries. These walks avoid the two preceding steps, though they are not true self-avoiding walks. Systems of difference equations are introduced to calculate the statistics of the restricted random walks. They yield that the fraction of loops is (2M - 2)/(2M + 1), the fraction of ties 3/(2M + 1), the average length of loops 2M - 0.5, the average length of ties 2/3M2 + 2/3M - 4/3, the average length of walks equals 3M - 3, the variance of the loop length 16/15M3 + O(M2), the variance of the tie length 28/45M4 + O(M3), and the variance of the walk length 2M3 + O(M2).

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.

A new high molecular weight immunoglobulin class from the carcharhine shark: implications for the properties of the primordial immunoglobulin.

All immunoglobulins and T-cell receptors throughout phylogeny share regions of highly conserved amino acid sequence. To identify possible primitive immunoglobulins and immunoglobulin-like molecules, we utilized 3' RACE (rapid amplification of cDNA ends) and a highly conserved constant region consensus amino acid sequence to isolate a new immunoglobulin class from the sandbar shark Carcharhinus plumbeus. The immunoglobulin, termed IgW, in its secreted form consists of 782 amino acids and is expressed in both the thymus and the spleen. The molecule overall most closely resembles mu chains of the skate and human and a new putative antigen binding molecule isolated from the nurse shark (NAR). The full-length IgW chain has a variable region resembling human and shark heavy-chain (VH) sequences and a novel joining segment containing the WGXGT motif characteristic of H chains. However, unlike any other H-chain-type molecule, it contains six constant (C) domains. The first C domain contains the cysteine residue characteristic of C mu1 that would allow dimerization with a light (L) chain. The fourth and sixth domains also contain comparable cysteines that would enable dimerization with other H chains or homodimerization. Comparison of the sequences of IgW V and C domains shows homology greater than that found in comparisons among VH and C mu or VL, or CL thereby suggesting that IgW may retain features of the primordial immunoglobulin in evolution.

Yeast artificial chromosome contigs reveal that distal variable-region genes reside at least 3 megabases from the joining regions in the murine immunoglobulin kappa locus.

The immunoglobulin kappa gene locus encodes 95% of the light chains of murine antibody molecules and is thought to contain up to 300 variable (V kappa)-region genes generally considered to comprise 20 families. To delineate the locus we have isolated 29 yeast artificial chromosome genomic clones that form two contigs, span > 3.5 megabases, and contain two known non-immunoglobulin kappa markers. Using PCR primers specific for 19 V kappa gene families and Southern analysis, we have refined the genetically defined order of these V kappa gene families. Of these, V kappa 2 maps at least 3.0 Mb from the joining (J kappa) region and appears to be the most distal V kappa gene segment.

Restricted and conserved T-cell repertoires involved in allorecognition of class II major histocompatibility complex.

The nature of the alloreactive T-cell response is not yet clearly understood. These strong cellular responses are thought to be the basis of allograft rejection and graft-vs.-host disease. The question of the extent of responding T-cell repertoires has so far been addressed by cellular cloning, often combined with molecular T-cell receptor (TCR) analysis. Here we present a broad repertoire analysis of primed responder cells from mixed lymphocyte cultures in which two different DR1/3 responders were stimulated with DR3/4 cells. Repertoire analysis was performed by TCR spectratyping, a method by which T cells are analyzed on the basis of the complementarity-determining region 3 length of different variable region (V) families. Strikingly, both responders showed very similar repertoires when the TCR V beta was used as a lineage marker. This was not seen when TCR V alpha was analyzed. A different pattern of TCR V beta was observed if the stimulating alloantigen was changed. This finding indicates that alloreactive T cells form a specific repertoire for each alloantigen. Since conservation appears to be linked to TCR V beta, the question of different roles of alpha and beta chains in allorecognition is raised.

Evolutionary consequences of food chain length in kelp forest communities.

Kelp forests are strongly influenced by macroinvertebrate grazing on fleshy macroalgae. In the North Pacific Ocean, sea otter predation on macroinvertebrates substantially reduces the intensity of herbivory on macroalgae. Temperate Australasia, in contrast, has no known predator of comparable influence. These ecological and biogeographic patterns led us to predict that (i) the intensity of herbivory should be greater in temperate Australasia than in the North Pacific Ocean; thus (ii) Australasian seaweeds have been under stronger selection to evolve chemical defenses and (iii) Australasian herbivores have been more strongly selected to tolerate these compounds. We tested these predictions first by measuring rates of algal tissue loss to herbivory at several locations in Australasian and North Pacific kelp forests. There were significant differences in grazing rates among sea otter-dominated locations in the North Pacific (0-2% day-1), Australasia (5-7% day-1), and a North Pacific location lacking sea otters (80% day-1). The expectations that chronically high rates of herbivory in Australasia have selected for high concentrations of defensive secondary metabolites (phlorotannins) in brown algae and increased tolerance of these defenses in the herbivores also were supported. Phlorotannin concentrations in kelps and fucoids from Australasia were, on average, 5-6 times higher than those in a comparable suite of North Pacific algae, confirming earlier findings. Furthermore, feeding rates of Australasian herbivores were largely unaffected by phlorotannins, regardless of the compounds' regional source. North Pacific herbivores, in contrast, were consistently deterred by phlorotannins from both Australasia and the North Pacific. These findings suggest that top-level consumers, acting through food chains of various lengths, can strongly influence the ecology and evolution of plantherbivore interactions.

Calibration of the length of a chain of single gold atoms

Using a scanning tunnelling microscope or mechanically controllable break junction it has been shown that it is possible to control the formation of a wire made of single gold atoms. In these experiments an interatomic distance between atoms in the chain of ∼3.6 Å was reported which is not consistent with recent theoretical calculations. Here, using precise calibration procedures for both techniques, we measure the length of the atomic chains. Based on the distance between the peaks observed in the chain length histogram we find the mean value of the interatomic distance before chain rupture to be 2.5±0.2 Å. This value agrees with the theoretical calculations for the bond length. The discrepancy with the previous experimental measurements was due to the presence of He gas, that was used to promote the thermal contact, and which affects the value of the work function that is commonly used to calibrate distances in scanning tunnelling microscopy and mechanically controllable break junctions at low temperatures.

The high-bias stability of monatomic chains

For the metals Au, Pt and Ir it is possible to form freely suspended monatomic chains between bulk electrodes. The atomic chains sustain very large current densities, but finally fail at high bias. We investigate the breaking mechanism, that involves current-induced heating of the atomic wires and electromigration forces. We find good agreement of the observations for Au based on models due to Todorov and co-workers. The high-bias breaking of atomic chains for Pt can also be described by the models, although here the parameters have not been obtained independently. In the limit of long chains the breaking voltage decreases inversely proportional to the length.

Local probe of fractional edge states of S=1 Heisenberg spin chains

Spin chains are among the simplest physical systems in which electron-electron interactions induce novel states of matter. Here we propose to combine atomic scale engineering and spectroscopic capabilities of state of the art scanning tunnel microscopy to probe the fractionalized edge states of individual atomic scale S=1 spin chains. These edge states arise from the topological order of the ground state in the Haldane phase. We also show that the Haldane gap and the spin-spin correlation length can be measured with the same technique.

Quantum theory of spin waves in finite chiral spin chains

We calculate the effect of spin waves on the properties of finite-size spin chains with a chiral spin ground state observed on biatomic Fe chains deposited on iridium(001). The system is described with a Heisenberg model supplemented with a Dzyaloshinskii-Moriya coupling and a uniaxial single ion anisotropy that presents a chiral spin ground state. Spin waves are studied using the Holstein-Primakoff boson representation of spin operators. Both the renormalized ground state and the elementary excitations are found by means of Bogoliubov transformation, as a function of the two variables that can be controlled experimentally, the applied magnetic field and the chain length. Three main results are found. First, because of the noncollinear nature of the classical ground state, there is a significant zero-point reduction of the ground-state magnetization of the spin spiral. Second, there is a critical external field from which the ground state changes from chiral spin ground state to collinear ferromagnetic order. The character of the two lowest-energy spin waves changes from edge modes to confined bulk modes over this critical field. Third, in the spin-spiral state, the spin-wave spectrum exhibits oscillatory behavior as function of the chain length with the same period of the spin helix.

Sequences of immunoglobulin chains : tabulation and analysis of amino acid sequences of precursors, V-regions, C-regions, J-chain and -microglobulins /

1976 ed. issued under title: Variable regions of immunoglobulin chains.

Diffusion modeling of percutaneous absorption kinetics: 3. Variable diffusion and partition coefficients, consequences for stratum corneum depth profiles and desorption kinetics

Stratum corneum (SC) desorption experiments have yielded higher calculated steady-state fluxes than those obtained by epidermal penetration studies. A possible explanation of this result is a variable diffusion or partition coefficient across the SC. We therefore developed the diffusion model for percutaneous penetration and desorption to study the effects of either a variable diffusion coefficient or variable partition coefficient in the SC over the diffusion path length. Steady-state flux, lag time, and mean desorption time were obtained from Laplace domain solutions. Numerical inversion of the Laplace domain solutions was used for simulations of solute concentration-distance and amount penetrated (desorbed)-time profiles. Diffusion and partition coefficients heterogeneity were examined using six different models. The effect of heterogeneity on predicted flux from desorption studies was compared with that obtained in permeation studies. Partition coefficient heterogeneity had a more profound effect on predicted fluxes than diffusion coefficient heterogeneity. Concentration-distance profiles show even larger dependence on heterogeneity, which is consistent with experimental tape-stripping data reported for clobetasol propionate and other solutes. The clobetasol propionate tape-stripping data were most consistent with the partition coefficient decreasing exponentially for half the SC and then becoming a constant for the remaining SC. (C) 2004 Wiley-Liss, Inc.