Probabilistic representation of gene regulatory networks

MOTIVATION: Recent experiments have established unambiguously that biological systems can have significant cell-to-cell variations in gene expression levels even in isogenic populations. Computational approaches to studying gene expression in cellular systems should capture such biological variations for a more realistic representation. RESULTS: In this paper, we present a new fully probabilistic approach to the modeling of gene regulatory networks that allows for fluctuations in the gene expression levels. The new algorithm uses a very simple representation for the genes, and accounts for the repression or induction of the genes and for the biological variations among isogenic populations simultaneously. Because of its simplicity, introduced algorithm is a very promising approach to model large-scale gene regulatory networks. We have tested the new algorithm on the synthetic gene network library bioengineered recently. The good agreement between the computed and the experimental results for this library of networks, and additional tests, demonstrate that the new algorithm is robust and very successful in explaining the experimental data. AVAILABILITY: The simulation software is available upon request. SUPPLEMENTARY INFORMATION: Supplementary material will be made available on the OUP server.     

An efficient method for the physical mapping of transgenes in barley using in situ hybridization.

The genetic transformation of crops by particle bombardment and Agrobacterium tumefaciens systems have the potential to complement conventional plant breeding programmes. However, before deployment, transgenic plants need to be characterized in detail, and physical mapping is an integral part of this process. Therefore, it is important to have a highly efficient method for transgene detection by fluorescence in situ hybridization (FISH). This study describes a new approach, which provides efficient control of probe length and labelling, both of which play an important role in in situ hybridization of transgenes. The approach is based on reducing the size of the plasmid prior to labelling by nick translation, rather than using the whole or linearized plasmid, or varying the amounts of DNaseI in the nick translation mixture. This provided much more efficient labelling of the probe, which yielded optimal hybridization. minimal fluorescent background, and accurate physical location of the transgene.     

Functional evaluation of domain-domain interactions and human protein interaction networks

MOTIVATION: Large amounts of protein and domain interaction data are being produced by experimental high-throughput techniques and computational approaches. To gain insight into the value of the provided data, we used our new similarity measure based on the Gene Ontology (GO) to evaluate the molecular functions and biological processes of interacting proteins or domains. The applied measure particularly addresses the frequent annotation of proteins or domains with multiple GO terms. RESULTS: Using our similarity measure, we compare predicted domain-domain and human protein-protein interactions with experimentally derived interactions. The results show that our similarity measure is of significant benefit in quality assessment and confidence ranking of domain and protein networks. We also derive useful confidence score thresholds for dividing domain interaction predictions into subsets of low and high confidence. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

OMWSA: detection of DNA repeats using moving window spectral analysis

Repetitive DNA sequences play paramount biological roles, such as gene variation and regulatory functions on gene expressions. Until now, detection of various kinds of DNA repeats accurately is still an open problem. In this article, we propose a new method and a visualization tool for detecting DNA repeats in a 2D plane of location and frequency by using optimized moving window spectral analysis. The spectrogram can display the general distribution of repetitive sequences while showing the repeat period, length and location without any prior knowledge. Experimental results demonstrate that our method is accurate and robust even under the condition of excessive mutating and interleaving. AVAILABILITY: Available on http://www.hy8.com/~tec/sw01/omwsa01.zip. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Influence of transgenesis on genome variability in cucumber lines with a thaumatin II gene

The development of new plant varieties by genetic modification aims at improving their features or introducing new qualities. However, concerns about the unintended effects of transgenes and negative environmental impact of genetically modified plants are an obstacle for the use of these plants in crops. To analyze the impact of transgenesis on plant genomes, we analyze three cucumber transgenic lines with an introduced thaumatin II gene. After genomes sequencing, we analyzed the transgene insertion site and performed variant prediction. As a result, we obtained similar number of variants for all analyzed lines (average of 4307 polymorphisms), with high abundance in one region of chromosome 4. According to SnpEff analysis, the presence of genomic variants generally does not influence the genome functionality, as less than 2% of polymorphisms have high impact. Moreover, analysis indicates that these changes were more likely induced by in vitro culture than by the transgenesis itself. The insertion site analysis shows that the region of transgene integration could cause changes in gene expression, by gene disruption or loss of promoter region continuity.Supplementary InformationThe online version contains supplementary material available at 10.1007/s12298-021-00990-8.     

Comparison of the accuracy of semiempirical and some DFT methods for predicting heats of formation

A comparison is made of the relative accuracy of some NDDO semiempirical methods and the DFT functionals LYP and PW91 using both double and triple zeta basis sets. The comparison is between the calculated heat of formation and that reported in the NIST database.Electronic Supplementary Material Parameters for the tailored method are available in the supplementary material; these are suitable for use with MOPAC2002. All raw data (experimental heats of formation, geometries, total energies and heats of formation for the various methods, etc.) are also provided in CAChe format. Supplementary material is available for this article if you access the article at . A link in the frame on the left on that page takes you directly to the supplementary material.     

