Contact-based sequence alignment

This paper introduces the novel method of contact-based protein sequence alignment, where structural information in the form of contact mutation probabilities is incorporated into an alignment routine using contact-mutation matrices (CAO: Contact Accepted mutatiOn). The contact-based alignment routine optimizes the score of matched contacts, which involves four (two per contact) instead of two residues per match in pairwise alignments. The first contact refers to a real side-chain contact in a template sequence with known structure, and the second contact is the equivalent putative contact of a homologous query sequence with unknown structure. An algorithm has been devised to perform a pairwise sequence alignment based on contact information. The contact scores were combined with PAM-type (Point Accepted Mutation) substitution scores after parameterization of gap penalties and score weights by means of a genetic algorithm. We show that owing to the structural information contained in the CAO matrices, significantly improved alignments of distantly related sequences can be obtained. This has allowed us to annotate eight putative Drosophila IGF sequences. Contact-based sequence alignment should therefore prove useful in comparative modelling and fold recognition.     

Immunocytochemical localization of Na+-HCO3- cotransporters and carbonic anhydrase dependence of fluid transport in corneal endothelial cells

In corneal endothelium, there is evidence for basolateral entry of HCO(3)(-) into corneal endothelial cells via Na(+)-HCO(3)(-) cotransporter (NBC) proteins and for net HCO(3)(-) flux from the basolateral to the apical side. However, how HCO(3)(-) exits the cells through the apical membrane is unclear. We determined that cultured corneal endothelial cells transport HCO(3)(-) similarly to fresh tissue. In addition, Cl(-) channel inhibitors decreased fluid transport by at most 16%, and inhibition of membrane-bound carbonic anhydrase IV by benzolamide or dextran-bound sulfonamide decreased fluid transport by at most 29%. Therefore, more than half of the fluid transport cannot be accounted for by anion transport through apical Cl(-) channels, CO(2) diffusion across the apical membrane, or a combination of these two mechanisms. However, immunocytochemistry using optical sectioning by confocal microscopy and cryosections revealed the presence of NBC transporters in both the basolateral and apical cell membranes of cultured bovine corneal endothelial cells and freshly isolated rabbit endothelia. This newly detected presence of an apical NBC transporter is consistent with its being the missing mechanism sought. We discuss discrepancies with other reports and provide a model that accounts for the experimental observations by assuming different stoichiometries of the NBC transport proteins at the basolateral and apical sides of the cells. Such functional differences might arise either from the expression of different isoforms or from regulatory factors affecting the stoichiometry of a single isoform.     

Multiple genetic processes result in heterogeneous rates of evolution within the major cluster disease resistance genes in lettuce

Resistance Gene Candidate2 (RGC2) genes belong to a large, highly duplicated family of nucleotide binding site-leucine rich repeat (NBS-LRR) encoding disease resistance genes located at a single locus in lettuce (Lactuca sativa). To investigate the genetic events occurring during the evolution of this locus, approximately 1.5- to 2-kb 3' fragments of 126 RGC2 genes from seven genotypes were sequenced from three species of Lactuca, and 107 additional RGC2 sequences were obtained from 40 wild accessions of Lactuca spp. The copy number of RGC2 genes varied from 12 to 32 per genome in the seven genotypes studied extensively. LRR number varied from 40 to 47; most of this variation had resulted from 13 events duplicating two to five LRRs because of unequal crossing-over within or between RGC2 genes at one of two recombination hot spots. Two types of RGC2 genes (Type I and Type II) were initially distinguished based on the pattern of sequence identities between their 3' regions. The existence of two types of RGC2 genes was further supported by intron similarities, the frequency of sequence exchange, and their prevalence in natural populations. Type I genes are extensive chimeras caused by frequent sequence exchanges. Frequent sequence exchanges between Type I genes homogenized intron sequences, but not coding sequences, and obscured allelic/orthologous relationships. Sequencing of Type I genes from additional wild accessions confirmed the high frequency of sequence exchange and the presence of numerous chimeric RGC2 genes in nature. Unlike Type I genes, Type II genes exhibited infrequent sequence exchange between paralogous sequences. Type II genes from different genotype/species within the genus Lactuca showed obvious allelic/orthologous relationships. Trans-specific polymorphism was observed for different groups of orthologs, suggesting balancing selection. Unequal crossover, insertion/deletion, and point mutation events were distributed unequally through the gene. Different evolutionary forces have impacted different parts of the LRR.     

