Multiple sequence alignment was performed against eight proteases from the Flaviviridae family using ClustalW to illustrate conserved domains. Two sets of prediction approaches were applied and the results compared. Firstly, secondary structure prediction was performed using available structure prediction servers. The second approach made use of the information on the secondary structures extracted from structure prediction servers, threading techniques and DSSP database of some of the templates used in the threading techniques. Consensus on the one-dimensional secondary structure of Den2 protease was obtained from each approach and evaluated against data from the recently crystallised Den2 NS2B/NS3 obtained from the Protein Data Bank (PDB). Results indicated the second approach to show higher accuracy compared to the use of prediction servers only. Thus, it is plausible that this approach is applicable to the initial stage of structural studies of proteins with low amino acid sequence homology against other available proteins in the PDB.
Phosphatidylinositol 4-phosphate 5-kinases (PIP5Ks) play diverse roles in the cellular biology of many organisms, including signal transduction, secretion and vesicular trafficking, and regulation of cytoskeleton assembly. Discovery of the PIP5K gene in Eimeria tenella may shed light on its role in the biology of this avian protozoan, and afford further understanding of the cell-host interaction, particularly during the invasion process. In this study, we report the identification of the PIP5K coding region in the genome sequence of Eimeria tenella using in silico gene prediction approaches. Prediction of the PIP5K coding sequence was confirmed by mapping the full-length cDNA sequence, generated via the Rapid Amplification of cDNA Ends (RACE) method, to the genomic sequence. The putative PIP5K gene of Eimeria tenella is located on the complementary strand of the E1080B12.b1 contig, and comprises 12 exons. Further analysis showed that the coding region spans from exon 1 to exon 7, with all exons obeying the adopted 'gt...ag' splicing rule of intronic sequences. Consensus of the hexameric 5' donor-splice site was deduced as GTRDBB... and the consensus for the 3' acceptor-splice sites as ...BHDYAG. The gene encodes a 252-amino acid residue protein. Domain search and protein fold recognition analyses provide compelling evidences that the deduced protein is a PIP5K.
To study the genetic diversity of re-emerging Plasmodium vivax in the Republic of Korea, nucleotide sequence variations at the merozoite surface protein-3alpha (PvMSP-3alpha) locus were analyzed using 24 re-emerging isolates and 4 isolates from imported cases. Compared with the well known Belem strain (Brazil), a large number of amino acid substitutions, deletions, and insertions were found at the locus of the isolates examined. The Korean isolates were divided into two allelic types; type I (15 isolates), similar to the Belem strain, and type II (9), similar to the Chess strain (New Guinea). Isolates from imported cases were classified into three types; type III (1 from Malaysia), similar to type B from western Thailand, type IV (1 each from Indonesia and India), and type V (1 from Pakistan), both being new types. Our results have shown that the MSP-3alpha locus of re-emerging Korean P. vivax is dimorphic with two allelic types coexisting in the endemic area.
Six papaya samples showing downward leaf curling were collected in Guangdong and Guangxi provinces, China. The result of TAS-ELISA showed they were all infected by geminiviruses. Comparison of partial DNA-A sequences reveals that these virus isolates can be classified into two groups. Group I includes isolates G2, G4, G5, G28 and G29 from Guangxi province, while isolate GD2 from Guangdong province belongs to Group II. The complete DNA-A sequence of G2 and GD2 were characterized. Sequence comparisons showed that the DNA-A of G2 and GD2 were most closely related to that of Ageratum yellow vein China virus- [Hn2] and Ageratum yellow vein virus , respectively, with 83.4 and 75.2% nucleotide sequence identity, while DNA-A sequence between G2 and GD2 had only 73.4% sequence identity. The molecular data suggests that G2 and GD2 are two distinct begomoviruses, for which the name Papaya leaf curl China virus (PaLCuCNV) for G2 and Papaya leaf curl Guangdong virus (PaLCuGDV) for GD2 are proposed. Comparison of individual encoded proteins showed the coat protein of G2 and GD2 shared highest amino acid sequence identity (97.7 and 94.2%, respectively) with that of Pepper leaf curl virus -[Malaysia] (PepLCV-[MY]), suggesting the CP of these viruses may have identical ancestor.
Association analysis was conducted in a core collection of 94 genotypes of Solanum pimpinellifolium to identify variations linked to salt tolerance traits (physiological and yield traits under salt stress) in four candidate genes viz., DREB1A, VP1.1, NHX1, and TIP. The candidate gene analysis covered a concatenated length of 4594 bp per individual and identified five SNP/Indels in DREB1A and VP1.1 genes explaining 17.0% to 25.8% phenotypic variation for various salt tolerance traits. Out of these five alleles, one at 297 bp in DREB1A had in-frame deletion of 6 bp (CTGCAT) or 12 bp (CTGCATCTGCAT), resulting in two alleles, viz., SpDREB1A_297_6 and SpDREB1A_297_12. These alleles individually or as haplotypes accounted for maximum phenotypic variance of about 25% for various salt tolerance traits. Design of markers for selection of the favorable alleles/haplotypes will hasten marker-assisted introgression of salt tolerance from S. pimpinellifolium into cultivated tomato.
Infection with Plasmodium knowlesi, a zoonotic primate malaria, is a growing human health problem in Southeast Asia. P. knowlesi is being used in malaria vaccine studies, and a number of proteins are being considered as candidate malaria vaccine antigens, including the Apical Membrane Antigen 1 (AMA1). In order to determine genetic diversity of the ama1 gene and to identify epitopes of AMA1 under strongest immune selection, the ama1 gene of 52 P. knowlesi isolates derived from human infections was sequenced. Sequence analysis of isolates from two geographically isolated regions in Sarawak showed that polymorphism in the protein is low compared to that of AMA1 of the major human malaria parasites, P. falciparum and P. vivax. Although the number of haplotypes was 27, the frequency of mutations at the majority of the polymorphic positions was low, and only six positions had a variance frequency higher than 10%. Only two positions had more than one alternative amino acid. Interestingly, three of the high-frequency polymorphic sites correspond to invariant sites in PfAMA1 or PvAMA1. Statistically significant differences in the quantity of three of the six high frequency mutations were observed between the two regions. These analyses suggest that the pkama1 gene is not under balancing selection, as observed for pfama1 and pvama1, and that the PkAMA1 protein is not a primary target for protective humoral immune responses in their reservoir macaque hosts, unlike PfAMA1 and PvAMA1 in humans. The low level of polymorphism justifies the development of a single allele PkAMA1-based vaccine.
The DNA of three biological variants, G1, Ic and G2, which originated from the same greenhouse isolate of rice tungro bacilliform virus (RTBV) at the International Rice Research Institute (IRRI), was cloned and sequenced. Comparison of the sequences revealed small differences in genome sizes. The variants were between 95 and 99% identical at the nucleotide and amino acid levels. Alignment of the three genome sequences with those of three published RTBV sequences (Phi-1, Phi-2 and Phi-3) revealed numerous nucleotide substitutions and some insertions and deletions. The published RTBV sequences originated from the same greenhouse isolate at IRRI 20, 11 and 9 years ago. All open reading frames (ORFs) and known functional domains were conserved across the six variants. The cysteine-rich region of ORF3 showed the greatest variation. When the six DNA sequences from IRRI were compared with that of an isolate from Malaysia (Serdang), similar changes were observed in the cysteine-rich region in addition to other nucleotide substitutions and deletions across the genome. The aligned nucleotide sequences of the IRRI variants and Serdang were used to analyse phylogenetic relationships by the bootstrapped parsimony, distance and maximum-likelihood methods. The isolates clustered in three groups: Serdang alone; Ic and G1; and Phi-1, Phi-2, Phi-3 and G2. The distribution of phylogenetically informative residues in the IRRI sequences shared with the Serdang sequence and the differing tree topologies for segments of the genome suggested that recombination, as well as substitutions and insertions or deletions, has played a role in the evolution of RTBV variants. The significance and implications of these evolutionary forces are discussed in comparison with badnaviruses and caulimoviruses.
Recently, a new paramyxovirus, now known as Nipah virus (NV), emerged in Malaysia and Singapore, causing fatal encephalitis in humans and a respiratory syndrome in pigs. Initial studies had indicated that NV is antigenically and genetically related to Hendra virus (HV). We generated the sequences of the N, P/C/V, M, F, and G genes of NV and compared these sequences with those of HV and other members of the family Paramyxoviridae. The intergenic regions of NV were identical to those of HV, and the gene start and stop sequences of NV were nearly identical to those of HV. The open reading frames (ORFs) for the V and C proteins within the P gene were found in NV, but the ORF encoding a potential short basic protein found in the P gene of HV was not conserved in NV. The N, P, C, V, M, F, and G ORFs in NV have nucleotide homologies ranging from 88% to 70% and predicted amino acid homologies ranging from 92% to 67% in comparison with HV. The predicted fusion cleavage sequence of the F protein of NV had a single amino acid substitution (K to R) in comparison with HV. Phylogenetic analysis demonstrated that although HV and NV are closely related, they are clearly distinct from any of the established genera within the Paramyxoviridae and should be considered a new genus.
Two hundred and forty nucleotides from the pre-M gene region of 10 Japanese encephalitis (JE) virus strains isolated in Malaysia in 1992 were sequenced and compared with the other JE virus strains from different geographic areas in Asia. Our JE virus strains belong to the largest genotypic group that includes strains isolated in temperate regions such as Japan, China, and Taiwan. Our Malaysian JE virus strains differed in 32 nucleotides (13.3%) from WTP/70/22 strain isolated from Malaysia in 1970, which belonged to another distinct genotypic group.
DNA mismatch repair proteins play an essential role in maintaining genomic integrity during replication and genetic recombination. We successfully isolated a full length MSH2 and partial MSH7 cDNAs from tomato, based on sequence similarity between MutS and plant MSH homologues. Semi-quantitative RT-PCR reveals higher levels of mRNA expression of both genes in young leaves and floral buds. Genetic mapping placed MSH2 and MSH7 on chromosomes 6 and 7, respectively, and indicates that these genes exist as single copies in the tomato genome. Analysis of protein sequences and phylogeny of the plant MSH gene family show that these proteins are evolutionarily conserved, and follow the classical model of asymmetric protein evolution. Genetic manipulation of the expression of these MSH genes in tomato will provide a potentially useful tool for modifying genetic recombination and hybrid fertility between wide crosses.
Presynaptic neurotoxins are one of the major components in Bungarus venom. Unlike other Bungarus species that have been studied, β-bungarotoxin has never been isolated from Bungarus fasciatus venom. It was hypothesized that the absence of β-bungarotoxin in this species was due to divergence during evolution prior to evolution of β-bungarotoxin. In this study, we have isolated a β-bungarotoxin isoform we named P-elapitoxin-Bf1a by using gel filtration, cation-exchange and reverse-phase chromatography from Malaysian B. fasciatus venom. The toxin consists of two heterogeneous subunits, subunit A and subunit B. LCMS/MS data showed that subunit A was homologous to acidic phospholipase A2 subunit A3 from Bungarus candidus and B. multicinctus venoms, whereas subunit B was homologous with subunit B1 from B. fasciatus venom that was previously detected by cDNA cloning. The toxin showed concentration- and time-dependent reduction of indirect-twitches without affecting contractile responses to ACh, CCh or KCl at the end of experiment in the chick biventer preparation. Toxin modification with 4-BPB inhibited the neurotoxic effect suggesting the importance of His-48. Tissue pre-incubation with monovalent B. fasciatus (BFAV) or neuro-polyvalent antivenom (NPV), at the recommended titer, was unable to inhibit the twitch reduction induced by the toxin. This study indicates that Malaysian B. fasciatus venom has a unique β-bungarotoxin isoform which was not neutralized by antivenoms. This suggests that there might be other presynaptic neurotoxins present in the venom and there is a variation in the enzymatic neurotoxin composition in venoms from different localities.
The IntFOLD server provides a unified resource for the automated prediction of: protein tertiary structures with built-in estimates of model accuracy (EMA), protein structural domain boundaries, natively unstructured or disordered regions in proteins, and protein-ligand interactions. The component methods have been independently evaluated via the successive blind CASP experiments and the continual CAMEO benchmarking project. The IntFOLD server has established its ranking as one of the best performing publicly available servers, based on independent official evaluation metrics. Here, we describe significant updates to the server back end, where we have focused on performance improvements in tertiary structure predictions, in terms of global 3D model quality and accuracy self-estimates (ASE), which we achieve using our newly improved ModFOLD7_rank algorithm. We also report on various upgrades to the front end including: a streamlined submission process, enhanced visualization of models, new confidence scores for ranking, and links for accessing all annotated model data. Furthermore, we now include an option for users to submit selected models for further refinement via convenient push buttons. The IntFOLD server is freely available at: http://www.reading.ac.uk/bioinf/IntFOLD/.
There is very little information on the capacity of freshwater carnivorous fish to biosynthesize highly unsaturated fatty acids (HUFA). The striped snakehead fish (Channa striata) is a carnivorous species cultured inland of several Southeast Asian countries due to its pharmaceutical properties in wound healing enhancement. We described here the full-length cDNA cloning of a striped snakehead fatty acid desaturase (fads), which is responsible for desaturation of unsaturated fatty acids in the HUFA biosynthesis. Bioinformatics analysis reveals a protein coding region with length of 445 amino acids containing all characteristic features of desaturase enzyme, including a cytochrome b5-domain with the heme-binding motif, two transmembrane domains and three histidine-rich regions. The striped snakehead fads amino acid sequence shares high similarity with known fads of other teleosts. The mRNA expression of striped snakehead fads also showed an ontogenic-related increase in expression in 0-20 days after hatch larva. Using ISH, we localized the presence of fads in larva brain, liver and intestinal tissues.
CATH (https://www.cathdb.info) identifies domains in protein structures from wwPDB and classifies these into evolutionary superfamilies, thereby providing structural and functional annotations. There are two levels: CATH-B, a daily snapshot of the latest domain structures and superfamily assignments, and CATH+, with additional derived data, such as predicted sequence domains, and functionally coherent sequence subsets (Functional Families or FunFams). The latest CATH+ release, version 4.3, significantly increases coverage of structural and sequence data, with an addition of 65,351 fully-classified domains structures (+15%), providing 500 238 structural domains, and 151 million predicted sequence domains (+59%) assigned to 5481 superfamilies. The FunFam generation pipeline has been re-engineered to cope with the increased influx of data. Three times more sequences are captured in FunFams, with a concomitant increase in functional purity, information content and structural coverage. FunFam expansion increases the structural annotations provided for experimental GO terms (+59%). We also present CATH-FunVar web-pages displaying variations in protein sequences and their proximity to known or predicted functional sites. We present two case studies (1) putative cancer drivers and (2) SARS-CoV-2 proteins. Finally, we have improved links to and from CATH including SCOP, InterPro, Aquaria and 2DProt.
1. The receptors for relaxin in the rat atria and uterus were investigated and compared by use of a series of synthetic and native relaxin analogues. The assays used were the positive chronotropic and inotropic effects in rat spontaneously beating, isolated right atrium and electrically driven left atrium and the relaxation of K+ precontracted uterine smooth muscle. 2. Relaxin analogues with an intact A- and B-chain were active in producing powerful chronotropic and inotropic effects in the rat isolated atria at nanomolar concentrations. Single-chain analogues and structural homologues of relaxin such as human insulin and sheep insulin-like growth factor I had no agonist action and did not antagonize the effect of the B29 form of human gene 2 relaxin. 3. Shortening the B-chain carboxyl terminal of human gene 1 (B2-29) relaxin to B2-26 reduced the activity of the peptide and removal of another 2 amino acid residues (B2-24) abolished the activity. This suggests that the B-chain length may be important for determination of the activity of relaxin. More detailed studies are needed to determine the effect of progressive amino acid removal on the structure and the bioactivity of relaxin. 4. Porcine prorelaxin was as active as porcine relaxin on a molar basis, suggesting that the presence of the intact C-peptide did not affect the binding of the prorelaxin to the receptor to produce functional responses. 5. Relaxin caused relaxation of uterine longitudinal and circular smooth muscle precontracted with 40 mM K+. The pEC50 values for human gene 2 and porcine relaxins were lower than those in the atrial assay, but rat relaxin had similar pEC50 values in both atrial and uterine assays. Rat relaxin was significantly less potent than either human gene 2 or porcine relaxin in the atrial assay, but in the uterine assay they were equipotent. The results suggest that the relaxin receptor or the signalling pathway in rat atria may differ from that in the uterus.
Beta-amyloid precursor protein cleavage enzyme 1 (BACE1) and beta-amyloid precursor protein cleavage enzyme 2 (BACE2), members of aspartyl protease family, are close homologues and have high similarity in their protein crystal structures. However, their enzymatic properties differ leading to disparate clinical consequences. In order to identify the residues that are responsible for such differences, we used evolutionary trace (ET) method to compare the amino acid conservation patterns of BACE1 and BACE2 in several mammalian species. We found that, in BACE1 and BACE2 structures, most of the ligand binding sites are conserved which indicate their enzymatic property of aspartyl protease family members. The other conserved residues are more or less randomly localized in other parts of the structures. Four group-specific residues were identified at the ligand binding site of BACE1 and BACE2. We postulated that these residues would be essential for selectivity of BACE1 and BACE2 biological functions and could be sites of interest for the design of selective inhibitors targeting either BACE1 or BACE2.
GalNAc-T1, a key candidate of GalNac-transferases genes family that is involved in mucin-type O-linked glycosylation pathway, is expressed in most biological tissues and cell types. Despite the reported association of GalNAc-T1 gene mutations with human disease susceptibility, the comprehensive computational analysis of coding, noncoding and regulatory SNPs, and their functional impacts on protein level, still remains unknown. Therefore, sequence- and structure-based computational tools were employed to screen the entire listed coding SNPs of GalNAc-T1 gene in order to identify and characterize them. Our concordant in silico analysis by SIFT, PolyPhen-2, PANTHER-cSNP, and SNPeffect tools, identified the potential nsSNPs (S143P, G258V, and Y414D variants) from 18 nsSNPs of GalNAc-T1. Additionally, 2 regulatory SNPs (rs72964406 and #x26; rs34304568) were also identified in GalNAc-T1 by using FastSNP tool. Using multiple computational approaches, we have systematically classified the functional mutations in regulatory and coding regions that can modify expression and function of GalNAc-T1 enzyme. These genetic variants can further assist in better understanding the wide range of disease susceptibility associated with the mucin-based cell signalling and pathogenic binding, and may help to develop novel therapeutic elements for associated diseases.
In this study, we have reported a full length of small heat shock protein 37 (designated MrHSP37) gene, identified from the transcriptome database of freshwater prawn Macrobrachium rosenbergii. The complete gene sequence of the MrHSP37 is 2,425 base pairs in length, and encodes 338 amino acids. MrHSP37 contains a long heat shock protein family profile in the amino acid sequence between 205 and 288. The mRNA expressions of MrHSP37 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). MrHSP37 is highly expressed in hepatopancreas and all the other tissues (walking leg, gills, muscle, stomach, haemocyte, intestine, pleopods, brain and eye stalk) of M. rosenbergii taken for analysis. The expression is strongly up-regulated after IHHNV challenge. To understand its biological activity, the recombinant MrHSP37 gene was constructed and expressed in Escherichia coli BL21 (DE3). The results of ATPase assay showed that the recombinant MrHSP37 protein exhibited apparent ATPase activity which increased with the concentration of the protein. And also the purified recombinant MrHSP37 protein was used for thermal aggregation assay (chaperone activity). It showed that the recombinant MrHSP37 protein is an active chaperone in this assay. Taken together, these results suggest that MrHSP37 is potentially involved in the immune responses against IHHNV challenge in M. rosenbergii.
Bats are increasingly being recognized as important reservoir hosts for a large number of viruses, some of them can be highly virulent when they infect human and livestock animals. Among the new bat zoonotic viruses discovered in recent years, several reoviruses (respiratory enteric orphan viruses) were found to be able to cause acute respiratory infections in humans, which included Melaka and Kampar viruses discovered in Malaysia, all of them belong to the genus Orthoreovirus, family Reoviridae. In this report, we describe the isolation of a highly related virus from an adult patient who suffered acute respiratory illness in Malaysia. Although there was no direct evidence of bat origin, epidemiological study indicated the potential exposure of the patient to bats before the onset of disease. The current study further demonstrates that spillover events of different strains of related orthoreoviruses from bats to humans are occurring on a regular basis, which calls for more intensive and systematic surveillances to fully assess the true public health impact of these newly discovered bat-borne zoonotic reoviruses.