A comprehensive software suite for the analysis of cDNAs

We have developed a comprehensive software suite for bioinformatics research of cDNAs; it is aimed at rapid characterization of the features of genes and the proteins they code. Methods implemented include the detection of translation initia- tion and termination signals, statistical analysis of codon usage, comparative study of amino acid composition, comparative modeling of the structures of product proteins, prediction of alternative splice forms, and metabolic pathway reconstruction. The software package is freely available under the GNU General Public License at http： / /www.g-language.org/ data/cdna/.     

The genome of the diploid anuran <Emphasis Type="Italic">Xenopus tropicalis</Emphasis> contains a novel array of sarcoplasmic myosin heavy chain genes expressed in larval muscle and larynx

The sarcomeric myosin heavy chain (MyHC) proteins are a family of molecular motors responsible for the transduction of chemical energy into mechanical work in striated muscle. The vertebrate genome contains multiple copies of the MyHC gene, and expression of different isoforms correlates with differences in the physiological properties of muscle fibers. Most MyHC isoforms are found in two arrays, one containing the "fast-twitch" skeletal muscle isoforms and the other the "slow-twitch" or cardiac isoforms. To extend our understanding of MyHC evolution, we have examined the genome of the anuran Xenopus tropicalis. The X. tropicalis genome includes15 full-length MyHC genes organized in seven genomic locations. One unique array of MyHC genes is similar to the mammalian fast-skeletal array, but is not found in amniotes. The isoforms in this array are expressed during larval stages and in muscles of the adult larynx. Duplication of the fast-skeletal MyHC array appears to have led to expression divergence of muscle proteins in the larval and adult stages of the anuran life cycle. A striking similarity of gene order between regions flanking X. tropicalis MyHC arrays and human arrays was evident; genomic organization of MyHC isoforms may thus be highly conserved across tetrapods.     

Computational identification of novel chitinase-like proteins in the Drosophila melanogaster genome

MOTIVATION: Multiple chitinases as well as lectins closely related to them have been characterized previously from many insect species and the corresponding genes/cDNAs have been cloned. However, the identification of the entire assortment of genes for chitinase family proteins and their differences in biochemical properties have not been carried out in any individual insect species. The completion of the entire DNA sequence of Drosophila melanogaster (fruit fly) genome and identification of open reading frames presents an opportunity to study the structures and functions of chitinase-like proteins, and also to identify new members of this family in DROSOPHILA: We are, therefore, interested in studying the functional genomics of chitinase-like gene families in insects. METHODS: We searched the Drosophila protein sequences database using fully characterized insect chitinase sequences and BLASTP software, identified all the putative chitinase-like proteins encoded in Drosophila genome, and predicted their structures using domain analysis tools. A phylogenetic analysis of the chitinase-like proteins from Drosophila and several other insect species was carried out. The structures of these chitinases were modeled using homology modeling software. RESULTS: Our analysis revealed the presence of 18 chitinase-like proteins in the Drosophila protein database. Among these are seven novel chitinase-like proteins that contain four signature amino acid sequences of chitinases belonging to family 18 glycosylhydrolases, including both acidic and hydrophobic amino acid residues critical for enzyme activity. All the proteins contain at least one catalytic domain with one having four catalytic domains. Phylogenetic analysis of chitinase-like proteins from Drosophila and other insects revealed an evolutionary relationship among all these proteins, which indicated gene duplication and domain shuffling to generate the observed diversity in the encoded proteins. Homology modeling showed that all the Drosophila chitinase-like proteins contain one or more catalytic domains with a (alpha/beta)8 barrel-like structure. Our results suggest that insects utilize multiple family 18 chitinolytic enzymes and also non-enzymatic chitinase-like proteins for degrading/remodeling/binding to chitin in different insect anatomical extracellular structures, such as the cuticle, peritrophic membrane, trachea and mouth parts during insect development, and possibly for other roles including chitin synthesis. AVAILABILITY: Perl program and supplementary material are available at http://www.ksu.edu/bioinformatics/supplementary.htm     

Voltage sensitive phosphoinositide phosphatases of Xenopus: their tissue distribution and voltage dependence

