Biomimetic in silico devices

We introduce biomimetic in silico devices, and means for validation along with methods for testing and refining them. The devices are constructed from adaptable software components designed to map logically to biological components at multiple levels of resolution. In this report we focus on the liver; the goal is to validate components that mimic features of the lobule (the hepatic primary functional unit) and dynamic aspects of liver behavior, structure, and function. An assembly of lobule-mimetic devices represents an in silico liver. We validate against outflow profiles for sucrose administered as a bolus to isolated, perfused rat livers. Acceptable in silico profiles are experimentally indistinguishable from those of the in situ referent. This new technology is intended to provide powerful Dew tools for challenging our understanding of how biological functional units function in vivo.

On systems thinking, systems biology and the in Silico Plant

The recent summary report of a Department of Energy Workshop on Plant Systems Biology (P.V. Minorsky [2003] Plant Physiol 132: 404-409) offered a welcomed advocacy for systems analysis as essential in understanding plant development, growth, and production. The goal of the Workshop was to consider methods for relating the results of molecular research to real-world challenges in plant production for increased food supplies, alternative energy sources, and environmental improvement. The rather surprising feature of this report, however, was that the Workshop largely overlooked the rich history of plant systems analysis extending over nearly 40 years (Sinclair and Seligman, 1996) that has considered exactly those challenges targeted by the Workshop. Past systems research has explored and incorporated biochemical and physiological knowledge into plant simulation models from a number of perspectives. The research has resulted in considerable understanding and insight about how to simulate plant systems and the relative contribution of various factors in influencing plant production. These past activities have contributed directly to research focused on solving the problems of increasing biomass production and crop yields. These modeling approaches are also now providing an avenue to enhance integration of molecular genetic technologies in plant improvement (Hammer et al., 2002).

In vitro and in silico analysis of signal peptides from the human blood fluke, Schistosoma mansoni

Proteins secreted by and anchored on the surfaces of parasites are in intimate contact with host tissues. The transcriptome of infective cercariae of the blood fluke, Schistosoma mansoni, was screened using signal sequence trap to isolate cDNAs encoding predicted proteins with an N-terminal signal peptide. Twenty cDNA fragments were identified, most of which contained predicted signal peptides or transmembrane regions, including a novel putative seven-transmembrane receptor and a membrane-associated mitogen-activated protein kinase. The developmental expression pattern within different life-cycle stages ranged from ubiquitous to a transcript that was highly upregulated in the cercaria. A bioinformatics-based comparison of 100 signal peptides from each of schistosomes, humans, a parasitic nematode and Escherichia coli showed that differences in the sequence composition of signal peptides, notably the residues flanking the predicted cleavage site, might account for the negative bias exhibited in the processing of schistosome signal peptides in mammalian cells. (c) 2005 Federation of European Microbiological Societies. Published by Elsevier B.V. All rights reserved.

Reconstruction and in silico analysis of the MAPK signaling pathways in the human blood fluke, Schistosoma japonicum

At present, little is known about signal transduction mechanisms in schistosomes, which cause the disease of schistosomiasis. The mitogen-activated protein kinase (MAPK) signaling pathways, which are evolutionarily conserved from yeast to Homo sapiens, play key roles in multiple cellular processes. Here, we reconstructed the hypothetical MAPK signaling pathways in Schistosoma japonicum and compared the schistosome pathways with those of model eukaryote species. We identified 60 homologous components in the S. japoncium MAPK signaling pathways. Among these, 27 were predicted to be full-length sequences. Phylogenetic analysis of these proteins confirmed the evolutionary conservation of the MAPK signaling pathways. Remarkably, we identified S. japonicum homologues of GTP-binding protein beta and alpha-I subunits in the yeast mating pathway, which might be involved in the regulation of different life stages and female sexual maturation processes as well in schistosomes. In addition, several pathway member genes, including ERK, JNK, Sja-DSP, MRAS and RAS, were determined through quantitative PCR analysis to be expressed in a stage-specific manner, with ERK, JNK and their inhibitor Sja-DSP markedly upregulated in adult female schistosomes. (c) 2006 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved.

Conserved modularity and potential for alternate splicing in mouse and human Slit genes

