Broken symmetry and strangeness of the semiconductor impurity band metal-insulator transition

The filamentary model of the metal-insulator transition in randomly doped semiconductor impurity bands is geometrically equivalent to similar models for continuous transitions in dilute antiferromagnets and even to the λ transition in liquid He, but the critical behaviors are different. The origin of these differences lies in two factors: quantum statistics and the presence of long range Coulomb forces on both sides of the transition in the electrical case. In the latter case, in addition to the main transition, there are two satellite transitions associated with disappearance of the filamentary structure in both insulating and metallic phases. These two satellite transitions were first identified by Fritzsche in 1958, and their physical origin is explained here in geometrical and topological terms that facilitate calculation of critical exponents.

Different Domains of the M-Band Protein Myomesin Are Involved in Myosin Binding and M-Band TargetingV⃞

Myomesin is a 185-kDa protein located in the M-band of striated muscle where it interacts with myosin and titin, possibly connecting thick filaments with the third filament system. By using expression of epitope-tagged myomesin fragments in cultured cardiomyocytes and biochemical binding assays, we could demonstrate that the M-band targeting activity and the myosin-binding site are located in different domains of the molecule. An N-terminal immunoglobulin-like domain is sufficient for targeting to the M-band, but solid-phase overlay assays between individual N-terminal domains and the thick filament protein myosin revealed that the unique head domain contains the myosin-binding site. When expressed in cardiomyocytes, the head domains of rat and chicken myomesin showed species-specific differences in their incorporation pattern. The head domain of rat myomesin localized to a central area within the A-band, whereas the head domain of chicken myomesin was diffusely distributed in the cytoplasm. We therefore conclude that the head domain of myomesin binds to myosin but that this affinity is not sufficient for the restriction of the domain to the M-band in vivo. Instead, the neighboring immunoglobulin-like domain is essential for the precise incorporation of myomesin into the M-band, possibly because of interaction with a yet unknown protein of the sarcomere.

Mapping of a Myosin-binding Domain and a Regulatory Phosphorylation Site in M-Protein, a Structural Protein of the Sarcomeric M Band

The myofibrils of cross-striated muscle fibers contain in their M bands cytoskeletal proteins whose main function seems to be the stabilization of the three-dimensional arrangement of thick filaments. We identified two immunoglobin domains (Mp2–Mp3) of M-protein as a site binding to the central region of light meromyosin. This binding is regulated in vitro by phosphorylation of a single serine residue (Ser76) in the immediately adjacent amino-terminal domain Mp1. M-protein phosphorylation by cAMP-dependent kinase A inhibits binding to myosin LMM. Transient transfection studies of cultured cells revealed that the myosin-binding site seems involved in the targeting of M-protein to its location in the myofibril. Using the same method, a second myofibril-binding site was uncovered in domains Mp9–Mp13. These results support the view that specific phosphorylation events could be also important for the control of sarcomeric M band formation and remodeling.

The estimation of statistical parameters for local alignment score distributions

The distribution of optimal local alignment scores of random sequences plays a vital role in evaluating the statistical significance of sequence alignments. These scores can be well described by an extreme-value distribution. The distribution’s parameters depend upon the scoring system employed and the random letter frequencies; in general they cannot be derived analytically, but must be estimated by curve fitting. For obtaining accurate parameter estimates, a form of the recently described ‘island’ method has several advantages. We describe this method in detail, and use it to investigate the functional dependence of these parameters on finite-length edge effects.

BAliBASE (Benchmark Alignment dataBASE): enhancements for repeats, transmembrane sequences and circular permutations

