Biomolecular characterization of Antheraea mylitta cocoonase: A secreted protease

Cocoonase is a protease secreted during the emergence of silk moths. In the present study cocoonase of Antheraea mylitta was collected, purified and secondary structure was determined using circular dichroism (CD) spectroscopy which revealed the presence of α-helix 4.3%, β-sheet 55%, turn 8% and random coil 32.7%. The thermal stability of cocoonase was studied using CD spectroscopy while the thermal property was observed using Differential Scanning Calorimetry (DSC). Furthermore, MALDI-TOF peptide mass fingerprinting (PMF) was performed for similar protein identification using the MASCOT server. Using casein as the substrate, the kinetic constants Km and Vmax were 13 × 103 mg/ml and 15.09 × 10-2 μg/mg.s1 respectively. The specific activity of cocoonase was observed to be maximum at temperature 40 °C, pH-8.0. The effect of heavy metals Hg2+, Cd2+, Co2+, Pb2+ showed inhibitory activity at higher concentrations, while few metals like Mn2+, Fe3+ enhanced the activity while the effect of Ca2+ was not much on the activity. Soybean trypsin inhibitor and PMSF showed an inhibitory effect on the activity of cocoonase. Additionally, antioxidant scavenging and fibrinolytic properties were also observed. Furthermore, the imperative information generated through the present study will serve to explore cocoonase for its prospective pharmaceutical applications.

Identification of cocoonase and cocoonase like protein using polyclonal antibody of Antheraea mylitta cocoonase

BACKGROUND

Cocoonase is a proteolytic enzyme released by silk moths during pupal adult emergence. Without damaging the silk fibroin, this enzyme dissolves the shell of the tasar cocoon by exclusively targeting the protein sericin. Prior to this study, there was no available antibody against Antheraea mylitta cocoonase to identify or screen out similar variants or cocoonase like protein.

RESULTS

In the present study, naturally secreted A. mylitta cocoonase was purified and used to immunize New Zealand white rabbits. The developed polyclonal antibody of cocoonase was purified and its specific interaction with cocoonase was determined using Indirect ELISA. The confirmation of its specificity and immuno-reactivity was evaluated by western blot using native cocoonase of tasar silkworm A. mylitta. The efficacy and specificity of the polyclonal antibody were further verified and confirmed by western blot which was performed to detect ten different ecotypes of A. mylitta cocoonase.

CONCLUSION

The developed antibody successfully detected the cocoonase of different ecotypes. Thus, in future this antibody can serve as one of the molecular detection method for cocoonase and cocoonase-like proteins.

Karanjin, A Promising Bioactive Compound Possessing Anti-cancer Activity against Experimental Model of Non-small Cell Lung Cancer Cells

AIMS

The aim of this study is to isolate the Millettia pinnata (Karanj) leaf extract for pure compound with anticancer properties and to study the molecular target of the isolates in non-small cell lung cancer cell lines.

BACKGROUND

In our earlier research Millettia pinnata leaf extract has demonstrated potential anticancer activities. Thus, in pursuit of the bioactive compounds, the most potential active extract from our previous study was purified. Furthermore, the anticancer properties of the isolated compound karanjin was studied and aimed for apoptosis and restraining growth.

METHODS

A novel method was developed through column chromatography for isolation and purification of the compound karanjin from leaf chloroform extract. The purified component was then characterised using FTIR, mass spectrometry, and NMR. An MTT-based cytotoxicity assay was used to analyse cell cytotoxicity, whereas fluorescence staining was used for apoptosis and reactive oxygen species inhibition quantification. Furthermore, the real-time PCR assay was used to determine the molecular mechanism of action in cells causing cytotoxicity induced by karanjin dosing.

RESULTS

The anticancer activity of karanjin in A549 cell line exhibited prominent activity revealing IC50 value of 4.85 μM. Conferring the predicted molecular pathway study, karanjin restrains the proliferation of cancer cells through apoptosis, which is controlled by extrinsic pathway proteins FAS/FADD/Caspases 8/3/9. Downregulation of KRAS and dependent gene expression also stopped cell proliferation.

CONCLUSION

Karanjin has been identified as a compound with potential effect in non-small cell lung cancer cells. Molecular mechanism for apoptosis and inhibition of reactive oxygen species induced through H2O2 were observed, concluding karanjin have medicinal and antioxidant properties.

Mutational analysis of flavonol synthase of M. pinnata towards enhancement of binding affinity: a computational approach

Millettia pinnata is an important medicinal plant that has been used as a treatment of various diseases due to presence of wide range of pharmacological properties. The plant contains quercetin, kaempferol, karanjin, pongaglabrone, kanjone, kanugin, gammatin, pongaglabol, and other bioflavonoids. Kaempferol is a natural flavonol that shows many pharmacological properties including anti-inflammatory, antioxidant, anticancer, and antidiabetic activities etc. The enzyme flavonol synthase (FLS, EC 1.14.20.6) catalyses the conversion of dihydroflavonols to flavonols, i.e. biosynthesis of kaempferol from dihydrokaempferol. The current work examined the binding affinity-based approach to improve the enzyme catalytic activity using computational methods. Sequential site-directed mutagenesis was used to create four mutants with the goal to increase hydrogen bonds and further improving the ligand (dihydrokaempferol) binding efficiency. Simulations were done to monitor the stability of the mutants followed by molecular docking to confirm interactions with ligand. For structure validation, various dynamic analysis like RMSD, RMSF, ROG, SASA, H-bond, PCA, DCCM, and FEL were performed, which predicts the stability of wild-type (WT) proteins and mutants. The Mutant_2 and Mutant_3 showed maximum H-bonding and better stability than other mutants and WT that proved higher affinity suggesting improved catalysis. Mutant_2 and Mutant_3 exhibited binding affinities of -7.6 and -8.2 kcal/mol, respectively for the ligand. The outcome of present study will provide significant improvement in synthesis of kaempferol and other plant-based flavonoids.Communicated by Ramaswamy H. Sarma.

