Identification and Characterization of New Members of Vacuolar H<Superscript>+</Superscript>-Pyrophosphatase Family from <Emphasis Type="Italic">Oryza sativa</Emphasis> Genome

The vacuolar H⁺-pyrophosphatase (V-PPase) is an electrogenic H⁺ pump localized in the plant vacuolar membrane. V-PPase from many species has been characterized previously and the corresponding genes/cDNAs have been cloned. Cloning of the V-PPase genes from many plant species has revealed conserved motifs that may correspond to catalytic sites. The completion of the entire DNA sequence of Oryza sativa (430 Mb) presented an opportunity to study the structure and function of V-PPase proteins, and also to identify new members of this family in Oryza sativa. Our analysis identified three novel V-PPase proteins in the Oryza sativa genome that contain functional domains typical of V-PPase. We have designated them as OVP3 to OVP5. The new predicted OVPs have chromosomal locations different from previously characterized V-PPases (OVP1 and OVP2) located on chromosome 6. They all contain three characteristic motifs of V-PPase and also a conserved motif [DE]YYTS, specific to type I V-PPases and involved in coupling PP_i hydrolysis to H⁺ translocation.     

The predictive performance of short-linear motif features in the prediction of calmodulin-binding proteins

Background

The prediction of calmodulin-binding (CaM-binding) proteins plays a very important role in the fields of biology and biochemistry, because the calmodulin protein binds and regulates a multitude of protein targets affecting different cellular processes. Computational methods that can accurately identify CaM-binding proteins and CaM-binding domains would accelerate research in calcium signaling and calmodulin function. Short-linear motifs (SLiMs), on the other hand, have been effectively used as features for analyzing protein-protein interactions, though their properties have not been utilized in the prediction of CaM-binding proteins.

Results

We propose a new method for the prediction of CaM-binding proteins based on both the total and average scores of known and new SLiMs in protein sequences using a new scoring method called sliding window scoring (SWS) as features for the prediction module. A dataset of 194 manually curated human CaM-binding proteins and 193 mitochondrial proteins have been obtained and used for testing the proposed model. The motif generation tool, Multiple EM for Motif Elucidation (MEME), has been used to obtain new motifs from each of the positive and negative datasets individually (the SM approach) and from the combined negative and positive datasets (the CM approach). Moreover, the wrapper criterion with random forest for feature selection (FS) has been applied followed by classification using different algorithms such as k-nearest neighbors (k-NN), support vector machines (SVM), naive Bayes (NB) and random forest (RF).

Conclusions

Our proposed method shows very good prediction results and demonstrates how information contained in SLiMs is highly relevant in predicting CaM-binding proteins. Further, three new CaM-binding motifs have been computationally selected and biologically validated in this study, and which can be used for predicting CaM-binding proteins.

     

Alterations in Digestive Enzymes of Three Freshwater Teleostean Fishes by Almix Herbicide: A Comparative Study

Three freshwater teleostean fishes viz., Anabas testudineus (Bloch), Heteropneustes fossilis (Bloch) and Oreochromis niloticus (Linnaeus) were exposed to almix (66.67 mg/l) herbicide for 30 days to investigate the activity of digestive enzymes (amylase, lipase and protease) in stomach, intestine and liver. Amylase activity showed significantly high (p < 0.05) in all the fishes compared to control value and highest activity was observed in liver of A. testudineus (721.99 %) and minimum in intestine of H. fossilis (195.37 %). Lipase activity was also significantly increased (p < 0.05) in all the tissues; but highest in intestine of O. niloticus (235.51 %) and minimum in intestine of A. testudineus (130.51 %). Protease activity also showed similar trends of enhancement; it was maximum in stomach of O. niloticus (362.69 %), whereas in liver of H. fossilis it was rather less (173.72 %). Increased activity of digestive enzymes resulting from tissue damage ultimately affected the fish health due to impairment of digestive physiology confirming the herbicidal contamination on fish species. The sensitivity to the almix herbicide was pronounced in the order of O. niloticus > A. testudineus > H. fossilis.     