BAliBASE is specifically designed to serve as an evaluation resource to address all the problems encountered when aligning complete sequences. The database contains high quality, manually constructed multiple sequence alignments together with detailed annotations. The alignments are all based on three-dimensional structural superpositions, with the exception of the transmembrane sequences. The first release provided sets of reference alignments dealing with the problems of high variability, unequal repartition and large N/C-terminal extensions and internal insertions. Here we describe version 2.0 of the database, which incorporates three new reference sets of alignments containing structural repeats, trans­membrane sequences and circular permutations to evaluate the accuracy of detection/prediction and alignment of these complex sequences. BAliBASE can be viewed at the web site http://www-igbmc.u-strasbg.fr/BioInfo/BAliBASE2/index.html or can be downloaded from ftp://ftp-igbmc.u-strasbg.fr/pub/BAliBASE2/.

A database and tools for 3-D protein structure comparison and alignment using the Combinatorial Extension (CE) algorithm

The database reported here is derived using the Combinatorial Extension (CE) algorithm which compares pairs of protein polypeptide chains and provides a list of structurally similar proteins along with their structure alignments. Using CE, structure–structure alignments can provide insights into biological function. When a protein of known function is shown to be structurally similar to a protein of unknown function, a relationship might be inferred; a relationship not necessarily detectable from sequence comparison alone. Establishing structure–structure relationships in this way is of great importance as we enter an era of structural genomics where there is a likelihood of an increasing number of structures with unknown functions being determined. Thus the CE database is an example of a useful tool in the annotation of protein structures of unknown function. Comparisons can be performed on the complete PDB or on a structurally representative subset of proteins. The source protein(s) can be from the PDB (updated monthly) or uploaded by the user. CE provides sequence alignments resulting from structural alignments and Cartesian coordinates for the aligned structures, which may be analyzed using the supplied Compare3D Java applet, or downloaded for further local analysis. Searches can be run from the CE web site, http://cl.sdsc.edu/ce.html, or the database and software downloaded from the site for local use.

PALI—a database of Phylogeny and ALIgnment of homologous protein structures

PALI (release 1.2) contains three-dimensional (3-D) structure-dependent sequence alignments as well as structure-based phylogenetic trees of homologous protein domains in various families. The data set of homologous protein structures has been derived by consulting the SCOP database (release 1.50) and the data set comprises 604 families of homologous proteins involving 2739 protein domain structures with each family made up of at least two members. Each member in a family has been structurally aligned with every other member in the same family (pairwise alignment) and all the members in the family are also aligned using simultaneous super­position (multiple alignment). The structural alignments are performed largely automatically, with manual interventions especially in the cases of distantly related proteins, using the program STAMP (version 4.2). Every family is also associated with two dendrograms, calculated using PHYLIP (version 3.5), one based on a structural dissimilarity metric defined for every pairwise alignment and the other based on similarity of topologically equivalent residues. These dendrograms enable easy comparison of sequence and structure-based relationships among the members in a family. Structure-based alignments with the details of structural and sequence similarities, superposed coordinate sets and dendrograms can be accessed conveniently using a web interface. The database can be queried for protein pairs with sequence or structural similarities falling within a specified range. Thus PALI forms a useful resource to help in analysing the relationship between sequence and structure variation at a given level of sequence similarity. PALI also contains over 653 ‘orphans’ (single member families). Using the web interface involving PSI_BLAST and PHYLIP it is possible to associate the sequence of a new protein with one of the families in PALI and generate a phylogenetic tree combining the query sequence and proteins of known 3-D structure. The database with the web interfaced search and dendrogram generation tools can be accessed at http://pa uling.mbu.iisc.ernet.in/~pali.

STACK: Sequence Tag Alignment and Consensus Knowledgebase

STACK is a tool for detection and visualisation of expressed transcript variation in the context of developmental and pathological states. The datasystem organises and reconstructs human transcripts from available public data in the context of expression state. The expression state of a transcript can include developmental state, pathological association, site of expression and isoform of expressed transcript. STACK consensus transcripts are reconstructed from clusters that capture and reflect the growing evidence of transcript diversity. The comprehensive capture of transcript variants is achieved by the use of a novel clustering approach that is tolerant of sub-sequence diversity and does not rely on pairwise alignment. This is in contrast with other gene indexing projects. STACK is generated at least four times a year and represents the exhaustive processing of all publicly available human EST data extracted from GenBank. This processed information can be explored through 15 tissue-specific categories, a disease-related category and a whole-body index and is accessible via WWW at http://www.sanbi.ac.za/Dbases.html. STACK represents a broadly applicable resource, as it is the only reconstructed transcript database for which the tools for its generation are also broadly available (http://www.sanbi.ac.za/CODES).