Protein-protein docking and molecular dynamics studies of sericin and cocoonase of silkworm: an insight for cocoon softening

Cocoonase is known to digest the sericin protein that encapsulates the silkworm cocoon's fibroin protein. Silk fibroin and sericin are two types of proteins that make up silk, and accounts for around 20-30% of the overall cocoon weight. The aim of the study was to see the protein-protein interaction (PPI) and molecular dynamic study of sericin, cocoonase and protein-protein docked complex of silkworm by computational approaches. Here motif analysis, phylogenetic analysis, principal component analysis, root-mean-square deviation (RMSD), root mean square fluctuation, radius of gyration, structural and functional study of cocoonase and sericin as well as molecular docking study were carried out. The 33 amino acid residues of cocoonase shows interaction with 38 aa residues of sericin involving 4 disulphide bonds, 22 hydrogen bonds and 319 non-bonded contacts. The confirmational stability and flexibility of both the proteins as well as protein-protein complex were achieved at 70 ns of MD simulation study. RMSD-based data indicated that cocoonase is more stable than sericin and complex, and complex has a greater fluctuation with more compact (higher Rg) value than cocoonase and sericin, inferring higher conformational stability and flexibility of protein-protein complex than cocoonase and sericin. This study provides a new dimension for PPI study by computational approaches.Communicated by Ramaswamy H. Sarma.

Molecular Beacon Probe (MBP)-Based Real-Time PCR

In the advancement of molecular biology techniques, several probe-based techniques, like molecular beacon probe (MBP) assay, TaqMan probe, and minor groove binder (MGB) probe assay, have been reported to identify specific sequences through real-time polymerase chain reaction (PCR). All probe-based methods are more sensitive than the conventional PCR for the detection and quantification of target genes. MBP is a hydrolysis probe that emits fluorescence when getting the specific sequences on the gene. Here, we describe the application of MBP for the identification of the motif sequences present in the promoters of differentially expressed genes.

Computational models for prediction of protein-protein interaction in rice and Magnaporthe grisea

INTRODUCTION

Plant-microbe interactions play a vital role in the development of strategies to manage pathogen-induced destructive diseases that cause enormous crop losses every year. Rice blast is one of the severe diseases to rice Oryza sativa (O. sativa) due to Magnaporthe grisea (M. grisea) fungus. Protein-protein interaction (PPI) between rice and fungus plays a key role in causing rice blast disease.

METHODS

In this paper, four genomic information-based models such as (i) the interolog, (ii) the domain, (iii) the gene ontology, and (iv) the phylogenetic-based model are developed for predicting the interaction between O. sativa and M. grisea in a whole-genome scale.

RESULTS AND DISCUSSION

A total of 59,430 interacting pairs between 1,801 rice proteins and 135 blast fungus proteins are obtained from the four models. Furthermore, a machine learning model is developed to assess the predicted interactions. Using composition-based amino acid composition (AAC) and conjoint triad (CT) features, an accuracy of 88% and 89% is achieved, respectively. When tested on the experimental dataset, the CT feature provides the highest accuracy of 95%. Furthermore, the specificity of the model is verified with other pathogen-host datasets where less accuracy is obtained, which confirmed that the model is specific to O. sativa and M. grisea. Understanding the molecular processes behind rice resistance to blast fungus begins with the identification of PPIs, and these predicted PPIs will be useful for drug design in the plant science community.

Evaluating the role of trypsin in silk degumming: An in silico approach

The trypsin being universal enzyme forming family of proteases catalyzes the hydrolysis of proteins into amino acids and regenerates the serine hydroxyl an active site. The trypsin enzyme from D. saccharalis, uses sericin as its preferred substrate. Presence of catalytic triad (serine, aspartic acid and histidine) at the substrate binding site of this enzyme is very important for the catalytic activity. In the current study, the interacting mechanism between the substrate sericin protein and enzyme trypsin protein were explored by integrating various computational approaches including physico-chemical properties, biophysical properties, dynamics, gene ontology, molecular docking, protein - protein interactions, binding free energy calculation and structural motifs were studied. The evolutionary study performed by MEGA X showed that trypsin protein sequence (ALE15212.1) is closely related to cocoonase protein sequence (ADG26770.1) from Antheraea pernyi. 3-D models of trypsin and sericin proteins were predicted using I-TASSER and further validated by PROCHECK, and ProSAweb softwares. The predicted trypsin structure model was assigned E.C. no. 3.4.21.4 which refers hydrolytic mechanism. Gene Ontology predicted by QuickGO showed that trypsin has serine hydrolase activity (GO: 00017171), and part of proteolysis (GO: 0006508) as well as protein metabolic process (GO:0019538) actvity. Molecular docking studies between trypsin and sericin proteins were conducted by the HADDOCK 2.4 having best docked protein complex with Z-score - 1.9. 2D and 3D protein-protein interaction was performed with LIGPLOT⁺ and HAWKDOCK, PDBsum, respectively. The amino acid residues interacting across proteins interface are sericin_chain A representing "Ser133, Tyr214, Thr188, Thr243, Ser225, Ser151, Ser156, His294, Arg293, Gly296″ and trypsin_chain B "Lys120, Tyr246, Asn119, Glu239, Ser62, Tyr194, Ile197, Ser171, Tyr169, Gly170″. Based on our results trypsin shows similarity with cocoonase and presumably trypsin can be used as an alternative source in cocoon degumming.

In silico structural and functional characterization of Antheraea mylitta cocoonase

Background

