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.
The nucleotide sequences of the second internal transcribed spacer of rDNA were determined for adult worms of Necator americanus originating from Togo (Africa) and Sarawak (Malaysia). The length of the sequences of specimens from Togo (325 bp) were shorter than those from Sarawak (327 bp). There were six fixed genetic differences in the aligned sequences of N. americanus from Sarawak and Togo, excluding one or two polymorphic sites within the sequence of N. americanus from each geographical region. These findings suggest that there is either population variation in the sequence of N. americanus, or that N. americanus from the two countries may represent genetically distinct but morphologically similar (i.e. cryptic) species, however, comparison of the sequence differences among other hookworm species supports the latter conclusion.
Genome data mining of the Antarctic yeast, Glaciozyma antarctica PI12 revealed an expansin-like protein encoding sequence (GaEXLX1). The GaEXLX1 protein is 24.8 kDa with a high alkaline pI of 9.81. Homology modeling of GaEXLX1 showed complete D1 and D2 domains of a conventional expansin. The protein exhibited 36% sequence similarity to Clavibacter michiganensis EXLX1 (PDB: 4JCW). Subsequently, a recombinant GaEXLX1 protein was produced using Escherichia coli expression system. Incubation with Avicel, filter paper and cotton fiber showed that the protein can disrupt the surface of crystalline and pure cellulose, suggesting a cell wall modification activity usually exhibited by expansin-like proteins. Binding assays displayed that GaEXLX1 can bind to polymeric substrates, including those postulated to be present in the sea ice ecosystem such as crab chitin and moss lichenan. GaEXLX1 may assist in the recognition and loosening of these substrates in the sea ice prior to hydrolysis by other extracellular enzymes. Similar loosening mechanism to classical expansin-like protein has been postulated for this psychrophilic protein based on several conserved residues of GaEXLX1 involved in binding interaction identified by docking analyses.
Chitinases play a vital role during the pathogenic invasion and immunosuppression in various organisms including invertebrates and vertebrates. In this study, we have investigated the participation of MrChit-3 (Macrobrachium rosenbergii Chitinase-3) during host-pathogenic interaction in freshwater prawn, M. rosenbergii. Quantitative real-time PCR analysis showed that the expression of MrChit-3 was up-regulated during bacterial, viral and laminarin challenge. Moreover, to understand the antimicrobial role of the GH18 domain, a putative membrane-targeting antimicrobial peptide (MrVG) was identified from the GH18 domain region of the protein and it was chemically synthesized. Physico-chemical features of the GH18 derived antimicrobial peptide (AMP) was assessed by various in silico tools and the antimicrobial property of the peptide was confirmed from in vitro studies. The membrane targeting mechanism of the peptide was determined by flow cytometry (FACS) and scanning electron microscope (SEM) analysis. Interestingly, the peptide was able to inhibit the growth of a chitinolytic fungal pathogen, Aspergillus niger, which was isolated from the shells of M. rosenbergii. The toxicity studies such as hemolysis activity on human blood erythrocytes and cell viability assay with primary kidney cells, HEK293 of MrVG revealed that the peptide was not involved in inducing any toxicity.
Studies related to the engineering of calcium binding sites of CGTase are limited. The calcium binding regions that are known for thermostability function were subjected to site-directed mutagenesis in this study. The starting gene-protein is a variant of CGTase Bacillus sp. G1, reported earlier and denoted as "parent CGTase" herein. Four CGTase variants (S182G, S182E, N132R and N28R) were constructed. The two variants with a mutation at residue 182, located adjacent to the Ca-I site and the active site cleft, possessed an enhanced thermostability characteristic. The activity half-life of variant S182G at 60 °C was increased to 94 min, while the parent CGTase was only 22 min. This improvement may be attributed to the formation of a shorter α-helix and the alleviation of unfavorable steric strains by glycine at the corresponding region. For the variant S182E, an extra ionic interaction at the A/B domain interface increased the half-life to 31 min, yet it reduced CGTase activity. The introduction of an ionic interaction at the Ca-I site via the mutation N132R disrupted CGTase catalytic activity. Conversely, the variant N28R, which has an additional ionic interaction at the Ca-II site, displayed increased cyclization activity. However, thermostability was not affected.
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.
Species-specific simple sequence repeat (SSR) markers are favored for genetic studies and marker-assisted selection (MAS) breeding for oil palm genetic improvement. This report characterizes 20 SSR markers from an Elaeis oleifera genomic library (gSSR). Characterization of the repeat type in 2000 sequences revealed a high percentage of di-nucleotides (63.6%), followed by tri-nucleotides (24.2%). Primer pairs were successfully designed for 394 of the E. oleifera gSSRs. Subsequent analysis showed the ability of the 20 selected E. oleifera gSSR markers to reveal genetic diversity in the genus Elaeis. The average Polymorphism Information Content (PIC) value for the SSRs was 0.402, with the tri-repeats showing the highest average PIC (0.626). Low values of observed heterozygosity (H(o)) (0.164) and highly positive fixation indices (F(is)) in the E. oleifera germplasm collection, compared to the E. guineensis, indicated an excess of homozygosity in E. oleifera. The transferability of the markers to closely related palms, Elaeis guineensis, Cocos nucifera and ornamental palms is also reported. Sequencing the amplicons of three selected E. oleifera gSSRs across both species and palm taxa revealed variations in the repeat-units. The study showed the potential of E. oleifera gSSR markers to reveal genetic diversity in the genus Elaeis. The markers are also a valuable genetic resource for studying E. oleifera and other genus in the Arecaceae family.
A thermophilic lipolytic bacterium identified as Bacillus sp. L2 via 16S rDNA was previously isolated from a hot spring in Perak, Malaysia. Bacillus sp. L2 was confirmed to be in Group 5 of bacterial classification, a phylogenically and phenotypically coherent group of thermophilic bacilli displaying very high similarity among their 16S rRNA sequences (98.5-99.2%). Polymerase chain reaction (PCR) cloning of L2 lipase gene was conducted by using five different primers. Sequence analysis of the L2 lipase gene revealed an open reading frame (ORF) of 1251 bp that codes for 417 amino acids. The signal peptides consist of 28 amino acids. The mature protein is made of 388 amino acid residues. Recombinant lipase was successfully overexpressed with a 178-fold increase in activity compared to crude native L2 lipase. The recombinant L2 lipase (43.2 kDa) was purified to homogeneity in a single chromatography step. The purified lipase was found to be reactive at a temperature range of 55-80 °C and at a pH of 6-10. The L2 lipase had a melting temperature (Tm) of 59.04 °C when analyzed by circular dichroism (CD) spectroscopy studies. The optimum activity was found to be at 70 °C and pH 9. Lipase L2 was strongly inhibited by ethylenediaminetetraacetic acid (EDTA) (100%), whereas phenylmethylsulfonyl fluoride (PMSF), pepstatin-A, 2-mercaptoethanol and dithiothreitol (DTT) inhibited the enzyme by over 40%. The CD spectra of secondary structure analysis showed that the L2 lipase structure contained 38.6% α-helices, 2.2% ß-strands, 23.6% turns and 35.6% random conformations.
Klebsiella pneumoniae causes neonatal sepsis and nosocomial infections. One of the strains, K. pneumoniae MGH 78578, shows high level of resistance to multiple microbial agents. In this study, domain family, amino acid sequence and topology analyses were performed on one of its hypothetical protein, YggG (KPN_03358). Structural bioinformatics approaches were used to predict the structure and functionality of YggG protein. The open reading frame (ORF) of yggG, which was a putative metalloprotease gene, was also cloned, expressed and characterized. The ORF was PCR amplified from K. pneumoniae MGH 78578 genomic DNA and cloned into a pET14-b vector for heterologous expression in Escherichia coli. The purified YggG protein was subsequently assayed for casein hydrolysis under different conditions. This protein was classified as peptidase M48 family and subclan gluzincin. It was predicted to contain one transmembrane domain by TMpred. Optimal protein expression was achieved by induction with 0.6 mM isopropyl thiogalactoside (IPTG) at 25 °C for six hours. YggG was purified as soluble protein and confirmed to be proteolytically active under the presence of 1.25 mM zinc acetate and showed optimum activity at 37 °C and pH 7.4. We confirmed for the first time that the yggG gene product is a zinc-dependent metalloprotease.
In recent years, studies on psychrophilic lipases have become an emerging area of research in the field of enzymology. The study described here focuses on the cold-adapted organic solvent tolerant lipase strain Pseudomonas sp. LSK25 isolated from Signy Station, South Orkney Islands, maritime Antarctic. Strain LSK25 lipase was successfully cloned, sequenced, and over-expressed in an Escherichia coli system. Sequence analysis revealed that the lipase gene of Pseudomonas sp. LSK25 consists of 1432 bp, lacks an N-terminal signal peptide and encodes a mature protein consisting of 476 amino acids. The recombinant LSK25 lipase was purified by single-step purification using Ni-Sepharose affinity chromatography and had a molecular mass of approximately 65 kDa. The final recovery and purification fold were 44% and 1.3, respectively. The LSK25 lipase was optimally active at 30 °C and at pH 6. Stable lipolytic activity was reported between temperatures of 5⁻30 °C and at pH 6⁻8. A significant enhancement of lipolytic activity was observed in the presence of Ca2+ ions, the organic lipids of rice bran oil and coconut oil, a synthetic C12 ester and a wide range of water immiscible organic solvents. Overall, lipase strain LSK25 is a potentially desirable candidate for biotechnological application, due to its stability at low temperatures, across a range of pH and in organic solvents.
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.
Bacteria of the genus Komagataeibacter are described to be the most noteworthy for having several of its species being efficient and strong cellulose producers. The 16S ribosomal RNA (rRNA) gene analysis is often used for the identification and taxonomic classification of these bacteria species. In order to observe the phylogenetic relationship among Komagataeibacter sp., twelve sequences of the 16S rRNA gene with three sequences each for species namely Komagataeibacter europaeus, Komagataeibacter hansenii, Komagataeibacter intermedius and Komagataeibacter xylinus were retrieved from NCBI GenBank database. The sequences were aligned and analysed using PAUP, OrthoANI and BLAST, followed by the phylogenetic tree construction using a Maximum Likelihood method. The parsimony character diagnostic analysis showed very few numbers of parsimony- informative characters present in the aligned sequences which is only 1.5% of the total characters. The inferred phylogenetic relationships demonstrated the unexpected positioning of K. xylinus (GQ240638: Gluconacetobacter xylinus strain) and K. xylinus (KC11853: G. xylinus strain) into the clades of K. europaeus and K. hansenii respectively. The also very low bootstrap values of the branch points linking the K. europaeus species indicated low support for the produced topologies. The findings of this study indicate that more phylogenies information can be attained by increasing the taxon sampling. In addition, more robust molecular data are needed to infer the phylogenetic relationships between the Komagataeibacter species more accurately.
Four ADP-glucose pyrophosphorylase cDNA clones were isolated from mature leaves and pith of sago palm by the polymerase chain reaction (PCR) technique. Three of them (agpp10, agpp12 and agpl19) encoded the AGP large subunit, while the fourth clone (agpl1) encoded the small subunit. agpp10 and agpp12 were isolated from pith, agpl19 was isolated from mature leaves, while agpl1 from both tissues. In addition, a full-length cDNA of agpl1 was successfully isolated from a cDNA library of mature leaves by a PCR-based screening technique. Semi-quantitative analysis suggests that agpp10 and agpp12 were detectable only in pith, agpl19 only in leaves, while agpl1 was expressed in both leaves and pith tissues.
Protein structure alignment and comparisons that are based on an alphabetical demonstration of protein structure are more simple to run with faster evaluation processes; thus, their accuracy is not as reliable as three-dimension (3D)-based tools. As a 1D method candidate, TS-AMIR used the alphabetic demonstration of secondary-structure elements (SSE) of proteins and compared the assigned letters to each SSE using the [Formula: see text]-gram method. Although the results were comparable to those obtained via geometrical methods, the SSE length and accuracy of adjacency between SSEs were not considered in the comparison process. Therefore, to obtain further information on accuracy of adjacency between SSE vectors, the new approach of assigning text to vectors was adopted according to the spherical coordinate system in the present study. Moreover, dynamic programming was applied in order to account for the length of SSE vectors. Five common datasets were selected for method evaluation. The first three datasets were small, but difficult to align, and the remaining two datasets were used to compare the capability of the proposed method with that of other methods on a large protein dataset. The results showed that the proposed method, as a text-based alignment approach, obtained results comparable to both 1D and 3D methods. It outperformed 1D methods in terms of accuracy and 3D methods in terms of runtime.
A genetic similarity algorithm is introduced in this study to find a group of semantically similar Gene Ontology terms. The genetic similarity algorithm combines semantic similarity measure algorithm with parallel genetic algorithm. The semantic similarity measure algorithm is used to compute the similitude strength between the Gene Ontology terms. Then, the parallel genetic algorithm is employed to perform batch retrieval and to accelerate the search in large search space of the Gene Ontology graph. The genetic similarity algorithm is implemented in the Gene Ontology browser named basic UTMGO to overcome the weaknesses of the existing Gene Ontology browsers which use a conventional approach based on keyword matching. To show the applicability of the basic UTMGO, we extend its structure to develop a Gene Ontology -based protein sequence annotation tool named extended UTMGO. The objective of developing the extended UTMGO is to provide a simple and practical tool that is capable of producing better results and requires a reasonable amount of running time with low computing cost specifically for offline usage. The computational results and comparison with other related tools are presented to show the effectiveness of the proposed algorithm and tools.
A Cucumber green mottle mosaic virus (CGMMV) was used to present a truncated dengue virus type 2 envelope (E) protein binding region from amino acids 379 to 423 (EB4). The EB4 gene was inserted at the terminal end of the CGMMV coat protein (CP) open reading frame (ORF). Read-through sequences of TMV or CGMMV, CAA-UAG-CAA-UUA, or AAA-UAG-CAA-UUA were, respectively, inserted in between the CP and the EB4 genes. The chimeric clones, pRT, pRG, and pCG+FSRTRE, were transcribed into full-length capped recombinant CGMMV transcripts. Only constructs with the wild-type CGMMV read-through sequence yielded infectious viruses following infection of host plant, muskmelon (Cucumis melo) leaves. The ratio of modified to unmodified CP for the read-through expression clone developed was also found to be approximately 1:1, higher than what has been previously reported. It was also observed that infectivity was not affected by differences in pI between the chimera and its wild counterpart. Analysis of recombinant viruses after 21-days-postinculation (dpi) revealed that deletions occurred resulting in partial reversions of the viral population to near wild type and suggesting that this would be the limiting harvest period for obtaining true to type recombinants with this construct.
Fungi of the Trichoderma species are valued industrial enzymes in support of the 'zero-waste' technology to convert agro-industrial biomass into valuable products, i.e. nanocellulose (NC). In this study, an in silico approach using substrate docking and molecular dynamic (MD) simulation was used to predict the order of which the multilayers of cellulosic polymers, i.e. lignin, hemicellulose and cellulose in oil palm leaves (OPL) are degraded by fungal enzymes, endocellulase and exocellulase. The study aimed to establish the catalytic tendencies of the enzymes to optimally degrade the cellulosic components of OPL for high yield production of NC. Energy minimized endocellulase and exocellulase models revealed satisfactory scores of PROCHECK (90.0% and 91.2%), Verify3D (97.23% and 98.85%) and ERRAT (95.24% and 91.00%) assessments. Active site prediction by blind docking, COACH meta-server and multiple sequence alignment indicated the catalytic triads for endocellulase and exocellulase were Ser116-His205-Glu249 and Ser382-Arg124-Asp385, respectively. Binding energy of endocellulase docked with hemicellulose (-6.0 kcal mol-1) was the most favourable followed by lignin (-5.6 kcal mol-1) and cellulose (-4.4 kcal mol-1). Exocellulase, contrarily, bonded favorably with lignin (-8.7 kcal mol-1), closely followed by cellulose (-8.5 kcal mol-1) and hemicellulose (-8.4 kcal mol-1). MDs simulations showed that interactions of complexes, endocellulase-hemicellulose and the exocellulase-cellulose being the most stable. Thus, the findings of the study successfully identified the specific actions of sugar-acting enzymes for NC production. Communicated by Ramaswamy H. Sarma.
DARPP-32 (dopamine and adenosine 3', 5'-monophosphate-regulated phosphoprotein of 32 kDa), which belongs to PPP1R1 gene family, is known to act as an important integrator in dopamine-mediated neurotransmission via the inhibition of protein phosphatase-1 (PP1). Besides its neuronal roles, this protein also behaves as a key player in pathological and pharmacological aspects. Use of bioinformatics and phylogenetics approaches to further characterize the molecular features of DARPP-32 can guide future works. Predicted phosphorylation sites on DARPP-32 show conservation across vertebrates. Phylogenetics analysis indicates evolutionary strata of phosphorylation site acquisition at the C-terminus, suggesting functional expansion of DARPP-32, where more diverse signalling cues may involve in regulating DARPP-32 in inhibiting PP1 activity. Moreover, both phylogenetics and synteny analyses suggest de novo origination of PPP1R1 gene family via chromosomal rearrangement and exonization.
The polyhydroxyalkanoate (PHA) producing capability of four bacterial strains isolated from Antarctica was reported in a previous study. This study analyzed the PHA synthase genes and the PHA-associated gene clusters from the two antarctic Pseudomonas isolates (UMAB-08 and UMAB-40) and the two antarctic Janthinobacterium isolates (UMAB-56 and UMAB-60) through whole-genome sequence analysis. The Pseudomonas isolates were found to carry PHA synthase genes which fall into two different PHA gene clusters, namely Class I and Class II, which are involved in the biosynthesis of short-chain-length-PHA (SCL-PHA) and medium-chain-length-PHA (MCL-PHA), respectively. On the other hand, the Janthinobacterium isolates carry a Class I and an uncharacterized putative PHA synthase genes. No other gene involved in PHA synthesis was detected in close proximity to the uncharacterized putative PHA synthase gene in the Janthinobacterium isolates, therefore it falls into a separate clade from the ordinary Class I, II, III and IV clades of PHA synthase (PhaC) phylogenetic tree. Multiple sequence alignment showed that the uncharacterized putative PHA synthase gene contains all the highly conserved amino acid residues and the proposed catalytic triad of PHA synthase. PHA biosynthesis and in vitro PhaC enzymatic assay results showed that this uncharacterized putative PHA synthase from Janthinobacterium sp. UMAB-60 is funtional. This report adds new knowledge to the PHA synthase database as we describe scarce information of PHA synthase genes and PHA-associated gene clusters from the antarctic bacterial isolates (extreme and geographically isolated environment) and comparing with those from non-antarctic PHA-producing bacteria.
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.