In vitro and in silico annotation of conserved and nonconserved microRNAs in the genome of the marsupial Monodelphis domestica

The gray short-tailed opossum, Monodelphis domestica, is the world's most widely utilized marsupial model for biomedical research. Recent completion of the initial M. domestica genome assembly offers the first opportunity to examine genome-wide phenomena in a marsupial. Using in silico methods, we have mapped 124 conserved microRNAs (miRNAs) to 94 loci in the M. domestica genome. In addition, using RNA pooled from 5 tissues, we cloned 85 miRNAs. Seventy-two of these are conserved miRNAs that we had mapped in silico. The additional 13 are nonconserved candidate miRNAs in 11 loci. Nine of these 13 are also found in the wallaby (Macropus eugenii) genome. Two of the candidate miRNA clones, located on the X chromosome, are part of a cluster containing a total of 24 potential miRNAs spanning more than 100 kb.     

In silico identification and characterization of stress and virulence associated repeats in Salmonella

So much genomic similarities yet causing different diseases, is like a paradox in Salmonella biology. Repeat is one of the probes that can explain such differences. Here, a comparative genomics approach is followed to identify and characterize repeats that might play role in adaptation and pathogenesis. Repeats are non-randomly distributed in the genomes except few typhoid causing strains. Perfect long repeats are rare compare to polymorphic ones and both are statistically consistent. Significant differences in repeat densities in stress related genes manifest its probable participation in survival and virulence. 573 and 1053 repeat loci have been identified which are exclusively associated with stress and virulent genes respectively. In Salmonella Typhi, an octameric VNTR locus is found in between acrD and yffB genes having more than 25 perfect copies across Salmonella Typhi but possesses only single copy in other serovars. This repeat can be used as a diagnostic probe for typhoid.     

In silico identification and biochemical characterization of novel inhibitors of NQO1

From in silico docking and COMPARE analysis, novel inhibitors of human NAD(P)H quinone oxidoreductase (NQO1) have been identified from the NCI compound database, the most potent of which has an observed IC₅₀ of 0.7 μM. The inhibitors exhibit a diverse range of scaffolds. The ability of docking calculations to predict experimentally determined binding affinities for NQO1 is discussed, considering the influence of target flexibility and scoring function.     

In silico cloning of evolutionarily conserved mouse F-LANa]

Differentially expressed genes between normal liver and hepatocellular carcinomas were investigated using differential display. In previous study, human F-LANa was identified as a differentially expressed gene, up-regulated in hepatocellular carcinoma. In this study, we developed an in silico cloning approach to rapidly and accurately characterize the mouse ortholog of the human F-LANa. Mouse F-LANa encodes a 239 aa protein exhibiting 97.9% similarity to the human ortholog gene. Homology analysis was carried out in various species and showed that F-LANa was evolutionarily conserved from yeast to human. Based on the alignment results, phylogenetic tree was established here.     

In silico characterization of thermostable lipases

Thermostable lipases are of high priority for industrial applications as they are endowed with the capability of carrying out diversified reactions at elevated temperatures. Extremophiles are their potential source. Sequence and structure annotation of thermostable lipases can elucidate evolution of lipases from their mesophilic counterparts with enhanced thermostability hence better industrial potential. Sequence analysis highlighted the conserved residues in bacterial and fungal thermostable lipases. Higher frequency of AXXXA motif and poly Ala residues in lid domain of thermostable Bacillus lipases were distinguishing characteristics. Comparison of amino acid composition among thermostable and mesostable lipases brought into light the role of neutral, charged and aromatic amino acid residues in enhancement of thermostability. Structural annotation of thermostable lipases with that of mesostable lipases revealed some striking features which are increment of gamma turns in thermostable lipases; being first time reported in our paper, longer beta strands, lesser beta-branched residues in helices, increase in charged-neutral hydrogen bonding pair, hydrophobic-hydrophobic contact and differences in the N-cap and C-cap residues of the α helices. Conclusively, it can be stated that subtle changes in the arrangement of amino acid residues in the tertiary structure of lipases contributes to enhanced thermostability.     

In silico identification of functional protein interfaces