Cocoonase is a serine protease present in sericigenous insects and majorly involved in dissolving of sericin protein allowing moth to escape. Cocoon structure is made up of sericin protein which holds fibroin filaments together. Cocoonase enzyme hydrolyzes sericin protein without harming the fibroin. However, until date, no detailed characterization of cocoonase enzyme and its presence in wild silk moth Antheraea mylitta has been carried out. Therefore, current study aimed for detailed characterization of amplified cocoonase enzyme, secondary and tertiary structure prediction, sequence and structural alignment, phylogenetic analysis, and computational validation. Several computational tools such as ProtParam, Iterative Threading Assembly Refinement (I-TASSER), PROCHECK, SAVES v6.0, TM-align, Molecular Evolutionary Genetics Analysis (MEGA) X, and Figtree were employed for characterization of cocoonase protein.

Results

The present study elucidates about the isolation of RNA, cDNA preparation, PCR amplification, and in silico characterization of cocoonase from Antheraea mylitta. Here, total RNA was isolated from head region of A. mylitta, and gene-specific primers were designed using Primer3 followed by PCR-based amplification and sequencing. The newly constructed 377-bp length sequence of cocoonase was subjected to in silico characterization. In silico study of A. mylitta cocoonase showed 26% similarity to A. pernyi strain Qing-6 cocoonase using Blastp and belongs to member of chymotrypsin-like serine protease superfamily. From phylogenetic study, it was found that A. mylitta cocoonase sequence is closely related to A. pernyi cocoonase sequence.

Conclusions

The present study revealed about the detailed in silico characterization of cocoonase gene and encoded protein obtained from A. mylitta head region. The results obtained infer the presence of cocoonase enzyme in the wild silkworm A. mylitta and can be used for cocoon degumming which will be a valuable and cost-effective strategy in silk industry.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43141-022-00367-8.

ACC deaminase producing rhizobacterium Enterobacter cloacae ZNP-4 enhance abiotic stress tolerance in wheat plant

Plant growth promoting rhizobacterium (PGPR) designated as ZNP-4, isolated from the rhizosphere of Ziziphus nummularia, was identified as Enterobacter cloacae following 16S rRNA sequence analysis. The isolated strain exhibited various plant growth promoting (PGP) traits. The 1-aminocyclopropane-1-carboxylic acid deaminase (ACCD) activity was evaluated under diverse physiological conditions that could be useful for minimizing the abiotic stress-induced inhibitory effects on wheat plants. The strain showed resistance to salt (NaCl) and metal (ZnSO₄) stress. The effect of E. cloacae ZNP-4 on the augmentation of plant growth was studied under salinity stress of 150 mM (T1 treatment) & 200 mM (T2 treatment) NaCl. The inoculation of strain ZNP-4 significantly improved the various growth parameters of wheat plant such as shoot length (41%), root length (31%), fresh weight (28%), dry weight (29%), photosynthetic pigments chlorophyll a (62%) and chlorophyll b (34%). Additionally, the strain was found to be efficient for minimizing the imposed Zn stress in terms of improving plant growth, biomass and photosynthetic pigments in pots containing different levels of metal stress of 150 mg kg^-1 (treatment T1) and 250 mg kg^-1 (treatment T2). Isolate ZNP-4 also improved the proline content and decreased malondialdehyde (MDA) level under both salinity and metal stress, therefore maintaining the membrane integrity. Furthermore, bacterial inoculation increased the activities of antioxidative enzymes such as superoxide dismutase (SOD), catalase (CAT), and peroxidase (POX). The positive effects of PGPR occurred concurrently with the decrease in abiotic stress-induced reactive oxygen species (ROS) molecules such as hydrogen peroxide (H₂O₂) and superoxide (O₂^-) contents. Overall, the observed results indicate that use of bacteria with such beneficial traits could be used as bio-fertilizers for many crops growing under stress conditions.

The molecular docking and molecular dynamics study of flavonol synthase and flavonoid 3'-monooxygenase enzymes involved for the enrichment of kaempferol

Kaempferol is a natural flavonol that shows many pharmacological properties including anti-inflammatory, antioxidant, anticancer, antidiabetic activities etc. It has been reported in many vegetables, fruits, herbs and medicinal plants. The enzyme flavonol synthase (FLS, EC 1.14.20.6) catalyses the conversion of dihydroflavonols to flavonols. Whereas flavonoid 3'-monooxygenase (F3'H, EC 1.14.14.82) catalyses the hydroxylation of dihydroflavonol, and flavonol. FLS is involved in the synthesis of the kaempferol whereas F3'H causes degradation of kaempferol. The present study aimed to analyse the binding affinity, stability and activating activity of enzyme FLS as well as inhibitory activity of enzyme F3'H involved in the enrichment of the kaempferol using the in-silico approaches. Computational study for physico-chemical properties, conserved domain identification, 3-D structure prediction and its validation, conservation analysis, molecular docking followed by molecular dynamics analysis of FLS and F3'H, protein-activator (FLS-LIG Complex) and protein-inhibitor (F3'H-LIG Complex) complexes have been performed. Other structural analyses like root mean square fluctuation (RMSF), root mean square deviation (RMSD), surface area solvent accessibility (SASA), radius of gyration (Rg), hydrogen bond analysis, principal component analysis (PCA), Poisson-Boltzmann analysis (MM_PBSA) and the dynamic cross correlation map (DCCM) analysis to explore the structural, functional and thermodynamic stability of the proteins and the complexes were also studied. The molecular docking result showed that FLS binds strongly with the activator ascorbate (CID _54670067) while F3'H binds with the inhibitor ketoconazole (CID_456201). The most powerful inhibitor (ketoconazole for F3'H) and activator (ascorbate for FLS) is determined by computing the thermodynamic binding free energy through MM_PBSA analysis. The current work provides wide-ranging structural and functional information about FLS and F3'H enzymes showing detailed molecular mechanism of kaempferol biosynthesis and its degradation and hence kaempferol enrichment. Finding of the present work opens up new possibilities for future research towards enrichment of kaempferol by using activator (ascorbate) for FLS and inhibitor (ketoconazole) for F3'H as well as for its large-scale production using in vitro approaches.Communicated by Ramaswamy H. Sarma.