Allelic loss on chromosome band 18p11.3 occurs early and reveals heterogeneity in breast cancer progression

We examined the stage specificity and heterogeneity of 18p11 alterations in a series of tumors representing 96 microdissected samples. Significant loss of heterozygosity (LOH) (63%) was found, with 56% occurring early in ductal carcinoma in situ. Although most cases indicated LOH was clonally inherited, heterogeneity for 18p LOH occurred in 27% of tumors. When compared with other LOH data, 18p LOH was found in conjunction with allelic deletion on 3p, 9p, 17p and 17q, while 13q, 16q, and 11p were less frequently associated. These analyses suggest chromosome 18p11 alteration is a common and early event in breast disease.

ParAlign: a parallel sequence alignment algorithm for rapid and sensitive database searches

There is a need for faster and more sensitive algorithms for sequence similarity searching in view of the rapidly increasing amounts of genomic sequence data available. Parallel processing capabilities in the form of the single instruction, multiple data (SIMD) technology are now available in common microprocessors and enable a single microprocessor to perform many operations in parallel. The ParAlign algorithm has been specifically designed to take advantage of this technology. The new algorithm initially exploits parallelism to perform a very rapid computation of the exact optimal ungapped alignment score for all diagonals in the alignment matrix. Then, a novel heuristic is employed to compute an approximate score of a gapped alignment by combining the scores of several diagonals. This approximate score is used to select the most interesting database sequences for a subsequent Smith–Waterman alignment, which is also parallelised. The resulting method represents a substantial improvement compared to existing heuristics. The sensitivity and specificity of ParAlign was found to be as good as Smith–Waterman implementations when the same method for computing the statistical significance of the matches was used. In terms of speed, only the significantly less sensitive NCBI BLAST 2 program was found to outperform the new approach. Online searches are available at http://dna.uio.no/search/

Extending the Microtubule/Microfibril Paradigm1 : Cellulose Synthesis Is Required for Normal Cortical Microtubule Alignment in Elongating Cells

The cortical microtubule array provides spatial information to the cellulose-synthesizing machinery within the plasma membrane of elongating cells. Until now data indicated that information is transferred from organized cortical microtubules to the cellulose-synthesizing complex, which results in the deposition of ordered cellulosic walls. How cortical microtubules become aligned is unclear. The literature indicates that biophysical forces, transmitted by the organized cellulose component of the cell wall, provide a spatial cue to orient cortical microtubules. This hypothesis was tested on tobacco (Nicotiana tabacum L.) protoplasts and suspension-cultured cells treated with the cellulose synthesis inhibitor isoxaben. Isoxaben (0.25–2.5 μm) inhibited the synthesis of cellulose microfibrils (detected by staining with 1 μg mL−1 fluorescent dye and polarized birefringence), the cells failed to elongate, and the cortical microtubules failed to become organized. The affects of isoxaben were reversible, and after its removal microtubules reorganized and cells elongated. Isoxaben did not depolymerize microtubules in vivo or inhibit the polymerization of tubulin in vitro. These data are consistent with the hypothesis that cellulose microfibrils, and hence cell elongation, are involved in providing spatial cues for cortical microtubule organization. These results compel us to extend the microtubule/microfibril paradigm to include the bidirectional flow of information.

Multiple DNA and protein sequence alignment based on segment-to-segment comparison.