Proteins perform many of their biological roles through protein-protein, protein-DNA or protein-ligand interfaces. The identification of the amino acids comprising these interfaces often enhances our understanding of the biological function of the proteins. Many methods for the detection of functional interfaces have been developed, and large-scale analyses have provided assessments of their accuracy. Among them are those that consider the size of the protein interface, its amino acid composition and its physicochemical and geometrical properties. Other methods to this effect use statistical potential functions of pairwise interactions, and evolutionary information. The rationale of the evolutionary approach is that functional and structural constraints impose selective pressure; hence, biologically important interfaces often evolve at a slower pace than do other external regions of the protein. Recently, an algorithm, Rate4Site, and a web-server, ConSurf (http://consurf.tau.ac.il/), for the identification of functional interfaces based on the evolutionary relations among homologous proteins as reflected in phylogenetic trees, were developed in our laboratory. The explicit use of the tree topology and branch lengths makes the method remarkably accurate and sensitive. Here we demonstrate its potency in the identification of the functional interfaces of a hypothetical protein, the structure of which was determined as part of the international structural genomics effort. Finally, we propose to combine complementary procedures, in order to enhance the overall performance of methods for the identification of functional interfaces in proteins.     

In silico identification and characterization of the MAPK family members of unicellular model eukaryote Tetrahymena thermophila

The biological function and evolutionary diversity of the mitogen-activated protein kinase (MAPK) family have mostly been studied in fungi, animals and plants, with very limited information from lower eukaryotes. This study aimed to describe the MAPKs of unicellular Tetrahymena thermophila. Eight members of the T. thermophila MAPK (TtMPK) gene family, in addition to previously reported TtMPK1, TtMPK2 and TtMPK3, were identified bioinformatically using a T. thermophila genome database. Phylogenetic analysis assigned the TtMPKs into two major groups, ERK1/2-like (TtMPK1, 2, 3, 5, 6, 7, 8, and 9) as stress-responsive MAPKs for biotic and abiotic stresses, and ERK7/8-like (TtMPK4, 10, and 11) as cell-cycle-associated protein kinases for biotic factors. Semi-quantitative RT-PCR analysis of the TtMPKs showed high mRNA expression at 30 °C; however, only TtMPK5 and TtMPK6 showed high expression at 37 °C. Osmotic shock by 100 mM NaCl only increased the expression of TtMPK2, whereas 20 mM NaCl reduced the expression of all MPKs to almost zero. The results suggested that T. thermophila MAPKs are among the closest representatives of the ancestors of the eukaryotic MAPK family. Although no functional characterization of MPKs was performed, this study is the first report of the genome-wide MAPK family in T. thermophila.     

In silico identification of putative metal binding motifs

Metal ion binding domains are found in proteins that mediate transport, buffering or detoxification of metal ions. In this study, we have performed an in silico analysis of metal binding proteins and have identified putative metal binding motifs for the ions of cadmium, cobalt, zinc, arsenic, mercury, magnesium, manganese, molybdenum and nickel. A pattern search against the UniProtKB/Swiss-Prot and UniProtKB/TrEMBL databases yielded true positives in each case showing the high-specificity of the motifs. Motifs were also validated against PDB structures and site directed mutagenesis studies.     

In silico identification and computational characterization of endogenous small interfering RNAs from diverse grapevine tissues and stages

Small interfering RNAs (siRNAs) are effectors of regulatory pathways underlying plant development, metabolism, and stress- and nutrient-signaling regulatory networks. The endogenous siRNAs are generally not conserved between plants; consequently, it is necessary and important to identify and characterize siRNAs from various plants. To address the nature and functions of siRNAs, and understand the biological roles of the huge siRNA population in grapevine (Vitis vinifera L.). The high-throughput sequencing technology was used to identify a large set of putative endogenous siRNAs from six grapevine tissues/organs. Subsequently, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed to classify the target genes of siRNA. In total, 520,519 candidate siRNAs were identified and their expression profiles exhibited typical temporal characters during grapevine development. In addition, we identified two grapevine trans-acting siRNA (TAS) gene homologs (VvTAS3 and VvTAS4) and the derived trans-acting siRNAs (tasiRNAs) that could target grapevine auxin response factor (ARF) and myeloblastosis (MYB) genes. Furthermore, the GO and KEGG analysis of target genes showed that most of them covered a broad range of functional categories, especially involving in disease-resistance process. The large-scale and completely genome-wide level identification and characterization of grapevine endogenous siRNAs from the diverse tissues by high throughput technology revealed the nature and functions of siRNAs in grapevine.