The l-2-haloacid dehalogenases (EC 3.8.1.2) specifically cleave carbon-halogen bonds in the L-isomers of halogenated organic acids. These enzymes have potential applications for the bioremediation and synthesis of various industrial products. One such enzyme is DehL, the l-2-haloacid dehalogenase from Rhizobium sp. RC1, which converts the L-isomers of 2-halocarboxylic acids into the corresponding D-hydroxycarboxylic acids. However, its catalytic mechanism has not been delineated, and to enhance its efficiency and utility for environmental and industrial applications, knowledge of its catalytic mechanism, which includes identification of its catalytic residues, is required. Using ab initio fragment molecular orbital calculations, molecular mechanics Poisson-Boltzmann surface area calculations, and classical molecular dynamic simulation of a three-dimensional model of DehL-l-2-chloropropionic acid complex, we predicted the catalytic residues of DehL and propose its catalytic mechanism. We found that when Asp13, Thr17, Met48, Arg51, and His184 were individually replaced with an alanine in silico, a significant decrease in the free energy of binding for the DehL-l-2-chloropropionic acid model complex was seen, indicating the involvement of these residues in catalysis and/or structural integrity of the active site. Furthermore, strong inter-fragment interaction energies calculated for Asp13 and L-2-chloropropionic acid, and for a water molecule and His184, and maintenance of the distances between atoms in the aforementioned pairs during the molecular dynamics run suggest that Asp13 acts as the nucleophile and His184 activates the water involved in DehL catalysis. The results of this study should be important for the rational design of a DehL mutant with improved catalytic efficiency.
Hepcidin is a principal regulator of iron homeostasis and its dysregulation has been recognised as a causative factor in cancers and iron disorders. The strategy of manipulating the presence of hepcidin peptide has been used for cancer treatment. However, this has demonstrated poor efficiency and has been short-lived in patients. Many studies reported using minihepcidin therapy as an alternative way to treat hepcidin dysregulation, but this was only applied to non-cancer patients. Highly conserved fish hepcidin protein, HepTH1-5, was investigated to determine its potential use in developing a hepcidin replacement for human hepcidin (Hepc25) and as a therapeutic agent by targeting the tumour suppressor protein, p53, through structure-function analysis. The authors found that HepTH1-5 is stably bound to ferroportin, compared to Hepc25, by triggering the ferroportin internalisation via Lys42 and Lys270 ubiquitination, in a similar manner to the Hepc25 activity. Moreover, the residues Ile24 and Gly24, along with copper and zinc ligands, interacted with similar residues, Lys24 and Asp1 of Hepc25, respectively, showing that those molecules are crucial to the hepcidin replacement strategy. HepTH1-5 interacts with p53 and activates its function through phosphorylation. This finding shows that HepTH1-5 might be involved in the apoptosis signalling pathway upon a DNA damage response. This study will be very helpful for understanding the mechanism of the hepcidin replacement and providing insights into the HepTH1-5 peptide as a new target for hepcidin and cancer therapeutics.Communicated by Ramaswamy H. Sarma.
Many studies reported that the activation of tumour suppressor protein, p53 induced the human hepcidin expression. However, its expression decreased when p53 was silenced in human hepatoma cells. Contrary to Tilapia hepcidin TH1-5, HepTH1-5 was previously reported to trigger the p53 activation through the molecular docking approach. The INhibitor of Growth (ING) family members are also shown to directly interact with p53 and promote cell cycle arrest, senescence, apoptosis and participate in DNA replication and DNA damage responses to suppress the tumour initiation and progression. However, the interrelation between INGs and HepTH1-5 remains unknown. Therefore, this study aims to identify the mechanism and their protein interactions using in silico approaches. The finding revealed that HepTH1-5 and its ligands had interacted mostly on hotspot residues of ING proteins which involved in histone modifications via acetylation, phosphorylation, and methylation. This proves that HepTH1-5 might implicate in an apoptosis signalling pathway and preserve the protein structure and function of INGs by reducing the perturbation of histone binding upon oxidative stress response. This study would provide theoretical guidance for the design and experimental studies to decipher the role of HepTH1-5 as a potential therapeutic agent for cancer therapy. Communicated by Ramaswamy H. Sarma.
Worldwide, breast cancer is the leading type of cancer among women. Overexpression of various prognostic indicators, including nuclear receptors, is linked to breast cancer features. To date, no effective drug has been discovered to block the proliferation of breast cancer cells. This study has been designed to discover target-based small molecular-like natural drug candidates that have anti-cancer potential without causing any serious side effects. A comprehensive substrate-based drug design was carried out to discover the potential plant compounds against the target breast cancer biomarkers including phytochemicals screening, active site identification, molecular docking, pharmacokinetic (PK) properties prediction, toxicity prediction, and molecular dynamics (MD) simulation approaches. Twenty plant compounds extracted from the rambutan (Nephelium lappaceum) were obtained from PubChem Database; and screened against the breast cancer biomarkers including estrogen receptor (ER), progesterone receptor (PR), and androgen receptor (AR). The best docking interaction was chosen based on the higher binding affinity. Analyzing the pharmacokinetic properties and toxicity prediction results indicated that the fifteen selected plant compounds have good potency without toxicity and are safe for humans. Four phytochemicals with a higher binding affinity were chosen for each breast cancer biomarker to study their stability in interaction with the target proteins using MD simulation. Among the above compounds, Ellagic acid showed the high binding affinity against all three breast cancer biomarkers.Communicated by Ramaswamy H. Sarma.
Pathologically relevant behaviors of Vibrio, such as the expression of virulence factors, biofilm production, and swarming motility, have been shown to be controlled by quorum sensing. The autoinducer-2 quorum sensing receptor protein LuxP is one of the target proteins for drug development to suppress the virulence of Vibrio. Here, we reported the potential molecular interaction of fatty acids identified in vibriosis-resistant grouper with LuxP. Fatty acid, 4-oxodocosahexaenoic acid (4R8) showed significant binding affinity toward LuxP (-6.0 kcal/mol) based on molecular docking analysis. The dynamic behavior of the protein-ligand complex was illustrated by molecular dynamic simulations. The fluctuation of the protein backbone, the stability of ligand binding, and hydrogen bond interactions were assessed, suggesting 4R8 possesses potential interaction with LuxP, which was supported by the low binding free energy (-29.144 kJ/mol) calculated using the molecular mechanics Poisson-Boltzmann surface area.
Here, we present a novel psychrophilic β-glucanase from Glaciozyma antarctica PI12 yeast that has been structurally modeled and analyzed in detail. To our knowledge, this is the first attempt to model a psychrophilic laminarinase from yeast. Because of the low sequence identity (<40%), a threading method was applied to predict a 3D structure of the enzyme using the MODELLER9v12 program. The results of a comparative study using other mesophilic, thermophilic, and hyperthermophilic laminarinases indicated several amino acid substitutions on the surface of psychrophilic laminarinase that totally increased the flexibility of its structure for efficient catalytic reactions at low temperatures. Whereas several structural factors in the overall structure can explain the weak thermal stability, this research suggests that the psychrophilic adaptation and catalytic activity at low temperatures were achieved through existence of longer loops and shorter or broken helices and strands, an increase in the number of aromatic and hydrophobic residues, a reduction in the number of hydrogen bonds and salt bridges, a higher total solvent accessible surface area, and an increase in the exposure of the hydrophobic side chains to the solvent. The results of comparative molecular dynamics simulation and principal component analysis confirmed the above strategies adopted by psychrophilic laminarinase to increase its catalytic efficiency and structural flexibility to be active at cold temperature.
Halogenated compounds are recalcitrant environmental pollutants prevalent in agricultural fields, waste waters and industrial by-products, but they can be degraded by dehalogenase-containing microbes. Notably, 2-haloalkanoic acid dehalogenases are employed to resolve optically active chloropropionates, as exemplified by the d-specific dehalogenase from Rhizobium sp. RCI (DehD), which acts on d-2-chloropropionate but not on its l-enantiomer. The catalytic residues of this dehalogenase responsible for its affinity toward d-2-chloropropionate have not been experimentally determined, although its three-dimensional crystal structure has been solved. For this study, we performed in silico docking and molecular dynamic simulations of complexes formed by this dehalogenase and d- or l-2-chloropropionate. Arg134 of the enzyme plays the key role in the stereospecific binding and Arg16 is in a position that would allow it to activate a water molecule for hydrolytic attack on the d-2-chloropropionate chiral carbon for release of the halide ion to yield l-2-hydroxypropionate. We propose that within the DehD active site, the NH group of Arg134 can form a hydrogen bond with the carboxylate of d-2-chloropropionate with a strength of ∼4 kcal/mol that may act as an acid-base catalyst, whereas, when l-2-chloropropionate is present, this bond cannot be formed. The significance of the present work is vital for rational design of this dehalogenase in order to confirm the involvement of Arg16 and Arg134 residues implicated in hydrolysis and binding of d-2-chloropropionate in the active site of d-specific dehalogenase from Rhizobium sp. RC1.
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.
Population surveys and species recognition for roosting bats are either based on capture, sight or optical-mechanical count methods. However, these methods are intrusive, are tedious and, at best, provide only statistical estimations. Here, we demonstrated the successful use of a terrestrial Light Detection and Ranging (LIDAR) laser scanner for remotely identifying and determining the exact population of roosting bats in caves. LIDAR accurately captured the 3D features of the roosting bats and their spatial distribution patterns in minimal light. The high-resolution model of the cave enabled an exact count of the visibly differentiated Hipposideros larvatus and their roosting pattern within the 3D topology of the cave. We anticipate that the development of LIDAR will open up new research possibilities by allowing researchers to study roosting behaviour within the topographical context of a cave's internal surface, thus facilitating rigorous quantitative characterisations of cave roosting behaviour.
Currently, there is no three-dimensional structure of D-specific dehalogenase (DehD) in the protein database. We modeled DehD using ab initio technique, performed molecular dynamics (MD) simulation and docking of D-2-chloropropionate (D-2CP), D-2-bromopropionate (D-2BP), monochloroacetate (MCA), monobromoacetate (MBA), 2,2-dichloropropionate (2,2-DCP), d,l-2,3-dichloropropionate (d,l-2,3-DCP), and 3-chloropropionate (3-CP) into the DehD active site. The sequences of DehD and D-2-haloacid dehalogenase (HadD) from Pseudomonas putida AJ1 have 15% sequence similarity. The model had 80% of the amino acid residues in the most favored region when compared to the crystal structure of DehI from Pseudomonas putida PP3. Docking analysis revealed that Arg107, Arg134 and Tyr135 interacted with D-2CP, and Glu20 activated the water molecule for hydrolytic dehalogenation. Single residue substitutions at 25-30 °C showed that polar residues of DehD were stable when substituted with nonpolar residues and showed a decrease in activity within the same temperature range. The molecular dynamics simulation of DehD and its variants showed that in R134A variant, Arg107 interacted with D-2CP, while in Y135A, Gln221 and Arg231 interacted with D-2CP. It is our emphatic belief that the new model will be useful for the rational design of DehDs with enhanced potentials.
Recombinant protein overexpression, an important biotechnological process, is ruled by complex biological rules which are mostly unknown, is in need of an intelligent algorithm so as to avoid resource-intensive lab-based trial and error experiments in order to determine the expression level of the recombinant protein. The purpose of this study is to propose a predictive model to estimate the level of recombinant protein overexpression for the first time in the literature using a machine learning approach based on the sequence, expression vector, and expression host. The expression host was confined to Escherichia coli which is the most popular bacterial host to overexpress recombinant proteins. To provide a handle to the problem, the overexpression level was categorized as low, medium and high. A set of features which were likely to affect the overexpression level was generated based on the known facts (e.g. gene length) and knowledge gathered from related literature. Then, a representative sub-set of features generated in the previous objective was determined using feature selection techniques. Finally a predictive model was developed using random forest classifier which was able to adequately classify the multi-class imbalanced small dataset constructed. The result showed that the predictive model provided a promising accuracy of 80% on average, in estimating the overexpression level of a recombinant protein.
The D-2-haloacid dehalogenase of D-specific dehalogenase (DehD) from Rhizobium sp. RC1 catalyses the hydrolytic dehalogenation of D-haloalkanoic acids, inverting the substrate-product configuration and thereby forming the corresponding L-hydroxyalkanoic acids. Our investigations were focused on DehD mutants: R134A and Y135A. We examined the possible interactions between these mutants with haloalkanoic acids and characterized the key catalytic residues in the wild-type dehalogenase, to design dehalogenase enzyme(s) with improved potential for dehalogenation of a wider range of substrates. Three natural substrates of wild-type DehD, specifically, monochloroacetate, monobromoacetate and D,L-2,3-dichloropropionate, and eight other non-natural haloalkanoic acids substrates of DehD, namely, L-2-chloropropionate; L-2-bromopropionate; 2,2-dichloropropionate; dichloroacetate; dibromoacetate; trichloroacetate; tribromoacetate; and 3-chloropropionate, were docked into the active site of the DehD mutants R134A and Y135A, which produced altered catalytic functions. The mutants interacted strongly with substrates that wild-type DehD does not interact with or degrade. The interaction was particularly enhanced with 3-chloropropionate, in addition to monobromoacetate, monochloroacetate and D,L-2,3-dichloropropionate. In summary, DehD variants R134A and Y135A demonstrated increased propensity for binding haloalkanoic acid and were non-stereospecific towards halogenated substrates. The improved characteristics in these mutants suggest that their functionality could be further exploited and harnessed in bioremediations and biotechnological applications.
To date, the genus Parvularcula consists of 6 species and no potential application of this genus was reported. Current study presents the genome sequence of Parvularcula flava strain NH6-79 T and its cellulolytic enzyme analysis. The assembled draft genome of strain NH6-79 T consists of 9 contigs and 7 scaffolds with 3.68 Mbp in size and GC content of 59.87%. From a total of 3,465 genes predicted, 96 of them are annotated as glycoside hydrolases (GHs). Within these GHs, 20 encoded genes are related to cellulosic biomass degradation, including 12 endoglucanases (5 GH10, 4 GH5, and 3 GH51), 2 exoglucanases (GH9) and 6 β-glucosidases (GH3). In addition, highest relative enzyme activities (endoglucanase, exoglucanase, and β-glucosidase) were observed at 27th hour when the strain was cultured in the carboxymethyl cellulose/Avicel®-containing medium for 45 h. The combination of genome analysis with experimental studies indicated the ability of strain NH6-79 T to produce extracellular endoglucanase, exoglucanase, and β-glucosidase. These findings suggest the potential of Parvularcula flava strain NH6-79 T in cellulose-containing biomass degradation and that the strain could be used in cellulosic biorefining process.
Dehalogenases continue to garner interest of the scientific community due to their potential applications in bioremediation of halogen-contaminated environment and in synthesis of various industrially relevant products. Example of such enzymes is DehL, an L-2-haloacid dehalogenase (EC 3.8.1.2) from Rhizobium sp. RC1 that catalyses the specific cleavage of halide ion from L-2-halocarboxylic acids to produce the corresponding D-2-hydroxycarboxylic acids. Recently, the catalytic residues of DehL have been identified and its catalytic mechanism has been fully elucidated. However, the enantiospecificity determinants of the enzyme remain unclear. This information alongside a well-defined catalytic mechanism are required for rational engineering of DehL for substrate enantiospecificity. Therefore, using quantum mechanics/molecular mechanics and molecular mechanics Poisson-Boltzmann surface area calculations, the current study theoretically investigated the molecular basis of DehL enantiospecificity. The study found that R51L mutation cancelled out the dehalogenation activity of DehL towards it natural substrate, L-2-chloropropionate. The M48R mutation, however introduced a new activity towards D-2-chloropropionate, conveying the possibility of inverting the enantiospecificity of DehL from L-to d-enantiomer with a minimum of two simultaneous mutations. The findings presented here will play important role in the rational design of DehL dehalogenase for improving substrate utility.
The Anoxybacillus sp. SK 3-4, previously isolated from a hot spring, was screened for its heavy metals resistance (Al(3+), Mn(2+), Cu(2+), Co(2+), Zn(2+), and Ni(2+)) and the strain was found to be most resistant to aluminum. Significant growth of the strain was observed when it was grown in medium containing aluminum (200 mg L(-1)-800 mg L(-1)) with relative growth rates ranging between 77% and 100%. A gene encoding the aluminum resistance protein (accession number: WP_021095658.1) was found in genome of strain SK 3-4, which revealed high sequence identity (>95%) to its homologues from Anoxybacillus species. Sequence comparisons with two functionally characterized aluminum resistance proteins, namely G2alt and ALU1-P, showed 97% and 81% of sequence identity, respectively. Four putative metal binding sites were detected in SK 3-4 aluminum resistance protein and G2alt at same amino acid residue positions of 186, 195, 198, and 201. Strain SK 3-4 was found to be able to remove aluminum from aqueous solution. This study demonstrated that Anoxybacillus sp. SK 3-4 could be applied in the treatment of aluminum contaminated wastewater.
Anoxybacillus sp. SK 3-4 is a Gram-positive, rod-shaped bacterium and a member of family Bacillaceae. We had previously reported that the strain is an aluminum resistant thermophilic bacterium. This is the first report to provide a detailed analysis of the global transcriptional response of Anoxybacillus when the cells were exposed to 600 mg L(-1) of aluminum. The transcriptome was sequenced using Illumina MiSeq sequencer. Total of 708 genes were differentially expressed (fold change >2.00) with 316 genes were up-regulated while 347 genes were down-regulated, in comparing to control with no aluminum added in the culture. Based on Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, the majority of genes encoding for cell metabolism such as glycolysis, sulfur metabolism, cysteine and methionine metabolism were up-regulated; while most of the gene associated with tricarboxylic acid cycle (TCA cycle) and valine, leucine and isoleucine metabolism were down-regulated. In addition, a significant number of the genes encoding ABC transporters, metal ions transporters, and some stress response proteins were also differentially expressed following aluminum exposure. The findings provide further insight and help us to understand on the resistance of Anoxybacillus sp. SK 3-4 toward aluminium.
Beta-glucosidase (BGL) is an important industrial enzyme for food, waste and biofuel processing. Jeotgalibacillus is an understudied halophilic genus, and no beta-glucosidase from this genus has been reported. A novel beta-glucosidase gene (1344 bp) from J. malaysiensis DSM 28777T was cloned and expressed in E. coli. The recombinant protein, referred to as BglD5, consists of a total 447 amino acids. BglD5 purified using a Ni-NTA column has an apparent molecular mass of 52 kDa. It achieved the highest activity at pH 7 and 65 °C. The activity and stability were increased when CaCl2 was supplemented to the enzyme. The enzyme efficiently hydrolyzed salicin and (1 → 4)-beta-glycosidic linkages such as in cellobiose, cellotriose, cellotetraose, cellopentose, and cellohexanose. Similar to many BGLs, BglD5 was not active towards polysaccharides such as Avicel, carboxymethyl cellulose, Sigmacell cellulose 101, alpha-cellulose and xylan. When BglD5 blended with Cellic® Ctec2, the total sugars saccharified from oil palm empty fruit bunches (OPEFB) was enhanced by 4.5%. Based on sequence signatures and tree analyses, BglD5 belongs to the Glycoside Hydrolase family 1. This enzyme is a novel beta-glucosidase attributable to its relatively low sequence similarity with currently known beta-glucosidases, where the closest characterized enzyme is the DT-Bgl from Anoxybacillus sp. DT3-1.
Rhodothermaceae bacterium RA is a halo-thermophile isolated from a saline hot spring. Previously, the genome of this bacterium was sequenced using a HiSeq 2500 platform culminating in 91 contigs. In this report, we report on the resequencing of its complete genome using a PacBio RSII platform. The genome has a GC content of 68.3%, is 4,653,222 bp in size, and encodes 3711 genes. We are interested in understanding the carbohydrate metabolic pathway, in particular the lignocellulosic biomass degradation pathway. Strain RA harbors 57 glycosyl hydrolase (GH) genes that are affiliated with 30 families. The bacterium consists of cellulose-acting (GH 3, 5, 9, and 44) and hemicellulose-acting enzymes (GH 3, 10, and 43). A crude cell-free extract of the bacterium exhibited endoglucanase, xylanase, β-glucosidase, and β-xylosidase activities. The complete genome information coupled with biochemical assays confirms that strain RA is able to degrade cellulose and xylan. Therefore, strain RA is another excellent member of family Rhodothermaceae as a repository of novel and thermostable cellulolytic and hemicellulolytic enzymes.