BioPP: a tool for web-publication of biological networks

<p>Abstract</p> <p>Background</p> <p>Cellular processes depend on the function of intracellular molecular networks. The curation of the literature relevant to specific biological pathways is important for many theoretical and experimental research teams and communities. No current tool supports web publication or hosting of user-developed large scale annotated pathway diagrams. Sharing via web publication is needed to allow real-time access to the current literature pathway knowledgebase, both privately within a research team or publicly among the outside research community. Web publication also facilitates team and/or community input into the curation process while allowing centralized control of the curation and validation process. We have developed new tool to address these needs. Biological Pathway Publisher (BioPP) is a software suite for converting CellDesigner Systems Biology Markup Language (CD-SBML) formatted pathways into a web viewable format. The BioPP suite is available for private use and for depositing knowledgebases into a newly created public repository.</p> <p>Results</p> <p>BioPP suite is a web-based application that allows pathway knowledgebases stored in CD-SBML to be web published with an easily navigated user interface. The BioPP suite consists of four interrelated elements: a pathway publisher, an upload web-interface, a pathway repository for user-deposited knowledgebases and a pathway navigator. Users have the option to convert their CD-SBML files to HTML for restricted use or to allow their knowledgebase to be web-accessible to the scientific community. All entities in all knowledgebases in the repository are linked to public database entries as well as to a newly created public wiki which provides a discussion forum.</p> <p>Conclusion</p> <p>BioPP tools and the public repository facilitate sharing of pathway knowledgebases and interactive curation for research teams and scientific communities. BioPP suite is accessible at <url>http://tsb.mssm.edu/pathwayPublisher/broadcast/</url></p

Method for multiplex cellular detection of mRNAs using quantum dot fluorescent in situ hybridization

The photostability and narrow emission spectra of non-organic quantum dot fluorophores (QDs) make them desirable candidates for fluorescent in situ hybridization (FISH) to study the expression of specific mRNA transcripts. We developed a novel method for direct QD labeling of modified oligonucleotide probes through streptavidin and biotin interactions, as well as protocols for their use in multiple-label FISH. We validated this technique in mouse brainstem sections. The subcellular localization of the vesicular monoamine transporter (Vmat2) mRNA corresponds when using probes labeled with two different QDs in the same hybridization. We developed protocols for combined direct QD FISH and QD immunohistochemical labeling within the same neurons as well as for simultaneous study of the subcellular distribution of multiple mRNA targets. We demonstrated increased sensitivity of FISH using QDs in comparison with organic fluorophores. These techniques gave excellent histological results both for multiplex FISH and combined FISH and immunohistochemistry. This approach can facilitate the ultrasensitive simultaneous study of multiple mRNA and protein markers in tissue culture and histological section

Elucidation of molecular kinetic schemes from macroscopic traces using system identification

Overall cellular responses to biologically-relevant stimuli are mediated by networks of simpler lower-level processes. Although information about some of these processes can now be obtained by visualizing and recording events at the molecular level, this is still possible only in especially favorable cases. Therefore the development of methods to extract the dynamics and relationships between the different lower-level (microscopic) processes from the overall (macroscopic) response remains a crucial challenge in the understanding of many aspects of physiology. Here we have devised a hybrid computational-analytical method to accomplish this task, the SYStems-based MOLecular kinetic scheme Extractor (SYSMOLE). SYSMOLE utilizes system-identification input-output analysis to obtain a transfer function between the stimulus and the overall cellular response in the Laplace-transformed domain. It then derives a Markov-chain state molecular kinetic scheme uniquely associated with the transfer function by means of a classification procedure and an analytical step that imposes general biological constraints. We first tested SYSMOLE with synthetic data and evaluated its performance in terms of its rate of convergence to the correct molecular kinetic scheme and its robustness to noise. We then examined its performance on real experimental traces by analyzing macroscopic calcium-current traces elicited by membrane depolarization. SYSMOLE derived the correct, previously known molecular kinetic scheme describing the activation and inactivation of the underlying calcium channels and correctly identified the accepted mechanism of action of nifedipine, a calcium-channel blocker clinically used in patients with cardiovascular disease. Finally, we applied SYSMOLE to study the pharmacology of a new class of glutamate antipsychotic drugs and their crosstalk mechanism through a heteromeric complex of G protein-coupled receptors. Our results indicate that our methodology can be successfully applied to accurately derive molecular kinetic schemes from experimental macroscopic traces, and we anticipate that it may be useful in the study of a wide variety of biological systems

MOCCASIN: converting MATLAB ODE models to SBML

MATLAB is popular in biological research for creating and simulating models that use ordinary differential equations (ODEs). However, sharing or using these models outside of MATLAB is often problematic. A community standard such as Systems Biology Markup Language (SBML) can serve as a neutral exchange format, but translating models from MATLAB to SBML can be challenging—especially for legacy models not written with translation in mind. We developed MOCCASIN (Model ODE Converter for Creating Automated SBML INteroperability) to help. MOCCASIN can convert ODE-based MATLAB models of biochemical reaction networks into the SBML format

Reconstruction of regulatory networks through temporal enrichment profiling and its application to H1N1 influenza viral infection