Study on cocoonase, sericin, and degumming of silk cocoon: computational and experimental

Background

Cocoonase is a proteolytic enzyme that helps in dissolving the silk cocoon shell and exit of silk moth. Chemicals like anhydrous Na₂CO₃, Marseille soap, soda, ethylene diamine and tartaric acid-based degumming of silk cocoon shell have been in practice. During this process, solubility of sericin protein increased resulting in the release of sericin from the fibroin protein of the silk. However, this process diminishes natural color and softness of the silk. Cocoonase enzyme digests the sericin protein of silk at the anterior portion of the cocoon without disturbing the silk fibroin. However, no thorough characterization of cocoonase and sericin protein as well as imaging analysis of chemical- and enzyme-treated silk sheets has been carried out so far. Therefore, present study aimed for detailed characterization of cocoonase and sericin proteins, phylogenetic analysis, secondary and tertiary structure prediction, and computational validation as well as their interaction with other proteins. Further, identification of tasar silkworm (Antheraea mylitta) pupa stage for cocoonase collection, its purification and effect on silk sheet degumming, scanning electron microscope (SEM)-based comparison of chemical- and enzyme-treated cocoon sheets, and its optical coherence tomography (OCT)-based imaging analysis have been investigated. Various computational tools like Molecular Evolutionary Genetics Analysis (MEGA) X and Figtree, Iterative Threading Assembly Refinement (I-TASSER), self-optimized predicted method with alignment (SOPMA), PROCHECK, University of California, San Francisco (UCSF) Chimera, and Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) were used for characterization of cocoonase and sericin proteins. Sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), protein purification using Sephadex G 25-column, degumming of cocoon sheet using cocoonase enzyme and chemical Na₂CO₃, and SEM and OCT analysis of degummed cocoon sheet were performed.

Results

Predicted normalized B-factors of cocoonase and sericin with respect to α and β regions showed that these regions are structurally more stable in cocoonase while less stable in sericin. Conserved domain analysis revealed that B. mori cocoonase contains a trypsin-like serine protease with active site range 45 to 180 query sequences while substrate binding site from 175 to 200 query sequences. SDS-PAGE analysis of cocoonase indicated its molecular weight of 25–26 kDa. Na₂CO₃ treatment showed more degumming effect (i.e., cocoon sheet weight loss) as compared to degumming with cocoonase. However, cocoonase-treated silk cocoon sheet holds the natural color of tasar silk, smoothness, and luster compared with the cocoon sheet treated with Na₂CO₃. SEM-based analysis showed the noticeable variation on the surface of silk fiber treated with cocoonase and Na₂CO₃. OCT analysis also exemplified the variations in the cross-sectional view of the cocoonase and Na₂CO₃-treated silk sheets.

Conclusions

Present study enlightens on the detailed characteristics of cocoonase and sericin proteins, comparative degumming activity, and image analysis of cocoonase enzyme and Na₂CO₃ chemical-treated silk sheets. Obtained findings illustrated about use of cocoonase enzyme in the degumming of silk cocoon at larger scale that will be a boon to the silk industry.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43141-021-00125-2.

Method development for optimised green synthesis of gold nanoparticles from Millettia pinnata and their activity in non-small cell lung cancer cell lines

Green synthesis of gold nanoparticles (GNPs) has received substantial attention, because nanoparticles are produced in an eco-friendly way using biomolecules present in plant extracts in a single step reaction. This research article highlights GNPs obtained using shade-dried leaf extracts of Millettia pinnata (L.) with aqueous auric chloride (HAuCl₄) at ambient temperature. In the present study, GNPs with average particle size 37 nm in size were fabricated. Furthermore, the synthesis method to obtain stable and monodispersed GNPs was advanced by optimising enzyme concentration 100 μg/ml, pH 5.4, substrate concentration 0.45 mM and 12 h time of reaction. The confirmation of GNPs formation and characterisation was followed by UV-vis-absorption spectroscopy, dynamic light scattering (DLS), and zeta potential (ZP) for the analysis of shape, size, and stability, respectively. TEM images and powder XRD revealed the GNPs synthesis of spherical-shaped nanoparticles in the face-centred cubic arrangement. Cytotoxicity of GNPs was studied against A549 lung cancer cells with IC₅₀ 14.76 μg/ml and found lower as compared to doxorubicin IC₅₀ 11.23 μg/ml but significant enough to be used as a vehicle GNPs produced using green source can be used as significant therapeutic agents and drug delivery carriers.

Anticancer activity of plant leaves extract collected from a tribal region of India