Acid/amide bonding for anthranilic acid derivatives: crystal structures of [W(X)Cl3(HO2CC6H4NH-2)] (X = O, NPh)

The complexes [W(X)Cl₃(HO₂CC₆H₄NH-2)] [X = O (1), NPh (2)] have been obtained by reaction of either [WOCl₄] or [W(NPh)Cl₄(Et₂O)] with anthranilic acid {1,2-(NH₂)(CO₂H)C₆H₄}, respectively. The X-ray crystal structures reveal pseudo-octahedral metal centres, each with a mer-arrangement of chlorines and a chelating acid/amide ligand derived from anthranilic acid. The acid group of this chelate ligand is trans to either the oxo or organoimido functionality.     

Accurate recognition of cis-regulatory motifs with the correct lengths in prokaryotic genomes

We present a new computational method for solving a classical problem, the identification problem of cis-regulatory motifs in a given set of promoter sequences, based on one key new idea. Instead of scoring candidate motifs individually like in all the existing motif-finding programs, our method scores groups of candidate motifs with similar sequences, called motif closures, using a P-value, which has substantially improved the prediction reliability over the existing methods. Our new P-value scoring scheme is sequence length independent, hence allowing direct comparisons among predicted motifs with different lengths on the same footing. We have implemented this method as a Motif Recognition Computer (MREC) program, and have extensively tested MREC on both simulated and biological data from prokaryotic genomes. Our test results indicate that MREC can accurately pick out the actual motif with the correct length as the best scoring candidate for the vast majority of the cases in our test set. We compared our prediction results with two motif-finding programs Cosmo and MEME, and found that MREC outperforms both programs across all the test cases by a large margin. The MREC program is available at http://csbl.bmb.uga.edu/~bingqiang/MREC1/.     

Multiple Motif Scanning to Identify Methyltransferases from the Yeast Proteome