InFiRe -- a novel computational method for the identification of insertion sites in transposon mutagenized bacterial genomes

MOTIVATION: InFiRe, Insertion Finder via Restriction digest, is a novel software tool that allows for the computational identification of transposon insertion sites in known bacterial genome sequences after transposon mutagenesis experiments. The approach is based on the fact that restriction endonuclease digestions of bacterial DNA yield a unique pattern of DNA fragments with defined sizes. Transposon insertion changes the size of the hosting DNA fragment by a known number of base pairs. The exact size of this fragment can be determined by Southern blot hybridization. Subsequently, the position of insertion can be identified with computational analysis. The outlined method provides a solid basis for the establishment of a new high-throughput technology. AVAILABILITY AND IMPLEMENTATION: The software is freely available on our web server at www.infire.tu-bs.de. The algorithm was implemented in the statistical programming language R. For the most flexible use, InFiRe is provided in two different versions. A web interface offers the convenient use in a web browser. In addition, the software and source code is freely available for download as R-packages on our website. CONTACT: m.steinert@tu-bs.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

<Emphasis Type="Italic">Drosophila</Emphasis>-based <Emphasis Type="Italic">in vivo</Emphasis> assay for the validation of inhibitors of the epidermal growth factor receptor/Ras pathway

Overexpression of epidermal growth factor receptor (EGFR) is a common phenomenon observed in most cancers. Clinical treatment of such cancer involves the use of chemotherapeutic agents such as gefitinib and erlotinib which are inhibitors of tyrosine kinase (TK). These small molecules bind to the ATP-binding sites of the TK domain of EGFR. Our in silico analysis suggests that the TK domains of Drosophila and human EGFR are highly conserved. We therefore employed the Drosophila system to validate the in silico observations made with two important anticancer drugs. Since a large number of mutant flies are available, it was possible to investigate the various components of the EGFR/Ras/Raf/ MAPK pathways and the phosphorylation status of diphosphorylated extracellular signal-regulated kinase (dp-ERK1/2). These studies confirm the binding of the anilinoquinazolines to the Drosophila EGFR protein and modulation of its activity. Thus, Drosophila appears to be a robust and simple model system for screening newer anticancer drugs that act as TK inhibitors (TKIs). Electronic Supplementary Material  Supplementary material is available for this article at and is accessible for authorized users. Supplementary figures and tables pertaining to this article are available on the Journal of Biosciences Website at http://www.ias.ac.in/jbiosci/dec2008/pp731-742-suppl.pdf     

CRONOS: the cross-reference navigation server

Summary: Cross-mapping of gene and protein identifiers betweendifferent databases is a tedious and time-consuming task. Toovercome this, we developed CRONOS, a cross-reference serverthat contains entries from five mammalian organisms presentedby major gene and protein information resources. Sequence similarityanalysis of the mapped entries shows that the cross-referencesare highly accurate. In total, up to 18 different identifiertypes can be used for identification of cross-references. Thequality of the mapping could be improved substantially by exclusionof ambiguous gene and protein names which were manually validated.Organism-specific lists of ambiguous terms, which are valuablefor a variety of bioinformatics applications like text miningare available for download. Availability: CRONOS is freely available to non-commercial usersat http://mips.gsf.de/genre/proj/cronos/index.html, web servicesare available at http://mips.gsf.de/CronosWSService/CronosWS?wsdl. Contact: brigitte.waegele{at}helmholtz-muenchen.de Supplementary information: Supplementary data are availableat Bioinformatics online. The online Supplementary Materialcontains all figures and tables referenced by this article. Associate Editor: Martin Bishop     

Functional genomics via multiscale analysis: application to gene expression and ChIP-on-chip data

We present a fast, versatile and adaptive-multiscale algorithm for analyzing a wide-variety of DNA microarray data. Its primary application is in normalization of array data as well as subsequent identification of 'enriched targets', e.g. differentially expressed genes in expression profiling arrays and enriched sites in ChIP-on-chip experimental data. We show how to accommodate the unique characteristics of ChIP-on-chip data, where the set of 'enriched targets' is large, asymmetric and whose proportion to the whole data varies locally. SUPPLEMENTARY INFORMATION: Supplementary figures, related preprint, free software as well as our raw DNA microarray data with PCR validations are available at http://www.math.umn.edu/~lerman/supp/bioinfo06 as well as Bioinformatics online.     

Optimal spliced alignment of homologous cDNA to a genomic DNA template