BACKGROUND: H1N1 influenza viruses were responsible for the 1918 pandemic that caused millions of deaths worldwide and the 2009 pandemic that caused approximately twenty thousand deaths. The cellular response to such virus infections involves extensive genetic reprogramming resulting in an antiviral state that is critical to infection control. Identifying the underlying transcriptional network driving these changes, and how this program is altered by virally-encoded immune antagonists, is a fundamental challenge in systems immunology. RESULTS: Genome-wide gene expression patterns were measured in human monocyte-derived dendritic cells (DCs) infected in vitro with seasonal H1N1 influenza A/New Caledonia/20/1999. To provide a mechanistic explanation for the timing of gene expression changes over the first 12 hours post-infection, we developed a statistically rigorous enrichment approach integrating genome-wide expression kinetics and time-dependent promoter analysis. Our approach, TIme-Dependent Activity Linker (TIDAL), generates a regulatory network that connects transcription factors associated with each temporal phase of the response into a coherent linked cascade. TIDAL infers 12 transcription factors and 32 regulatory connections that drive the antiviral response to influenza. To demonstrate the generality of this approach, TIDAL was also used to generate a network for the DC response to measles infection. The software implementation of TIDAL is freely available at http://tsb.mssm.edu/primeportal/?q=tidal_prog. CONCLUSIONS: We apply TIDAL to reconstruct the transcriptional programs activated in monocyte-derived human dendritic cells in response to influenza and measles infections. The application of this time-centric network reconstruction method in each case produces a single transcriptional cascade that recapitulates the known biology of the response with high precision and recall, in addition to identifying potentially novel antiviral factors. The ability to reconstruct antiviral networks with TIDAL enables comparative analysis of antiviral responses, such as the differences between pandemic and seasonal influenza infections

Time and length scales of autocrine signals in three dimensions

A model of autocrine signaling in cultures of suspended cells is developed on the basis of the effective medium approximation. The fraction of autocrine ligands, the mean and distribution of distances traveled by paracrine ligands before binding, as well as the mean and distribution of the ligand lifetime are derived. Interferon signaling by dendritic immune cells is considered as an illustration.Comment: 15 page

Functional microdomains in G-protein-coupled receptors: the conserved arginine-cage motif in the gonadotropin-releasing hormone receptor

An Arg present in the third transmembrane domain of all rhodopsin-like G-protein-coupled receptors is required for efficient signal transduction. Mutation of this Arg in the gonadotropin-releasing hormone receptor to Gln, His, or Lys abolished or severely impaired agonist-stimulated inositol phosphate generation, consistent with Arg having a role in receptor activation. To investigate the contribution of the surrounding structural domain in the actions of the conserved Arg, an integrated microdomain modeling and mutagenesis approach has been utilized. Two conserved residues that constrain the Arg side chain to a limited number of conformations have been identified. In the inactive wild-type receptor, the Arg side chain is proposed to form an ionic interaction with Asp3.49(138). Experimental results for the Asp3. 49(138) --> Asn mutant receptor show a modestly enhanced receptor efficiency, consistent with the hypothesis that weakening the Asp3. 49(138)-Arg3.50(139) interaction by protonation of the Asp or by the mutation to Asn favors activation. With activation, the Asp3. 49(138)-Arg3.50(139) ionic bond would break, and the unrestrained Arg would be prevented from orienting itself toward the water phase by a steric clash with Ile3.54(143). The mutation Ile3.54(143) --> Ala, which eliminates this clash in simulations, causes a marked reduction in measured receptor signaling efficiency, implying that solvation of Arg3.50(139) prevents it from functioning in the activation of the receptor. These data are consistent with residues Asp3.49(138) and Ile3.54(143) forming a structural motif, which helps position Arg in its appropriate inactive and active receptor conformations

The Functional Microdomain in Transmembrane Helices 2 and 7 Regulates Expression, Activation, and Coupling Pathways of the Gonadotropin-releasing Hormone Receptor

Structural microdomains of G protein-coupled receptors (GPCRs) consist of spatially related side chains that mediate discrete functions. The conserved helix 2/helix 7 microdomain was identified because the gonadotropin-releasing hormone (GnRH) receptor appears to have interchanged the Asp(2.50) and Asn(7.49) residues which are conserved in transmembrane helices 2 and 7 of rhodopsin-like GPCRs. We now demonstrate that different side chains of this microdomain contribute specifically to receptor expression, heterotrimeric G protein-, and small G protein-mediated signaling. An Asn residue is required in position 2.50(87) for expression of the GnRH receptor at the cell surface, most likely through an interaction with the conserved Asn(1.50(53)) residue, which we also find is required for receptor expression. Most GPCRs require an Asp side chain at either the helix 2 or helix 7 locus of the microdomain for coupling to heterotrimeric G proteins, but the GnRH receptor has transferred the requirement for an acidic residue from helix 2 to 7. However, the presence of Asp at the helix 7 locus precludes small G protein-dependent coupling to phospholipase D. These results implicate specific components of the helix 2/helix 7 microdomain in receptor expression and in determining the ability of the receptor to adopt distinct activated conformations that are optimal for interaction with heterotrimeric and small G proteins