A new program (Multiple Motif Scanning) was developed to scan the Saccharomyces cerevisiae proteome for Class I S-adenosylmethionine-dependent methyltransferases. Conserved Motifs I, Post I, II, and III were identified and expanded in known methyltransferases by primary sequence and secondary structural analysis through hidden Markov model profiling of both a yeast reference database and a reference database of methyltransferases with solved three-dimensional structures. The roles of the conserved amino acids in the four motifs of the methyltransferase structure and function were then analyzed to expand the previously defined motifs. Fisher-based negative log statistical matrix sets were developed from the prevalence of amino acids in the motifs. Multiple Motif Scanning is able to scan the proteome and score different combinations of the top fitting sequences for each motif. In addition, the program takes into account the conserved number of amino acids between the motifs. The output of the program is a ranked list of proteins that can be used to identify new methyltransferases and to reevaluate the assignment of previously identified putative methyltransferases. The Multiple Motif Scanning program can be used to develop a putative list of enzymes for any type of protein that has one or more motifs conserved at variable spacings and is freely available (www.chem.ucla.edu/files/MotifSetup.Zip). Finally hidden Markov model profile clustering analysis was used to subgroup Class I methyltransferases into groups that reflect their methyl-accepting substrate specificity.Enzymes that catalyze the transfer of a methyl group from S-adenosylmethionine to protein, nucleic acid, lipid, and small molecule substrates are widely distributed in nature and function in a variety of biological pathways including metabolic regulation, gene expression, the repair of aging biomolecules, and biosynthesis (1). There are several classes of AdoMet¹ dependent methyltransferases. Class I enzymes are the most abundant and share a common three-dimensional structural core that includes a seven-strand twisted β sheet (2–5). It has been estimated that about 0.6–1.6% of genes in organisms ranging from Escherichia coli, Saccharomyces cerevisiae, Caenorhabditis elegans, Drosophila melanogaster, Arabidopsis thaliana, and humans encode Class I methyltransferases (6). These results suggest that there are some 50 species in yeast and some 300 species in humans. However, most of these assignments have not been confirmed, and only a relatively small fraction of them have been functionally identified. Previous bioinformatics studies on Class I methyltransferases have taken advantage of the fact that the structure of the AdoMet binding site is conserved in the primary sequences of four short signature motifs designated Motifs I, Post I, II, and III (6–8). These motifs are present with variable, but often conserved, spacing in the primary sequences. Initial identifications of Class I methyltransferases in the yeast proteome were based on searches with a “Motifs in Protein Data Bases” program with individual motifs to generate a list of putative methyltransferases (8). More recently, Katz et al. (6) performed a search using the Motif Alignment and Search Tool (MAST) with Multiple Em for Motif Elicitation (MEME)-generated matrices combining Motifs I/Post I as well as PSI-BLAST searches to identify methyltransferases in both yeast and a variety of other organisms. However, the success of these methods was limited by difficulties in identifying these motifs in the known methyltransferases. In fact, it was not possible to take advantage of the information content of Motifs II and III in the latter study (6).In recent years, two developments have provided new windows to improve the screening of proteomes to identify the complete cast of Class I methyltransferases in various organisms. In the first place, three-dimensional structures are now known for a large number of Class I methyltransferases. Because the motifs are closely linked to structural features (2–4), their identification is unambiguous. Secondly secondary structure prediction algorithms can now be used as independent confirmations of the structure-linked sequence motifs in methyltransferases whose structures have not been determined. We wanted to take advantage of the improved motif identification by developing advanced software for searching proteomes for multiple motifs at varying, but partially conserved, spacing. We have now used secondary structure prediction to obtain the motifs of a reference group of 32 known yeast methyltransferases and a search of the Research Collaboratory for Structural Bioinformatics Protein Data Bank to identify the motifs of a second reference group of 33 distinct types of Class I methyltransferases with three-dimensional structures. We describe a method for generating search matrices for each of these reference groups and a program “Multiple Motif Scanning” to score these matrices in the yeast proteome. This approach not only identified new methyltransferases but allowed us to reject some of the previously described candidate methyltransferases. Additionally we used HMM profile clustering analysis to extract more information about the possible substrates of these putative methyltransferases (9).     

Microbial rhamnolipid production: a critical re-evaluation of published data and suggested future publication criteria

High production cost and potential pathogenicity of Pseudomonas aeruginosa, commonly used for rhamnolipid synthesis, have led to extensive research for safer producing strains and cost-effective production methods. This has resulted in numerous research publications claiming new non-pathogenic producing strains and novel production techniques many of which are unfortunately without proper characterisation of product and/or producing strain/s. Genes responsible for rhamnolipid production have only been confirmed in P. aeruginosa, Burkholderia thailandensis and Burkholderia pseudomallei. Comparing yields in different publications is also generally unreliable especially when different methodologies were used for rhamnolipid quantification. After reviewing the literature in this area, we strongly feel that numerous research outputs have insufficient evidence to support claims of rhamnolipid-producing strains and/or yields. We therefore recommend that standards should be set for reporting new rhamnolipid-producing strains and production yields. These should include (1) molecular and bioinformatic tools to fully characterise new microbial isolates and confirm the presence of the rhamnolipid rhl genes for all bacterial strains, (2) using gravimetric methods to quantify crude yields and (3) use of a calibrated method (high-performance liquid chromatography or ultra-performance liquid chromatography) for absolute quantitative yield determination.     

Isolation,characterization and analysis of the agglutinative activity of a lectin from <Emphasis Type="Italic">Crotalaria spectabilis</Emphasis>

