The purification and characterization of phosphonopyruvate hydrolase, a novel carbon-phosphorus bond cleavage enzyme from Variovorax sp. Pal2.

Phosphonopyruvate hydrolase, a novel bacterial carbon-phosphorus bond cleavage enzyme, was purified to homogeneity by a series of chromatographic steps from cell extracts of a newly isolated environmental strain of Variovorax sp. Pal2. The enzyme was inducible in the presence of phosphonoalanine or phosphonopyruvate; unusually, its expression was independent of the phosphate status of the cell. The native enzyme had a molecular mass of 63 kDa with a subunit mass of 31.2 kDa. Activity of purified phosphonopyruvate hydrolase was Co2+-dependent and showed a pH optimum of 6.7–7.0. The enzyme had a Km of 0.53 mM for its sole substrate, phosphonopyruvate, and was inhibited by the analogues phosphonoformic acid, 3-phosphonopropionic acid, and hydroxymethylphosphonic acid. The nucleotide sequence of the phosphonopyruvate hydrolase structural gene indicated that it is a member of the phosphoenolpyruvate phosphomutase/isocitrate lyase superfamily with 41% identity at the amino acid level to the carbon-to-phosphorus bond-forming enzyme phosphoenolpyruvate phosphomutase from Tetrahymena pyriformis. Thus its apparently ancient evolutionary origins differ from those of each of the two carbon-phosphorus hydrolases that have been reported previously; phosphonoacetaldehyde hydrolase is a member of the haloacetate dehalogenase family, whereas phosphonoacetate hydrolase belongs to the alkaline phosphatase superfamily of zinc-dependent hydrolases. Phosphonopyruvate hydrolase is likely to be of considerable significance in global phosphorus cycling, because phosphonopyruvate is known to be a key intermediate in the formation of all naturally occurring compounds that contain the carbon-phosphorus bond.

Structural influence of gene networks on their inference: analysis of C3NET

Background: The availability of large-scale high-throughput data possesses considerable challenges toward their functional analysis. For this reason gene network inference methods gained considerable interest. However, our current knowledge, especially about the influence of the structure of a gene network on its inference, is limited.

Promoter region of the Escherichia coli O7-specific lipopolysaccharide gene cluster:structural and functional characterization of an upstream untranslated mRNA sequence

We report the identification of the promoter region of the Escherichia coli O7-specific lipopolysaccharide (LPS) gene cluster (wbEcO7). Typical -10 and -35 sequences were found to be located in the intervening region between galF and rlmB, the first gene of the wbEcO7 cluster. Data from RNase protection experiments revealed the existence of an untranslated leader mRNA segment of 173 bp, including the JUMPStart and two ops sequences. We characterized the structure of this leader mRNA by using the program Mfold and a combination of nested and internal deletions transcriptionally fused to a promoterless lac operon. Our results indicated that the leader mRNA may fold into a series of complex stem-loop structures, one of which includes the JUMPStart element. We have also found that one of the ops sequences resides on the predicted stem and the other resides on the loop region, and we confirmed that these sequences are essential for the RfaH-mediated regulation of the O polysaccharide cluster. A very similar stem-loop structure could be predicted in the promoter region of the LPS core operon encoding the waaQGPSBIJYZK genes. We observed another predicted stem-loop, located immediately downstream from the wbEcO7 transcription initiation site, which appeared to be involved in premature termination of transcription. This putative stem-loop is common to many other O polysaccharide gene clusters but is not present in core oligosaccharide genes. wbEcO7-lac transcriptional fusions in single copy numbers were also used to determine the effects of various environmental cues in the transcriptional regulation of O polysaccharide synthesis. No effects were detected with temperature, osmolarity, Mg2+ concentration, and drugs inducing changes in DNA supercoiling. We therefore conclude that the wbEcO7 promoter activity may be constitutive and that regulation takes place at the level of elongation of the mRNA in a RfaH-mediated manner.

Harnessing the complexity of gene expression data from cancer: from single gene to structural pathway methods

