Budheswar Dehury

In silico identification and characterization of conserved miRNAs and their target genes in sweet potato (Ipomoea batatas L.) Expressed Sequence Tags (ESTs)

The endogenous small non-coding micro RNAs (miRNAs), which are typically ~21–24 nt nucleotides, play a crucial role in regulating the intrinsic normal growth of cells and development of the plants as well as in maintaining the integrity of genomes. These small non-coding RNAs function as the universal specificity factors in post-transcriptional gene silencing. Discovering miRNAs, identifying their targets, and further inferring miRNA functions is a routine process to understand normal biological processes of miRNAs and their roles in the development of plants. Comparative genomics based approach using expressed sequence tags (EST) and genome survey sequences (GSS) offer a cost-effective platform for identification and characterization of miRNAs and their target genes in plants. Despite the fact that sweet potato (Ipomoea batatas L.) is an important staple food source for poor small farmers throughout the world, the role of miRNA in various developmental processes remains largely unknown. In this paper, we report the computational identification of miRNAs and their target genes in sweet potato from their ESTs. Using comparative genomics-based approach, 8 potential miRNA candidates belonging to miR168, miR2911, and miR156 families were identified from 23 406 ESTs in sweet potato. A total of 42 target genes were predicted and their probable functions were illustrated. Most of the newly identified miRNAs target transcription factors as well as genes involved in plant growth and development, signal transduction, metabolism, defense, and stress response. The identification of miRNAs and their targets is expected to accelerate the pace of miRNA discovery, leading to an improved understanding of the role of miRNA in development and physiology of sweet potato, as well as stress response.

Computational identification and characterization of conserved miRNAs and their target genes in garlic (Allium sativum L.) expressed sequence tags.

The endogenous small non-coding functional microRNAs (miRNAs) are short in size, range from ~21 to 24 nucleotides in length, play a pivotal role in gene expression in plants and animals by silencing genes either by destructing or blocking of translation of homologous mRNA. Although various high-throughput, time consuming and expensive techniques like forward genetics and direct cloning are employed to detect miRNAs in plants but comparative genomics complemented with novel bioinformatic tools pave the way for efficient and cost-effective identification of miRNAs through homologous sequence search with previously known miRNAs. In this study, an attempt was made to identify and characterize conserved miRNAs in garlic expressed sequence tags (ESTs) through computational means. For identification of novel miRNAs in garlic, a total 3227 known mature miRNAs of plant kingdom Viridiplantae were searched for homology against 21,637 EST sequences resulting in identification of 6 potential miRNA candidates belonging to 6 different miRNA families. The psRNATarget server predicted 33 potential target genes and their probable functions for the six identified miRNA families in garlic. Most of the garlic miRNA target genes seem to encode transcription factors as well as genes involved in stress response, metabolism, plant growth and development. The results from the present study will shed more light on the understanding of molecular mechanisms of miRNA in garlic which may aid in the development of novel and precise techniques to understand some post-transcriptional gene silencing mechanism in response to stress tolerance.

Mining for SSRs and FDMs from expressed sequence tags of Camellia sinensis

Simple Sequence Repeats (SSRs) developed from Expressed Sequence Tags (ESTs), known as EST-SSRs are most widely used and potentially valuable source of gene based markers for their high levels of crosstaxon portability, rapid and less expensive development. The EST sequence information in the publicly available databases is increasing in a faster rate. The emerging computational approach provides a better alternative process of development of SSR markers from the ESTs than the conventional methods. In the present study, 12,851 EST sequences of Camellia sinensis, downloaded from National Center for Biotechnology Information (NCBI) were mined for the development of Microsatellites. 6148 (4779 singletons and 1369 contigs) non redundant EST sequences were found after preprocessing and assembly of these sequences using various computational tools. Out of total 3822.68 kb sequence examined, 1636 (26.61%) EST sequences containing 2371 SSRs were detected with a density of 1 SSR/1.61 kb leading to development of 245 primer pairs. These mined EST-SSR markers will help further in the study of variability, mapping, evolutionary relationship in Camellia sinensis. In addition, these developed SSRs can also be applied for various studies across species.