The vertebrate Slit gene family currently consists of three members;Slit1,Slit2 and Slit3. Each gene encodes a protein containing multiple epidermal growth factor and leucine rich repeat motifs, which are likely to have importance in cell-cell interactions. In this study, we sought to fully define and characterise the vertebrate Slit gene family. Using long distance PCR coupled with in silico mapping, we determined the genomic structure of all three Slit genes in mouse and man. Analysis of EST and genomic databases revealed no evidence of further Slit family members in either organism. All three Slit genes were encoded by 36 (Slit3) or 37 (Slit1 and Slit2) exons covering at least 143 kb or 183 kb of mouse or human genomic DNA respectively. Two additional potential leucine-rich repeat encoding exons were identified within intron 12 of Slit2. These could be inserted in frame, suggesting that alternate splicing may occur in Slit2 A search for STS sequences within human Slit3 anchored this gene to D5S2075 at the 5' end (exon 4) and SGC32449 within the 3' UTR, suggesting that Slit3 may cover greater than 693 kb. The genomic structure of all Slit genes demonstrated considerable modularity in the placement of exon-intron boundaries such that individual leucine-rich repeat motifs were encoded by individual 72 by exons. This further implies the potential generation of multiple Slit protein isoforms varying in their number of repeat units. cDNA library screening and EST database searching verified that such alternate splicing does occur.

Structural basis and prediction of substrate specificity in protein serine/threonine kinases

The large number of protein kinases makes it impractical to determine their specificities and substrates experimentally. Using the available crystal structures, molecular modeling, and sequence analyses of kinases and substrates, we developed a set of rules governing the binding of a heptapeptide substrate motif (surrounding the phosphorylation site) to the kinase and implemented these rules in a web-interfaced program for automated prediction of optimal substrate peptides, taking only the amino acid sequence of a protein kinase as input. We show the utility of the method by analyzing yeast cell cycle control and DNA damage checkpoint pathways. Our method is the only available predictive method generally applicable for identifying possible substrate proteins for protein serine/threonine kinases and helps in silico construction of signaling pathways. The accuracy of prediction is comparable to the accuracy of data from systematic large-scale experimental approaches.

Intracellular compartmentation in planctomycetes

The phylum Planctomycetes of the domain Bacteria consists of budding, peptidoglycan-less organisms important for understanding the origins of complex cell organization. Their significance for cell biology lies in their possession of intracellular membrane compartmentation. All planctomycetes share a unique cell plan, in which the cell cytoplasm is divided into compartments by one or more membranes, including a major cell compartment containing the nucleoid. Of special significance is Gemmata obscuriglobus, in which the nucleoid is enveloped in two membranes to form a nuclear body that is analogous to the structure of a eukaryotic nucleus. Planctomycete compartmentation may have functional physiological roles, as in the case of anaerobic ammonium-oxidizing anammox planctomycetes, in which the anammoxosome harbors specialized enzymes and is wrapped in an envelope possessing unique ladderane lipids. Organisms in phyla other than the phylum Planctomycetes may possess compartmentation similar to that of some planctomycetes, as in the case of members of the phylum Poribacteria from marine sponges.

Identification and analysis of novel genes expressed in the mouse embryonic facial primordia

Craniofacial anomalies are a common feature of human congenital dysmorphology syndromes, suggesting that genes expressed in the developing face are likely to play a wider role in embryonic development. To facilitate the identification of genes involved in embryogenesis, we previously constructed an enriched cDNA library by subtracting adult mouse liver cDNA from that of embryonic day (E)10.5 mouse pharyngeal arch cDNA. From this library, 273 unique clones were sequenced and known proteins binned into functional categories in order to assess enrichment of the library (1). We have now selected 31 novel and poorly characterised genes from this library and present bioinformatic analysis to predict proteins encoded by these genes, and to detect evolutionary conservation. Of these genes 61% (19/31) showed restricted expression in the developing embryo, and a subset of these was chosen for further in silico characterisation as well as experimental determination of subcellular localisation based on transient transfection of predicted full-length coding sequences into mammalian cell lines. Where a human orthologue of these genes was detected, chromosomal localisation was determined relative to known loci for human congenital disease.

Potencies of human immunodeficiency virus protease inhibitors in vitro against Plasmodium falciparum and in vivo against murine malaria

Parasite resistance to antimalarial drugs is a serious threat to human health, and novel agents that act on enzymes essential for parasite metabolism, such as proteases, are attractive targets for drug development. Recent studies have shown that clinically utilized human immunodeficiency virus (HIV) protease inhibitors can inhibit the in vitro growth of Plasmodium falciparum at or below concentrations found in human plasma after oral drug administration. The most potent in vitro antimalarial effects have been obtained for parasites treated with saquinavir, ritonavir, or lopinavir, findings confirmed in this study for a genetically distinct P. falciparum line (3D7). To investigate the potential in vivo activity of antiretroviral protease inhibitors (ARPIs) against malaria, we examined the effect of ARPI combinations in a murine model of malaria. In mice infected with Plasmodium chabaudi AS and treated orally with ritonavir-saquinavir or ritonavir-lopinavir, a delay in patency and a significant attenuation of parasitemia were observed. Using modeling and ligand docking studies we examined putative ligand binding sites of ARPIs in aspartyl proteases of P. falciparum (plasmepsins II and IV) and P. chabaudi (plasmepsin) and found that these in silico analyses support the antimalarial activity hypothesized to be mediated through inhibition of these enzymes. In addition, in vitro enzyme assays demonstrated that P. falciparum plasmepsins II and IV are both inhibited by the ARPIs saquinavir, ritonavir, and lopinavir. The combined results suggest that ARPIs have useful antimalarial activity that may be especially relevant in geographical regions where HIV and P. falciparum infections are both endemic.