The goal of this research was to explore the preliminary anticancer properties of five plants namely Calotropis procera, Moringa oleifera, Millettia pinnata, Basela alba and Euphorbia neriifolia available in Jharkhand which is used for the medicinal purpose by local tribes. In the present study, plant leaves from five species were collected, dried and extracted with solvents of increasing polarity, followed by assessment of their cytotoxicity in A549 non-small-cell lung cancer cells. In the antimicrobial assay, the methanol extract of the M. pinnata leaves exhibited comparatively higher zone of inhibition of 0.7 ± 0.20 cm against a Salmonella typhi culture than the other extracts. M. pinnata leaves extract also displayed the maximum percentage inhibition in the DPPH, 83.97 ± 0.01 FRAP, 193.14 ± 3.01 mM assays. Furthermore, the cytotoxicity of the chloroform (37.45 ± 1.04) and ethyl acetate extracts (34.20 ± 0.81) of M. pinnata against A549 cells was found relatively higher with respect to another extract. In contrast, a study with the L132 normal epithelial lung cell line revealed less toxicity from the chloroform extract (0.33 ± 0.19) compared to the ethyl acetate extract (6.65 ± 0.59). Based on these findings, phytochemical investigation on chloroform and ethyl acetate extract of M. pinnata was performed using UPLC-ESI-MS/MS analysis revealing the presence of β-sitosterol, lanceolatin B, karanjin, and stigmasterol. Congruently, a complete phytochemical and cytotoxic investigation of the M. pinnata extract constituents might infer the potency of this extract/s as anticancer, antioxidant and antimicrobial agents.

Gene network modules associated with abiotic stress response in tolerant rice genotypes identified by transcriptome meta-analysis

Abiotic stress tolerance is a complex trait regulated by multiple genes and gene networks in plants. A range of abiotic stresses are known to limit rice productivity. Meta-transcriptomics has emerged as a powerful approach to decipher stress-associated molecular network in model crops. However, retaining specificity of gene expression in tolerant and susceptible genotypes during meta-transcriptome analysis is important for understanding genotype-dependent stress tolerance mechanisms. Addressing this aspect, we describe here "abiotic stress tolerant" (ASTR) genes and networks specifically and differentially expressing in tolerant rice genotypes in response to different abiotic stress conditions. We identified 6,956 ASTR genes, key hub regulatory genes, transcription factors, and functional modules having significant association with abiotic stress-related ontologies and cis-motifs. Out of the 6956 ASTR genes, 73 were co-located within the boundary of previously identified abiotic stress trait-related quantitative trait loci. Functional annotation of 14 uncharacterized ASTR genes is proposed using multiple computational methods. Around 65% of the top ASTR genes were found to be differentially expressed in at least one of the tolerant genotypes under different stress conditions (cold, salt, drought, or heat) from publicly available RNAseq data comparison. The candidate ASTR genes specifically associated with tolerance could be utilized for engineering rice and possibly other crops for broad-spectrum tolerance to abiotic stresses.

Identification of GCC-box and TCC-box motifs in the promoters of differentially expressed genes in rice (Oryza sativa L.): Experimental and computational approaches

The transcription factor selectively binds with the cis-regulatory elements of the promoter and regulates the differential expression of genes. In this study, we aimed to identify and validate the presence of GCC-box and TCC-box motifs in the promoters of upregulated differentially expressed genes (UR-DEGs) and downregulated differentially expressed genes (DR-DEGs) under anoxia using molecular beacon probe (MBP) based real-time PCR. The GCC-box motif was detected in UR-DEGs (DnaJ and 60S ribosomal protein L7 genes), whereas, the TCC-box was detected in DR-DEGs (DnaK and CPuORF11 genes). In addition, the mechanism of interaction of AP2/EREBP family transcription factor (LOC_Os03g22170) with GCC-box promoter motif present in DnaJ gene (LOC_Os06g09560) and 60S ribosomal protein L7 gene (LOC_Os08g42920); and TCC-box promoter motif of DnaK gene (LOC_Os02g48110) and CPuORF11 gene (LOC_Os02g01240) were explored using molecular dynamics (MD) simulations analysis including binding free energy calculations, principal component analyses, and free energy landscapes. The binding free energy analysis revealed that AP2/EREBP model residues such as Arg68, Arg72, Arg83, Lys87, and Arg90 were commonly involved in the formation of hydrogen bonds with GCC and TCC-box promoter motifs, suggesting that these residues are critical for strong interaction. The movement of the entire protein bound to DNA was restricted, confirming the stability of the complex. This study provides comprehensive binding information and a more detailed view of the dynamic interaction between proteins and DNA.

Lac dye as a potential anti-neoplastic agent

OBJECTIVE

Certain anthraquinone compounds are effectively used for treatment of cancer. The present study describes the inhibitory activity of lac dye, a mixture of polyhydroxy anthraquinone compounds (Laccaic acid A, B, C, D and E), produced by the Indian lac insect, Kerria lacca (Kerr).

MATERIALS AND METHODS

In vitro testing for anticancer activity of lac dye was done at four concentrations (10, 20, 40, 80 μg/ml) on ten human malignant cell lines including six human leukemia cell lines, U973, Raji, K562, Jurkat, MOLT 4 and HL60 using SRB assay.

RESULTS AND CONCLUSION

Both crimson and yellow lac dye inhibits proliferation of Human leukemia cell lines, Raji, U937, K562, HL60 and Jurkat, therefore can be considered as a potential anticancer agent for leukemia.

In Silico Characterization and Analysis of RTBP1 and NgTRF1 Protein Through MD Simulation and Molecular Docking: A Comparative Study

Gaining access to sequence and structure information of telomere-binding proteins helps in understanding the essential biological processes involve in conserved sequence-specific interaction between DNA and the proteins. Rice telomere-binding protein (RTBP1) and Nicotiana glutinosa telomere repeat binding factor (NgTRF1) are helix-turn-helix motif type of proteins that plays role in telomeric DNA protection and length regulation. Both the proteins share same type of domain, but till now there is very less communication on the in silico studies of these complete proteins. Here we intend to do a comparative study between two proteins through modeling of the complete proteins, physiochemical characterization, MD simulation and DNA-protein docking. I-TASSER and CLC protein work bench was performed to find out the protein 3D structure as well as the different parameters to characterize the proteins. MD simulation was completed by GROMOS forcefield of GROMACS for 10 ns of time stretch. The simulated 3D structures were docked with template DNA (3D DNA modeled through 3D-DART) of TTTAGGG conserved sequence motif using HADDOCK Web server. By digging up all the facts about the proteins, it was revealed that around 120 amino acids in the tail part were showing a good sequence similarity between the proteins. Molecular modeling, sequence characterization and secondary structure prediction also indicate the similarity between the protein's structure and sequence. The result of MD simulation highlights on the RMSD, RMSF, Rg, PCA and energy plots which also conveys the similar type of motional behavior between them. The best complex formation for both the proteins in docking result also indicates for the first interaction site which is mainly the helix3 region of the DNA-binding domain. The overall computational analysis reveals that RTBP1 and NgTRF1 proteins display good amount of similarity in their physicochemical properties, structure, dynamics and binding mode.