In silico analyses of superoxide dismutases (SODs) of rice (Oryza sativa L.)

Superoxide dismutases (SODs), members of the metalloenzymes family are most effective intracellular enzymatic antioxidant in aerobic organisms. These enzymes provide the first line of defense in plants against the toxic effects of elevated levels of reactive oxygen species (ROS) generated during various environmental stresses. The availability of high-throughput computational tools has provided better opportunities to characterize the protein features and determine their function. In the present study an attempt was made to gain an insight into the structure and evolution of subunits of SODs (Cu-Zn, Mn and Fe SODs) of rice. The 3-Dimensional structures of SODs were modeled based on available X-ray crystal structures and further validated. The primary sequence, secondary and tertiary structure analysis revealed Mn and Fe SOD to be structurally homologous while Cu-Zn SOD is un-related to either of them. Comparative structural study also revealed former two were dominated by α-helices followed by β-strands in contrast; Cu-Zn SOD dominated by β-strands. Molecular phylogeny indicated a common evolutionary origin of Mn and Fe SOD while Cu-Zn SOD may have evolved separately.

Mining the proteome of Leishmania donovani for the development of novel MHC class I restricted epitope for the control of visceral leishmaniasis

Although, the precise host defence mechanism(s) is not completely understood, T cell‐mediated immune responses is believed to play a pivotal role in controlling parasite infection. Here we target the stage dependent over expressed gene. Here, the consensus based computational approach was adopted for the screening of potential major histocompatibility complex class I restricted epitopes. Based on the computational analysis and previously published report, a set 19 antigenic proteins derived from Leishmania donovani were screened for further characterization as vaccine candidates. A total of 49 epitopes were predicted, which revealed a comprehensive binding affinity to the 40 different MHC class I supertypes. Based on the population coverage and HLA cross presentation, nine highly promiscuous epitopes such as LTYDDVWTV (P1), FLFPQRTAL(P2), FLFSNGAVV (P3), YIYNFGIRV (P4), YMTAAFAAL (P5), KLLRPFAPL (P6), FMLGWIVTI (P7), SLFERNKRV (P8), and SVWNRIFTL (P9) which have either a high or an intermediate TAP binding affinity were selected for further analysis. Theoretical population coverage analysis of polytope vaccine (P1–P9) revealed more than 92% population. Stimulation with the cocktail of peptide revealed a proliferative CD8+ T cell response and increased IFN‐γ production. An upregulated NF‐κB activity is thought to be play a pivotal role in T cell proliferation against the selected peptide. The Th1‐type cytokine profile (presence of IFN‐γ and absence of IL‐10) suggests the potentiality of the cocktail of epitope as a subunit vaccine against leishmaniasis. However, the efficiency of these epitopes to trigger other Th1 cytokines and chemokines in a humanized mice model could explore its plausibility as a vaccine candidate. J. Cell. Biochem. 119: 378–391, 2018.

Structural Models of Zebrafish (Danio rerio) NOD1 and NOD2 NACHT Domains Suggest Differential ATP Binding Orientations: Insights from Computational Modeling, Docking and Molecular Dynamics Simulations