Lectins are proteins with ability to recognize specific carbohydrates. These are present in virtually all organisms and have increasing applications in biotechnology. Here, our aim was to purify lectins from seeds of Crotalaria spectabilis Roth and determine their agglutinative ability. In this study, 45 g of seeds were milled, their proteins were precipitated by acetone or ammonium sulfate and purified by exclusion and ion-exchange chromatography. An isolated lectin was submitted to tests for hemagglutination and inhibition of hemagglutinating activity by carbohydrates as well as tests for its response to chelating and reducing agents. Our results show that the apparent molecular weight (as determined by SDS-PAGE) of the lectin is 30 kDa, and the tests for inhibition of erythrocytes’ agglutinative activity by sugars were positive for d-galactose and N-acetyl-d-galactosamine. Data obtained with the chelating agent EDTA demonstrated the presence of divalent cations in the protein structure. However, the reducing agent 2-mercaptoethanol was unable to inhibit the protein’s bioactivity. The lectin agglutinated the blood groups A, B, AB and O, as well as bacterial lineages from the species Leptospira interrogans and Leptospira biflexa, indicating a prospective application in the diagnosis and treatment of leptospirosis.     

Experimental evolution and genome sequencing reveal variation in levels of clonal interference in large populations of bacteriophage φX174

Background

The first step of GPI anchor biosynthesis is catalyzed by PIG-A, an enzyme that transfers N -acetylglucosamine from UDP- N -acetylglucosamine to phosphatidylinositol. This protein is present in all eukaryotic organisms ranging from protozoa to higher mammals, as part of a larger complex of five to six 'accessory' proteins whose individual roles in the glycosyltransferase reaction are as yet unclear. The PIG-A gene has been shown to be an essential gene in various eukaryotes. In humans, mutations in the protein have been associated with paroxysomal noctural hemoglobuinuria. The corresponding PIG-A gene has also been recently identified in the genome of many archaeabacteria although genes of the accessory proteins have not been discovered in them. The present study explores the evolution of PIG-A and the phylogenetic relationship between this protein and other glycosyltransferases.

Results

In this paper we show that out of the twelve conserved motifs identified by us eleven are exclusively present in PIG-A and, therefore, can be used as markers to identify PIG-A from newly sequenced genomes. Three of these motifs are absent in the primitive eukaryote, G. lamblia. Sequence analyses show that seven of these conserved motifs are present in prokaryote and archaeal counterparts in rudimentary forms and can be used to differentiate PIG-A proteins from glycosyltransferases. Using partial least square regression analysis and data involving presence or absence of motifs in a range of PIG-A and glycosyltransferases we show that (i) PIG-A may have evolved from prokaryotic glycosyltransferases and lipopolysaccharide synthases, members of the GT4 family of glycosyltransferases and (ii) it is possible to uniquely classify PIG-A proteins versus glycosyltransferases.

Conclusion

Besides identifying unique motifs and showing that PIG-A protein from G. lamblia and some putative PIG-A proteins from archaebacteria are evolutionarily closer to glycosyltransferases, these studies provide a new method for identification and classification of PIG-A proteins.     

Identification of a Class of Chromatin Boundary Elements