Voltage-sensitive phosphatases (VSPs) are unique proteins in which membrane potential controls enzyme activity. They are comprised of the voltage sensor domain of an ion channel coupled to a lipid phosphatase specific for phosphoinositides, and for ascidian and zebrafish VSPs, the phosphatase activity has been found to be activated by membrane depolarization. The physiological functions of these proteins are unknown, but their expression in testis and embryos suggests a role in fertilization or development. Here we investigate the expression pattern and voltage dependence of VSPs in two frog species, Xenopus laevis and Xenopus tropicalis, that are well suited for experimental studies of these possible functions. X. laevis has two VSP genes (Xl-VSP1 and Xl-VSP2), whereas X. tropicalis has only one gene (Xt-VSP). The highest expression of these genes was observed in testis, ovary, liver, and kidney. Our results show that while Xl-VSP2 activates only at positive membrane potentials outside of the physiological range, Xl-VSP1 and Xt-VSP phosphatase activity is regulated in the voltage range that regulates sperm-egg fusion at fertilization.     

Three members of the human pyruvate dehydrogenase kinase gene family are direct targets of the peroxisome proliferator-activated receptor beta/delta

The nuclear receptors peroxisome proliferator-activated receptors (PPARs) are known for their critical role in the metabolic syndrome. Here, we show that they are direct regulators of the family of pyruvate dehydrogenase kinase (PDK) genes, whose products act as metabolic homeostats in sensing hunger and satiety levels in key metabolic tissues by modulating the activity of the pyruvate dehydrogenase complex. Mis-regulation of this tightly controlled network may lead to hyperglycemia. In human embryonal kidney cells we found the mRNA expression of PDK2, PDK3 and PDK4 to be under direct primary control of PPAR ligands, and in normal mouse kidney tissue Pdk2 and Pdk4 are PPAR targets. Both, treatment of HEK cells with PPARbeta/delta-specific siRNA and the genetic disruption of the Pparbeta/delta gene in mouse fibroblasts resulted in reduced expression of Pdk genes and abolition of induction by PPARbeta/delta ligands. These findings suggest that PPARbeta/delta is a key regulator of PDK genes, in particular the PDK4/Pdk4 gene. In silico analysis of the human PDK genes revealed two candidate PPAR response elements in the PDK2 gene, five in the PDK3 gene and two in the PDK4 gene, but none in the PDK1 gene. For seven of these sites we could demonstrate both PPARbeta/delta ligand responsiveness in context of their chromatin region and simultaneous association of PPARbeta/delta with its functional partner proteins, such as retinoidXreceptor, co-activator and mediator proteins and phosphorylated RNA polymerase II. In conclusion, PDK2, PDK3 and PDK4 are primary PPARbeta/delta target genes in humans underlining the importance of the receptor in the control of metabolism.     

基于SwissPdbViewer和MATLAB对东北虎Sox蛋白三级结构建模及分析

SOX蛋白具有一个与DNA特异结合的高保守HMG-box结构域。为研究东北虎SOX蛋白三级结构的分子机理,利用MATLAB的Bioinformatics工具从GenBank中下载东北虎SOX蛋白序列信息,以三级结构已知的SOX2为模板,联合SwissPdbViewer与MATLAB,采用同源建模方法对SOX蛋白HMG-box进行建模、预测;利用MATLAB的Visualization Tool分析预测结果的三维结构。结果显示PtSox蛋白的HMG-box由3个α-螺旋和2个loop区构成;热稳定性分析表明PtSox蛋白loop区的热力学结构不稳定;表面静电分布显示出PtSox蛋白C-端的中间有一个可能与其它小分子或蛋白质的相互作用位点的N/C腔,上述空间结构可能与其活性与功能的调控有关。     

Molecular modeling approaches for determining gene function: application to a putative poly-A binding protein from Leishmania amazonensis (LaPABP)

The great expansion in the number of genome sequencing projects has revealed the importance of computational methods to speed up the characterization of unknown genes. These studies have been improved by the use of three dimensional information from the predicted proteins generated by molecular modeling techniques. In this work, we disclose the structure-function relationship of a gene product from Leishmania amazonensis by applying molecular modeling and bioinformatics techniques. The analyzed sequence encodes a 159 amino acids polypeptide (estimated 18 kDa) and was denoted LaPABP for its high homology with poly-A binding proteins from trypanosomatids. The domain structure, clustering analysis and a three dimensional model of LaPABP, basically obtained by homology modeling on the structure of the human poly-A binding protein, are described. Based on the analysis of the electrostatic potential mapped on the model's surface and conservation of intramolecular contacts responsible for folding stabilization we hypothesize that this protein may have less avidity to RNA than it's L. major counterpart but still account for a significant functional activity in the parasite. The model obtained will help in the design of mutagenesis experiments aimed to elucidate the mechanism of gene expression in trypanosomatids and serve as a starting point for its exploration as a potential source of targets for a rational chemotherapy.     