Deterministic extinction effect of parasites on host populations

 Experimental studies have shown that parasites can reduce host density and even drive host population to extinction. Conventional mathematical models for parasite-host interactions, while can address the host density reduction scenario, fail to explain such deterministic extinction phenomena. In order to understand the parasite induced host extinction, Ebert et al. (2000) formulated a plausible but ad hoc epidemiological microparasite model and its stochastic variation. The deterministic model, resembles a simple SI type model, predicts the existence of a globally attractive positive steady state. Their simulation of the stochastic model indicates that extinction of host is a likely outcome in some parameter regions. A careful examination of their ad hoc model suggests an alternative and plausible model assumption. With this modification, we show that the revised parasite-host model can exhibit the observed parasite induced host extinction. This finding strengthens and complements that of Ebert et al. (2000), since all continuous models are likely break down when all population densities are small. This extinction dynamics resembles that of ratio-dependent predator-prey models. We report here a complete global study of the revised parasite-host model. Biological implications and limitations of our findings are also presented. Received: 30 October 2001 / Revised version: 11 February 2002 / Published online: 17 October 2002 Work is partially supported by NSF grant DMS-0077790 Mathematics Subject Classification (2000): 34C25, 34C35, 92D25. Keywords or phrases: Microparasite model – Ratio-dependent predator-prey model – Host extinction – Global stability – Biological control     

Coordinate expression of fibulin-5/DANCE and elastin during lung injury repair

Fibulin-5, previously known as DANCE and EVEC, is a secreted extracellular matrix protein that functions as a scaffold for elastin fiber assembly and as a ligand for integrins alphavbeta3, alphavbeta5, and alpha9beta1. Fibulin-5 is developmentally regulated in the lung, and lung air space enlargement develops in mice deficient in fibulin-5. Fibulin-5 is also induced in adult lung following lung injury by hyperoxia. To further examine the role of fibulin-5 during repair of lung injury, we assessed fibulin-5 expression during elastase-induced emphysema in C57/b mice. Mice were treated with either saline or elastase via the trachea, and the lung was examined 20 days after treatment. Fibulin-5 mRNA was induced almost fourfold, whereas elastin mRNA was minimally elevated. Immunohistochemistry studies showed that fibulin-5 was induced in cells within the alveolar wall following elastase treatment. Western analysis demonstrates that fibulin-5 was strongly expressed in isolated primary lung interstitial fibroblasts. Fibulin-5 protein was localized to the fibroblast cell layer in culture, and brief elastase treatment degraded the protein. Intact fibulin-5 did not accumulate in the culture media. Treatment of fibroblasts with the proinflammatory cytokine interleukin-1beta abolished fibulin-5 mRNA expression. Our results indicate that fibulin-5 is coordinately expressed and regulated with elastin in lung fibroblasts and may serve a key role during lung injury and repair.     

Amino acid encoding schemes from protein structure alignments: multi-dimensional vectors to describe residue types

Bioinformatic software has used various numerical encoding schemes to describe amino acid sequences. Orthogonal encoding, employing 20 numbers to describe the amino acid type of one protein residue, is often used with artificial neural network (ANN) models. However, this can increase the model complexity, thus leading to difficulty in implementation and poor performance. Here, we use ANNs to derive encoding schemes for the amino acid types from protein three-dimensional structure alignments. Each of the 20 amino acid types is characterized with a few real numbers. Our schemes are tested on the simulation of amino acid substitution matrices. These simplified schemes outperform the orthogonal encoding on small data sets. Using one of these encoding schemes, we generate a colouring scheme for the amino acids in which comparable amino acids are in similar colours. We expect it to be useful for visual inspection and manual editing of protein multiple sequence alignments.     