Boundary elements are thought to define the ends of functionally independent domains of genetic activity. An assay for boundary activity based on this concept measures the ability to insulate a bracketed, chromosomally integrated reporter gene from position effects. Despite their presumed importance, the few examples identified to date apparently do not share sequence motifs or DNA binding proteins. The Drosophila protein BEAF binds the scs′ boundary element of the 87A7 hsp70 locus and roughly half of polytene chromosome interband loci. To see if these sites represent a class of boundary elements that have BEAF in common, we have isolated and studied several genomic BEAF binding sites as candidate boundary elements (cBEs). BEAF binds with high affinity to clustered, variably arranged CGATA motifs present in these cBEs. No other sequence homologies were found. Two cBEs were tested and found to confer position-independent expression on a mini-white reporter gene in transgenic flies. Furthermore, point mutations in CGATA motifs that eliminate binding by BEAF also eliminate the ability to confer position-independent expression. Taken together, these findings suggest that clustered CGATA motifs are a hallmark of a BEAF-utilizing class of boundary elements found at many loci. This is the first example of a class of boundary elements that share a sequence motif and a binding protein.Chromatin appears to be partitioned into chromosomal domains that are operationally defined by bracketing DNA regions called boundary elements or insulators (10; see reference 34 for a review). Boundary elements are presumably necessary to curtail the potentially promiscuous behavior of enhancers, limiting their action to the domain in which they reside. The biological activity of a boundary element is experimentally measured by either position-independent expression or enhancer-blocking assays. If this view of chromosomal organization is correct, boundary elements play a very important functional role. Yet only a few examples have been identified, and each is so far a unique case, as they do not appear to have notable sequence homologies or to have binding activities in common.The best-characterized boundary elements are the scs and scs′ regions found to bracket the 87A7 hsp70 heat shock puff of Drosophila melanogaster polytene chromosomes (33) and a 340-bp fragment from the gypsy retrotransposon (11). The scs/scs′ and the gypsy-derived elements have a boundary function in both of the assays mentioned above. They confer position-independent expression on a bracketed reporter gene by insulating the transgene from both activating and repressive effects at the site of chromosomal integration, and they block communication between a specific enhancer and promoter when interposed (20, 21, 31). It is important to note that boundary elements do not inactivate promoters or enhancers; they only block communication when interposed (2, 3, 21, 32). For instance, if an enhancer and boundary element are located between two divergently transcribed promoters, the enhancer cannot activate the promoter with the intervening boundary element but can activate the other promoter. Thus, the positional functioning of boundary elements is distinct from the bidirectional repressive effect of silencer elements.The boundary activity of the gypsy-derived element is known to be mediated by the binding of the zinc finger protein su(Hw) to its reiterated binding sites (31). The su(Hw) protein has been studied in some detail, and regions involved in DNA binding, enhancer blocking, and interactions with mod(mdg4) have been identified (8, 13, 22). Interactions between the mod(mdg4) gene product and the su(Hw) protein are necessary for boundary function (9). In addition to loss of enhancer blocking, it has been suggested that some mod(mdg4) mutations lead to an unmasked activity that represses certain promoters (3).To address the boundary activity of scs′ at a biochemical level, we previously characterized two cDNAs encoding the related scs′ boundary element-associated factors BEAF-32A and -32B (14, 38). The BEAF activity in Drosophila nuclear extracts appears to be composed predominantly of trimers of one 32A and two 32B subunits. Interactions between BEAF subunits results in cooperative binding to the three CGATA motifs of the high-affinity binding site in scs′ which, in turn, facilitates binding to the lower-affinity binding site located some 200 bp away (14).Evidence of a role for BEAF in boundary activity derives from an enhancer-blocking assay in Drosophila D1 cells: seven tandem copies of a 48-bp oligonucleotide containing the scs′ high-affinity binding site had enhancer-blocking activity (although less than that obtained by using scs′), while point mutations that eliminated BEAF binding further reduced this activity (38). We immunolocalized BEAF to numerous interbands and puff boundaries on polytene chromosomes, suggesting the existence of a common class of boundary elements in Drosophila and that the band-interband structure of polytene chromosomes could be related to the localization of boundary elements.In this study, we isolated some of these genomic BEAF binding sites and used transgenic flies to demonstrate that the newly isolated sequences tested represent boundary elements. The only homology found between these candidate boundary elements (cBEs) and scs′ are clusters of CGATA motifs. Despite the varied spacing and orientations of the motifs in the different clusters, BEAF interacts with all of the clusters. We also used transgenic flies to directly establish the functional importance of BEAF binding sites by mutagenesis of CGATA motifs. This strongly indicates that the hundreds of BEAF binding sites in the Drosophila genome represent an abundant class of boundary elements, providing the first example of a class of binding elements that share a sequence motif and a binding protein.     

Discovery of protein phosphorylation motifs through exploratory data analysis

Background

