Conserved water-mediated H-bonding dynamics of catalytic Asn 175 in plant thiol protease

The role of invariant water molecules in the activity of plant cysteine protease is ubiquitous in nature. On analysing the 11 different Protein DataBank (PDB) structures of plant thiol proteases, the two invariant water molecules W I and W2 (W220 and W222 in the template 1PPN structure) were observed to form H-bonds with the Ob atom of Asn 175. Extensive energy minimization and molecular dynamics simulation studies up to 2 ns on all the PDB and solvated structures clearly revealed the involvement of the H-bonding association of the two water molecules in fixing the orientation of the asparagine residue of the catalytic triad. From this study, it is suggested that H-bonding of the water molecule at the W1 invariant site better stabilizes the Asn residue at the active site of the catalytic triad.

sRNAscanner: A Computational Tool for Intergenic Small RNA Detection in Bacterial Genomes

Background:Bacterial non-coding small RNAs (sRNAs) have attracted considerable attention due to their ubiquitous nature and contribution to numerous cellular processes including survival, adaptation and pathogenesis. Existing computational approaches for identifying bacterial sRNAs demonstrate varying levels of success and there remains considerable room for improvement. Methodology/Principal Findings: Here we have proposed a transcriptional signal-based computational method to identify intergenic sRNA transcriptional units (TUs) in completely sequenced bacterial genomes. Our sRNAscanner tool uses position weight matrices derived from experimentally defined E. coli K-12 MG1655 sRNA promoter and rho-independent terminator signals to identify intergenic sRNA TUs through sliding window based genome scans. Analysis of genomes representative of twelve species suggested that sRNAscanner demonstrated equivalent sensitivity to sRNAPredict2, the best performing bioinformatics tool available presently. However, each algorithm yielded substantial numbers of known and uncharacterized hits that were unique to one or the other tool only. sRNAscanner identified 118 novel putative intergenic sRNA genes in Salmonella enterica Typhimurium LT2, none of which were flagged by sRNAPredict2. Candidate sRNA locations were compared with available deep sequencing libraries derived from Hfq-co-immunoprecipitated RNA purified from a second Typhimurium strain (Sittka et al. (2008) PLoS Genetics 4: e1000163). Sixteen potential novel sRNAs computationally predicted and detected in deep sequencing libraries were selected for experimental validation by Northern analysis using total RNA isolated from bacteria grown under eleven different growth conditions. RNA bands of expected sizes were detected in Northern blots for six of the examined candidates. Furthermore, the 5'-ends of these six Northern-supported sRNA candidates were successfully mapped using 5'-RACE analysis. Conclusions/Significance: We have developed, computationally examined and experimentally validated the sRNAscanner algorithm. Data derived from this study has successfully identified six novel S. Typhimurium sRNA genes. In addition, the computational specificity analysis we have undertaken suggests that similar to 40% of sRNAscanner hits with high cumulative sum of scores represent genuine, undiscovered sRNA genes. Collectively, these data strongly support the utility of sRNAscanner and offer a glimpse of its potential to reveal large numbers of sRNA genes that have to date defied identification. sRNAscanner is available from: http://bicmku.in:8081/sRNAscanner or http://cluster.physics.iisc.ernet.in/sRNAscanner/.

Reducing the babel in plant volatile communication: using the forest to see the trees

