Protein structure prediction from amino acid sequence has been one of the most challenging aspects in computational structural biology despite significant progress in recent years showed by critical assessment of protein structure prediction (CASP) experiments. When experimentally determined structures are unavailable, the predictive structures may serve as starting points to study a protein. If the target protein consists of homologous region, high-resolution (typically <1.5 Å) model can be built via comparative modelling. However, when confronted with low sequence similarity of the target protein (also known as twilight-zone protein, sequence identity with available templates is less than 30%), the protein structure prediction has to be initiated from scratch. Traditionally, twilight-zone proteins can be predicted via threading or ab initio method. Based on the current trend, combination of different methods brings an improved success in the prediction of twilight-zone proteins. In this mini review, the methods, progresses and challenges for the prediction of twilight-zone proteins were discussed.
Angiostrongylus mackerrasae is a parasitic nematode of rats found in Australia. When first reported, it was referred to as A. cantonensis. Recent molecular studies, including the mitochondrial genome, indicate that it is highly similar to A. cantonensis. These studies did not include A. malaysiensis, another member of the A. cantonensis species complex, for comparison. The present study examined the genetic distance and phylogenetic relationship between the component taxa (A. cantonensis, A. mackerrasae and A. malaysiensis) of the A. cantonensis species complex, based on the 12 protein-coding genes (PCGs) of their mitochondrial genome. Both the nucleotide and amino acid sequences were analysed. Angiostrongylus mackerrasae and A. cantonensis are members of the same genetic lineage and both are genetically distinct from A. malaysiensis. The genetic distance based on concatenated nucleotide sequences of 12 mt-PCGs between A. mackerrasae and A. cantonensis from Thailand is p = 1.73%, while that between the Thai and Chinese taxa of A. cantonensis is p = 3.52%; the genetic distance between A. mackerrasae and A. cantonensis from China is p = 3.70%. The results indicate that A. mackerrasae and A. cantonensis belong to the same genetic lineage, and that A. mackerrasae may be conspecific with A. cantonensis. It remains to be resolved whether A. mackerrasae is conspecific with A. cantonensis or undergoing incipient speciation.
Filamentous phage random peptide libraries were used to identify the epitopes of Burkholderia pseudomallei protease by panning against IgG polyclonal sera that exhibited protease neutralizing properties. The isolated fusion peptides presented a consensus peptide sequence, TKSMALSG, which closely resembles part of the active site sequence, 435GTSMATPHVAG445, of B. pseudomallei serine metalloprotease. By comparing the consensus sequence, TKSMALSG, with the predicted three-dimensional molecular model of B. pseudomallei serine metalloprotease, it appears that the potential antibody binding epitope was buried within the molecule. This active site was conformational whereby one continuous sub-region (SMA) was located between two discontinuous sub-regions, supplied by the flanking residues in the same polypeptide. All phages selected from the biopanning with IgG polyclonal sera showed good binding towards the polyclonal antibodies when compared to the negative control. In addition, these peptide-bearing phages showed competitive inhibition of B. pseudomallei serine metalloprotease binding to the polyclonal IgG.
Interests in Acinetobacter haemolyticus lipases are showing an increasing trend concomitant with growth of the enzyme industry and the widening search for novel enzymes and applications. Here, we present a structural model that reveals the key catalytic residues of lipase KV1 from A. haemolyticus. Homology modeling of the lipase structure was based on the structure of a carboxylesterase from the archaeon Archaeoglobus fulgidus as the template, which has a sequence that is 58% identical to that of lipase KV1. The lipase KV1 model is comprised of a single compact domain consisting of seven parallel and one anti-parallel β-strand surrounded by nine α-helices. Three structurally conserved active-site residues, Ser165, Asp259, and His289, and a tunnel through which substrates access the binding site were identified. Docking of the substrates tributyrin and palmitic acid into the pH 8 modeled lipase KV1 active sites revealed an aromatic platform responsible for the substrate recognition and preference toward tributyrin. The resulting binding modes from the docking simulation correlated well with the experimentally determined hydrolysis pattern, for which pH 8 and tributyrin being the optimum pH and preferred substrate. The results reported herein provide useful insights into future structure-based tailoring of lipase KV1 to modulate its catalytic activity.
Owing to their ability to maintain a thermodynamically stable fold at extremely high temperatures, thermophilic proteins (TTPs) play a critical role in basic research and a variety of applications in the food industry. As a result, the development of computation models for rapidly and accurately identifying novel TTPs from a large number of uncharacterized protein sequences is desirable. In spite of existing computational models that have already been developed for characterizing thermophilic proteins, their performance and interpretability remain unsatisfactory. We present a novel sequence-based thermophilic protein predictor, termed SCMTPP, for improving model predictability and interpretability. First, an up-to-date and high-quality dataset consisting of 1853 TPPs and 3233 non-TPPs was compiled from published literature. Second, the SCMTPP predictor was created by combining the scoring card method (SCM) with estimated propensity scores of g-gap dipeptides. Benchmarking experiments revealed that SCMTPP had a cross-validation accuracy of 0.883, which was comparable to that of a support vector machine-based predictor (0.906-0.910) and 2-17% higher than that of commonly used machine learning models. Furthermore, SCMTPP outperformed the state-of-the-art approach (ThermoPred) on the independent test dataset, with accuracy and MCC of 0.865 and 0.731, respectively. Finally, the SCMTPP-derived propensity scores were used to elucidate the critical physicochemical properties for protein thermostability enhancement. In terms of interpretability and generalizability, comparative results showed that SCMTPP was effective for identifying and characterizing TPPs. We had implemented the proposed predictor as a user-friendly online web server at http://pmlabstack.pythonanywhere.com/SCMTPP in order to allow easy access to the model. SCMTPP is expected to be a powerful tool for facilitating community-wide efforts to identify TPPs on a large scale and guiding experimental characterization of TPPs.
Information about Plasmodium malariae is scanty worldwide due to its "benign" nature and low infection rates. Consequently, studies on the genetic polymorphisms of P. malariae are lacking. Here, we report genetic polymorphisms of 28 P. malariae circumsporozoite protein (Pmcsp) isolates from Malaysia which were compared with those in other regions in Asia as well as those from Africa. Phylogenetic analysis revealed that most Malaysian P. malariae isolates clustered together but independently from other Asian isolates. Low nucleotide diversity was observed in Pmcsp non-repeat regions in contrast to high nucleotide diversity observed in non-repeat regions of Plasmodium knowlesi CSP gene, the current major cause of malaria in Malaysia. This study contributes to the characterisation of naturally occurring polymorphisms in the P. malariae CSP gene.
The aim of this study was to determine whether mutations could occur in the dengue virus genome following three subpassages of the virus in a mosquito cell line. This was done because sources of virus isolates used for sequencing studies are usually maintained in cell lines rather than in patients' sera. Therefore it must be assured that no mutation occurred during the passaging. For this purpose, sequencing was carried out using the polymerase chain reaction (PCR) products of the envelope/non-structural protein 1 junction region (280 nucleotides) of dengue type 3 virus. Sequence data were compared between the virus from a patient's serum against the virus subpassaged three times in the C6/36 cell line. We found that the sequence data of the virus from serum was identical to the virus that was subpassaged three times in C6/36 cell line.
Human leukocyte antigen (HLA) is a key genetic factor conferring risk of systemic lupus erythematosus (SLE), but precise independent localization of HLA effects is extremely challenging. As a result, the contribution of specific HLA alleles and amino-acid residues to the overall risk of SLE and to risk of specific autoantibodies are far from completely understood. Here, we dissected (a) overall SLE association signals across HLA, (b) HLA-peptide interaction, and (c) residue-autoantibody association. Classical alleles, SNPs, and amino-acid residues of eight HLA genes were imputed across 4,915 SLE cases and 13,513 controls from Eastern Asia. We performed association followed by conditional analysis across HLA, assessing both overall SLE risk and risk of autoantibody production. DR15 alleles HLA-DRB1*15:01 (P = 1.4x10-27, odds ratio (OR) = 1.57) and HLA-DQB1*06:02 (P = 7.4x10-23, OR = 1.55) formed the most significant haplotype (OR = 2.33). Conditioned protein-residue signals were stronger than allele signals and mapped predominantly to HLA-DRB1 residue 13 (P = 2.2x10-75) and its proxy position 11 (P = 1.1x10-67), followed by HLA-DRB1-37 (P = 4.5x10-24). After conditioning on HLA-DRB1, novel associations at HLA-A-70 (P = 1.4x10-8), HLA-DPB1-35 (P = 9.0x10-16), HLA-DQB1-37 (P = 2.7x10-14), and HLA-B-9 (P = 6.5x10-15) emerged. Together, these seven residues increased the proportion of explained heritability due to HLA to 2.6%. Risk residues for both overall disease and hallmark autoantibodies (i.e., nRNP: DRB1-11, P = 2.0x10-14; DRB1-13, P = 2.9x10-13; DRB1-30, P = 3.9x10-14) localized to the peptide-binding groove of HLA-DRB1. Enrichment for specific amino-acid characteristics in the peptide-binding groove correlated with overall SLE risk and with autoantibody presence. Risk residues were in primarily negatively charged side-chains, in contrast with rheumatoid arthritis. We identified novel SLE signals in HLA Class I loci (HLA-A, HLA-B), and localized primary Class II signals to five residues in HLA-DRB1, HLA-DPB1, and HLA-DQB1. These findings provide insights about the mechanisms by which the risk residues interact with each other to produce autoantibodies and are involved in SLE pathophysiology.
The Rab GTPase family plays a vital role in several plant physiological processes including fruit ripening. Fruit softening during ripening involves trafficking of cell wall polymers and enzymes between cellular compartments. Mango, an economically important fruit crop, is known for its delicious taste, exotic flavour and nutritional value. So far, there is a paucity of information on the mango Rab GTPase family. In this study, 23 genes encoding Rab proteins were identified in mango by a comprehensive in silico approach. Sequence alignment and similarity tree analysis with the model plant Arabidopsis as a reference enabled the bona fide assignment of the deduced mango proteins to classify into eight subfamilies. Expression analysis by RNA-Sequencing (RNA-Seq) showed that the Rab genes were differentially expressed in ripe and unripe mangoes suggesting the involvement of vesicle trafficking during ripening. Interaction analysis showed that the proteins involved in vesicle trafficking and cell wall softening were interconnected providing further evidence of the involvement of the Rab GTPases in fruit softening. Correlation analyses showed a significant relationship between the expression level of the RabA3 and RabA4 genes and fruit firmness at the unripe stage of the mango varieties suggesting that the differences in gene expression level might be associated with the contrasting firmness of these varieties. This study will not only provide new insights into the complexity of the ripening-regulated molecular mechanism but also facilitate the identification of potential Rab GTPases to address excessive fruit softening.
Lignosus rhinocerotis (Tiger milk mushroom) is an important folk medicine for indigenous peoples in Southeast Asia. We previously reported its de novo assembled 34.3 Mb genome encoding a repertoire of proteins including a putative bioactive fungal immunomodulatory protein. Here we report the cDNA of this new member (FIP-Lrh) with a homology range of 54-64% to FIPs from other mushroom species, the closest is with FIP-glu (LZ-8) (64%) from Ganoderma lucidum. The FIP-Lrh of 112 amino acids (12.59 kDa) has a relatively hydrophobic N-terminal. Its predicted 3-dimensional model has identical folding patterns to FIP-fve and contains a partially conserved and more positively charged carbohydrates binding pocket. Docking predictions of FIP-Lrh on 14 glycans commonly found on cellular surfaces showed the best binding energy of -3.98 kcal/mol to N-acetylgalactosamine and N-acetylglucosamine. Overexpression of a 14.9 kDa soluble 6xHisFIP-Lrh was achieved in pET-28a(+)/BL21 and the purified recombinant protein was sequence verified by LC-MS/MS (QTOF) analysis. The ability to haemagglutinate both mouse and human blood at concentration ≥0.34 μM, further demonstrated its lectin nature. In addition, the cytotoxic effect of 6xHisFIP-Lrh on MCF-7, HeLa and A549 cancer cell lines was detected at IC50 of 0.34 μM, 0.58 μM and 0.60 μM, respectively.
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.
The mantled floral phenotype of oil palm (Elaeis guineensis) affects somatic embryogenesis-derived individuals and is morphologically similar to mutants defective in the B-class MADS-box genes. This somaclonal variation has been previously demonstrated to be associated to a significant deficit in genome-wide DNA methylation. In order to elucidate the possible role of DNA methylation in the transcriptional regulation of EgDEF1, the APETALA3 ortholog of oil palm, we studied this epigenetic mark within the gene in parallel with transcript accumulation in both normal and mantled developing inflorescences. We also examined the methylation and expression of two neighboring retrotransposons that might interfere with EgDEF1 regulation. We show that the EgDEF1 gene is essentially unmethylated and that its methylation pattern does not change with the floral phenotype whereas expression is dramatically different, ruling out a direct implication of DNA methylation in the regulation of this gene. Also, we find that both the gypsy element inserted within an intron of the EgDEF1 gene and the copia element located upstream from the promoter are heavily methylated and show little or no expression. Interestingly, we identify a shorter, alternative transcript produced by EgDEF1 and characterize its accumulation with respect to its full-length counterpart. We demonstrate that, depending on the floral phenotype, the respective proportions of these two transcripts change differently during inflorescence development. We discuss the possible phenotypical consequences of this alternative splicing and the new questions it raises in the search for the molecular mechanisms underlying the mantled phenotype in the oil palm.
Glutathione transferases (GST) were purified from locally isolated bacteria, Acinetobacter calcoaceticus Y1, by glutathione-affinity chromatography and anion exchange, and their substrate specificities were investigated. SDS-polyacrylamide gel electrophoresis revealed that the purified GST resolved into a single band with a molecular weight (MW) of 23 kDa. 2-dimensional (2-D) gel electrophoresis showed the presence of two isoforms, GST1 (pI 4.5) and GST2 (pI 6.2) with identical MW. GST1 was reactive towards ethacrynic acid, hydrogen peroxide, 1-chloro-2,4-dinitrobenzene, and trans,trans-hepta-2,4-dienal while GST2 was active towards all substrates except hydrogen peroxide. This demonstrated that GST1 possessed peroxidase activity which was absent in GST2. This study also showed that only GST2 was able to conjugate GSH to isoproturon, a herbicide. GST1 and GST2 were suggested to be similar to F0KLY9 (putative glutathione S-transferase) and F0KKB0 (glutathione S-transferase III) of Acinetobacter calcoaceticus strain PHEA-2, respectively.
The structure of a novel psychrophilic β-mannanase enzyme from Glaciozyma antarctica PI12 yeast has been modelled and analysed in detail. To our knowledge, this is the first attempt to model a psychrophilic β-mannanase from yeast. To this end, a 3D structure of the enzyme was first predicted using a threading method because of the low sequence identity (<30%) using MODELLER9v12 and simulated using GROMACS at varying low temperatures for structure refinement. Comparisons with mesophilic and thermophilic mannanases revealed that the psychrophilic mannanase contains longer loops and shorter helices, increases in the number of aromatic and hydrophobic residues, reductions in the number of hydrogen bonds and salt bridges and numerous amino acid substitutions on the surface that increased the flexibility and its efficiency for catalytic reactions at low temperatures.
The monkey malaria parasite Plasmodium knowlesi is now recognized as the fifth species of Plasmodium that can cause human malaria. Like the region II of the Duffy binding protein of P. vivax (PvDBPII), the region II of the P. knowlesi Duffy binding protein (PkDBPαII) plays an essential role in the parasite's invasion into the host's erythrocyte. Numerous polymorphism studies have been carried out on PvDBPII, but none has been reported on PkDBPαII. In this study, the genetic diversity, haplotyes and allele groups of PkDBPαII of P. knowlesi clinical isolates from Peninsular Malaysia were investigated.
Angiotensin I-converting enzyme (ACE) inhibitors derived from foods are valuable auxiliaries to agents such as captopril. Eight highly functional ACE inhibitory peptides from the mushroom, Agaricus bisporus, were identified by LC-MS/MS. Among these peptides, the most potent ACE inhibitory activity was exhibited by AHEPVK, RIGLF and PSSNK with IC₅₀ values of 63, 116 and 129 μM, respectively. These peptides exhibited high ACE inhibitory activity after gastrointestinal digestion. Lineweaver-Burk plots suggested that AHEPVK and RIGLF act as competitive inhibitors against ACE, whereas PSSNK acts as a non-competitive inhibitor. Mushrooms can be a good component of dietary supplement due to their readily available source and, in addition, they rarely cause food allergy. Compared to ACE inhibitory peptides isolated from other edible mushrooms, AHEPVK, RIGLF and PSSNK have lower IC₅₀ values. Therefore, these peptides may serve as an ideal ingredient in the production of antihypertensive supplements.
The structural comparison of proteins is a vital step in structural biology that is used to predict and analyse a new unknown protein function. Although a number of different techniques have been explored, the study to develop new alternative methods is still an active research area. The present paper introduces a text modelling-based technique for the structural comparison of proteins. The method models the secondary and tertiary structure of proteins in two linear sequences and then applies them to the comparison of two structures. The technique used for pairwise comparison of the sequences has been adopted from computational linguistics and its well-known techniques for analysing and quantifying textual sequences. To this end, an n-gram modelling technique is used to capture regularities between sequences, and then, the cross-entropy concept is employed to measure their similarities. Several experiments are conducted to evaluate the performance of the method and compare it with other commonly used programs. The assessments for information retrieval evaluation demonstrate that the technique has a high running speed, which is similar to other linear encoding methods, such as 3D-BLAST, SARST, and TS-AMIR, whereas its accuracy is comparable to CE and TM-align, which are high accuracy comparison tools. Accordingly, the results demonstrate that the algorithm has high efficiency compared with other state-of-the-art methods.
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.
Bioassay and immunohistochemical studies have detected the presence of prosaposin in the central nervous system (CNS) of mammals. Here, first time, we have determined the partial cDNA sequence of pigeon prosaposin and mapped the distribution of its mRNA in the pigeon CNS. The predicted amino acid sequence of pigeon prosaposin showed 93 and 60% identity to chicken and human prosaposin, respectively. In situ hybridization, autoradiograms showed that the prosaposin mRNA expression was found in the olfactory bulb, prepiriform cortex, Wulst, mesopallium, nidopallium, hippocampal formation, thalamus, tuberis nucleus, pre-tectal nucleus, nucleus mesencephalicus lateralis, pars dorsalis, nucleus isthmi, pars parvocellularis and magnocellularis, Edinger-Westphal nucleus, optic tectum, cerebellar cortex and nuclei, vestibular nuclei and gray matter of the spinal cord. These results suggest that the cDNA sequence of pigeon prosaposin is comparable to other vertebrates, and the general distribution pattern of prosaposin mRNA resembles those are found in mammals.
The structure of psychrophilic chitinase (CHI II) from Glaciozyma antarctica PI12 has yet to be studied in detail. Due to its low sequence identity (<30 %), the structural prediction of CHI II is a challenge. A 3D model of CHI II was built by first using a threading approach to search for a suitable template and to generate an optimum target-template alignment, followed by model building using MODELLER9v7. Analysis of the catalytic insertion domain structure in CHI II revealed an increase in the number of aromatic residues and longer loops compared to mesophilic and thermophilic chitinases. A molecular dynamics simulation was used to examine the stability of the CHI II structure at 273, 288 and 300 K. Structural analysis of the substrate-binding cleft revealed a few exposed aromatic residues. Substitutions of certain amino acids in the surface and loop regions of CHI II conferred an increased flexibility to the enzyme, allowing for an adaptation to cold temperatures. A substrate binding comparison of CHI II with the mesophilic chitinase from Coccidioides immitis, 1D2K, suggested that the psychrophilic adaptation and catalytic activity at low temperatures were achieved through a reduction in the number of salt bridges, fewer hydrogen bonds and an increase in the exposure of the hydrophobic side chains to the solvent.