Nucleotide-binding oligomerization domain-containing protein 1 (NOD1) and NOD2 are cytosolic pattern recognition receptors playing pivotal roles in innate immune signaling. NOD1 and NOD2 recognize bacterial peptidoglycan derivatives iE-DAP and MDP, respectively and undergoes conformational alternation and ATP-dependent self-oligomerization of NACHT domain followed by downstream signaling. Lack of structural adequacy of NACHT domain confines our understanding about the NOD-mediated signaling mechanism. Here, we predicted the structure of NACHT domain of both NOD1 and NOD2 from model organism zebrafish (Danio rerio) using computational methods. Our study highlighted the differential ATP binding modes in NOD1 and NOD2. In NOD1, γ-phosphate of ATP faced toward the central nucleotide binding cavity like NLRC4, whereas in NOD2 the cavity was occupied by adenine moiety. The conserved ‘Lysine’ at Walker A formed hydrogen bonds (H-bonds) and Aspartic acid (Walker B) formed electrostatic interaction with ATP. At Sensor 1, Arg328 of NOD1 exhibited an H-bond with ATP, whereas corresponding Arg404 of NOD2 did not. ‘Proline’ of GxP motif (Pro386 of NOD1 and Pro464 of NOD2) interacted with adenine moiety and His511 at Sensor 2 of NOD1 interacted with γ-phosphate group of ATP. In contrast, His579 of NOD2 interacted with the adenine moiety having a relatively inverted orientation. Our findings are well supplemented with the molecular interaction of ATP with NLRC4, and consistent with mutagenesis data reported for human, which indicates evolutionary shared NOD signaling mechanism. Together, this study provides novel insights into ATP binding mechanism, and highlights the differential ATP binding modes in zebrafish NOD1 and NOD2.

Structural and dynamic investigation of bovine folate receptor alpha (FOLR1), and role of ultra-high temperature processing on conformational and thermodynamic characteristics of FOLR1–folate complex

The folate receptor alpha (FOLR1) present in milk has widely been studied to investigate the effects of pasteurization, ultra-high temperature (UHT) processing and fermentation on net folate concentration. However, the folate binding mechanism with FOLR1, and effect of temperature on FOLR1-folate complex is poorly explored till now in bovine milk which is a chief resource of folate. Despite of enormous importance of folic acid and the routine intake of bovine milk, folic acid deficiency diseases are common in human race. To understand the folate deficiency in milk after processing, in absence of experimental structure, 3D model of bovine FOLR1 (bvFOLR1) was built followed by 40ns molecular dynamics (MD) simulation. The folate and its derivatives binding sites in bvFOLR1 were anticipated by molecular docking using AutoDock 4.2. Essential MD studies suggested the presence of a longer signal peptide (22 residues) and a short propeptide (7 residues) at the C-terminus that may cleaved during post-translational modification. MD analysis of bvFOLR1-folate complex at 298, 323, 353, 373 and 408K followed by binding energy (BE) calculation showed maximum binding affinity at ∼353K. However, at 373K and UHT (408K), the folate BE is significantly decreased with substantial conformational alteration. Heating at UHT followed by cooling within 298-408K range demoed no structural reformation with temperature reduction, and the folate was displaced from the active site. This study presented the disintegration of folate from bvFOLR1 during high temperature processing and revealed a lower folate concentration in UHT milk and dairy products.

Molecular recognition of avirulence protein (avrxa5) by eukaryotic transcription factor xa5 of rice (Oryza sativa L.): Insights from molecular dynamics simulations

The avirulence gene avrxa5 of bacterial blight pathogen Xanthomonas oryzae pv. oryzae (Xoo) recognized by the resistant rice lines having corresponding resistance (xa5) gene in a gene-for-gene manner. We used a combinatorial approach involving protein-protein docking, molecular dynamics (MD) simulations and binding free energy calculations to gain novel insights into the gene-for-gene mechanism that governs the direct interaction of R-Avr protein. From the best three binding poses predicted by molecular docking, MD simulations were performed to explore the dynamic binding mechanism of xa5 and avrxa5. Molecular Mechanics/Poisson Boltzmann Surface Area (MM/PBSA) techniques were employed to calculate the binding free energy and to uncover the thriving force behind the molecular recognition of avrxa5 by eukaryotic transcription factor xa5. Binding free energy analysis revealed van der Waals term as the most constructive component that favors the xa5 and avrxa5 interaction. In addition, hydrogen bonds (H-bonds) and essential electrostatic interactions analysis highlighted amino acid residues Lys54/Asp870, Lys56/Ala868, Lys56/Ala866, Lys56/Glu871, Ile59/His862, Gly61/Phe858, His62/Arg841, His62/Leu856, Ser101/Ala872 and Ser105/Asp870 plays pivotal role for the energetically stability of the R-Avr complex. Insights gained from the present study are expected to unveil the molecular mechanisms that define the transcriptional activator mediated transcriptome modification in host plants.