While plants of a single species emit a diversity of volatile organic compounds (VOCs) to attract or repel interacting organisms, these specific messages may be lost in the midst of the hundreds of VOCs produced by sympatric plants of different species, many of which may have no signal content. Receivers must be able to reduce the babel or noise in these VOCs in order to correctly identify the message. For chemical ecologists faced with vast amounts of data on volatile signatures of plants in different ecological contexts, it is imperative to employ accurate methods of classifying messages, so that suitable bioassays may then be designed to understand message content. We demonstrate the utility of `Random Forests' (RF), a machine-learning algorithm, for the task of classifying volatile signatures and choosing the minimum set of volatiles for accurate discrimination, using datam from sympatric Ficus species as a case study. We demonstrate the advantages of RF over conventional classification methods such as principal component analysis (PCA), as well as data-mining algorithms such as support vector machines (SVM), diagonal linear discriminant analysis (DLDA) and k-nearest neighbour (KNN) analysis. We show why a tree-building method such as RF, which is increasingly being used by the bioinformatics, food technology and medical community, is particularly advantageous for the study of plant communication using volatiles, dealing, as it must, with abundant noise.

DNA Free Energy-Based Promoter Prediction and Comparative Analysis of Arabidopsis and Rice Genomes

The cis-regulatory regions on DNA serve as binding sites for proteins such as transcription factors and RNA polymerase. The combinatorial interaction of these proteins plays a crucial role in transcription initiation, which is an important point of control in the regulation of gene expression. We present here an analysis of the performance of an in silico method for predicting cis-regulatory regions in the plant genomes of Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) on the basis of free energy of DNA melting. For protein-coding genes, we achieve recall and precision of 96% and 42% for Arabidopsis and 97% and 31% for rice, respectively. For noncoding RNA genes, the program gives recall and precision of 94% and 75% for Arabidopsis and 95% and 90% for rice, respectively. Moreover, 96% of the false-positive predictions were located in noncoding regions of primary transcripts, out of which 20% were found in the first intron alone, indicating possible regulatory roles. The predictions for orthologous genes from the two genomes showed a good correlation with respect to prediction scores and promoter organization. Comparison of our results with an existing program for promoter prediction in plant genomes indicates that our method shows improved prediction capability.

Hybrid and Rogue Kinases Encoded in the Genomes of Model Eukaryotes

The highly modular nature of protein kinases generates diverse functional roles mediated by evolutionary events such as domain recombination, insertion and deletion of domains. Usually domain architecture of a kinase is related to the subfamily to which the kinase catalytic domain belongs. However outlier kinases with unusual domain architectures serve in the expansion of the functional space of the protein kinase family. For example, Src kinases are made-up of SH2 and SH3 domains in addition to the kinase catalytic domain. A kinase which lacks these two domains but retains sequence characteristics within the kinase catalytic domain is an outlier that is likely to have modes of regulation different from classical src kinases. This study defines two types of outlier kinases: hybrids and rogues depending on the nature of domain recombination. Hybrid kinases are those where the catalytic kinase domain belongs to a kinase subfamily but the domain architecture is typical of another kinase subfamily. Rogue kinases are those with kinase catalytic domain characteristic of a kinase subfamily but the domain architecture is typical of neither that subfamily nor any other kinase subfamily. This report provides a consolidated set of such hybrid and rogue kinases gleaned from six eukaryotic genomes-S. cerevisiae, D. melanogaster, C. elegans, M. musculus, T. rubripes and H. sapiens-and discusses their functions. The presence of such kinases necessitates a revisiting of the classification scheme of the protein kinase family using full length sequences apart from classical classification using solely the sequences of kinase catalytic domains. The study of these kinases provides a good insight in engineering signalling pathways for a desired output. Lastly, identification of hybrids and rogues in pathogenic protozoa such as P. falciparum sheds light on possible strategies in host-pathogen interactions.

Typical and atypical domain combinations in human protein kinases: functions, disease causing mutations and conservation in other primates

Ser/Thr and Tyr protein kinases orchestrate many signalling pathways and hence loss in this balance leads to many disease phenotypes. Due to their high abundance, diversity and importance, efforts have been made in the past to classify kinases and annotate their functions at both gross and fine levels. These kinases are conventionally classified into subfamilies based on the sequences of catalytic domains. Usually the domain architecture of a full-length kinase is consistent with the subfamily classification made based on the sequence of kinase domain. Important contributions of modular domains to the overall function of the kinase are well known. Recently occurrence of two kinds of outlier kinases-''Hybrid'' and ``Rogue'' has been reported. These show considerable deviations in their domain architectures from the typical domain architecture known for the classical kinase subfamilies. This article provides an overview of the different subfamilies of human kinases and the role of non-kinase domains in functions and diseases. Importantly this article provides analysis of hybrid and rogue kinases encoded in the human genome and highlights their conservation in closely related primate species. These kinases are examples of elegant rewiring to bring about subtle functional differences compared to canonical variants.

Evolution, Homology Conservation, and Identification of Unique Sequence Signatures in GH19 Family Chitinases

The discovery of GH (Glycoside Hydrolase) 19 chitinases in Streptomyces sp. raises the possibility of the presence of these proteins in other bacterial species, since they were initially thought to be confined to higher plants. The present study mainly concentrates on the phylogenetic distribution and homology conservation in GH19 family chitinases. Extensive database searches are performed to identify the presence of GH19 family chitinases in the three major super kingdoms of life. Multiple sequence alignment of all the identified GH19 chitinase family members resulted in the identification of globally conserved residues. We further identified conserved sequence motifs across the major sub groups within the family. Estimation of evolutionary distance between the various bacterial and plant chitinases are carried out to better understand the pattern of evolution. Our study also supports the horizontal gene transfer theory, which states that GH19 chitinase genes are transferred from higher plants to bacteria. Further, the present study sheds light on the phylogenetic distribution and identifies unique sequence signatures that define GH19 chitinase family of proteins. The identified motifs could be used as markers to delineate uncharacterized GH19 family chitinases. The estimation of evolutionary distance between chitinase identified in plants and bacteria shows that the flowering plants are more related to chitinase in actinobacteria than that of identified in purple bacteria. We propose a model to elucidate the natural history of GH19 family chitinases.

Hyper-activation of the TCP4 transcription factor in Arabidopsis thaliana accelerates multiple aspects of plant maturation

Plant organs are initiated as primordial outgrowths, and require controlled cell division and differentiation to achieve their final size and shape. Superimposed on this is another developmental program that orchestrates the switch from vegetative to reproductive to senescence stages in the life cycle. These require sequential function of heterochronic regulators. Little is known regarding the coordination between organ and organismal growth in plants. The TCP gene family encodes transcription factors that control diverse developmental traits, and a subgroup of class II TCP genes regulate leaf morphogenesis. Absence of these genes results in large, crinkly leaves due to excess division, mainly at margins. It has been suggested that these class II TCPs modulate the spatio-temporal control of differentiation in a growing leaf, rather than regulating cell proliferation per se. However, the link between class II TCP action and cell growth has not been established. As loss-of-function mutants of individual TCP genes in Arabidopsis are not very informative due to gene redundancy, we generated a transgenic line that expressed a hyper-activated form of TCP4 in its endogenous expression domain. This resulted in premature onset of maturation and decreased cell proliferation, leading to much smaller leaves, with cup-shaped lamina in extreme cases. Further, the transgenic line initiated leaves faster than wild-type and underwent precocious reproductive maturation due to a shortened adult vegetative phase. Early senescence and severe fertility defects were also observed. Thus, hyper-activation of TCP4 revealed its role in determining the timing of crucial developmental events, both at the organ and organism level.

Characterization of structural and free energy properties of promoters associated with Primary and Operon TSS in Helicobacter pylori genome and their orthologs

Promoter regions in the genomes of all domains of life show similar trends in several structural properties such as stability, bendability, curvature, etc. In current study we analysed the stability and bendability of various classes of promoter regions (based on the recent identification of different classes of transcription start sites) of Helicobacter pylori 26695 strain. It is found that primary TSS and operon-associated TSS promoters show significantly strong features in their promoter regions. DNA free-energy-based promoter prediction tool PromPredict was used to annotate promoters of different classes, and very high recall values (similar to 80%) are obtained for primary TSS. Orthologous genes from other strains of H. pylori show conservation of structural properties in promoter regions as well as coding regions. PromPredict annotates promoters of orthologous genes with very high recall and precision.

A multi-level multi-scale approach to study essential genes in Mycobacterium tuberculosis

Background: The set of indispensable genes that are required by an organism to grow and sustain life are termed as essential genes. There is a strong interest in identification of the set of essential genes, particularly in pathogens, not only for a better understanding of the pathogen biology, but also for identifying drug targets and the minimal gene set for the organism. Essentiality is inherently a systems property and requires consideration of the system as a whole for their identification. The available experimental approaches capture some aspects but each method comes with its own limitations. Moreover, they do not explain the basis for essentiality in most cases. A powerful prediction method to recognize this gene pool including rationalization of the known essential genes in a given organism would be very useful. Here we describe a multi-level multi-scale approach to identify the essential gene pool in a deadly pathogen, Mycobacterium tuberculosis. Results: The multi-level workflow analyses the bacterial cell by studying (a) genome-wide gene expression profiles to identify the set of genes which show consistent and significant levels of expression in multiple samples of the same condition, (b) indispensability for growth by using gene expression integrated flux balance analysis of a genome-scale metabolic model, (c) importance for maintaining the integrity and flow in a protein-protein interaction network and (d) evolutionary conservation in a set of genomes of the same ecological niche. In the gene pool identified, the functional basis for essentiality has been addressed by studying residue level conservation and the sub-structure at the ligand binding pockets, from which essential amino acid residues in that pocket have also been identified. 283 genes were identified as essential genes with high-confidence. An agreement of about 73.5% is observed with that obtained from the experimental transposon mutagenesis technique. A large proportion of the identified genes belong to the class of intermediary metabolism and respiration. Conclusions: The multi-scale, multi-level approach described can be generally applied to other pathogens as well. The essential gene pool identified form a basis for designing experiments to probe their finer functional roles and also serve as a ready shortlist for identifying drug targets.

Contribution of forest fire ash and plant litter decay on stream dissolved composition in a sub-humid tropical watershed (Mule Hole, Southern India)

The current understanding of wildfire effects on water chemistry is limited by the quantification of the elemental dissolution rates from ash and element release rate from the plant litter, as well as quantification of the specific ash contribution to stream water chemistry. The main objective of the study was to provide such knowledge through combination of experimental modelling, field data and end-member mixing analysis (EMMA) of wildfire impact on a watershed scale. The study concerns watershed effects of fire in the Indian subcontinent, a region that is typically not well represented in the fire science literature. In plant litter ash, major elements are either hosted in readily-soluble phases (K, Mg) such as salts, carbonates and oxides or in less-soluble carrier-phases (Si, Ca) such as amorphous silica, quartz and calcite. Accordingly, elemental release rates, inferred from ash leaching experiments in batch reactor, indicated that the element release into solution followed the order K > Mg > Na > Si > Ca. Experiments on plant litter leaching in mixed-flow reactor indicated two dissolution regimes: rapid, over the week and slower over the month. The mean dissolution rates at steady-state (R-ss) indicated that the release of major elements from plant litter followed the order Ca > Si > Cl > Mg > K > Na. R-ss for Si and Ca for tree leaves and herbaceous species are similar to those reported for boreal and European tree species and are higher than that from the dissolution of soil clay minerals. This identifies tropical plant litters as important source of Si and Ca for tropical surface waters. In the wildfire-impacted year 2004, the EMMA indicated that the streamflow composition (Ca, K, Mg, Na, Si, Cl) was controlled by four main sources: rainwater, throughfall, ash leaching and soil solution. The influence of the ash end-member was maximal early in the rainy season (the two first storm events) and decreased later in the rainy season, when the stream was dominated by the throughfall end-member. The contribution of plant litter decay to the streamwater composition for a year not impacted by wildfire is significant with estimated solute fluxes originating from this decay greatly exceed, for most major elements, the annual elemental dissolved fluxes at the Mule Hole watershed outlet. This highlighted the importance of solute retention and vegetation back uptake processes within the soil profile. Overall, the fire increased the mobility and export of major elements from the soils to the stream. It also shifted the vegetation-related contribution to the elemental fluxes at the watershed outlet from long-term (seasonal) to short-term (daily to monthly). (C) 2014 Elsevier B.V. All rights reserved.

Distribution, Demography, and Conservation of Lion-tailed Macaques (Macaca silenus) in the Anamalai Hills Landscape, Western Ghats, India

The status of the endemic and endangered lion-tailed macaque (Macaca silenus) has not been properly assessed in several regions of the Western Ghats of southern India. We conducted a study in Parambikulam Forest Reserve in the state of Kerala to determine the distribution, demography, and status of lion-tailed macaques. We laid 5km(2) grid cells on the map of the study area (644km(2)) and made four replicated walks in each grid cell using GPS. We gathered data on lion-tailed macaque group locations, demography, and site covariates including trail length, duration of walk, proportion of evergreen forest, height of tallest trees, and human disturbance index. We also performed occupancy modeling using PRESENCE ver. 3.0. We estimated a minimum of 17 groups of macaques in these hills. Low detection and occupancy probabilities indicated a low density of lion-tailed macaques in the study area. Height of the tallest trees correlated positively whereas human disturbance and proportion of evergreen forest correlated negatively with occupancy in grid cells. We also used data from earlier studies carried out in the surrounding Anamalai Tiger Reserve and Nelliyampathy Hills to discuss the conservation status in the large Anamalai Hills Landscape. This landscape harbors an estimated population of 1108 individuals of lion-tailed macaques, which is about one third of the entire estimated wild population of this species. A conservation plan for this landscape could be used as a model for conservation in other regions of the Western Ghats.

Distribution, relative abundance, and conservation status of Asian elephants in Karnataka, southern India

Karnataka state in southern India supports a globally significant and the country's largest population of the Asian elephant Elephas maximus. A reliable map of Asian elephant distribution and measures of spatial variation in their abundance, both vital needs for conservation and management action, are unavailable not only in Karnataka, but across its global range. Here, we use various data gathered between 2000 and 2015 to map the distribution of elephants in Karnataka at the scale of the smallest forest management unit, the `beat', while also presenting data on elephant dung density for a subset of `elephant beats.' Elephants occurred in 972 out of 2855 forest beats of Karnataka. Sixty percent of these 972 beats and 55% of the forest habitat lay outside notified protected areas (PM), and included lands designated for agricultural production and human dwelling. While median elephant dung density inside protected areas was nearly thrice as much as outside, elephants routinely occurred in or used habitats outside PM where human density, land fraction under cultivation, and the interface between human-dominated areas and forests were greater. Based on our data, it is clear that India's framework for elephant conservation which legally protects the species wherever it occurs, but protects only some of its habitats while being appropriate in furthering their conservation within PM, seriously falters in situations where elephants reside in and/or seasonally use areas outside PAs. Attempts to further elephant conservation in production and dwelling areas have extracted high costs in human, elephant, material and monetary terms in Karnataka. In such settings, conservation planning exercises are necessary to determine where the needs of elephants or humans must take priority over the other, and to achieve that in a manner that is based not only on reliable scientific data but also on a process of public reasoning. (C) 2015 Elsevier Ltd. All rights reserved.