MOTIVATION: Supplementary cDNA or EST evidence is often decisive for discriminating between alternative gene predictions derived from computational sequence inspection by any of a number of requisite programs. Without additional experimental effort, this approach must rely on the occurrence of cognate ESTs for the gene under consideration in available, generally incomplete, EST collections for the given species. In some cases, particular exon assignments can be supported by sequence matching even if the cDNA or EST is produced from non-cognate genomic DNA, including different loci of a gene family or homologous loci from different species. However, marginally significant sequence matching alone can also be misleading. We sought to develop an algorithm that would simultaneously score for predicted intrinsic splice site strength and sequence matching between the genomic DNA template and a related cDNA or EST. In this case, weakly predicted splice sites may be chosen for the optimal scoring spliced alignment on the basis of surrounding sequence matching. Strongly predicted splice sites will enter the optimal spliced alignment even without strong sequence matching. RESULTS: We designed a novel algorithm that produces the optimal spliced alignment of a genomic DNA with a cDNA or EST based on scoring for both sequence matching and intrinsic splice site strength. By example, we demonstrate that this combined approach appears to improve gene prediction accuracy compared with current methods that rely only on either search by content and signal or on sequence similarity. AVAILABILITY: The algorithm is available as a C subroutine and is implemented in the SplicePredictor and GeneSeqer programs. The source code is available via anonymous ftp from ftp. zmdb.iastate.edu. Both programs are also implemented as a Web service at http://gremlin1.zool.iastate.edu/cgi-bin/s p.cgiand http://gremlin1.zool.iastate.edu/cgi-bin/g s.cgi, respectively. CONTACT: vbrendel@iastate.edu     

The hypolipidemic compound cetaben induces changes in Golgi morphology and vesicle movement

The human hepatoma cell line HepG2 was used to study the effect of cetaben, a non-fibrate hypolipidemic drug, on cell morphology and vesicle distribution. Cetaben treatment correlated with a fragmentation and/or condensation of Golgi cisternae and the appearance of large electron-lucent vesicles. The Golgi apparatus, demonstrated, for example, by fluorescence-lectin histochemistry, was fragmented after cetaben treatment. The lectin-positive remnants were dispersed throughout the cytoplasm, but with a preference for being transported to tips of cells. However, microtubules and the intermediate filaments as well as the actin microfilaments were unchanged after cetaben treatment indicating that changes in Golgi morphology are not caused by alterations in the cytoskeleton. Cetaben decreases the cholesterol content due to inhibition of cholesterol biosynthesis. Changes in the intracellular cholesterol content are known to influence the intracellular vesicle distribution and are most likely responsible for cetaben-induced Golgi alterations, as depletion of cellular cholesterol by starvation or lovastatin and/or cyclodextrin treatment resulted in a similar redistribution of Golgi-derived wheat germ agglutinin vesicles. These lectin-stained vesicles colocalized with lysosomal marker proteins such as Limp-1 and Lamp-2, but not with the early endosomal markers Rab5 and EEA1. Upon removal of cetaben the lectin- and Limp-1/Lamp-2-costained vesicles dissociated and were transported back to the perinuclear region. Thus, cetaben-induced changes such as fragmentation of the Golgi apparatus and the dispersion of lysosomes away from their juxtanuclear location were reversible.Electronic Supplementary Material Supplementary material is available in the online version of this article at http://dx.doi.org/10.1007/s00418-004-0689-6     

Classification of wheat low-molecular-weight glutenin subunit genes and its chromosome assignment by developing LMW-GS group-specific primers

On the basis of sequence analysis, 69 known low-molecular-weight glutenin subunit (LMW-GS) genes were experimentally classified into nine groups by the deduced amino acid sequence of the highly conserved N-terminal domain. To clarify the chromosomal locations of these groups, 11 specific primer sets were designed to carry out polymerase chain reactions (PCR) with the genomic DNA of group 1 ditelosomic lines of Chinese Spring, among which nine primer sets proved to be LMW-GS group-specific. Each group of LMW-GS genes was specifically assigned on a single chromosome arm and hence to a specific locus. Therefore, these results provided the possibility to predict the chromosome location of a new LMW-GS gene based on its deduced N-terminal sequence. The validity of the classification was confirmed by the amplifications in 27 diploid wheat and Aegilops accessions. The length polymorphisms of LMW-GS genes of groups 1 and 2, and groups 3 and 4.1 were detected in diploid A-genome and S-genome accessions, respectively. The diploid wheat and Aegilops species could be used as valuable resources of novel allele variations of LMW-GS gene in the improvement of wheat quality. The nine LMW-GS group-specific primer sets could be utilized to select specific allele variations of LMW-GS genes in the marker-assisted breeding. Electronic Supplementary Material Supplementary material is available for this article at Hai Long and Yu-Ming Wei are the two authors who have contributed equally to this paper     

