Identification of phosphorylation sites is an important step in the function study and drug design of proteins. In recent years, there have been increasing applications of the computational method in the identification of phosphorylation sites because of its low cost and high speed. Most of the currently available methods focus on using local information around potential phosphorylation sites for prediction and do not take the global information of the protein sequence into consideration. Here, we demonstrated that the global information of protein sequences may be also critical for phosphorylation site prediction. In this paper, a new deep neural network model, called DeepPSP, was proposed for the prediction of protein phosphorylation sites. In the DeepPSP model, two parallel modules were introduced to extract both local and global features from protein sequences. Two squeeze-and-excitation blocks and one bidirectional long short-term memory block were introduced into each module to capture effective representations of the sequences. Comparative studies were carried out to evaluate the performance of DeepPSP, and four other prediction methods using public data sets The F1-score, area under receiver operating characteristic curves (AUROC), and area under precision-recall curves (AUPRC) of DeepPSP were found to be 0.4819, 0.82, and 0.50, respectively, for S/T general site prediction and 0.4206, 0.73, and 0.39, respectively, for Y general site prediction. Compared with the MusiteDeep method, the F1-score, AUROC, and AUPRC of DeepPSP were found to increase by 8.6, 2.5, and 8.7%, respectively, for S/T general site prediction and by 20.6, 5.8, and 18.2%, respectively, for Y general site prediction. Among the tested methods, the developed DeepPSP method was also found to produce best results for different kinase-specific site predictions including CDK, mitogen-activated protein kinase, CAMK, AGC, and CMGC. Taken together, the developed DeepPSP method may offer a more accurate phosphorylation site prediction by including global information. It may serve as an alternative model with better performance and interpretability for protein phosphorylation site prediction.
Genome sequences of members of a potential fourth rhinovirus (RV) species, provisionally denoted as rhinovirus A clade D, from patients with acute respiratory infection were determined. Bayesian coalescent analysis estimated that clade D emerged around the 1940s and diverged further around 2006-2007 into two distinctive sublineages (RV-A8-like and RV-A45-like) that harbored unique "clade-defining" substitutions. Similarity plots and bootscan mapping revealed a recombination breakpoint located in the 5'-UTR region of members of the RV-A8-like sublineage. Phylogenetic reconstruction revealed the distribution of clade D viruses in the Asia Pacific region and in Europe, underlining its worldwide distribution.
Vaccination is an essential component in controlling infectious bursal disease (IBD), however, there is a lack of information on the genetic characteristics of a recent infectious bursal disease virus (IBDV) that was isolated from IBD vaccinated commercial flocks in Malaysia. The present study investigated 11 IBDV isolates that were isolated from commercial poultry farms. The isolates were detected using reverse transcription-polymerase chain reaction (RT-PCR) targeting the hypervariable region (HVR) of VP2. Based on the HVR sequences, five isolates (IBS536/2017, IBS624/2017, UPM766/2018, UPM1056/2018, and UPM1432/2019) were selected for whole-genome sequencing using the MiSeq platform. The nucleotide and amino acid (aa) sequences were compared with the previously characterized IBDV strains. Deduced aa sequences of VP2HVR revealed seven isolates with 94-99% aa identity to very virulent strains (genogroup 3), two isolates with 97-100% aa identity to variant strains (genogroup 2), and two strains with 100% identity to the vaccine strain (genogroup 1) of IBDV. The phylogenetic analysis also showed that the isolates formed clusters with the respective genogroups. The characteristic motifs 222T, 249K, 286I, and 318D are typical of the variant strain and were observed for UPM1219/2019 and UPM1432/2019. In comparison, very virulent residues such as 222A, 249Q, 286T, and 318G were found for the vvIBDV, except for the UPM1056/2018 strain with a A222T substitution. In addition, the isolate has aa substitutions such as D213N, G254D, S315T, S317R, and A321E that are not commonly found in previously reported vvIBDV strains. Unlike the other vvIBDVs characterized in this study, UPM766/2018 lacks the MLSL aa residues in VP5. The aa tripeptides 145/146/147 (TDN) of VP1 were conserved for the vvIBDV, while a different motif, NED, was observed for the Malaysian variant strain. The phylogenetic tree showed that the IBDV variant clustered with the American and Chinese variant viruses and are highly comparable to the novel Chinese variants, with 99.9% identity. Based on the sequences and phylogenetic analyses, this is the first identification of an IBDV variant being reported in Malaysia. Further research is required to determine the pathogenicity of the IBDV variant and the protective efficacy of the current IBD vaccines being used against the virus.
The study was conducted to determine anti-tyrosinase and antioxidant activities of the extracted collagen hydrolysate (CH) derived from Malaysian jellyfish, Rhopilema hispidum. Collagen was extracted using 1:1 (w:v) 0.1 M NaOH solution at temperature 25 °C for 48 hr followed by treatment of 1:2 (w:v) distilled water for another 24 hr and freeze-dried. The extracted collagen was hydrolyzed using papain at optimum temperature, pH and enzyme/substrate ratio [E/S] of 60 °C, 7.0 and 1:50, respectively. CH was found to exhibit tyrosinase inhibitory activity, DPPH radical scavenging and metal ion-chelating assays up to 64, 28, and 83%, respectively, after 8 hr of hydrolysis process. The molecular weight of CH was found <10 kDa consisting of mainly Gly (19.219%), Glu (10.428%), and Arg (8.848%). The UV-visible spectrum analysis showed a major and minor peak at 218 and 276 nm, accordingly. The FTIR spectroscopy confirmed the amide groups in CH. The SEM images demonstrated spongy and porous structure of CH. In the cytotoxicity study, CH has no cytotoxicity against mouse embryonic 3T3 fibroblast cell line with IC50 value >500 µg/ml. Results revealed that the CH generated from this study has a potential to be developed as active ingredient in cosmeceutical application.
Dengue virus (DV) infection is one of the main public health concerns, affecting approximately 390 million people worldwide, as reported by the World Health Organization. Yet, there is no antiviral treatment for DV infection. Therefore, the development of potent and nontoxic anti-DV, as a complement for the existing treatment strategies, is urgently needed. Herein, we investigate a series of small peptides inhibitors of DV antiviral activity targeting the entry process as the promising strategy to block DV infection. The peptides were designed based on our previously reported peptide sequence, DN58opt (TWWCFYFCRRHHPFWFFYRHN), to identify minimal effective inhibitory sequence through molecular docking and dynamics studies. The in silico designed peptides were synthesized using conventional Fmoc solid-phase peptide synthesis chemistry, purified by RP-HPLC and characterized using LCMS. Later, they were screened for their antiviral activity. One of the peptides, AC 001, was able to reduce about 40% of DV plaque formation. This observation correlates well with the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) analysis - AC 001 showed the most favorable binding affinity through 60 ns simulations. Pairwise residue decomposition analysis has revealed four key residues that contributed to the binding of these peptides into the DV2 E protein pocket. This work identifies the minimal peptide sequence required to inhibit DV replication and explains the behavior observed on an atomic level using computational study.Communicated by Ramaswamy H. Sarma.
The alginate lyase AlyQ from Persicobacter sp. CCB-QB2 is a three-domained enzyme with a carbohydrate-binding module (CBM) from family 32. The CBM32 domain, AlyQB, binds enzymatically cleaved but not intact alginate. Co-crystallisation of AlyQB with the cleaved alginate reveals that it binds to the 4,5-unsaturated mannuronic acid of the non-reducing end. The binding pocket contains a conserved R248 that interacts with the sugar's carboxyl group, as well as an invariant W303 that stacks against the unsaturated pyranose ring. Targeting specifically the non-reducing end is more efficient than the reducing end since the latter consists of a mixture of mannuronic acid and guluronic acid. AlyQB also seems unable to bind these two saturated sugars as they contain OH groups that will clash with the pocket. Docking analysis of YeCBM32, which binds oligogalacturonic acid, shows that the stacking of the pyranose ring is shifted in order to accommodate the sugar's axial C1-OH, and its R69 is accordingly elevated to bind the sugar's carboxyl group. Unlike AlyQB, YeCBM32's binding pocket is able to accommodate both saturated and unsaturated galacturonic acid.
Previously, a hypothetical protein (HP) termed Bleg1_2437 (currently named Bleg1_2478) from Bacillus lehensis G1 was discovered to be an evolutionary divergent B3 subclass metallo-β-lactamase (MBL). Due to the scarcity of clinical inhibitors for B3 MBLs and the divergent nature of Bleg1_2478, this study aimed to design and characterise peptides as inhibitors against Bleg1_2478. Through in silico docking, RSWPWH and SSWWDR peptides with comparable binding energy to ampicillin were obtained. In vitro assay results showed RSWPWH and SSWWDR inhibited the activity of Bleg1_2478 by 50% at concentrations as low as 0.90 µM and 0.50 µM, respectively. At 10 µM of RSWPWH and 20 µM of SSWWDR, the activity of Bleg1_2478 was almost completely inhibited. Isothermal titration calorimetry (ITC) analyses showed slightly improved binding properties of the peptides compared to ampicillin. Docked peptide-protein complexes revealed that RSWPWH bound near the vicinity of the Bleg1_2478 active site while SSWWDR bound at the center of the active site itself. We postulate that the peptides caused the inhibition of Bleg1_2478 by reducing or blocking the accessibility of its active site from ampicillin, thus hampering its catalytic function.
We report two novel RNA viruses from yellow crazy ants, (Anoplolepis gracilipes) detected using next-generation sequencing. The complete genome sequences of the two viruses were 10,662 and 8,238 nucleotides in length, respectively, with both possessing two open reading frames with three conserved protein domains. The genome organization is characteristic of members of the genus Triatovirus in the family Dicistroviridae. The two novel viruses were tentatively named "Anoplolepis gracilipes virus 1" and "Anoplolepis gracilipes virus 2" (AgrV-1 and AgrV-2). Phylogenetic analyses based on amino acid sequences of the non-structural polyprotein (ORF1) suggest that the two viruses are triatovirus-like viruses. This is the first report on the discovery of novel triatovirus-like viruses in yellow crazy ants with a description of their genome structure (two ORFs and conserved domains of RNA helicase, RNA-dependent RNA polymerase, and capsid protein), complete sequences, and viral prevalence across the Asia-Pacific region.
The occurrences of multiple drug-resistant strains have been relentlessly increasing in recent years. The aquaculture industry has encountered major disease outbreaks and crucially affected by this situation. The usage of non-specific chemicals and antibiotics expedites the stimulation of resistant strains. Triggering the natural defense mechanism would provide an effective and safest way of protecting the host system. Hence, we have investigated the innate immune function of serine/threonine-protein kinase (STPK) in Macrobrachium rosenbergii (Mr). The in-silico protein analysis resulted in the identification of cationic antimicrobial peptide, MrSL-19, with interesting properties from STPK of M. rosenbergii. Antimicrobial assay, FACS and SEM analysis demonstrated that the peptide potentially inhibits Staphylococcus aureus by interacting with its membrane. The toxic study on MrSL-19 demonstrated that the peptide is not toxic against HEK293 cells as well as human erythrocytes. This investigation showed the significant innate immune property of an efficient cationic antimicrobial peptide, MrSL-19 of STPK from M. rosenbergii.
Little is known about the genetic features of Nipah virus (NiV) associated with virulence and transmission. Herein, phylogenetic and genetic analyses for all available NiV strains revealed sequence variations between the two genetic lineages of NiV with pathogenic differences, as well as among different strains within Bangladesh lineage. A total of 143 conserved amino acid differences, distributed among viral nucleocapsid (N), phosphoprotein (P), matrix protein (M), fusion protein (F) and glycoprotein (G), were revealed. Structural modeling revealed one key substitution (S3554N) in the viral G protein that might mediate a 12-amino-acid structural change from a loop into a β sheet. Multiple key amino acids substitutions in viral G protein were observed, which may alter viral fitness and transmissibility from bats to humans.
Inherited cutis laxa, or inelastic, sagging skin is a genetic condition of premature and generalised connective tissue ageing, affecting various elastic components of the extracellular matrix. Several cutis laxa syndromes are inborn errors of metabolism and lead to severe neurological symptoms. In a patient with cutis laxa, a choreoathetoid movement disorder, dysmorphic features and intellectual disability we performed exome sequencing to elucidate the underlying genetic defect. We identified the amino acid substitution R275W in phosphatidylinositol 4-kinase type IIα, caused by a homozygous missense mutation in the PI4K2A gene. We used lipidomics, complexome profiling and functional studies to measure phosphatidylinositol 4-phosphate synthesis in the patient and evaluated PI4K2A deficient mice to define a novel metabolic disorder. The R275W residue, located on the surface of the protein, is involved in forming electrostatic interactions with the membrane. The catalytic activity of PI4K2A in patient fibroblasts was severely reduced and lipid mass spectrometry showed that particular acyl-chain pools of PI4P and PI(4,5)P2 were decreased. Phosphoinositide lipids play a major role in intracellular signalling and trafficking and regulate the balance between proliferation and apoptosis. Phosphatidylinositol 4-kinases such as PI4K2A mediate the first step in the main metabolic pathway that generates PI4P, PI(4,5)P2 and PI(3,4,5)P3 . Although neurologic involvement is common, cutis laxa has not been reported previously in metabolic defects affecting signalling. Here we describe a patient with a complex neurological phenotype, premature ageing and a mutation in PI4K2A, illustrating the importance of this enzyme in the generation of inositol lipids with particular acylation characteristics.
COVID-19 caused by the novel SARS-CoV-2 has been declared a pandemic by the WHO is causing havoc across the entire world. As of May end, about 6 million people have been affected, and 367 166 have died from COVID-19. Recent studies suggest that the SARS-CoV-2 genome shares about 80% similarity with the SARS-CoV-1 while their protein RNA dependent RNA polymerase (RdRp) shares 96% sequence similarity. Remdesivir, an RdRp inhibitor, exhibited potent activity against SARS-CoV-2 in vitro. 3-Chymotrypsin like protease (also known as Mpro) and papain-like protease, have emerged as the potential therapeutic targets for drug discovery against coronaviruses owing to their crucial role in viral entry and host-cell invasion. Crystal structures of therapeutically important SARS-CoV-2 target proteins, namely, RdRp, Mpro, endoribonuclease Nsp15/NendoU and receptor binding domain of CoV-2 spike protein has been resolved, which have facilitated the structure-based design and discovery of new inhibitors. Furthermore, studies have indicated that the spike proteins of SARS-CoV-2 use the Angiotensin Converting Enzyme-2 (ACE-2) receptor for its attachment similar to SARS-CoV-1, which is followed by priming of spike protein by Transmembrane protease serine 2 (TMPRSS2) which can be targeted by a proven inhibitor of TMPRSS2, camostat. The current treatment strategy includes repurposing of existing drugs that were found to be effective against other RNA viruses like SARS, MERS, and Ebola. This review presents a critical analysis of druggable targets of SARS CoV-2, new drug discovery, development, and treatment opportunities for COVID-19.
Endo-β-1,3-glucanase from alkalophilic bacterium, Bacillus lehensis G1 (Blg32) composed of 284 amino acids with a predicted molecular mass of 31.6 kDa is expressed in Escherichia coli and purified to homogeneity. Herein, Blg32 characteristics, substrates and product specificity as well as structural traits that might be involved in the production of sugar molecules are analysed. This enzyme functions optimally at the temperature of 70 °C, pH value of 8.0 with its catalytic activity strongly enhanced by Mn2+. Remarkably, the purified enzyme is highly stable in high temperature and alkaline conditions. It exhibits the highest activity on laminarin (376.73 U/mg) followed by curdlan and yeast β-glucan. Blg32 activity increased by 62% towards soluble substrate (laminarin) compared to insoluble substrate (curdlan). Hydrolytic products of laminarin were oligosaccharides with degree of polymerisation (DP) of 1 to 5 with the main product being laminaritriose (DP3). This suggests that the active site of Blg32 could recognise up to five glucose units. High concentration of Blg32 mainly produces glucose whilst low concentration of Blg32 yields oligosaccharides with different DP (predominantly DP3). A theoretical structural model of Blg32 was constructed and structural analysis revealed that Trp156 is involved in multiple hydrophobic stacking interactions. The amino acid was predicted to participate in substrate recognition and binding. It was also exhibited that catalytic groove of Blg32 has a narrow angle, thus limiting the substrate binding reaction. All these properties and knowledge of the subsites are suggested to be related to the possible mode of action of how Blg32 produces glucooligosaccharides.
Sarcocystis scandentiborneensis sp. nov. was discovered in histological sections of striated musculature of treeshrews (Tupaia minor, T. tana) from Northern Borneo. Sarcocysts were cigar-shaped, 102 μm-545 μm long, and on average 53 μm in diameter. The striated cyst wall varied in thickness (2-10 μm), depending on whether the finger-like, villous protrusions (VP) were bent. Ultrastructurally, sarcocysts were similar to wall type 12 but basal microtubules extended into VPs that tapered off with a unique U-shaped, electron-dense apical structure. In phylogenetic trees of the nuclear 18S rRNA gene, S. scandentiborneensis formed a distinct branch within a monophyletic subclade of Sarcocystis spp. with (colubrid) snake-rodent life cycle. We mapped all intraspecific (two haplotypes) and interspecific nucleotide substitutions to the secondary structure of the 18S rRNA gene: in both cases, the highest variability occurred within helices V2 and V4 but intraspecific variability mostly related to transitions, while transition/transversion ratios between S. scandentiborneensis, S. zuoi, and S. clethrionomyelaphis were skewed towards transversions. Lack of relevant sequences restricted phylogenetic analysis of the mitochondrial Cytochrome C oxidase subunit I (COI) gene to include only one species of Sarcocystis recovered from a snake host (S. pantherophisi) with which the new species formed a sister relationship. We confirm the presence of the functionally important elements of the COI barcode amino acid sequence of S. scandentiborneensis, whereby the frequency of functionally important amino acids (Alanine, Serine) was markedly different to other taxa of the Sarcocystidae. We regard S. scandentiborneensis a new species, highlighting that structurally or functionally important aspects of the 18S rRNA and COI could expand their utility for delineation of species. We also address the question why treeshrews, believed to be close to primates, carry a parasite that is genetically close to a Sarcocystis lineage preferably developing in the Rodentia as intermediate hosts.
The ongoing outbreak of viral pneumonia in China and across the world is associated with a new coronavirus, SARS-CoV-21. This outbreak has been tentatively associated with a seafood market in Wuhan, China, where the sale of wild animals may be the source of zoonotic infection2. Although bats are probable reservoir hosts for SARS-CoV-2, the identity of any intermediate host that may have facilitated transfer to humans is unknown. Here we report the identification of SARS-CoV-2-related coronaviruses in Malayan pangolins (Manis javanica) seized in anti-smuggling operations in southern China. Metagenomic sequencing identified pangolin-associated coronaviruses that belong to two sub-lineages of SARS-CoV-2-related coronaviruses, including one that exhibits strong similarity in the receptor-binding domain to SARS-CoV-2. The discovery of multiple lineages of pangolin coronavirus and their similarity to SARS-CoV-2 suggests that pangolins should be considered as possible hosts in the emergence of new coronaviruses and should be removed from wet markets to prevent zoonotic transmission.
The aim of the present study was to analyze the sequence of the VP1 gene in enterovirus 71 (EV71) isolates and to explore their genetic evolution, so as to provide a scientific basis for the clinical prevention and treatment of hand, foot and mouth disease. The fecal samples of 590 patients with suspected hand, foot and mouth disease treated at Yan'an Hospital (Kunming, China) between January 2015 and December 2016 were collected and EV71 nucleic acid was detected by fluorescence PCR. The viral RNA of EV71-positive samples was extracted, the VP1 gene was amplified by PCR and the products were sequenced. The VP1 gene sequence was analyzed using DNAMAN and MEGA (version 4.0) software and homologous modeling was performed using Pymol software. A total of 50 EV71-positive samples were identified and the detection rate was 8.47% (50/590 cases). All of the 50 EV71 strains were of the C4 subtype. The genetic distance between the strains detected in the present study and EV71 strains detected in Beijing, Anhui and Malaysia was 0.01-0.03, while that between the strains detected in the present study and Australian strains was 2.11. Homologous modeling indicated that the amino acid sequence of the VP1 gene of the detected strains had a H144Y mutation. There was no significant genetic variation in the EV71 strain within the 2-year period. In conclusion, the EV71 strains detected in the present study was similar to that detected in Beijing, Anhui and Malaysia but different to that from Australia. A point mutation was present in the amino acid sequence of the VP1 gene.
Plasmodium vivax is the most widespread and difficult to treat cause of human malaria. The development of vaccines against the blood stages of P. vivax remains a key objective for the control and elimination of vivax malaria. Erythrocyte binding-like (EBL) protein family members such as Duffy binding protein (PvDBP) are of critical importance to erythrocyte invasion and have been the major target for vivax malaria vaccine development. In this study, we focus on another member of EBL protein family, P. vivax erythrocyte binding protein (PvEBP). PvEBP was first identified in Cambodian (C127) field isolates and has subsequently been showed its preferences for binding reticulocytes which is directly inhibited by antibodies. We analysed PvEBP sequence from 316 vivax clinical isolates from eight countries including China (n = 4), Ethiopia (n = 24), Malaysia (n = 53), Myanmar (n = 10), Papua New Guinea (n = 16), Republic of Korea (n = 10), Thailand (n = 174), and Vietnam (n = 25). PvEBP gene exhibited four different phenotypic clusters based on the insertion/deletion (indels) variation. PvEBP-RII (179-479 aa.) showed highest polymorphism similar to other EBL family proteins in various Plasmodium species. Whereas even though PvEBP-RIII-V (480-690 aa.) was the most conserved domain, that showed strong neutral selection pressure for gene purifying with significant population expansion. Antigenicity of both of PvEBP-RII (16.1%) and PvEBP-RIII-V (21.5%) domains were comparatively lower than other P. vivax antigen which expected antigens associated with merozoite invasion. Total IgG recognition level of PvEBP-RII was stronger than PvEBP-RIII-V domain, whereas total IgG inducing level was stronger in PvEBP-RIII-V domain. These results suggest that PvEBP-RII is mainly recognized by natural IgG for innate protection, whereas PvEBP-RIII-V stimulates IgG production activity by B-cell for acquired immunity. Overall, the low antigenicity of both regions in patients with vivax malaria likely reflects genetic polymorphism for strong positive selection in PvEBP-RII and purifying selection in PvEBP-RIII-V domain. These observations pose challenging questions to the selection of EBP and point out the importance of immune pressure and polymorphism required for inclusion of PvEBP as a vaccine candidate.
Phytol (PHY), a chlorophyll-derived diterpenoid, exhibits numerous pharmacological properties, including antioxidant, antimicrobial, and anticancer activities. This study evaluates the anti-diarrheal effect of phytol (PHY) along with its possible mechanism of action through in-vivo and in-silico models. The effect of PHY was investigated on castor oil-induced diarrhea in Swiss mice by using prazosin, propranolol, loperamide, and nifedipine as standards with or without PHY. PHY at 50 mg/kg (p.o.) and all other standards exhibit significant (p < 0.05) anti-diarrheal effect in mice. The effect was prominent in the loperamide and propranolol groups. PHY co-treated with prazosin and propranolol was found to increase in latent periods along with a significant reduction in diarrheal section during the observation period than other individual or combined groups. Furthermore, molecular docking studies also suggested that PHY showed better interactions with the α- and β-adrenergic receptors, especially with α-ADR1a and β-ADR1. In the former case, PHY showed interaction with hydroxyl group of Ser192 at a distance of 2.91Å, while in the latter it showed hydrogen bond interactions with Thr170 and Lys297 with a distance of 2.65 and 2.72Å, respectively. PHY exerted significant anti-diarrheal effect in Swiss mice, possibly through blocking α- and β-adrenergic receptors.
β-Glucosidases and β-mannosidases hydrolyze substrates that differ only in the epimer of the nonreducing terminal sugar moiety, but most such enzymes show a strong preference for one activity or the other. Rice Os3BGlu7 and Os7BGlu26 β-glycosidases show a less strong preference, but Os3BGlu7 and Os7BGlu26 prefer glucosides and mannosides, respectively. Previous studies of crystal structures with glucoimidazole (GIm) and mannoimidazole (MIm) complexes and metadynamic simulations suggested that Os7BGlu26 hydrolyzes mannosides via the B2,5 transition state (TS) conformation preferred for mannosides and glucosides via their preferred 4H3/4E TS conformation. However, MIm is weakly bound by both enzymes. In the present study, we found that MIm was not bound in the active site of crystallized Os3BGlu7, but GIm was tightly bound in the -1 subsite in a 4H3/4E conformation via hydrogen bonds with the surrounding residues. One-microsecond molecular dynamics simulations showed that GIm was stably bound in the Os3BGlu7 active site and the glycone-binding site with little distortion. In contrast, MIm initialized in the B2,5 conformation rapidly relaxed to a E3/4H3 conformation and moved out into a position in the entrance of the active site, where it bound more stably despite making fewer interactions. The lack of MIm binding in the glycone site in protein crystals and simulations implies that the energy required to distort MIm to the B2,5 conformation for optimal active site residue interactions is sufficient to offset the energy of those interactions in Os3BGlu7. This balance between distortion and binding energy may also provide a rationale for glucosidase versus mannosidase specificity in plant β-glycosidases.
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.