Multiple roles of phosphoinositide-specific phospholipase C isozymes

Phosphoinositide-specific phospholipase C is an effector molecule in the signal transduction process. It generates two second messengers, inositol-1,4,5-trisphosphate and diacylglycerol from phosphatidylinositol 4,5-bisphosphate. Currently, thirteen mammal PLC isozymes have been identified, and they are divided into six groups: PLC-beta, -gamma, -delta, -epsilon, -zeta and -eta. Sequence analysis studies demonstrated that each isozyme has more than one alternative splicing variant. PLC isozymes contain the X and Y domains that are responsible for catalytic activity. Several other domains including the PH domain, the C2 domain and EF hand motifs are involved in various biological functions of PLC isozymes as signaling proteins. The distribution of PLC isozymes is tissue and organ specific. Recent studies on isolated cells and knockout mice depleted of PLC isozymes have revealed their distinct phenotypes. Given the specificity in distribution and cellular localization, it is clear that each PLC isozyme bears a unique function in the modulation of physiological responses. In this review, we discuss the structural organization, enzymatic properties and molecular diversity of PLC splicing variants and study functional and physiological roles of each isozyme.     

Cytochrome c: gene structure, homology and ancestral relationships.

In this paper, the author notes the recommended definition of the word "homology" (i.e., indicating an ancestral relationship) and the recommended stipulation that "evidence for homology should be explicitly laid out". The postulated homology for somatic and testes-specific isozymes of cytochrome c is then examined, using recent data obtained from the study of cytochrome c genes. Consideration is also given to some newer findings of molecular biology and possibilities are considered for various types of change in the genome of an organism. Possible roles of introns, pseudogenes and multigene families are considered. The relationship of testes-specific cytochrome c to somatic cytochrome c is carefully considered from data obtained in experimental studies of genes of these two isozymes. If one assumes that these isozymes arose as a consequence of a gene duplication, data from rat and mouse genes indicate that the testes-specific isozyme has incorporated more amino acid changes than the somatic isozyme since the time of their divergence. However, when the 15 amino acid differences (testes-specific vs. somatic isozyme) are considered, there is virtually no similarity in these 15 positions of the testes-specific isozyme with any of the hypothetical ancestral sequences of the somatic isozyme. Nucleotide differences in cytochrome c genes have been evaluated by comparing genes for the two rodent cytochrome c isozymes to cytochrome c genes of fruit flies, chickens and humans. Comparisons of nucleotide substitution rates in genes for the two cytochrome c isozymes in rodents confirm the conclusions from amino acid sequence comparisons; namely, that more rapid nucleotide changes have occurred in the testes-specific cytochrome c gene, than in the somatic cytochrome c gene. Possible explanations for these findings are considered.     

Structural features that govern enzymatic activity in carbonic anhydrase from a low-temperature adapted fish, Chionodraco hamatus

The carbonic anhydrase (CA) family of zinc metalloenzymes includes many known isozymes that have different subcellular distributions. The study described here focuses on identification of the structural features that define low-temperature adaptation in a Chionodraco hamatus protein, both for the reaction center, at an atomic level, and for the tertiary structure of the protein. To this aim, an x-ray absorption near-edge spectroscopy/Minuit x-ray absorption near-edge spectroscopy analysis of the reaction center was undertaken for both a structurally characterized human CAII and CA of C. hamatus. Higher structural levels were analyzed by sequence comparison and homology modeling. To establish whether the structural insights acquired in fish CAs are general, theoretical models were generated by homology modeling for three temperate-climate-adapted fish CAs. The measured structural differences between the two proteins are discussed in terms of the differences in the electrostatic potential between human CAII and CA of C. hamatus. We conclude that modulation of the interaction between the catalytic water molecule and the zinc ion could depend on the effect of the electrostatic potential distribution.     

A Xenopus tropicalis oligonucleotide microarray works across species using RNA from Xenopus laevis