In this paper, a new way to think about, and to construct, pairwise as well as multiple alignments of DNA and protein sequences is proposed. Rather than forcing alignments to either align single residues or to introduce gaps by defining an alignment as a path running right from the source up to the sink in the associated dot-matrix diagram, we propose to consider alignments as consistent equivalence relations defined on the set of all positions occurring in all sequences under consideration. We also propose constructing alignments from whole segments exhibiting highly significant overall similarity rather than by aligning individual residues. Consequently, we present an alignment algorithm that (i) is based on segment-to-segment comparison instead of the commonly used residue-to-residue comparison and which (ii) avoids the well-known difficulties concerning the choice of appropriate gap penalties: gaps are not treated explicity, but remain as those parts of the sequences that do not belong to any of the aligned segments. Finally, we discuss the application of our algorithm to two test examples and compare it with commonly used alignment methods. As a first example, we aligned a set of 11 DNA sequences coding for functional helix-loop-helix proteins. Though the sequences show only low overall similarity, our program correctly aligned all of the 11 functional sites, which was a unique result among the methods tested. As a by-product, the reading frames of the sequences were identified. Next, we aligned a set of ribonuclease H proteins and compared our results with alignments produced by other programs as reported by McClure et al. [McClure, M. A., Vasi, T. K. & Fitch, W. M. (1994) Mol. Biol. Evol. 11, 571-592]. Our program was one of the best scoring programs. However, in contrast to other methods, our protein alignments are independent of user-defined parameters.

Functional cell surface expression of the anion transport domain of human red cell band 3 (AE1) in the yeast Saccharomyces cerevisiae.

We expressed the 52-kDa integral membrane domain (B3mem) of the human erythrocyte anion transporter (band 3; AE1) in a protease-deficient strain of the yeast Saccharomyces cerevisiae under the control of the inducible GAL10-CYC1 promoter. Immunoblots of total protein from transformed yeast cells confirmed that the B3mem polypeptide was overexpressed shortly after induction with galactose. Cell surface expression of the functional anion transporter was detected by using a simple transport assay to measure stilbene disulfonate-inhibitable chloride influx into intact yeast cells. The B3mem polypeptide was recycled and degraded by the cells with a half-life of approximately 1-3 hr, which led to a steady-state level of expression in exponentially growing cultures. Our data suggest that 5-10% of total B3mem is functionally active at the cell surface at any one time and that overexpression of this anion transport protein does not interfere with cell growth or survival. This is one of only a few reports of the functional expression of a plasma membrane transport protein in the plasma membrane of yeast cells and to our knowledge is the first report of red cell band 3-mediated anion transport at the plasma membrane of cDNA-transformed cells. The cell surface expression system we describe will provide a simple means for future study of the functional properties of band 3 by using site-directed mutagenesis.

Gene recognition via spliced sequence alignment.

Gene recognition is one of the most important problems in computational molecular biology. Previous attempts to solve this problem were based on statistics, and applications of combinatorial methods for gene recognition were almost unexplored. Recent advances in large-scale cDNA sequencing open a way toward a new approach to gene recognition that uses previously sequenced genes as a clue for recognition of newly sequenced genes. This paper describes a spliced alignment algorithm and software tool that explores all possible exon assemblies in polynomial time and finds the multiexon structure with the best fit to a related protein. Unlike other existing methods, the algorithm successfully recognizes genes even in the case of short exons or exons with unusual codon usage; we also report correct assemblies for genes with more than 10 exons. On a test sample of human genes with known mammalian relatives, the average correlation between the predicted and actual proteins was 99%. The algorithm correctly reconstructed 87% of genes and the rare discrepancies between the predicted and real exon-intron structures were caused either by short (less than 5 amino acids) initial/terminal exons or by alternative splicing. Moreover, the algorithm predicts human genes reasonably well when the homologous protein is nonvertebrate or even prokaryotic. The surprisingly good performance of the method was confirmed by extensive simulations: in particular, with target proteins at 160 accepted point mutations (PAM) (25% similarity), the correlation between the predicted and actual genes was still as high as 95%.