Reserve nutrient substance accumulation in Brassica rapa L. seeds in microgravity conditions (STS-87)

The results of study of Brassica embryo differentiation and reserve nutrient substance accumulation in the seeds were represented. Near resemblance of the spaceflight and around control embryo development was revealed. Different character of the reserve substance accumulation was noted, despite of the morphologic similarity in seeds produced in spaceflight and on the ground. It allows to consider spaceflight embryos morphologically more younger compared to the ground control.     

Conserved codon composition of ribosomal protein coding genes in Escherichia coli,Mycobacterium tuberculosis and Saccharomyces cerevisiae: lessons from supervised machine learning in functional genomics

Genomics projects have resulted in a flood of sequence data. Functional annotation currently relies almost exclusively on inter-species sequence comparison and is restricted in cases of limited data from related species and widely divergent sequences with no known homologs. Here, we demonstrate that codon composition, a fusion of codon usage bias and amino acid composition signals, can accurately discriminate, in the absence of sequence homology information, cytoplasmic ribosomal protein genes from all other genes of known function in Saccharomyces cerevisiae, Escherichia coli and Mycobacterium tuberculosis using an implementation of support vector machines, SVM(light). Analysis of these codon composition signals is instructive in determining features that confer individuality to ribosomal protein genes. Each of the sets of positively charged, negatively charged and small hydrophobic residues, as well as codon bias, contribute to their distinctive codon composition profile. The representation of all these signals is sensitively detected, combined and augmented by the SVMs to perform an accurate classification. Of special mention is an obvious outlier, yeast gene RPL22B, highly homologous to RPL22A but employing very different codon usage, perhaps indicating a non-ribosomal function. Finally, we propose that codon composition be used in combination with other attributes in gene/protein classification by supervised machine learning algorithms.     

Proteomics approach to identify wound-response related proteins from rice leaf sheath

In order to avoid the complex conditions of the intact plant for simple analysis of proteins in wound-response stress, we used the detached rice leaf sheath which is a very active part of the rice seedling. Proteins were extracted from rice leaf sheath at 0, 12, 24, 48 h after cutting and separated by two-dimensional (2-D) polyacrylamide gel electrophoresis. Changes in differentially displayed proteins were found in leaf sheaths after cutting in the 0-48 h time course. Ten proteins were up-regulated, while 19 proteins were down-regulated compared with those on the four 2-D gels. Among them, 14 proteins were analyzed by N-terminal, or internal amino acid sequence. The clear functions of nine proteins could be identified. Six proteins did not yield amino acid sequence information due to their blocked N-termini. Furthermore, 11 proteins were determined by matrix-assisted laser desorption/ionization-time of flight mass spectrometry, and identified protein database matching. It was shown that the down-regulated proteins were calreticulin (nos. 5, 6), histone H1 (no. 15) and hemoglobin (no. 17), putative peroxidase (no. 19); the up-regulated proteins were Bowman-Birk trypsin inhibitor (no. 23), putative receptor-like protein kinase (nos. 24, 25), calmodulin-related protein (no. 26), small subunit of ribulose-1,5-bisphosphate carboxylase/oxygenase (no. 27), mannose-binding rice lectin (nos. 28, 29). Among all the above proteins, four (nos. 23, 24, 25, 26) have been confirmed to be wound-response proteins. The others cannot be excluded as also being related to wound-responses, such as the signal transduction-related proteins (nos. 5, 6), photosynthesis-related protein (no. 27), and stress-response proteins (nos. 19, 28, 29). This is the first time protein changes in response to wounding in rice leaf sheath have been shown.