Annotation of Stress-Responsive Candidate Genes in Peanut ESTs

Peanut (Arachis hypogaea L.) is an internationally important crop for human consumption as a good source of protein and vegetable oil. Peanut is widely cultivated around the world in tropical, subtropical and warm temperate climate. Because of its huge genome size (2.8 Gb) and unsequenced genome, studies on genomics and genetic modification of peanut are less as compared to other model crops. As peanut can be cultivated in arid and semiarid regions, its growth is drastically affected by various stresses that reduce the yield. Therefore, study on stress-responsive genes and its regulation is very much important. Here we report about the identification and annotation of some stress-responsive candidate genes using peanut expressed sequence tags (ESTs). The selection of genes was based on the publically available expression data. Due to good expression data and lack of available literature in peanut, some of the stress-responsive genes were screened. Individual EST of the said group was further searched in peanut ESTs (1,78,490 whole EST sequences) using computational approach. Various tools like VecScreen, RepeatMasker, EST trimmer, DNA Baser and Wise2 were being used for stress-responsive gene identification and annotation. Research progress made toward contig assembly, determination of biological function of genes, and prediction of domain as well as 3D structure for related protein are included.

In silico characterization and analysis of RTBP1 and NgTRF1 protein through MD simulation and molecular docking - A comparative study

Gaining access to sequence and structure information of telomere binding proteins helps in understanding the essential biological processes involve in conserved sequence specific interaction between DNA and the proteins. Rice telomere binding protein (RTBP1) and Nicotiana glutinosa telomere repeat binding factor (NgTRF1) are helix turn helix motif type of proteins that plays role in telomeric DNA protection and length regulation. Both the proteins share same type of domain but till now there is very less communication on the in silico studies of these complete proteins.Here we intend to do a comparative study between two proteins through modeling of the complete proteins, physiochemical characterization, MD simulation and DNA-protein docking. I-TASSER and CLC protein work bench was performed to find out the protein 3D structure as well as the different parameters to characterize the proteins. MD simulation was completed by GROMOS forcefield of GROMACS for 10 ns of time stretch. The simulated 3D structures were docked with template DNA (3D DNA modeled through 3D-DART) of TTTAGGG conserved sequence motif using HADDOCK web server.Digging up all the facts about the proteins it was reveled that around 120 amino acids in the tail part was showing a good sequence similarity between the proteins. Molecular modeling, sequence characterization and secondary structure prediction also indicates the similarity between the protein's structure and sequence. The result of MD simulation highlights on the RMSD, RMSF, Rg, PCA and Energy plots which also conveys the similar type of motional behavior between them. The best complex formation for both the proteins in docking result also indicates for the first interaction site which is mainly the helix3 region of the DNA binding domain. The overall computational analysis reveals that RTBP1 and NgTRF1 proteins display good amount of similarity in their physicochemical properties, structure, dynamics and binding mode.

Annotation of stress responsive candidate genes in peanut ESTs

Peanut (Arachis hypogaea L.) is an internationally important crop for human consumption as a good source of protein and vegetable oil. Peanut is widely cultivated around the world in tropical, sub-tropical and warm temperate climate. Because of its huge genome size (2.8 Gb) and unsequenced genome, studies on genomics and genetic modification of peanut are less as compared to other model crops. As peanut can be cultivated in arid and semi-arid regions, and its growth is drastically affected by various stresses that reduces the yield. Therefore, study on stress responsive genes and its regulation are very much important. Here we report about the identification and annotation of some stress responsive candidate genes using peanut Expressed Sequences Tags (ESTs). The selection of genes was based on the publically available expression data. Due to good expression data and lack of available literature in peanut some of the stress responsive genes were screened. Individual EST of the said group were further searched in peanut ESTs (1, 78,490 whole EST sequences) using computational approach. Various tools like Vec-Screen, Repeat Masker, EST Trimmer, DNA Baser and WISE2 were being used for stress responsive gene identification and annotation. Research progress made towards contigs assembly, determination of biological function of genes, and prediction of domain as well as 3D structure for related protein are included.

Identification of new aquaporin genes and single nucleotide polymorphism in bread wheat

Major facilitators of water movement through plant cell membranes include aquaporin proteins. Wheat is among the largest and most important cereal crops worldwide; however, unlike other model plants such as rice, maize and Arabidopsis, little has been reported on wheat major intrinsic proteins (MIPs). This study presents a comprehensive computational identification of 349 new wheat expressed sequence tags (ESTs), encoding 13 wheat aquaporin genes. Identified aquaporins consist of 6 plasma membrane intrinsic proteins (PIP) and 1 TIP showing high sequence similarity with rice aquaporins. We also identified 4 NOD26-like intrinsic proteins (NIP) and 2 SIP members that showed more divergence. Further, expression analysis of the aquaporin genes using the available EST information in UniGene revealed their transcripts were differentially regulated in various stress- and tissue-specific libraries. Allele specific Polymerase chain reaction (PCR) primers based on single nucleotide polymorphism (SNP) were designed using PIP as the target gene and validated on a core set of Indian wheat genotypes. A 3D theoretical model of the wheat aquaporin protein was built by homology modeling and could prove to be useful in the further functional characterization of this protein. Collectively with expression and bioinformatics analysis, our results support the idea that the genes identified in this study signify an important genetic resource providing potential targets to modify the water use properties of wheat.