High-dimensional gene expression data provide a rich source of information because they capture the expression level of genes in dynamic states that reflect the biological functioning of a cell. For this reason, such data are suitable to reveal systems related properties inside a cell, e.g., in order to elucidate molecular mechanisms of complex diseases like breast or prostate cancer. However, this is not only strongly dependent on the sample size and the correlation structure of a data set, but also on the statistical hypotheses tested. Many different approaches have been developed over the years to analyze gene expression data to (I) identify changes in single genes, (II) identify changes in gene sets or pathways, and (III) identify changes in the correlation structure in pathways. In this paper, we review statistical methods for all three types of approaches, including subtypes, in the context of cancer data and provide links to software implementations and tools and address also the general problem of multiple hypotheses testing. Further, we provide recommendations for the selection of such analysis methods.

Partial structure of the phylloxin gene from the giant monkey frog, Phyllomedusa bicolor: parallel cloning of precursor cDNA and genomic DNA from lyophilized skin secretion

Phylloxin is a novel prototype antimicrobial peptide from the skin of Phyllomedusa bicolor. Here, we describe parallel identification and sequencing of phylloxin precursor transcript (mRNA) and partial gene structure (genomic DNA) from the same sample of lyophilized skin secretion using our recently-described cloning technique. The open-reading frame of the phylloxin precursor was identical in nucleotide sequence to that previously reported and alignment with the nucleotide sequence derived from genomic DNA indicated the presence of a 175 bp intron located in a near identical position to that found in the dermaseptins. The highly-conserved structural organization of skin secretion peptide genes in P. bicolor can thus be extended to include that encoding phylloxin (plx). These data further reinforce our assertion that application of the described methodology can provide robust genomic/transcriptomic/peptidomic data without the need for specimen sacrifice.

Morphology and localization of interstitial cells in the guinea pig bladder: structural relationships with smooth muscle and neurons.

PURPOSE: In the current study we examined the location of interstitial cell of Cajal (ICC)-like cells in the guinea pig bladder wall and studied their structural interactions with nerves and smooth muscle cells. MATERIALS AND METHODS: Whole mount samples and cryosections of bladder tissue were labeled with primary and fluorescent secondary antibodies, and imaged using confocal and multiphoton microscopy. RESULTS: Kit positive ICC-like cells were located below the urothelium, in the lamina propria region and throughout the detrusor. In the suburothelium they had a stellate morphology and appeared to network. They made connections with nerves, as shown by double labeling experiments with anti-kit and anti-protein gene product 9.5. A network of vimentin positive cells was also found, of which many but not all were kit positive. In the detrusor kit positive cells were most often seen at the edge of smooth muscle bundles. They were elongated with lateral branches, running in parallel with the bundles and closely associated with intramural nerves. Another population of kit positive cells was seen in the detrusor between muscle bundles. These cells had a more stellate-like morphology and made connections with each other. Kit positive cells were seen tracking nerve bundles and close to intramural ganglia. Vimentin positive cells were present in the detrusor, of which some were also kit positive. CONCLUSIONS: There are several populations of ICC-like cells throughout the guinea pig bladder wall. They differ in morphology and orientation but all make connections with intramural nerves and in the detrusor they are closely associated with smooth muscle cells.

Structure-activity studies on position 14 of human alpha-calcitonin gene-related peptide

A structure-activity study was performed to examine the role of position 14 of human alpha-calcitonin gene-related peptide (h-alpha-CGRP) in activating the CGRP receptor. Interestingly, position 14 of h-alpha-CGRP contains a glycyl residue and is part of an alpha-helix spanning residues 8-18. Analogues [Ala(14)]-h-alpha-CGRP, [Aib(14)]-h-alpha-CGRP, [Asp(14)]-h-alpha-CGRP, [Asn(14)]-h-alpha-CGRP, and [Pro(14)]-h-alpha-CGRP were synthesized by solid phase peptide methodology and purified by RP-HPLC. Secondary structure was measured by circular dichroism spectroscopy. Agonist activities were determined as the analogues' ability to stimulate amylase secretion from guinea pig pancreatic acini and to relax precontracted porcine coronary arteries. Analogues [Ala(1)4]-h-alpha-CGRP, [Aib(14)]-h-alpha-CGRP, [Asp(14)]-h-alpha-CGRP, and [Asn(14)]-h-alpha-CGRP, all containing residues with a high helical propensity in position 14, were potent full agonists compared to h-alpha-CGRP in both tissues. Interestingly, replacement of Gly(14) of h-alpha-CGRP with these residues did not substantially increase the helical content of these analogues. [Pro(14)]-h-alpha-CGRP, predictably, has significantly lower helical content and is a 20-fold less potent agonist on coronary artery, known to contain CGRP-1 receptor subtypes, and an antagonist on pancreatic acini, known to contain CGRP-2 receptor subtypes. In conclusion, the residue in position 14 plays a structural role in stabilizing the alpha-helix spanning residues 8-18. The alpha-helix is crucial for maintaining highly potent agonist effects of h-alpha-CGRP at CGRP receptors. The wide variety of functional groups that can be tolerated in position 14 with no substantial modification of agonist effects suggests the residue in this position is not in contact with the CGRP receptor. [Pro(14)]-h-alpha-CGRP may be a useful pharmacological tool to distinguish between CGRP-1 and CGRP-2 receptor subtypes.