Structural bioinformatics of DNA: a web-based tool for the analysis of molecular dynamics results and structure prediction

We report here the release of a web-based tool (MDDNA) to study and model the fine structural details of DNA on the basis of data extracted from a set of molecular dynamics (MD) trajectories of DNA sequences involving all the unique tetranucleotides. The dynamic web interface can be employed to analyze the first neighbor sequence context effects on the 10 unique dinucleotide steps of DNA. Functionality is included to build all atom models of any user-defined sequence based on the MD results. The backend of this interface is a relational database storing the conformational details of DNA obtained in 39 different MD simulation trajectories comprising all the 136 unique tetranucleotide steps. Examples of the use of this data to predict DNA structures are included. Availability: http://humphry.chem.wesleyan.edu:8080/MDDNA. Supplementary information: Supplementary data including color figures are available at Bioinformatics online.     

A high-throughput processing service for retention time alignment of complex proteomics and metabolomics LC-MS data

Warp2D is a novel time alignment approach, which uses the overlapping peak volume of the reference and sample peak lists to correct misleading peak shifts. Here, we present an easy-to-use web interface for high-throughput Warp2D batch processing time alignment service using the Dutch Life Science Grid, reducing processing time from days to hours. This service provides the warping function, the sample chromatogram peak list with adjusted retention times and normalized quality scores based on the sum of overlapping peak volume of all peaks. Heat maps before and after time alignment are created from the arithmetic mean of the sum of overlapping peak area rearranged with hierarchical clustering, allowing the quality control of the time alignment procedure. Taverna workflow and command line tool are provided for remote processing of local user data. AVAILABILITY: online data processing service is available at http://www.nbpp.nl/warp2d.html. Taverna workflow is available at myExperiment with title '2D Time Alignment-Webservice and Workflow' at http://www.myexperiment.org/workflows/1283.html. Command line tool is available at http://www.nbpp.nl/Warp2D_commandline.zip. CONTACT: p.l.horvatovich@rug.nl SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Computational approaches for modeling human intestinal absorption and permeability

Human intestinal absorption (HIA) is an important roadblock in the formulation of new drug substances. Computational models are needed for the rapid estimation of this property. The measurements are determined via in vivo experiments or in vitro permeability studies. We present several computational models that are able to predict the absorption of drugs by the human intestine and the permeability through human Caco-2 cells. The training and prediction sets were derived from literature sources and carefully examined to eliminate compounds that are actively transported. We compare our results to models derived by other methods and find that the statistical quality is similar. We believe that models derived from both sources of experimental data would provide greater consistency in predictions. The performance of several QSPR models that we investigated to predict outside the training set for either experimental property clearly indicates that caution should be exercised while applying any of the models for quantitative predictions. However, we are able to show that the qualitative predictions can be obtained with close to a 70% success rate.Electronic Supplementary Material Supplementary material is available for this article at .Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday.     

AM1* parameters for aluminum, silicon, titanium and zirconium

Our extension of the AM1 semiempirical molecular orbital technique, AM1*, has been parameterized for the elements Al, Si, Ti and Zr. The basis sets for all four metals contain a set of d-orbitals. Thus, AM1* parameters are now available for H, C, N, O and F (which use the original AM1 parameters), Al, Si, P, S, Cl, Ti, Mo and Zr. Special attention was paid to reproducing homolytic and heterolytic bond-dissociation energies correctly. Such bond-energy data help to avoid eccentricities in the parameterization caused by inaccurate experimental heats of formation. The performance and typical errors of AM1* for the newly parameterized elements are discussed. Generally, the new method performs less well than established techniques for heats of formation but considerably better for the heats of reaction. Electronic Supplementary Material Supplementary material is available for this article at     

Generation of trangenic Xenopus laevis using the Sleeping Beauty transposon system

Using the Sleeping Beauty (SB) transposon system, we have developed a simple method for the generation of Xenopus laevis transgenic lines. The transgenesis protocol is based on the co-injection of the SB transposase mRNA and a GFP-reporter transposon into one-cell stage embryos. Transposase-dependent reporter gene expression was observed in cell clones and in hemi-transgenic animals. We determined an optimal ratio of transposase mRNA versus transposon-carrying plasmid DNA that enhanced the proportion of hemi-transgenic tadpoles. The transgene is integrated into the genome and may be transmitted to the F1 offspring depending on the germline mosaicism. Although the transposase is necessary for efficient generation of transgenic Xenopus, the integration of the transgene occurred by an non-canonical transposition process. This was observed for two transgenic lines analysed. The transposon-based technique leads to a high transgenesis rate and is simple to handle. For these reasons, it could present an attractive alternative to the classical Restriction Enzyme Mediated Integration (REMI) procedure.Electronic Supplementary Material Supplementary material is available to authorised users in the online version of this article at .