SVM based model generation for binding site prediction on helix turn helix motif type of transcription factors in eukaryotes

Support vector machine is a class of machine learning algorithms which uses a set of related supervised learning methods for classification and regression. Nowadays this method is vividly applied to many detection problems related with secondary structure, tumor cell and binding residue prediction. In this work, support vector machines (SVMs) have been trained on 90 sequences of transcription factors with HTH motif. Four sequence features were used as attribute for the prediction of interaction site in HTH motif. A web page was also developed so that user can easily enter the protein sequence and receive the output as interaction site predicted or not predicted. The generated model shows a very high amount of accuracy, sensitivity and specificity which proves to be a good model for the selected case.

Identification of new stress-induced microRNA and their targets in wheat using computational approach

MicroRNAs (miRNAs) are a class of short endogenous non-coding small RNA molecules of about 18-22 nucleotides in length. Their main function is to downregulate gene expression in different manners like translational repression, mRNA cleavage and epigenetic modification. Computational predictions have raised the number of miRNAs in wheat significantly using an EST based approach. Hence, a combinatorial approach which is amalgamation of bioinformatics software and perl script was used to identify new miRNA to add to the growing database of wheat miRNA. Identification of miRNAs was initiated by mining the EST (Expressed Sequence Tags) database available at National Center for Biotechnology Information. In this investigation, 4677 mature microRNA sequences belonging to 50 miRNA families from different plant species were used to predict miRNA in wheat. A total of five abiotic stress-responsive new miRNAs were predicted and named Ta-miR5653, Ta-miR855, Ta-miR819k, Ta-miR3708 and Ta-miR5156. In addition, four previously identified miRNA, i.e., Ta-miR1122, miR1117, Ta-miR1134 and Ta-miR1133 were predicted in newly identified EST sequence and 14 potential target genes were subsequently predicted, most of which seems to encode ubiquitin carrier protein, serine/threonine protein kinase, 40S ribosomal protein, F-box/kelch-repeat protein, BTB/POZ domain-containing protein, transcription factors which are involved in growth, development, metabolism and stress response. Our result has increased the number of miRNAs in wheat, which should be useful for further investigation into the biological functions and evolution of miRNAs in wheat and other plant species.

Identification of conserved drought stress responsive gene-network across tissues and developmental stages in rice

Identification of genes that are coexpressed across various tissues and environmental stresses is biologically interesting, since they may play coordinated role in similar biological processes. Genes with correlated expression patterns can be best identified by using coexpression network analysis of transcriptome data. In the present study, we analyzed the temporal-spatial coordination of gene expression in root, leaf and panicle of rice under drought stress and constructed network using WGCNA and Cytoscape. Total of 2199 differentially expressed genes (DEGs) were identified in at least three or more tissues, wherein 88 genes have coordinated expression profile among all the six tissues under drought stress. These 88 highly coordinated genes were further subjected to module identification in the coexpression network. Based on chief topological properties we identified 18 hub genes such as ABC transporter, ATP-binding protein, dehydrin, protein phosphatase 2C, LTPL153 - Protease inhibitor, phosphatidylethanolaminebinding protein, lactose permease-related, NADP-dependent malic enzyme, etc. Motif enrichment analysis showed the presence of ABRE cis-elements in the promoters of > 62% of the coordinately expressed genes. Our results suggest that drought stress mediated upregulated gene expression was coordinated through an ABA-dependent signaling pathway across tissues, at least for the subset of genes identified in this study, while down regulation appears to be regulated by tissue specific pathways in rice.

Identification and annotation of abiotic stress responsive candidate genes in peanut ESTs

Peanut (Arachis hypogaea L.) ranks fifth among the world oil crops and is widely grown in India and neighbouring countries. Due to its large and unknown genome size, studies on genomics and genetic modification of peanut are still scanty as compared to other model crops like Arabidopsis, rice, cotton and soybean. Because of its favourable cultivation in semi-arid regions, study on abiotic stress responsive genes and its regulation in peanut is very much important. Therefore, we aim to identify and annotate the abiotic stress responsive candidate genes in peanut ESTs. Expression data of drought stress responsive corresponding genes and EST sequences were screened from dot blot experiments shown as heat maps and supplementary tables, respectively as reported by Govind et al. (2009). Some of the screened genes having no information about their ESTs in above mentioned supplementary tables were retrieved from NCBI. A phylogenetic analysis was performed to find a group of utmost similar ESTs for each selected gene. Individual EST of the said group were further searched in peanut ESTs (1,78,490 whole EST sequences) using stand alone BLAST. For the prediction as well as annotation of abiotic stress responsive selected genes, various tools (like Vec-Screen, Repeat Masker, EST-Trimmer, DNA Baser, WISE2 and I-TASSER) were used. Here we report the predicted result of Contigs, domain as well as 3D structure for HSP 17.3KDa protein, DnaJ protein and Type 2 Metallothionein protein.

Molecular Dynamics Simulation of Rap1 Myb-type domain in Saccharomyces cerevisiae