Microarrays have great potential for the study of developmental biology. As a model system Xenopus is well suited for making the most of this potential. However, Xenopus laevis has undergone a genome wide duplication meaning that most genes are represented by two paralogues. This causes a number of problems. Most importantly the presence of duplicated genes mean that a X. laevis microarray will have less or even half the coverage of a similar sized microarray from the closely related but diploid frog Xenopus tropicalis. However, to date, X. laevis is the most commonly used amphibian system for experimental embryology. Therefore, we have tested if a microarray based on sequences from X. tropicalis will work across species using RNA from X. laevis. We produced a pilot oligonucleotide microarray based on sequences from X. tropicalis. The microarray was used to identify genes whose expression levels changed during early X. tropicalis development. The same assay was then carried out using RNA from X. laevis. The cross species experiments gave similar results to those using X. tropicalis RNA. This was true at the whole microarray level and for individual genes, with most genes giving similar results using RNA from X. laevis and X. tropicalis. Furthermore, the overlap in genes identified between a X. laevis and a X. tropicalis set of experiments was only 12% less than the overlap between two sets of X. tropicalis experiments. Therefore researchers can work with X. laevis and still make use of the advantages offered by X. tropicalis microarrays.     

Dbmodeling

Genome sequencing efforts are providing us with complete genetic blueprints for hundreds of organisms. We are now faced with assigning, understanding, and modifying the functions of proteins encoded by these genomes. DBMODELING is a relational database of annotated comparative protein structure models and their metabolic pathway characterization, when identified. This procedure was applied to complete genomes such as Mycobacterium tuberculosis and Xylella fastidiosa. The main interest in the study of metabolic pathways is that some of these pathways are not present in humans, which makes them selective targets for drug design, decreasing the impact of drugs in humans. In the database, there are currently 1116 proteins from two genomes. It can be accessed by any researcher at http://www. biocristalografia.df.ibilce. unesp.br/tools/. This project confirms that homology modeling is a useful tool in structural bioinformatics and that it can be very valuable in annotating genome sequence information, contributing to structural and functional genomics, and analyzing protein-ligand docking.     

Combining structural genomics and enzymology: completing the picture in metabolic pathways and enzyme active sites

An important goal of structural genomics is to complete the structural analysis of all the enzymes in metabolic pathways and to understand the structural similarities and differences. A preliminary glimpse of this type of analysis was achieved before structural genomics efforts with the glycolytic pathway and efforts are underway for many other pathways, including that of catecholamine metabolism. Structural enzymology necessitates a complete structural characterization, even for highly homologous proteins (greater than 80% sequence homology), as every active site has distinct structural features and it is these active site differences that distinguish one enzyme from another. Short cuts with homology modeling cannot be taken with our current knowledge base. Each enzyme structure in a pathway needs to be determined, including structures containing bound substrates, cofactors, products and transition state analogs, in order to obtain a complete structural and functional understanding of pathway-related enzymes.     

Complementarity of conserved sequence elements present in 28S ribosomal RNA and in ribosomal protein genes of Xenopus laevis and Xenopus tropicalis

The sequence analysis of the L1 ribosomal protein (r-protein) gene of Xenopus laevis has revealed a strong homology in four out of the nine introns of the gene; this homology region spans 60 nucleotides (nt) with 80% homology [Loreni et al., EMBO J. 4 (1985) 3483-3488]. We have extended our analysis to X. tropicalis, a species which is closely related to X. laevis. Partial sequencing of the isolated L1 gene has revealed that these 60-nt homology regions are also present in at least two introns of the X. tropicalis L1 gene. Computer analysis has revealed that perfect nt sequence complementarity exists between 13 nt of this intron region and the 28S ribosomal RNA in a region which is conserved in all eukaryotes, suggesting a possible base-pairing interaction between these two sequences.     

Protein models: The Grand Challenge of protein docking