Variations in the levels of aminoacylation and modified nucleotide content between total tRNAs from chloroplasts and etioplasts in cucumber cotyledons

Total tRNAs isolated from chloroplasts and etioplasts of cucumber cotyledons were compared with respect toamino acid acceptance, isoacceptor distribution and extent of modification. Aminoacylation of the tRNAs with nine different amino acids studied indicated that the relative acceptor activities of chloroplast total tRNAs for four amino acids are significantly higher than etioplast total tRNAs. Two dimensional polyacrylamide gel electrophoresis(2D-PAGE) of chloroplast total tRNAs separated at least 32 spots, while approximately 41 spots were resolved from etioplast total tRNAs. Comparison of the reversed-phase chromatography (RPC-5) profiles of chloroplast and etioplast leucyl-, lysyl-, phenylalanyl-, and valyl-tRNA species showed no qualitative differences in the elution profiles. However, leucyl-, lysyl- and valyl-tRNA species showed quantitative differences in the relative amounts of the isoaccepting species present in chloroplasts and etioplasts. The analysis of modified nucleotides of total tRNAs from the two plastid types indicated that total tRNA from etioplasts was undermodified with respect to ribothymidine, isopentenyladenosine/hydroxy-isopentenyladenosine, 1 -methylguanosine and 2-o-methylguanosine. This indicates that illumination may cause de novo synthesis of chloroplast tRNAmodifying enzymes encoded for by nuclear genes leading to the formation of highly modified tRNAs in chloroplasts. Based on these results, we speculate that the observed decrease in levels of aminoacylation, variations in the relative amounts of certain isoacceptors, and differences in the electrophoretic mobilities of some extra tRNA spots in the etioplast total tRNAs as compared to chloroplast total tRNAs could be due to some partially undermodified etioplast tRNAs. Taken together, the data suggested that the light-induced transformation of etioplasts into chloroplasts is accompanied by increases in the relative levels of some functional chloroplast tRNAs by post transcriptional nucleotide modifications.