Telomere is a nucleoprotein complex that plays important role in stability and their maintenance and consists of random repeats of species specific motifs. In budding Saccharomyces cerevisiae, Repressor Activator Protein 1 (Rap1) is a sequence specific protein that involved in transcriptional regulation. Rap1 consist of three active domains like N-terminal BRCT-domain, DNA-binding domain and C-terminal RCT-domain. In this study the unknown 3D structure of Myb-type domain (having 61 residues) within DNAbinding domain was modeled by Modeller7, and verified using different online bioinformatics tools (ProCheck, WhatIf, Verify3D). Dynamics of Myb-type domain of Rap1was carried out through simulation studies using GROMACS software. Time dependent interactions among the molecules were analyzed by Root Mean Square Deviation (RMSD), Radius of Gyration (Rg) and Root Mean Square Fluctuation (RMSF) plots. Motional properties in reduced dimension were also performed by Principal Component Analysis (PCA). Result indicated that Rap1 interacts with DNA major groove through its Helix Turn Helix motifs. Helix 3 was rigid, less amount of fluctuation was found as it interacts with DNA major groove. Helix2 and N-terminal having considerable fluctuation in the time scale.

Identification of miRNAs in sorghum by using bioinformatics approach

MicroRNAs (miRNAs) regulate gene expression mainly by post-transcriptional gene silencing (PTGS) and in some cases by transcriptional genes silencing (TGS). miRNAs play critical roles in developmental processes, nutrient homeostasis, abiotic stress and pathogen responses of plants. In contrast to the large number of miRNAs predicted in cereal model plant rice, only 148 miRNAs were predicted in sorghum till date (miRBase release 17). This suggested that miRNAs identified in sorghum is far from saturation. Hence, we developed a bioinformatics pipeline using an in-house PERL script and publicly available structure prediction tools to identify miRNAs and their target genes from publically available Expressed Sequence Tags (EST) and Genomic Survey Sequence (GSS). About 1379 known and unique plant miRNAs from 33 different crops were used to predict new miRNAs in sorghum. We identified 31 new miRNAs belonging to 10 different miRNA families. We predicted 72 potential target genes for 31 miRNAs, and most of these target genes are predicted to be involved in plant growth and development.These newly identified miRNAs add to the growing database of miRNA and lay the foundation for further understanding of miRNA function in sorghum plant development.

In silico analysis of enzyme involved in enrichment of citronella oil

Citronella oil is one of the essential oils obtained from Cymbopogon sp. having medicinally important aromatic chemicals (like citronellal, citronellol, hydroxy-citronellol and geraniol) exhibiting insecticidal, anti-oxidant and anti-inflammatory effects. Geraniol Dehydrogenase (GDH) is responsible for the degradation of Citronella oil. Therefore, we aimed to generate 3D structure of GDH and a potent specific GDH inhibitor by homology modelling, virtual screening of ligand database and molecular docking. Inhibitor model indicated strong binding affinity to the binding pocket of GDH and varying affinity for different ligands. Obtained structures will open the possibility of testing new inhibitor families, in addition to new substituent for the already known lead structures.

Identification of single exon genes and their encoded proteins in rice (Oryza sativa L.) genome: an in silico approach

In this study, we predicted Single Exon Genes (SEGs) distributed in whole rice genome and their expressed proteins. Complete genome of rice was retrieved from TIGR. CDS annotation in the FEATURE (GenBank format) was used to predict SEGs sequences. Organelle gene sequences, pseudogenes, tRNA genes, rRNA genes and duplicated genes were eliminated through different bioinformatics tools. A sizeable number (8.1%) of SEGs in whole rice genome were detected. Predicted SEGs were further searched for their differential response under anoxia. Out of total detected SEGs, only 39.33% were anoxia responsive. Among the total detected anoxia-responsive SEG, only 23.48% encode the known proteins.

Molecular marker survey and expression analyses of the rice submergence-tolerance gene SUB1A

The major rice quantitative-trait locus Submergence 1 (Sub1) confers tolerance of submergence for about 2 weeks. To identify novel sources of tolerance, we have conducted a germplasm survey with allele-specific markers targeting SUB1A and SUB1C, two of the three transcription-factor genes within the Sub1 locus. The objective was to identify tolerant genotypes without the SUB1A gene or with the intolerant SUB1A-2 allele. The survey revealed that all tolerant genotypes possessed the tolerant Sub1 haplotype (SUB1A-1/SUB1C-1), whereas all accessions without the SUB1A gene were intolerant. Only the variety James Wee with the SUB1A-2 allele was moderately tolerant. However, some intolerant genotypes with the SUB1A-1 allele were identified and RT-PCR analyses were conducted to compare gene expression in tolerant and intolerant accessions. Initial analyses of leaf samples failed to reveal a clear association of SUB1A transcript abundance and tolerance. Temporal and spatial gene expression analyses subsequently showed that SUB1A expression in nodes and internodes associated best with tolerance across representative genotypes. In James Wee, transcript abundance was high in all tissues, suggesting that some level of tolerance might be conferred by high expression of the SUB1A-2 allele. To further assess tissue-specific expression, we have expressed the GUS reporter gene under the control of the SUB1A-1 promoter. The data revealed highly specific GUS expression at the base of the leaf sheath and in the leaf collar region. Specific expression in the growing part of rice leaves is well in agreement with the role of SUB1A in suppressing leaf elongation under submergence.

In silico analysis of motifs in promoters of differentially expressed genes in rice (Oryza sativa L.) under anoxia

The aim of this study was to characterise the molecular mechanisms of transcriptional regulation of Differentially Expressed Genes (DEGs) in rice coleoptiles under anoxia by identifying motifs that are common in the promoter region of co-regulated genes. Un-changed DEGs (<2 fold and >-2), up-regulated DEGs (>or=2 fold) and down-regulated DEGs (<or=-2 fold) were separated in three different data sets. Their gene promoters were extracted from eukaryotic promoter database. Statistically significant consensus promoter motifs were detected by in silico method. A significant variation in the number of promoter motifs, consensus promoter motif and their sequences between UR-DEGs and DR-DEGs were detected that might be responsible for their related expression.