Transcriptional network dynamics in macrophage activation

Transcriptional regulatory networks govern cell differentiation and the cellular response to external stimuli. However, mammalian model systems have not yet been accessible for network analysis. Here, we present a genome-wide network analysis of the transcriptional regulation underlying the mouse macrophage response to bacterial lipopolysaccharide (LPS). Key to uncovering the network structure is our combination of time-series cap analysis of gene expression with in silico prediction of transcription factor binding sites. By integrating microarray and qPCR time-series expression data with a promoter analysis, we find dynamic subnetworks that describe how signaling pathways change dynamically during the progress of the macrophage LPS response, thus defining regulatory modules characteristic of the inflammatory response. In particular, our integrative analysis enabled us to suggest novel roles for the transcription factors ATF-3 and NRF-2 during the inflammatory response. We believe that our system approach presented here is applicable to understanding cellular differentiation in higher eukaryotes. (c) 2006 Elsevier Inc. All rights reserved.

Estimating the potential refolding yield of recombinant proteins expressed as inclusion bodies

Recombinant protein production in bacteria is efficient except that insoluble inclusion bodies form when some gene sequences are expressed. Such proteins must undergo renaturation, which is an inefficient process due to protein aggregation on dilution from concentrated denaturant. In this study, the protein-protein interactions of eight distinct inclusion-body proteins are quantified, in different solution conditions, by measurement of protein second virial coefficients (SVCs). Protein solubility is shown to decrease as the SVC is reduced (i.e., as protein interactions become more attractive). Plots of SVC versus denaturant concentration demonstrate two clear groupings of proteins: a more aggregative group and a group having higher SVC and better solubility. A correlation of the measured SVC with protein molecular weight and hydropathicity, that is able to predict which group each of the eight proteins falls into, is presented. The inclusion of additives known to inhibit aggregation during renaturation improves solubility and increases the SVC of both protein groups. Furthermore, an estimate of maximum refolding yield (or solubility) using high-performance liquid chromatography was obtained for each protein tested, under different environmental conditions, enabling a relationship between yield and SVC to be demonstrated. Combined, the results enable an approximate estimation of the maximum refolding yield that is attainable for each of the eight proteins examined, under a selected chemical environment. Although the correlations must be tested with a far larger set of protein sequences, this work represents a significant move beyond empirical approaches for optimizing renaturation conditions. The approach moves toward the ideal of predicting maximum refolding yield using simple bioinformatic metrics that can be estimated from the gene sequence. Such a capability could potentially screen, in silico, those sequences suitable for expression in bacteria from those that must be expressed in more complex hosts. (C) 2004 Wiley Periodicals, Inc.

Trait physiology and crop modelling as a framework to link phenotypic complexity to underlying genetic systems

New tools derived from advances in molecular biology have not been widely adopted in plant breeding for complex traits because of the inability to connect information at gene level to the phenotype in a manner that is useful for selection. In this study, we explored whether physiological dissection and integrative modelling of complex traits could link phenotype complexity to underlying genetic systems in a way that enhanced the power of molecular breeding strategies. A crop and breeding system simulation study on sorghum, which involved variation in 4 key adaptive traits-phenology, osmotic adjustment, transpiration efficiency, stay-green-and a broad range of production environments in north-eastern Australia, was used. The full matrix of simulated phenotypes, which consisted of 547 location-season combinations and 4235 genotypic expression states, was analysed for genetic and environmental effects. The analysis was conducted in stages assuming gradually increased understanding of gene-to-phenotype relationships, which would arise from physiological dissection and modelling. It was found that environmental characterisation and physiological knowledge helped to explain and unravel gene and environment context dependencies in the data. Based on the analyses of gene effects, a range of marker-assisted selection breeding strategies was simulated. It was shown that the inclusion of knowledge resulting from trait physiology and modelling generated an enhanced rate of yield advance over cycles of selection. This occurred because the knowledge associated with component trait physiology and extrapolation to the target population of environments by modelling removed confounding effects associated with environment and gene context dependencies for the markers used. Developing and implementing this gene-to-phenotype capability in crop improvement requires enhanced attention to phenotyping, ecophysiological modelling, and validation studies to test the stability of candidate genetic regions.