Characterization of life processes at the molecular level requires structural details of protein–protein interactions (PPIs). The number of experimentally determined protein structures accounts only for a fraction of known proteins. This gap has to be bridged by modeling, typically using experimentally determined structures as templates to model related proteins. The fraction of experimentally determined PPI structures is even smaller than that for the individual proteins, due to a larger number of interactions than the number of individual proteins, and a greater difficulty of crystallizing protein–protein complexes. The approaches to structural modeling of PPI (docking) often have to rely on modeled structures of the interactors, especially in the case of large PPI networks. Structures of modeled proteins are typically less accurate than the ones determined by X‐ray crystallography or nuclear magnetic resonance. Thus the utility of approaches to dock these structures should be assessed by thorough benchmarking, specifically designed for protein models. To be credible, such benchmarking has to be based on carefully curated sets of structures with levels of distortion typical for modeled proteins. This article presents such a suite of models built for the benchmark set of the X‐ray structures from the Dockground resource ( http://dockground.bioinformatics.ku.edu ) by a combination of homology modeling and Nudged Elastic Band method. For each monomer, six models were generated with predefined C^α root mean square deviation from the native structure (1, 2, …, 6 Å). The sets and the accompanying data provide a comprehensive resource for the development of docking methodology for modeled proteins. Proteins 2014; 82:278–287. © 2013 Wiley Periodicals, Inc.     

Acid phosphatase isozymes secreted under phosphate-deficient conditions inPholiota nameko

We previously reported the purification of an acid phosphatase (APase52) secreted from the mycelia ofPholiota nameko under phosphate-deficient conditions. In the present study, two other isozymes (APase47 and APase48) were found and their structures were compared with that of APase52. Thirteen amino acid residues at theN-terminus of APase47 were completely identical with those of APase48 and had partial homology with those of APase52. The deglycosylation of proteins indicated that three APase isozymes differ in theN-linked oligosaccharide content. The protease-generated peptide maps of the APases differed from one another in the band pattern. These results suggest that the APases are the products of different genes.     

Expression of Xenopus tropicalis noggin1 and noggin2 in early development: two noggin genes in a tetrapod

We report the identification of two distinct noggin genes in the tetrapod Xenopus tropicalis. Noggin functions to antagonize BMP signaling in many developmental contexts, and much work has explored its role in early vertebrate development. We have identified two noggin genes in the tropical clawed frog, X. tropicalis, a diploid anuran which is being explored for its potential as a genetic model system for early vertebrate development. Here we report the cloning and characterization of the Xenopus tropicalis noggin1 and noggin2 genes, which have distinct expression domains in the early embryo with one overlapping domain in the anterior neural tissue. X. tropicalis noggin1 expression is very similar to that of noggin in Xenopus laevis, with expression beginning in the blastula organizer region and continuing through gastrulation and neurulation in the organizer and notochord. Later, it is also expressed in the anterior neural ridge and subsequent forebrain; noggin1 is also expressed in the pharyngeal arches after neural tube closure. At the tadpole stage expression is maintained in the dorsal neural tube and is present in the otic vesicle. However, the expression of noggin2 is much more similar to the expression of noggin2 in D. rerio with expression in the forebrain, hindbrain, and somites, but unlike D. rerio, X. tropicalis noggin2 is expressed in the heart by stage 28. This work presents the first example of a tetrapod with at least two noggin genes.     

低氧环境对三阴性乳腺癌细胞糖代谢途径的重编程

代谢改变是癌细胞的特征之一。研究表明,低氧会使癌细胞的糖代谢发生改变,但是更详细的分子机制仍有待进一步研究。本研究利用转录物组测序技术(RNA-sequencing,RNA-seq)和生物信息学分析发现,低氧导致BT549细胞中334个基因和MDA-MB-231细胞中215个基因在转录水平的表达改变。这些表达变化的基因多与糖代谢相关。进一步分析RNA-seq数据并应用Western 印迹、酶活性检测和代谢产物定量测定的结果显示,低氧通过升高BT549细胞中葡萄糖转运蛋白1(GLUT1)和MDA-MB-231细胞中GLUT1和GLUT3的表达以增加葡萄糖的摄入;低氧使催化糖的无氧氧化途径几乎全部反应的酶都至少有一种同工酶或酶蛋白亚基,以及调节酶6-磷酸果糖-2-激酶/果糖-2,6-二磷酸酶3(PFKFB3)和4(PFKFB4)同工酶的表达增加来促进了糖的无氧氧化;低氧还通过增加调节丙酮酸脱氢酶激酶1(PDK1)和3(PDK3)同工酶基因的表达,以及降低关键酶异柠檬酸脱氢酶3(IDH3)同工酶、琥珀酸脱氢酶B亚基和D亚基的表达来减少糖的有氧氧化途径进行;低氧可能还增加磷酸戊糖途径的关键酶葡糖-6-磷酸脱氢酶、糖原合成途径的关键酶糖原合酶GYS1同工酶的表达以促进这2条途径的进行,而对糖异生和糖原分解代谢途径酶基因的表达影响较小。生物信息学分析乳腺癌组织样本在线数据库中糖代谢途径酶基因在转录水平表达结果与细胞研究结果基本一致。总之,该文系统分析了低氧对糖代谢6条代谢途径中全部酶以及2种重要调节酶的影响,可见低氧会通过改变这些酶的同工酶或亚基的基因表达使糖代谢途径进行重编程,这对进一步认识低氧环境下癌细胞糖代谢的分子机制具有一定的意义。