We report a detailed structural analysis of the psychrophilic exo-β-1,3-glucanase (GaExg55) from Glaciozyma antarctica PI12. This study elucidates the structural basis of exo-1,3-β-1,3-glucanase from this psychrophilic yeast. The structural prediction of GaExg55 remains a challenge because of its low sequence identity (37 %). A 3D model was constructed for GaExg55. Threading approach was employed to determine a suitable template and generate optimal target-template alignment for establishing the model using MODELLER9v15. The primary sequence analysis of GaExg55 with other mesophilic exo-1,3-β-glucanases indicated that an increased flexibility conferred to the enzyme by a set of amino acids substitutions in the surface and loop regions of GaExg55, thereby facilitating its structure to cold adaptation. A comparison of GaExg55 with other mesophilic exo-β-1,3-glucanases proposed that the catalytic activity and structural flexibility at cold environment were attained through a reduced amount of hydrogen bonds and salt bridges, as well as an increased exposure of the hydrophobic side chains to the solvent. A molecular dynamics simulation was also performed using GROMACS software to evaluate the stability of the GaExg55 structure at varying low temperatures. The simulation result confirmed the above findings for cold adaptation of the psychrophilic GaExg55. Furthermore, the structural analysis of GaExg55 with large catalytic cleft and wide active site pocket confirmed the high activity of GaExg55 to hydrolyze polysaccharide substrates.
Malaria parasites, Plasmodium can infect a wide range of hosts including humans and rodents. There are two copies of mitogen activated protein kinases (MAPKs) in Plasmodium, namely MAPK1 and MAPK2. The MAPKs have been studied extensively in the human Plasmodium, P. falciparum. However, the MAPKs from other Plasmodium species have not been characterized and it is therefore the premise of presented study to characterize the MAPKs from other Plasmodium species-P. vivax, P knowlesi, P berghei, P chabaudi and P.yoelli using a series of publicly available bioinformatic tools. In silico data indicates that all Plasmodium MAPKs are nuclear-localized and contain both a nuclear localization signal (NLS) and a Leucine-rich nuclear export signal (NES). The activation motifs of TDY and TSH were found to be fully conserved in Plasmodium MAPK1 and MAPK2, respectively. The detailed manual inspection of a multiple sequence alignment (MSA) construct revealed a total of 17 amino acid stack patterns comprising of different amino acids present in MAPK1 and MAPK2 respectively, with respect to rodent and human Plasmodia. It is proposed that these amino acid stack patterns may be useful in explaining the disparity between rodent and human Plasmodium MAPKs.
Dimocarpus longan, commonly known as the longan, belongs to the family Sapindaceae, and is one of the most economically important fruits commonly cultivated in several regions in Asia. There are various cultivars of longan throughout the Thai-Malay peninsula region, but until now no phylogenetic analysis has been undertaken to determine the genetic relatedness of these cultivars. To address this issue, 6 loci, namely ITS2, matK, rbcL, trnH-psbA, trnL-I and trnL-trnF were amplified and sequenced from 40 individuals consisting of 26 longan cultivars 2 types of lychee and 8 herbarium samples. The sequencing results were used to construct a phylogenetic tree using the neighbor-joining (NJ), maximum likelihood (ML) and Bayesian inference (BI) criteria. The tree showed cryptic groups of D. longan from the Thailand-Malaysia region (Dimocarpus longan spp.). This is the first report of the genetic relationship of Dimocarpus based on multi-locus molecular markers and morphological characteristics. Multiple sequence alignments, phylogenetic trees and species delimitation support that Dimocarpus longan spp. longan var. obtusus and Dimocarpus longan spp. malesianus var. malesianus should be placed into a higher order and are two additional species in the genus Dimocarpus. Therefore these two species require nomenclatural changes as Dimocarpus malesianus and Dimocarpus obtusus, respectively.
The cell wall determines morphology and the environmental responses of plant cells. The primary cell wall (PCW) is produced during cell division and expansion, determining the cell shape and volume. After cell expansion, specific types of plant cells produce a lignified wall, known as a secondary cell wall (SCW). We functionally analyzed Group IIId Arabidopsis AP2/EREBP genes, namely ERF34, ERF35, ERF38, and ERF39, which are homologs of a rice ERF gene previously proposed to be related to SCW biosynthesis. Expression analysis revealed that these four genes are expressed in regions related to cell division and/or cell differentiation in seedlings (i.e., shoot apical meristems, the primordia of leaves and lateral roots, trichomes, and central cylinder of primary roots) and flowers (i.e., vascular tissues of floral organs and replums and/or valve margins of pistils). Overexpression of ERF genes significantly upregulated PCW-type, but not SCW-type, CESA genes encoding cellulose synthase catalytic subunits in Arabidopsis seedlings. Transient co-expression reporter analysis indicated that ERF35, ERF38, and ERF39 possess transcriptional activator activity, and that ERF34, ERF35, ERF38, and ERF39 upregulated the promoter activity of CESA1, a PCW-type CESA gene, through the DRECRTCOREAT elements, the core cis-acting elements known to be recognized by AP2/ERF proteins. Together, our findings show that Group IIId ERF genes are positive transcriptional regulators of PCW-type CESA genes in Arabidopsis and are possibly involved in modulating cellulose biosynthesis in response to developmental requirements and environmental stimuli.
In August 2008, 30% of tomato (Solanum lycopersicum) plants in plots in Lubbock County, Texas showed yellowing, lateral stem dieback, upward leaf curling, enlargement of stems, adventitious roots, and swollen nodes. Yellowing in leaves was similar to that seen with zebra chip disease (ZC) of potato that was confirmed in a potato field 112 km away in July 2008 and was associated with a 'Candidatus Liberibacter' species (1), similar to findings earlier in 2008 in New Zealand and California (2,3). Tissue from four symptomatic plants of cv. Spitfire and two of cv. Celebrity were collected and DNA was extracted from midribs and petioles with a FastDNA Spin Kit (Qbiogene, Inc., Carlsbad, CA,). PCR amplification was done with 16S rRNA gene primers OA2 and OI2c, which are specific for "Ca. Liberibacter solanacearum" from potato and tomato and amplify a 1.1-kb fragment of the 16S rRNA gene of this new species (1,3). Amplicons of 1.1 kb were obtained from all samples and these were sequenced in both orientations (McLab, San Francisco, CA). Sequences of the 16S rRNA gene were identical for both Spitfire and Celebrity and were submitted to the NCBI as GenBank Accession Nos. FJ939136 and FJ939137, respectively. On the basis of a BLAST search, sequence alignments revealed 99.9% identity with a new species of 'Ca. Liberibacter' from potato (EU884128 and EU884129) in Texas (1); 99.7% identity with the new species "Ca. Liberibacter solanacearum" described from potato and tomato (3) in New Zealand (EU849020 and EU834130, respectively) and from the potato psyllid Bactericera cockerelli in California (2) (EU812559, EU812556); 97% identity with 'Ca L. asiaticus' from citrus in Malaysia (EU224393) and 94% identity with both 'Ca. L. africanus' and 'Ca. L. americanus' from citrus (EU921620 and AY742824, respectively). A neighbor-joining cladogram constructed using the 16S rRNA gene fragments delineated four clusters corresponding to each species, and these sequences clustered with "Ca. L. solanacearum". A second PCR analysis was conducted with the CL514F/CL514R primer pair, which amplifies a sequence from the rplJ and rplL ribosomal protein genes of "Ca. L. solanacearum". The resulting 669-bp products were 100% identical to a sequence reported from tomato in Mexico (FJ498807). This sequence was submitted to NCBI (GU169328). ZC, a disease causing losses to the potato industry, is associated with a 'Candidatus Liberibacter' species (1-3) and was reported in Central America and Mexico in the 1990s, in Texas in 2000, and more recently in other states in the United States (4). In 2008, a "Ca. Liberibacter solanacearum" was detected on Capsicum annuum, S. betaceum, and Physalis peruviana in New Zealand (3). Several studies have shown that the potato psyllid, B. cockerelli, is a potential vector for this pathogen (2,4). To our knowledge, this is the first report of "Ca. Liberibacter solanacearum" in field tomatoes showing ZC-like foliar disease symptoms in the United States. References: (1). J. A. Abad et al. Plant Dis. 93:108, 2009 (2) A. K. Hansen et al. Appl. Environ. Microbiol. 74:5862, 2008. (3) L. W. Liefting et al. Plant Dis. 93:208, 2009. (4) G. A. Secor et al. Plant Dis. 93:574, 2009.
The Arracacia clade (Apiaceae, Apioideae) is a heterogeneous assemblage of 12 genera, comprising 111 known species distributed in high montane temperate and sub-alpine habitats of meso- and South America. Previous studies have indicated that the genera Arracacia, Coulterophytum, and Prionosciadium are polyphyletic, but for the most part relationships among the members of the clade are largely unknown. Initially, cladistic analyses of nrDNA ITS sequences were carried out on 212 accessions (122 taxa), representing 92 species of the Arracacia clade and outgroups from the closely-related páramo genera Cotopaxia, Niphogeton, and Perissocoeleum and members of the Perennial Endemic North American clade and its allies. Using the ITS results to inform sampling of a small subset of taxa, a pilot study examining the phylogenetic utility of 20 noncoding chloroplast loci was subsequently performed to identify those regions most useful at resolving relationships. A cost-benefit analysis determined that five loci (trnQ-5'rps16, trnD-trnT, rpl32-trnL, psbD-trnT, ndhA intron) would maximize resolution and branch support in the clade. Cladistic analyses of four of these loci (trnQ-5'rps16, trnD-trnT, rpl32-trnL, ndhA intron) and the ITS region, separately and combined, revealed that Arracacia, Coaxana, Coulterophytum, Prionosciadium, and Rhodosciadium are each polyphyletic and that Donnellsmithia and Myrrhidendron are each monophyletic. Although most relationships in the Arracacia clade and among the closely-related genera Cotopaxia, Niphogeton, and Perissocoeleum are poorly resolved and supported, ten groups are recognized for future revisionary studies. Polyploidy and rapid species radiation have likely confounded generic circumscriptions and interpretation of relationships.
BACKGROUND: Patient-Derived Tumour Xenografts (PDTXs) have emerged as the pre-clinical models that best represent clinical tumour diversity and intra-tumour heterogeneity. The molecular characterization of PDTXs using High-Throughput Sequencing (HTS) is essential; however, the presence of mouse stroma is challenging for HTS data analysis. Indeed, the high homology between the two genomes results in a proportion of mouse reads being mapped as human.
RESULTS: In this study we generated Whole Exome Sequencing (WES), Reduced Representation Bisulfite Sequencing (RRBS) and RNA sequencing (RNA-seq) data from samples with known mixtures of mouse and human DNA or RNA and from a cohort of human breast cancers and their derived PDTXs. We show that using an In silico Combined human-mouse Reference Genome (ICRG) for alignment discriminates between human and mouse reads with up to 99.9% accuracy and decreases the number of false positive somatic mutations caused by misalignment by >99.9%. We also derived a model to estimate the human DNA content in independent PDTX samples. For RNA-seq and RRBS data analysis, the use of the ICRG allows dissecting computationally the transcriptome and methylome of human tumour cells and mouse stroma. In a direct comparison with previously reported approaches, our method showed similar or higher accuracy while requiring significantly less computing time.
CONCLUSIONS: The computational pipeline we describe here is a valuable tool for the molecular analysis of PDTXs as well as any other mixture of DNA or RNA species.
Brugia malayi is a filarial nematode, which causes lymphatic filariasis in humans. In 1995, the disease has been identified by the World Health Organization (WHO) as one of the second leading causes of permanent and long-term disability and thus it is targeted for elimination by year 2020. Therefore, accurate filariasis diagnosis is important for management and elimination programs. A recombinant antigen (BmR1) from the Bm17DIII gene product was used for antibody-based filariasis diagnosis in "Brugia Rapid". However, the structure and dynamics of BmR1 protein is yet to be elucidated. Here we study the three dimensional structure and dynamics of BmR1 protein using comparative modeling, threading and ab initio protein structure prediction. The best predicted structure obtained via an ab initio method (Rosetta) was further refined and minimized. A total of 5 ns molecular dynamics simulation were performed to investigate the packing of the protein. Here we also identified three epitopes as potential antibody binding sites from the molecular dynamics average structure. The structure and epitopes obtained from this study can be used to design a binder specific against BmR1, thus aiding future development of antigen-based filariasis diagnostics to complement the current diagnostics.
High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are available for only a few non-microbial species1-4. To address this issue, the international Genome 10K (G10K) consortium5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling highly accurate and nearly complete reference genomes. Here we present lessons learned from generating assemblies for 16 species that represent six major vertebrate lineages. We confirm that long-read sequencing technologies are essential for maximizing genome quality, and that unresolved complex repeats and haplotype heterozygosity are major sources of assembly error when not handled correctly. Our assemblies correct substantial errors, add missing sequence in some of the best historical reference genomes, and reveal biological discoveries. These include the identification of many false gene duplications, increases in gene sizes, chromosome rearrangements that are specific to lineages, a repeated independent chromosome breakpoint in bat genomes, and a canonical GC-rich pattern in protein-coding genes and their regulatory regions. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an international effort to generate high-quality, complete reference genomes for all of the roughly 70,000 extant vertebrate species and to help to enable a new era of discovery across the life sciences.
The American mink (Neovison vison) is a semiaquatic species of mustelid native to North America. It's an important animal for the fur industry. Many efforts have been made to locate genes influencing fur quality and color, but this search has been impeded by the lack of a reference genome. Here we present the first draft genome of mink. In our study, two mink individuals were sequenced by Illumina sequencing with 797 Gb sequence generated. Assembly yielded 7,175 scaffolds with an N50 of 6.3 Mb and length of 2.4 Gb including gaps. Repeat sequences constitute around 31% of the genome, which is lower than for dog and cat genomes. The alignments of mink, ferret and dog genomes help to illustrate the chromosomes rearrangement. Gene annotation identified 21,053 protein-coding sequences present in mink genome. The reference genome's structure is consistent with the microsatellite-based genetic map. Mapping of well-studied genes known to be involved in coat quality and coat color, and previously located fur quality QTL provide new knowledge about putative candidate genes for fur traits. The draft genome shows great potential to facilitate genomic research towards improved breeding for high fur quality animals and strengthen our understanding on evolution of Carnivora.
Polyhydroxyalkanoates (PHA), are biodegradable polyesters derived from many microorganisms such as the pseudomonads. These polyesters are in great demand especially in the packaging industries, the medical line as well as the paint industries. The enzyme responsible in catalyzing the formation of PHA is PHA synthase. Due to the limited structural information, its functional properties including catalysis are lacking. Therefore, this study seeks to investigate the structural properties as well as its catalytic mechanism by predicting the three-dimensional (3D) model of the Type II Pseudomonas sp. USM 4-55 PHA synthase 1 (PhaC1P.sp USM 4-55).
Mycobacterium abscessus subspecies classification has important clinical implications. We used phylogenomic network and amino acid analyses to provide evidence for the separation of Mycobacterium bolletii and Mycobacterium massiliense into two distinct subspecies which can potentially be differentiated rapidly by their protein signatures.
An outbreak of coronavirus disease 2019 (COVID-19) caused by the 2019 novel coronavirus (SARS-CoV-2) began in the city of Wuhan in China and has widely spread worldwide. Currently, it is vital to explore potential intermediate hosts of SARS-CoV-2 to control COVID-19 spread. Therefore, we reinvestigated published data from pangolin lung samples from which SARS-CoV-like CoVs were detected by Liu et al. [1]. We found genomic and evolutionary evidence of the occurrence of a SARS-CoV-2-like CoV (named Pangolin-CoV) in dead Malayan pangolins. Pangolin-CoV is 91.02% and 90.55% identical to SARS-CoV-2 and BatCoV RaTG13, respectively, at the whole-genome level. Aside from RaTG13, Pangolin-CoV is the most closely related CoV to SARS-CoV-2. The S1 protein of Pangolin-CoV is much more closely related to SARS-CoV-2 than to RaTG13. Five key amino acid residues involved in the interaction with human ACE2 are completely consistent between Pangolin-CoV and SARS-CoV-2, but four amino acid mutations are present in RaTG13. Both Pangolin-CoV and RaTG13 lost the putative furin recognition sequence motif at S1/S2 cleavage site that can be observed in the SARS-CoV-2. Conclusively, this study suggests that pangolin species are a natural reservoir of SARS-CoV-2-like CoVs.
In this study, we reported a full length of catalase gene (designated as MrCat), identified from the transcriptome database of freshwater prawn Macrobrachium rosenbergii. The complete gene sequence of the MrCat is 2504 base pairs in length, and encodes 516 amino acids. The MrCat protein contains three domains such as catalase 1 (catalase proximal heme-ligand signature) at 350-358, catalase 2 (catalase proximal active site signature) at 60-76 and catalase 3 (catalase family profile) at 20-499. The mRNA expressions of MrCat in healthy and the infectious hypodermal and hematopoietic necrosis virus (IHHNV) challenged M. rosenbergii were examined using quantitative real time polymerase chain reaction (qRT-PCR). The MrCat is highly expressed in digestive tract and all the other tissues (walking leg, gills, muscle, hemocyte, hepatopancreas, pleopods, brain and eye stalk) of M. rosenbergii taken for analysis. The expression is strongly up-regulated in digestive tract after IHHNV challenge. To understand its biological activity, the recombinant MrCat gene was constructed and expressed in Escherichia coli BL21 (DE3). The recombinant MrCat existed in high thermal stability and broad spectrum of pH, which showed over 95% enzyme activity between pH 5 and 10.5, and was stable from 40 °C to 70 °C, and exhibited 85-100% enzyme activity from 30 °C to 40 °C.
This study reports the comprehensive comparative information of two different detoxification enzymes such as glutathione S-transferases (GSTs) delta and kappa from freshwater giant prawn Macrobrachium rosenbergii (designated as MrGSTD and MrGSTK) by investigating their in-silico characters and mRNA modulation against various biotic and abiotic oxidative stressors. The physico-chemical properties of these cDNA and their polypeptide structure were analyzed using various bioinformatics program. The analysis indicated the variation in size of the polypeptides, presence or absence of domains and motifs and structure. Homology and phylogenetic analysis revealed that MrGSTD shared maximum identity (83%) with crustaceans GST delta, whereas MrGSTK fell in arthropods GST kappa. It is interesting to note that MrGSTD and MrGSTK shared only 21% identity; it indicated their structural difference. Structural analysis indicated that MrGSTD to be canonical dimer like shape and MrGSTK appeared to be butterfly dimer like shape, in spite of four β-sheets being conserved in both GSTs. Tissue specific gene expression analysis showed that both MrGSTD and MrGSTK are highly expressed in immune organs such as haemocyte and hepatopancreas, respectively. To understand the role of mRNA modulation of MrGSTD and MrGSTK, the prawns were inducted with oxidative stressors such as bacteria (Vibrio harveyi), virus [white spot syndrome virus (WSSV)] and heavy metal, cadmium (Cd). The analysis revealed an interesting fact that both MrGSTD and MrGSTK showed higher (P
Babesia gibsoni is a tick-borne hemoprotozoan parasite of dogs that often causes fever and hemolytic illness. Detection of B. gibsoni has been predominantly reported in Asian countries, including Japan, Korea, Taiwan, Malaysia, Bangladesh and India. The present study shows the first molecular characterization of B. gibsoni detected from dogs in Bangladesh. Blood samples were collected on FTA® Elute cards from 50 stray dogs in Mymensingh District in Bangladesh. DNA eluted from the cards was subjected to nested PCR for the 18S rRNA gene of Babesia species. Approximately 800bp PCR products were detected in 15 of 50 dogs (30%). Based on restriction fragment length polymorphism (RFLP) and direct sequencing of the PCR products, all parasite isolates were identified as B. gibsoni. Furthermore, the BgTRAP (B. gibsoni thrombospondin-related adhesive protein) gene fragments were detected in 13 of 15 18S rRNA gene PCR positive blood samples. Phylogenetic analysis of the BgTRAP gene revealed that B. gibsoni parasites in Bangladesh formed a cluster, which was genetically different from other Asian B. gibsoni isolates. In addition, tandem repeat analysis of the BgTRAP gene clearly showed considerable genetic variation among Bangladeshi isolates. These results suggested that B. gibsoni parasites in a different genetic clade are endemic in dogs in Bangladesh. Further studies are required to elucidate the origin, distribution, vector and pathogenesis of B. gibsoni parasites circulating in dogs in Bangladesh.
Cytochrome P450s are a superfamily of heme monooxygenases which catalyze a wide range of biochemical reactions. The reactions involve the introduction of an oxygen atom into an inactivated carbon of a compound which is essential to produce an intermediate of a hydroxylated product. The diversity of chemical reactions catalyzed by cytochrome P450s has led to their increased demand in numerous industrial and biotechnology applications. A recent study showed that a gene sequence encoding a CYP was found in the genome of Bacillus lehensis G1, and this gene shared structural similarity with the bacterial vitamin D hydroxylase (Vdh) from Pseudonocardia autotrophica. The objectives of present study was to mine, for a novel CYP from a new isolate B. lehensis G1 alkaliphile and determine the biological properties and functionalities of CYP in this bacterium. Our study employed the usage of computational methods to search for the novel CYP from CYP structural databases to identify the conserved pattern, functional domain and sequence properties of the uncharacterized CYP from B. lehensis G1. A computational homology model of the protein's structure was generated and a docking analysis was performed to provide useful structural knowledge on the enzyme's possible substrate and their interaction. Sequence analysis indicated that the newly identified CYP, termed CYP107CB2, contained the fingerprint heme binding sequence motif FxxGxxxCxG at position 336-345 as well as other highly conserved motifs characteristic of cytochrome P450 proteins. Using docking studies, we identified Ser-79, Leu-81, Val-231, Val-279, Val-383, Ala-232, Thr-236 and Thr-283 as important active site residues capable of stabilizing interactions with several potential substrates, including vitamin D3, 25-hydroxyvitamin D3 and 1α-hydroxyvitamin D3, in which all substrates docked proximally to the enzyme's heme center. Biochemical analysis indicated that CYP107CB2 is a biologically active protein to produce 1α,25-dihydroxyvitamin D3 from 1α-hydroxyvitamin D3. Based on these results, we conclude that the novel CYP107CB2 identified from B. lehensis G1 is a putative vitamin D hydroxylase which is possibly capable of catalyzing the bioconversion of parental vitamin D3 to calcitriol, or related metabolic products.
beta-D-N-Acetylhexosaminidase, a family 20 glycosyl hydrolase, catalyzes the removal of β-1,4-linked N-acetylhexosamine residues from oligosaccharides and their conjugates. We constructed phylogenetic tree of β-hexosaminidases to analyze the evolutionary history and predicted functions of plant hexosaminidases. Phylogenetic analysis reveals the complex history of evolution of plant β-hexosaminidase that can be described by gene duplication events. The 3D structure of tomato β-hexosaminidase (β-Hex-Sl) was predicted by homology modeling using 1now as a template. Structural conformity studies of the best fit model showed that more than 98% of the residues lie inside the favoured and allowed regions where only 0.9% lie in the unfavourable region. Predicted 3D structure contains 531 amino acids residues with glycosyl hydrolase20b domain-I and glycosyl hydrolase20 superfamily domain-II including the (β/α)8 barrel in the central part. The α and β contents of the modeled structure were found to be 33.3% and 12.2%, respectively. Eleven amino acids were found to be involved in ligand-binding site; Asp(330) and Glu(331) could play important roles in enzyme-catalyzed reactions. The predicted model provides a structural framework that can act as a guide to develop a hypothesis for β-Hex-Sl mutagenesis experiments for exploring the functions of this class of enzymes in plant kingdom.
A key event in the domestication and breeding of the oil palm Elaeis guineensis was loss of the thick coconut-like shell surrounding the kernel. Modern E. guineensis has three fruit forms, dura (thick-shelled), pisifera (shell-less) and tenera (thin-shelled), a hybrid between dura and pisifera. The pisifera palm is usually female-sterile. The tenera palm yields far more oil than dura, and is the basis for commercial palm oil production in all of southeast Asia. Here we describe the mapping and identification of the SHELL gene responsible for the different fruit forms. Using homozygosity mapping by sequencing, we found two independent mutations in the DNA-binding domain of a homologue of the MADS-box gene SEEDSTICK (STK, also known as AGAMOUS-LIKE 11), which controls ovule identity and seed development in Arabidopsis. The SHELL gene is responsible for the tenera phenotype in both cultivated and wild palms from sub-Saharan Africa, and our findings provide a genetic explanation for the single gene hybrid vigour (or heterosis) attributed to SHELL, via heterodimerization. This gene mutation explains the single most important economic trait in oil palm, and has implications for the competing interests of global edible oil production, biofuels and rainforest conservation.
Klebsiella pneumoniae is a Gram-negative, cylindrical rod shaped opportunistic pathogen that is found in the environment as well as existing as a normal flora in mammalian mucosal surfaces such as the mouth, skin, and intestines. Clinically it is the most important member of the family of Enterobacteriaceae that causes neonatal sepsis and nosocomial infections. In this work, a combination of protein sequence analysis, structural modeling and molecular docking simulation approaches were employed to provide an understanding of the possible functions and characteristics of a hypothetical protein (KPN_02809) from K. pneumoniae MGH 78578. The computational analyses showed that this protein was a metalloprotease with zinc binding motif, HEXXH. To verify this result, a ypfJ gene which encodes for this hypothetical protein was cloned from K. pneumoniae MGH 78578 and the protein was overexpressed in Escherichia coli BL21 (DE3). The purified protein was about 32 kDa and showed maximum protease activity at 30 °C and pH 8.0. The enzyme activity was inhibited by metalloprotease inhibitors such as EDTA, 1,10-phenanthroline and reducing agent, 1,4-dithiothreitol (DTT). Each molecule of KPN_02809 protein was also shown to bind one zinc ion. Hence, for the first time, we experimentally confirmed that KPN_02809 is an active enzyme with zinc metalloprotease activity.