Bradykinin-related peptides, including a novel structural variant, (Val1)-bradykinin, from the skin secretion of Guenther's frog, Hylarana guentheri and their molecular precursors.

Multiple bradykinin-related peptides including a novel bradykinin structural variant, (Val1)-bradykinin, have been identified from the defensive skin secretion of Guenther's frog, Hylarana guentheri by a tandem mass spectrometry method. Subsequently, four different preprobradykinin cDNAs, which encoded multiple bradykinin copies and its structural variants, were consistently cloned from a skin derived cDNA library. These preprobradykinin cDNAs showed little structural similarity with mammalian kininogens and the kininogens from the skin of toads, but have regions that are highly conserved in the kininogens from another ranid frog, Odorrana schmackeri. Alignment of these preprobradykinins revealed that preprobradykinin 1, 2 and 3 may derive from a single gene by alternative exon splicing.

Expression of the phosphonoalanine-degradative gene cluster from Variovorax sp Pal2 is induced by growth on phosphonoalanine and phosphonopyruvate

The phosphonopyruvate hydrolase (PalA) found in Variovorax sp., Pal2, is a novel carbon-phosphorus bond cleavage enzyme, which is expressed even in the presence of high levels of phosphate, thus permitting phosphonopyruvate to be used as the sole carbon and energy source. Analysis of the regions adjacent to the palA gene revealed the presence of the five structural genes that constitute the 2-amino-3-phosphonopropionic acid (phosphonoalanine)-degradative operon. Reverse transcriptase-PCR (RT-PCR) experiments demonstrated that all five genes in the operon are transcribed as a single mRNA and that their transcription is induced by phosphonoalanine or phosphonopyruvate. Transcriptional fusions of the regulatory region of the phosphonoalanine degradative operon with the gfp gene were constructed. Expression analysis indicated that the presence of a LysR-type regulator (encoded by the palR gene) is essential for the transcription of the structural genes of the operon. Similar gene clusters were found in the sequenced genomes of six bacterial species from the Alpha-, Beta- and Gammaproteobacteria, and analysis of metagenomic libraries revealed that sequences related to palA are widely spread in the marine environment.

sscMap: An extensible Java application for connecting small-molecule drugs using gene-expression signatures

Background
Connectivity mapping is a process to recognize novel pharmacological and toxicological properties in small molecules by comparing their gene expression signatures with others in a database. A simple and robust method for connectivity mapping with increased specificity and sensitivity was recently developed, and its utility demonstrated using experimentally derived gene signatures.