Rediscovering Medicinal Plants' Potential with OMICS: Microsatellite Survey in Expressed Sequence Tags of Eleven Traditional Plants with Potent Antidiabetic Properties

Herbal medicines and traditionally used medicinal plants present an untapped potential for novel molecular target discovery using systems science and OMICS biotechnology driven strategies. Since up to 40% of the worlds poor people have no access to government health services, traditional and folk medicines are often the only therapeutics available to them. In this vein, North East (NE) India is recognized for its rich bioresources. As part of the Indo-Burma hotspot, it is regarded as an epicenter of biodiversity for several plants having myriad traditional uses, including medicinal use. However, the improvement of these valuable bioresources through molecular breeding strategies, for example, using genic microsatellites or Simple Sequence Repeats (SSRs) or Expressed Sequence Tags (ESTs)-derived SSRs has not been fully utilized in large scale to date. In this study, we identified a total of 47,700 microsatellites from 109,609 ESTs of 11 medicinal plants (pineapple, papaya, noyontara, bitter orange, bermuda brass, ratalu, barbados nut, mango, mulberry, lotus, and guduchi) having proven antidiabetic properties. A total of 58,159 primer pairs were designed for the non-redundant 8060 SSR-positive ESTs and putative functions were assigned to 4483 unique contigs. Among the identified microsatellites, excluding mononucleotide repeats, di-/trinucleotides are predominant, among which repeat motifs of AG/CT and AAG/CTT were most abundant. Similarity search of SSR containing ESTs and antidiabetic gene sequences revealed 11 microsatellites linked to antidiabetic genes in five plants. GO term enrichment analysis revealed a total of 80 enriched GO terms widely distributed in 53 biological processes, 17 molecular functions, and 10 cellular components associated with the 11 markers. The present study therefore provides concrete insights into the frequency and distribution of SSRs in important medicinal resources. The microsatellite markers reported here markedly add to the genetic stock for cross transferability in these plants and the literature on biomarkers and novel drug discovery for common chronic diseases such as diabetes.

A comparative proteomic approach to analyse structure, function and evolution of rice chitinases: a step towards increasing plant fungal resistance

Glycoside hydrolase family 19 chitinases (EC 3.2.1.14) widely distributed in plants, bacteria and viruses catalyse the hydrolysis of chitin and play a major role in plant defense mechanisms and development. Rice possesses several classes of chitinase, out of which a single structure of class I has been reported in PDB to date. In the present study an attempt was made to gain more insight into the structure, function and evolution of class I, II and IV chitinases of GH family 19 from rice. The three-dimensional structures of chitinases were modelled and validated based on available X-ray crystal structures. The structural study revealed that they are highly α-helical and bilobed in nature. These enzymes are single or multi domain and multi-functional in which chitin-binding domain (CBD) and catalytic domain (CatD) are present in class I and IV whereas class II lacks CBD. The CatD possesses a catalytic triad which is thought to be involved in catalytic process. Loop III, which is common in all three classes of chitinases, reflects that it may play a significant role in their function. Our study also confirms that the absence and presence of different loops in GH family 19 of rice may be responsible for various sized products. Molecular phylogeny revealed chitinases in monocotyledons and dicotyledons differed from each other forming two different clusters and may have evolved differentially. More structural study of this enzyme from different plants is required to enhance the knowledge of catalytic mechanism and substrate binding.