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.     

Ion transit pathways and gating in ClC chloride channels

ClC chloride channels possess a homodimeric structure in which each monomer contains an independent chloride ion pathway. ClC channel gating is regulated by chloride ion concentration, pH and voltage. Based on structural and physiological evidence, it has been proposed that a glutamate residue on the extracellular end of the selectivity filter acts as a fast gate. We utilized a new search algorithm that incorporates electrostatic information to explore the ion transit pathways through wild-type and mutant bacterial ClC channels. Examination of the chloride ion permeation pathways supports the importance of the glutamate residue in gating. An external chloride binding site previously postulated in physiological experiments is located near a conserved basic residue adjacent to the gate. In addition, access pathways are found for proton migration to the gate, enabling pH control at hyperpolarized membrane potentials. A chloride ion in the selectivity filter is required for the pH-dependent gating mechanism.     

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.     

Photosynthetic characteristics and tolerance to photo-oxidation of transgenic rice expressing C4 photosynthesis enzymes

The photosynthetic characteristics of four transgenic rice lines over-expressing rice NADP-malic enzyme (ME), and maize phosphoenolpyruvate carboxylase (PC), pyruvate,orthophosphate dikinase (PK), and PC+PK (CK) were investigated using outdoor-grown plants. Relative to untransformed wild-type (WT) rice, PC transgenic rice exhibited high PC activity (25-fold increase) and enhanced activity of carbonic anhydrase (more than two-fold increase), while the activity of ribulose-bisphosphate carboxylase/oxygenase (Rubisco) and its kinetic property were not significantly altered. The PC transgenic plants also showed a higher light intensity for saturation of photosynthesis, higher photosynthetic CO₂ uptake rate and carboxylation efficiency, and slightly reduced CO₂ compensation point. In addition, chlorophyll a fluorescence analysis indicates that PC transgenic plants are more tolerant to photo-oxidative stress, due to a higher capacity to quench excess light energy via photochemical and non-photochemical means. Furthermore, PC and CK transgenic rice produced 22–24% more grains than WT plants. Taken together, these results suggest that expression of maize C₄ photosynthesis enzymes in rice, a C₃ plant, can improve its photosynthetic capacity with enhanced tolerance to photo-oxidation. This revised version was published online in June 2006 with corrections to the Cover Date.     

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.