The need for efficient algorithms to uncover biologically relevant phosphorylation motifs has become very important with rapid expansion of the proteomic sequence database along with a plethora of new information on phosphorylation sites. Here we present a novel unsupervised method, called Motif Finder (in short, F-Motif) for identification of phosphorylation motifs. F-Motif uses clustering of sequence information represented by numerical features that exploit the statistical information hidden in some foreground data. Furthermore, these identified motifs are then filtered to find “actual” motifs with statistically significant motif scores.

Results and Discussion

We have applied F-Motif to several new and existing data sets and compared its performance with two well known state-of-the-art methods. In almost all cases F-Motif could identify all statistically significant motifs extracted by the state-of-the-art methods. More importantly, in addition to this, F-Motif uncovers several novel motifs. We have demonstrated using clues from the literature that most of these new motifs discovered by F-Motif are indeed novel. We have also found some interesting phenomena. For example, for CK2 kinase, the conserved sites appear only on the right side of S. However, for CDK kinase, the adjacent site on the right of S is conserved with residue P. In addition, three different encoding methods, including a novel position contrast matrix (PCM) and the simplest binary coding, are used and the ability of F-motif to discover motifs remains quite robust with respect to encoding schemes.

Conclusions

An iterative algorithm proposed here uses exploratory data analysis to discover motifs from phosphorylated data. The effectiveness of F-Motif has been demonstrated using several real data sets as well as using a synthetic data set. The method is quite general in nature and can be used to find other types of motifs also. We have also provided a server for F-Motif at http://f-motif.classcloud.org/, http://bio.classcloud.org/f-motif/ or http://ymu.classcloud.org/f-motif/.     

Herbicidal activity of the antibiotics geldanamycin and nigericin

Geldanamycin and nigericin, phytotoxic metabolites from a strain ofStreptomyces hygroscopicus, were tested for herbicidal activity and selectivity on a range of crop and weed species. In petri dish bioassays, geldanamycin reduced radicle growth of all species tested, whereas nigericin inhibited 7 of 10. The two compounds in mixture appeared to be additive rather than synergistic in effect. In assays with seeds and seedlings in field soil, geldanamycin showed significant preemergence activity on proso millet, barnyardgrass, garden cress, and giant foxtail. It had no postemergence herbicidal effect on any of the species tested. Nigericin had preemergence activity on garden cress and large crabgrass and postemergence activity on garden cress and velvetleaf. The postemergence effect of nigericin on velvetleaf was especially striking, with leaves showing symptoms of injury within 24–48 h of treatment. Doses as low as 0.3 kg/ha caused damage. The primary herbicidal effect of both compounds was slowing of seed germination or seedling growth, although some plants were killed, especially at higher rates of application. Herbicidal effects were most pronounced for 1 to 2 weeks after treatment and diminished thereafter.     

Development and characterization of 79 nuclear markers amplifying in viviparous and oviparous clades of the European common lizard

The European common lizard (Zootoca vivipara) is a widely distributed species across Europe and Asia exhibiting two reproductive modes (oviparity/viviparity), six major lineages and several sublineages. It has been used to tackle a large variety of research questions, nevertheless, few nuclear DNA sequence markers have been developed for this species. Here we developed 79 new nuclear DNA sequence markers using a clonation protocol. These markers were amplified in several oviparous and viviparous specimens including samples of all extant clades, to test the amplification success and their diversity. 49.4% of the markers were polymorphic and of those, 51.3% amplified in all and 94.9% amplified in 5–7 of the extant Z. vivipara clades. These new markers will be very useful for the study of the population structure, population dynamics, and micro/macro evolution of Z. vivipara. Cross-species amplification in four lizard species (Psammodromus edwardsianus, Podarcis muralis, Lacerta bilineata, and Takydromus sexlineatus) was positive in several of the markers, and six makers amplified in all five species. The large genetic distance between P. edwardsianus and Z. vivipara further suggests that these markers may as well be employed in many other species.