Results
This paper introduces sscMap (statistically significant connections' map), a Java application designed to undertake connectivity mapping tasks using the recently published method. The software is bundled with a default collection of reference gene-expression profiles based on the publicly available dataset from the Broad Institute Connectivity Map 02, which includes data from over 7000 Affymetrix microarrays, for over 1000 small-molecule compounds, and 6100 treatment instances in 5 human cell lines. In addition, the application allows users to add their custom collections of reference profiles and is applicable to a wide range of other 'omics technologies.

Conclusion
The utility of sscMap is two fold. First, it serves to make statistically significant connections between a user-supplied gene signature and the 6100 core reference profiles based on the Broad Institute expanded dataset. Second, it allows users to apply the same improved method to custom-built reference profiles which can be added to the database for future referencing. The software can be freely downloaded from http://purl.oclc.org/NET/sscMap

A simple and robust method for connecting small-molecule drugs using gene-expression signatures

Background

Interaction of a drug or chemical with a biological system can result in a gene-expression profile or signature characteristic of the event. Using a suitably robust algorithm these signatures can potentially be used to connect molecules with similar pharmacological or toxicological properties by gene expression profile. Lamb et al first proposed the Connectivity Map [Lamb et al (2006), Science 313, 1929–1935] to make successful connections among small molecules, genes, and diseases using genomic signatures.

Results

Here we have built on the principles of the Connectivity Map to present a simpler and more robust method for the construction of reference gene-expression profiles and for the connection scoring scheme, which importantly allows the valuation of statistical significance of all the connections observed. We tested the new method with two randomly generated gene signatures and three experimentally derived gene signatures (for HDAC inhibitors, estrogens, and immunosuppressive drugs, respectively). Our testing with this method indicates that it achieves a higher level of specificity and sensitivity and so advances the original method.

Conclusion

The method presented here not only offers more principled statistical procedures for testing connections, but more importantly it provides effective safeguard against false connections at the same time achieving increased sensitivity. With its robust performance, the method has potential use in the drug development pipeline for the early recognition of pharmacological and toxicological properties in chemicals and new drug candidates, and also more broadly in other 'omics sciences.

Inferring the conservative causal core of gene regulatory networks

Background
Inferring gene regulatory networks from large-scale expression data is an important problem that received much attention in recent years. These networks have the potential to gain insights into causal molecular interactions of biological processes. Hence, from a methodological point of view, reliable estimation methods based on observational data are needed to approach this problem practically.

Results
In this paper, we introduce a novel gene regulatory network inference (GRNI) algorithm, called C3NET. We compare C3NET with four well known methods, ARACNE, CLR, MRNET and RN, conducting in-depth numerical ensemble simulations and demonstrate also for biological expression data from E. coli that C3NET performs consistently better than the best known GRNI methods in the literature. In addition, it has also a low computational complexity. Since C3NET is based on estimates of mutual information values in conjunction with a maximization step, our numerical investigations demonstrate that our inference algorithm exploits causal structural information in the data efficiently.

Conclusions
For systems biology to succeed in the long run, it is of crucial importance to establish methods that extract large-scale gene networks from high-throughput data that reflect the underlying causal interactions among genes or gene products. Our method can contribute to this endeavor by demonstrating that an inference algorithm with a neat design permits not only a more intuitive and possibly biological interpretation of its working mechanism but can also result in superior results.

Functional and Structural Diversification of the Anguimorpha Lizard Venom System

Venom has only been recently discovered to be a basal trait of the Anguimorpha lizards. Consequently, very little is known about the timings of toxin recruitment events, venom protein molecular evolution, or even the relative physical diversifications of the venom system itself. A multidisciplinary approach was used to examine the evolution across the full taxonomical range of this similar to 130 million-year-old clade. Analysis of cDNA libraries revealed complex venom transcriptomes. Most notably, three new cardioactive peptide toxin types were discovered (celestoxin, cholecystokinin, and YY peptides). The latter two represent additional examples of convergent use of genes in toxic arsenals, both having previously been documented as components of frog skin defensive chemical secretions. Two other novel venom gland-overexpressed modified versions of other protein frameworks were also recovered from the libraries (epididymal secretory protein and ribonuclease). Lectin, hyaluronidase, and veficolin toxin types were sequenced for the first time from lizard venoms and shown to be homologous to the snake venom forms. In contrast, phylogenetic analyses demonstrated that the lizard natriuretic peptide toxins were recruited independently of the form in snake venoms. The de novo evolution of helokinestatin peptide toxin encoding do-mains within the lizard venom natriuretic gene was revealed to be exclusive to the helodermatid/anguid subclade. New isoforms were sequenced for cysteine-rich secretory protein, kallikrein, and phospholipase A 2 toxins. Venom gland morphological analysis revealed extensive evolutionary tinkering. Anguid glands are characterized by thin capsules and mixed glands, serous at the bottom of the lobule and mucous toward the apex. Twice, independently this arrangement was segregated into specialized serous protein-secreting glands with thick capsules with the mucous lobules now distinct (Heloderma and the Lanthanotus/Varanus clade). The results obtained highlight the importance of utilizing evolution-based search strategies for biodiscovery and emphasize the largely untapped drug design and development potential of lizard venoms. Molecular & Cellular Proteomics 9:2369-2390, 2010.